Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer 3 // * License and Disclaimer * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. 9 // * include a list of copyright holders. * 10 // * 10 // * * 11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. * 17 // * 17 // * * 18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************* 24 // ******************************************************************** 25 // 25 // 26 // 26 // >> 27 // GEANT4 tag $Name: geant4-08-03-patch-02 $ >> 28 // >> 29 // 27 // GEANT4 physics class: G4PhotoNuclearCrossSe 30 // GEANT4 physics class: G4PhotoNuclearCrossSection -- header file 28 // Created: M.V. Kossov, CERN/ITEP(Moscow), 10 31 // Created: M.V. Kossov, CERN/ITEP(Moscow), 10-OCT-01 29 // The last update: M.V. Kossov, CERN/ITEP (Mo 32 // The last update: M.V. Kossov, CERN/ITEP (Moscow) 17-May-02 >> 33 // 30 34 31 #ifndef G4PhotoNuclearCrossSection_h 35 #ifndef G4PhotoNuclearCrossSection_h 32 #define G4PhotoNuclearCrossSection_h 1 36 #define G4PhotoNuclearCrossSection_h 1 33 37 34 #include "G4VCrossSectionDataSet.hh" 38 #include "G4VCrossSectionDataSet.hh" >> 39 /////////#include "G4HadronCrossSections.hh" 35 #include "G4DynamicParticle.hh" 40 #include "G4DynamicParticle.hh" 36 #include "G4Element.hh" 41 #include "G4Element.hh" >> 42 //#include "G4QPDGCode.hh" 37 #include "G4ParticleTable.hh" 43 #include "G4ParticleTable.hh" 38 #include "G4NucleiProperties.hh" 44 #include "G4NucleiProperties.hh" 39 #include "G4NistManager.hh" << 45 #include "G4NucleiPropertiesTable.hh" 40 #include <vector> 46 #include <vector> 41 47 42 class G4PhotoNuclearCrossSection : public G4VC 48 class G4PhotoNuclearCrossSection : public G4VCrossSectionDataSet 43 { 49 { 44 public: 50 public: 45 << 46 G4PhotoNuclearCrossSection(); << 47 ~G4PhotoNuclearCrossSection() override; << 48 << 49 static const char* Default_Name() {return << 50 << 51 void CrossSectionDescription(std::ostream& << 52 << 53 G4bool IsIsoApplicable(const G4DynamicPart << 54 const G4Element* el << 55 const G4Material* mat = nullptr) over << 56 << 57 G4bool IsElementApplicable(const G4Dynamic << 58 const G4Materia << 59 << 60 G4double GetIsoCrossSection(const G4Dynami << 61 G4int Z, G4int << 62 const G4Isotop << 63 const G4Elemen << 64 const G4Materi << 65 << 66 G4double GetElementCrossSection(const G4Dy << 67 const G4Ma << 68 << 69 G4double ComputeElementXSection(G4double e << 70 << 71 G4double ComputeIsoXSection(G4double energ << 72 << 73 G4PhotoNuclearCrossSection& operator= << 74 (const G4PhotoNuclearCrossSection& right) << 75 G4PhotoNuclearCrossSection(const G4PhotoNu << 76 51 >> 52 G4PhotoNuclearCrossSection() // Constructor @@?? >> 53 { >> 54 //theHadronCrossSections = G4HadronCrossSections::Instance(); >> 55 } >> 56 >> 57 ~G4PhotoNuclearCrossSection() {} >> 58 >> 59 >> 60 G4bool IsApplicable(const G4DynamicParticle* particle, const G4Element* ) >> 61 { >> 62 return IsZAApplicable(particle, 0, 0); >> 63 } >> 64 >> 65 G4bool IsZAApplicable(const G4DynamicParticle* particle, >> 66 G4double /*ZZ*/, G4double /*AA*/) >> 67 { >> 68 G4bool result = false; >> 69 if( particle->GetDefinition()->GetPDGEncoding()==22) result = true; >> 70 return result; >> 71 } >> 72 >> 73 >> 74 G4double GetCrossSection(const G4DynamicParticle* particle, >> 75 const G4Element* element, G4double temp = 0.); >> 76 >> 77 G4double GetIsoZACrossSection(const G4DynamicParticle* particle, >> 78 G4double ZZ, G4double AA, >> 79 G4double /*aTemperature*/); >> 80 >> 81 >> 82 void BuildPhysicsTable(const G4ParticleDefinition&) {} >> 83 >> 84 void DumpPhysicsTable(const G4ParticleDefinition&) {} >> 85 >> 86 private: >> 87 G4int GetFunctions(G4double a, G4double* y, G4double* z); >> 88 //G4double LinearFit(G4double X, G4int N, const G4double* XN, const G4double* YN); >> 89 G4double EquLinearFit(G4double X, G4int N,const G4double X0,const G4double XD, const G4double* Y); >> 90 G4double ThresholdEnergy(G4int Z, G4int N); >> 91 >> 92 // Body 77 private: 93 private: 78 << 94 static G4int lastN; // The last N of calculated nucleus 79 G4int GetFunctions(G4double a, G4double* y << 95 static G4int lastZ; // The last Z of calculated nucleus 80 G4double EquLinearFit(G4double X, G4int N, << 96 static G4double lastSig; // Last value of the Cross Section 81 const G4double XD, c << 97 static G4double* lastGDR; // Pointer to the last array of GDR cross sections 82 G4double ThresholdEnergy(G4int Z, G4int N) << 98 static G4double* lastHEN; // Pointer to the last array of HEn cross sections 83 << 99 static G4double lastE; // Last used in the cross section Energy 84 G4int lastZ = 0; // The last Z << 100 static G4double lastTH; // Last value of the Energy Threshold (A-dependent) 85 G4double lastSig = 0.0; // Last value << 101 static G4double lastSP; // Last value of the ShadowingPomeron (A-dependent) 86 G4double* lastGDR = nullptr; // Pointer to << 102 87 G4double* lastHEN = nullptr; // Pointer to << 103 //G4HadronCrossSections* theHadronCrossSections; 88 G4double lastE = 0.0; // Last used << 89 G4double lastTH = 0.0; // Last value << 90 G4double lastSP = 0.0; // Last value << 91 << 92 // Vector of pointers to the GDRPhotonucle << 93 std::vector <G4double*> GDR; << 94 << 95 // store deuteron, triton, He3 XS << 96 G4double* deuteron_GDR = nullptr; << 97 G4double* deuteron_HR = nullptr; << 98 G4double deuteron_TH = 0.0; << 99 G4double deuteron_SP = 0.0; << 100 G4double* triton_GDR = nullptr; << 101 G4double* triton_HR = nullptr; << 102 G4double triton_TH = 0.0; << 103 G4double triton_SP = 0.0; << 104 G4double* he3_GDR = nullptr; << 105 G4double* he3_HR = nullptr; << 106 G4double he3_TH = 0.0; << 107 G4double he3_SP = 0.0; << 108 << 109 // Vector of pointers to the HighEnPhotonu << 110 std::vector <G4double*> HEN; << 111 << 112 std::vector <G4double> spA; // shadowing << 113 std::vector <G4double> eTH; // energy t << 114 << 115 G4NistManager* nistmngr; << 116 << 117 G4double mNeut; << 118 G4double mProt; << 119 }; 104 }; >> 105 >> 106 // Gives the threshold energy for different nuclei (min of p- and n-threshold) >> 107 inline G4double G4PhotoNuclearCrossSection::ThresholdEnergy(G4int Z, G4int N) >> 108 { >> 109 // CHIPS - Direct GEANT >> 110 //static const G4double mNeut = G4QPDGCode(2112).GetMass(); >> 111 //static const G4double mProt = G4QPDGCode(2212).GetMass(); >> 112 static const G4double mNeut = G4NucleiProperties::GetNuclearMass(1,0); >> 113 static const G4double mProt = G4NucleiProperties::GetNuclearMass(1,1); >> 114 // --------- >> 115 static const G4double infEn = 9.e27; >> 116 >> 117 G4int A=Z+N; >> 118 if(A<1) return infEn; >> 119 else if(A==1) return 134.9766; // Pi0 threshold for the nucleon >> 120 // CHIPS - Direct GEANT >> 121 //G4double mT= G4QPDGCode(111).GetNuclMass(Z,N,0); >> 122 G4double mT= 0.; >> 123 if(G4NucleiPropertiesTable::IsInTable(Z,A)) mT=G4NucleiProperties::GetNuclearMass(A,Z); >> 124 else >> 125 { >> 126 G4cerr<<"G4PhotoNucCrossSect.hh::ThreshEn:Z="<<Z<<",A="<<A<<" element isn't in G4NucPr"<<G4endl; >> 127 return 0.; // If it is not in the Table of Stable Nuclei, then the Threshold=inf >> 128 } >> 129 // --------- >> 130 G4double mP= infEn; >> 131 //if(Z) mP= G4QPDGCode(111).GetNuclMass(Z-1,N,0); >> 132 >> 133 if(Z && G4NucleiPropertiesTable::IsInTable(Z-1,A-1)) >> 134 { >> 135 mP = G4NucleiProperties::GetNuclearMass(A-1,Z-1); >> 136 } >> 137 else >> 138 { >> 139 G4cerr << "G4PhotoNucCrossSect.hh::ThrEn:Z=" << Z-1 << ",A=" >> 140 << A-1 << " element isn't in G4NucP" << G4endl; >> 141 } >> 142 G4double mN= infEn; >> 143 //if(N) mN= G4QPDGCode(111).GetNuclMass(Z,N-1,0); >> 144 if(N&&G4NucleiPropertiesTable::IsInTable(Z,A-1)) mN=G4NucleiProperties::GetNuclearMass(A-1,Z); >> 145 else >> 146 { >> 147 G4cerr<<"G4PhotoNucCrossSect.hh::ThreshEn:Z="<<Z<<",A="<<A-1<<" element isn't in G4NuP"<<G4endl; >> 148 } >> 149 G4double dP= mP+mProt-mT; >> 150 G4double dN= mN+mNeut-mT; >> 151 if(dP<dN)dN=dP; >> 152 return dN; >> 153 } 120 154 121 #endif 155 #endif 122 156