Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // * License and Disclaimer 4 // * 5 // * The Geant4 software is copyright of th 6 // * the Geant4 Collaboration. It is provided 7 // * conditions of the Geant4 Software License 8 // * LICENSE and available at http://cern.ch/ 9 // * include a list of copyright holders. 10 // * 11 // * Neither the authors of this software syst 12 // * institutes,nor the agencies providing fin 13 // * work make any representation or warran 14 // * regarding this software system or assum 15 // * use. Please see the license in the file 16 // * for the full disclaimer and the limitatio 17 // * 18 // * This code implementation is the result 19 // * technical work of the GEANT4 collaboratio 20 // * By using, copying, modifying or distri 21 // * any work based on the software) you ag 22 // * use in resulting scientific publicati 23 // * acceptance of all terms of the Geant4 Sof 24 // ******************************************* 25 // 26 // ------------------------------------------- 27 // 28 // GEANT4 Class file 29 // 30 // 31 // File name: G4BGGPionInelasticXS 32 // 33 // Author: Vladimir Ivanchenko 34 // 35 // Creation date: 01.10.2003 36 // Modifications: 37 // 38 // ------------------------------------------- 39 // 40 41 #include "G4BGGPionInelasticXS.hh" 42 #include "G4SystemOfUnits.hh" 43 #include "G4ComponentGGHadronNucleusXsc.hh" 44 #include "G4UPiNuclearCrossSection.hh" 45 #include "G4HadronNucleonXsc.hh" 46 #include "G4NuclearRadii.hh" 47 48 #include "G4Proton.hh" 49 #include "G4PionPlus.hh" 50 #include "G4PionMinus.hh" 51 #include "G4NistManager.hh" 52 #include "G4Pow.hh" 53 #include "G4HadronicParameters.hh" 54 55 G4double G4BGGPionInelasticXS::theGlauberFacPi 56 G4double G4BGGPionInelasticXS::theGlauberFacPi 57 G4double G4BGGPionInelasticXS::theLowEPiPlus[9 58 G4double G4BGGPionInelasticXS::theLowEPiMinus[ 59 G4int G4BGGPionInelasticXS::theA[93] = {0}; 60 61 G4BGGPionInelasticXS::G4BGGPionInelasticXS(con 62 : G4VCrossSectionDataSet("BarashenkovGlauberG 63 { 64 verboseLevel = 0; 65 fGlauberEnergy = 91.*CLHEP::GeV; 66 fLowEnergy = 20.*CLHEP::MeV; 67 fLowestEnergy = 1.*CLHEP::MeV; 68 SetMinKinEnergy(0.0); 69 SetMaxKinEnergy( G4HadronicParameters::Insta 70 71 fPion = new G4UPiNuclearCrossSection(); 72 fGlauber = new G4ComponentGGHadronNucleusXsc 73 fHadron = new G4HadronNucleonXsc(); 74 75 fG4pow = G4Pow::GetInstance(); 76 77 theProton = G4Proton::Proton(); 78 thePiPlus = G4PionPlus::PionPlus(); 79 isPiplus = (p == thePiPlus); 80 SetForAllAtomsAndEnergies(true); 81 82 if (0 == theA[0]) { Initialise(); } 83 } 84 85 //....oooOO0OOooo........oooOO0OOooo........oo 86 87 G4BGGPionInelasticXS::~G4BGGPionInelasticXS() 88 { 89 delete fHadron; 90 } 91 92 //....oooOO0OOooo........oooOO0OOooo........oo 93 94 G4bool 95 G4BGGPionInelasticXS::IsElementApplicable(cons 96 cons 97 { 98 return true; 99 } 100 101 //....oooOO0OOooo........oooOO0OOooo........oo 102 103 G4bool G4BGGPionInelasticXS::IsIsoApplicable(c 104 G 105 c 106 c 107 { 108 return (1 == Z); 109 } 110 111 //....oooOO0OOooo........oooOO0OOooo........oo 112 113 G4double 114 G4BGGPionInelasticXS::GetElementCrossSection(c 115 G 116 { 117 // this method should be called only for Z > 118 119 G4double cross = 0.0; 120 G4double ekin = std::max(dp->GetKineticEnerg 121 G4int Z = std::min(ZZ, 92); 122 123 if(1 == Z) { 124 cross = 1.0115*GetIsoCrossSection(dp,1,1); 125 } else if(ekin < fLowEnergy) { 126 cross = (isPiplus) ? theLowEPiPlus[Z]*Coul 127 : theLowEPiMinus[Z]*FactorPiMinus(ekin); 128 } else if(ekin > fGlauberEnergy) { 129 cross = (isPiplus) ? theGlauberFacPiPlus[Z 130 cross *= fGlauber->GetInelasticGlauberGrib 131 } else { 132 cross = fPion->GetInelasticCrossSection(dp 133 } 134 #ifdef G4VERBOSE 135 if(verboseLevel > 1) { 136 G4cout << "G4BGGPionInelasticXS::GetCrossS 137 << dp->GetDefinition()->GetParticle 138 << " Ekin(GeV)= " << dp->GetKineti 139 << " in nucleus Z= " << Z << " A= 140 << " XS(b)= " << cross/barn 141 << G4endl; 142 } 143 #endif 144 return cross; 145 } 146 147 //....oooOO0OOooo........oooOO0OOooo........oo 148 149 G4double 150 G4BGGPionInelasticXS::GetIsoCrossSection(const 151 G4int 152 const 153 const 154 const 155 { 156 // this method should be called only for Z = 157 fHadron->HadronNucleonXscNS(dp->GetDefinitio 158 dp->GetKineticEn 159 G4double cross = A*fHadron->GetInelasticHadr 160 161 #ifdef G4VERBOSE 162 if(verboseLevel > 1) { 163 G4cout << "G4BGGPionInelasticXS::GetCrossS 164 << dp->GetDefinition()->GetParticle 165 << " Ekin(GeV)= " << dp->GetKineti 166 << " in nucleus Z=1 A=" << A 167 << " XS(b)= " << cross/barn 168 << G4endl; 169 } 170 #endif 171 return cross; 172 } 173 174 //....oooOO0OOooo........oooOO0OOooo........oo 175 176 void G4BGGPionInelasticXS::BuildPhysicsTable(c 177 { 178 if(verboseLevel > 1) { 179 G4cout << "G4BGGPionInelasticXS::BuildPhys 180 << p.GetParticleName() << G4endl; 181 } 182 if(&p == G4PionPlus::PionPlus() || &p == G4P 183 isPiplus = (&p == G4PionPlus::PionPlus()); 184 } else { 185 G4ExceptionDescription ed; 186 ed << "This BGG cross section is applicabl 187 << p.GetParticleName() << G4endl; 188 G4Exception("G4BGGPionInelasticXS::BuildPh 189 FatalException, ed); 190 } 191 } 192 193 //....oooOO0OOooo........oooOO0OOooo........oo 194 195 void G4BGGPionInelasticXS::Initialise() 196 { 197 theA[0] = theA[1] = 1; 198 G4ThreeVector mom(0.0,0.0,1.0); 199 G4DynamicParticle dp(thePiPlus, mom, fGlaube 200 201 G4NistManager* nist = G4NistManager::Instanc 202 G4double csup, csdn; 203 204 for (G4int iz=2; iz<93; ++iz) { 205 G4int A = G4lrint(nist->GetAtomicMassAmu(i 206 theA[iz] = A; 207 208 csup = fGlauber->GetInelasticGlauberGribov 209 csdn = fPion->GetInelasticCrossSection(&dp 210 theGlauberFacPiPlus[iz] = csdn/csup; 211 } 212 213 dp.SetDefinition(G4PionMinus::PionMinus()); 214 for (G4int iz=2; iz<93; ++iz) { 215 csup = fGlauber->GetInelasticGlauberGribov 216 csdn = fPion->GetInelasticCrossSection(&dp 217 theGlauberFacPiMinus[iz] = csdn/csup; 218 219 if(verboseLevel > 1) { 220 G4cout << "Z= " << iz << " A= " << the 221 << " factorPiPlus= " << theGlauberFacPi 222 << " factorPiMinus= " << theGlauberFacP 223 << G4endl; 224 } 225 } 226 227 theLowEPiPlus[1] = theLowEPiMinus[1]= 1.0; 228 dp.SetDefinition(thePiPlus); 229 dp.SetKineticEnergy(fLowEnergy); 230 for (G4int iz=2; iz<93; ++iz) { 231 theLowEPiPlus[iz] = fPion->GetInelasticCro 232 /CoulombFactorPiPlus(fLowEnergy, iz); 233 } 234 235 dp.SetDefinition(G4PionMinus::PionMinus()); 236 for (G4int iz=2; iz<93; ++iz) { 237 theLowEPiMinus[iz] = fPion->GetInelasticCr 238 /FactorPiMinus(fLowEnergy); 239 240 if (verboseLevel > 1) { 241 G4cout << "Z= " << iz << " A= " << the 242 << " LowEtorPiPlus= " << theLowEPiPlus[ 243 << " LowEtorPiMinus= " << theLowEPiMinu 244 << G4endl; 245 } 246 } 247 } 248 249 //....oooOO0OOooo........oooOO0OOooo........oo 250 251 G4double G4BGGPionInelasticXS::CoulombFactorPi 252 { 253 return (kinEnergy > 0.0) ? 254 G4NuclearRadii::CoulombFactor(Z, theA[Z], 255 } 256 257 //....oooOO0OOooo........oooOO0OOooo........oo 258 259 G4double G4BGGPionInelasticXS::FactorPiMinus(G 260 { 261 return 1.0/std::sqrt(kinEnergy); 262 } 263 264 //....oooOO0OOooo........oooOO0OOooo........oo 265 266 void 267 G4BGGPionInelasticXS::CrossSectionDescription( 268 { 269 outFile << "The Barashenkov-Glauber-Gribov c 270 << "pion scattering from nuclei at a 271 << "parameterization is used below 9 272 << "parameterization is used above 9 273 } 274 275 //....oooOO0OOooo........oooOO0OOooo........oo 276