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