00001 // 00002 // ******************************************************************** 00003 // * License and Disclaimer * 00004 // * * 00005 // * The Geant4 software is copyright of the Copyright Holders of * 00006 // * the Geant4 Collaboration. It is provided under the terms and * 00007 // * conditions of the Geant4 Software License, included in the file * 00008 // * LICENSE and available at http://cern.ch/geant4/license . These * 00009 // * include a list of copyright holders. * 00010 // * * 00011 // * Neither the authors of this software system, nor their employing * 00012 // * institutes,nor the agencies providing financial support for this * 00013 // * work make any representation or warranty, express or implied, * 00014 // * regarding this software system or assume any liability for its * 00015 // * use. Please see the license in the file LICENSE and URL above * 00016 // * for the full disclaimer and the limitation of liability. * 00017 // * * 00018 // * This code implementation is the result of the scientific and * 00019 // * technical work of the GEANT4 collaboration. * 00020 // * By using, copying, modifying or distributing the software (or * 00021 // * any work based on the software) you agree to acknowledge its * 00022 // * use in resulting scientific publications, and indicate your * 00023 // * acceptance of all terms of the Geant4 Software license. * 00024 // ******************************************************************** 00025 // 00026 #ifndef G4LENDCrossSection_h 00027 #define G4LENDCrossSection_h 1 00028 00029 // Class Description 00030 // Cross Sections for a LEND (Low Energy Nuclear Data) 00031 // LEND is Geant4 interface for GIDI (General Interaction Data Interface) 00032 // which gives a discription of nuclear and atomic reactions, such as 00033 // Binary collision cross sections 00034 // Particle number multiplicity distributions of reaction products 00035 // Energy and angular distributions of reaction products 00036 // Derived calculational constants 00037 // GIDI is developped at Lawrence Livermore National Laboratory 00038 // Class Description - End 00039 00040 // 071025 First implementation done by T. Koi (SLAC/SCCS) 00041 // 101118 Name modifications for release T. Koi (SLAC/PPA) 00042 00043 #include "G4LENDHeader.hh" 00044 00045 #include "G4LENDManager.hh" 00046 #include "G4LENDUsedTarget.hh" 00047 00048 #include "G4VCrossSectionDataSet.hh" 00049 #include "G4DynamicParticle.hh" 00050 #include "G4Element.hh" 00051 #include "G4Material.hh" 00052 #include "G4ParticleDefinition.hh" 00053 #include "G4PhysicsTable.hh" 00054 00055 #include <map> 00056 00057 class G4LENDCrossSection : public G4VCrossSectionDataSet 00058 { 00059 00060 // 00061 // 00062 // G4bool IsIsoApplicable(const G4DynamicParticle*, G4int Z, G4int A, 00063 // const G4Element* elm = 0, 00064 // const G4Material* mat = 0); 00065 // 00066 // G4double GetIsoCrossSection(const G4DynamicParticle*, G4int Z, G4int A, 00067 // const G4Element* elm = 0, 00068 // const G4Material* mat = 0); 00069 // 00070 00071 public: 00072 00073 G4LENDCrossSection(const G4String name = "" ); 00074 00075 ~G4LENDCrossSection(); 00076 00077 //G4bool IsApplicable( const G4DynamicParticle* , const G4Element* ); 00078 //TK110811 00079 G4bool IsIsoApplicable( const G4DynamicParticle* , G4int /*Z*/ , G4int /*A*/ , 00080 const G4Element* , const G4Material* ); 00081 G4double GetIsoCrossSection( const G4DynamicParticle*, G4int /*Z*/, G4int /*A*/ , 00082 const G4Isotope* , const G4Element* , const G4Material* ); 00083 00084 void BuildPhysicsTable( const G4ParticleDefinition& ); 00085 00086 void DumpPhysicsTable( const G4ParticleDefinition& ); 00087 00088 00089 //TK110810 00090 //G4double GetCrossSection( const G4DynamicParticle* , const G4Element* , G4double aT ); 00091 //G4double GetCrossSection(const G4DynamicParticle*, G4int , const G4Material* ); 00092 00093 00094 void ChangeDefaultEvaluation( G4String name_tmp ){ default_evaluation = name_tmp; }; 00095 void AllowNaturalAbundanceTarget(){ allow_nat = true; }; 00096 void AllowAnyCandidateTarget(){ allow_any = true; }; 00097 00098 //Hadronic Framework still does not handle isotope in GPIL 00099 //G4VDiscreteProcess::PostStepGetPhysicalInteractionLenght() 00100 // G4HadronicProcess::GetMeanFreePath() 00101 // G4double GetCrossSection(const G4DynamicParticle*, const G4Material*) 00102 // G4double GetCrossSection(const G4DynamicParticle*, const G4Element*, G4double aTemperature); 00103 00104 //TK110810 00105 //G4bool IsIsoApplicable( const G4DynamicParticle* , G4int /*ZZ*/, G4int /*AA*/) { return true; } 00106 //G4bool IsZandAApplicable( const G4DynamicParticle* , G4int /*ZZ*/, G4int /*AA*/, const G4Element* , const G4Material* ) { return true; } 00107 //TK110810 00108 //G4double GetZandACrossSection(const G4DynamicParticle* , G4int /*Z*/, G4int /*A*/, G4double aTemperature); 00109 //G4double GetZandACrossSection(const G4DynamicParticle* , G4int /*Z*/, G4int /*A*/, const G4Material* mat); 00110 //TK110810 00111 //G4double GetIsoCrossSection( const G4DynamicParticle* , const G4Isotope* , G4double ); 00112 // G4double GetIsoCrossSection( const G4DynamicParticle* , const G4Isotope* , const G4Material* mat ); 00113 00114 private: 00115 00116 std::map< G4int , G4LENDUsedTarget* > usedTarget_map; 00117 00118 G4String default_evaluation; 00119 G4bool allow_nat; 00120 G4bool allow_any; 00121 00122 G4LENDManager* lend_manager; 00123 void recreate_used_target_map(); 00124 00125 protected : 00126 //G4String name; 00127 G4ParticleDefinition* proj; 00128 void create_used_target_map(); 00129 virtual G4double getLENDCrossSection( G4GIDI_target* , G4double , G4double ) { return 0.0; }; 00130 // elow ehigh xs_elow xs_ehigh ke (<elow) 00131 G4double GetUltraLowEnergyExtrapolatedXS( G4double , G4double , G4double , G4double , G4double ); 00132 }; 00133 #endif