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 // $Id: G4BGGNucleonInelasticXS.cc 79192 2014-02-20 10:07:18Z gcosmo $ >> 27 // 26 // ------------------------------------------- 28 // ------------------------------------------------------------------- 27 // 29 // 28 // GEANT4 Class file 30 // GEANT4 Class file 29 // 31 // 30 // 32 // 31 // File name: G4BGGNucleonInelasticXS 33 // File name: G4BGGNucleonInelasticXS 32 // 34 // 33 // Author: Vladimir Ivanchenko 35 // Author: Vladimir Ivanchenko 34 // 36 // 35 // Creation date: 13.03.2007 37 // Creation date: 13.03.2007 36 // Modifications: 38 // Modifications: 37 // 39 // 38 // 40 // 39 // ------------------------------------------- 41 // ------------------------------------------------------------------- 40 // 42 // 41 43 42 #include "G4BGGNucleonInelasticXS.hh" 44 #include "G4BGGNucleonInelasticXS.hh" 43 #include "G4SystemOfUnits.hh" 45 #include "G4SystemOfUnits.hh" 44 #include "G4ComponentGGHadronNucleusXsc.hh" << 46 #include "G4GlauberGribovCrossSection.hh" 45 #include "G4NucleonNuclearCrossSection.hh" 47 #include "G4NucleonNuclearCrossSection.hh" 46 #include "G4HadronNucleonXsc.hh" 48 #include "G4HadronNucleonXsc.hh" 47 #include "G4ComponentSAIDTotalXS.hh" 49 #include "G4ComponentSAIDTotalXS.hh" 48 #include "G4Proton.hh" 50 #include "G4Proton.hh" 49 #include "G4Neutron.hh" 51 #include "G4Neutron.hh" 50 #include "G4NistManager.hh" 52 #include "G4NistManager.hh" 51 #include "G4Material.hh" 53 #include "G4Material.hh" 52 #include "G4Element.hh" 54 #include "G4Element.hh" 53 #include "G4Isotope.hh" 55 #include "G4Isotope.hh" 54 #include "G4Log.hh" 56 #include "G4Log.hh" 55 #include "G4Exp.hh" 57 #include "G4Exp.hh" 56 #include "G4NuclearRadii.hh" << 57 58 58 //....oooOO0OOooo........oooOO0OOooo........oo << 59 #include "G4CrossSectionDataSetRegistry.hh" 59 60 60 namespace << 61 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 61 { << 62 const G4double llog10 = G4Log(10.); << 63 } << 64 62 65 G4double G4BGGNucleonInelasticXS::theGlauberFa << 63 const G4double llog10 = G4Log(10.); 66 G4double G4BGGNucleonInelasticXS::theCoulombFa << 67 G4double G4BGGNucleonInelasticXS::theGlauberFa << 68 G4double G4BGGNucleonInelasticXS::theCoulombFa << 69 G4int G4BGGNucleonInelasticXS::theA[93] = {0}; << 70 64 71 G4BGGNucleonInelasticXS::G4BGGNucleonInelastic 65 G4BGGNucleonInelasticXS::G4BGGNucleonInelasticXS(const G4ParticleDefinition* p) 72 : G4VCrossSectionDataSet("BarashenkovGlauberG << 66 : G4VCrossSectionDataSet("Barashenkov-Glauber") 73 { 67 { 74 verboseLevel = 0; 68 verboseLevel = 0; 75 fGlauberEnergy = 91.*CLHEP::GeV; << 69 fGlauberEnergy = 91.*GeV; 76 fLowEnergy = 14.*CLHEP::MeV; << 70 fLowEnergy = 14.*MeV; 77 << 71 fHighEnergy = 5.*GeV; 78 fNucleon = new G4NucleonNuclearCrossSection( << 72 fSAIDHighEnergyLimit = 1.3*GeV; 79 fGlauber = new G4ComponentGGHadronNucleusXsc << 73 fLowestXSection = millibarn; 80 fHadron = new G4HadronNucleonXsc(); << 74 for (G4int i = 0; i < 93; ++i) { >> 75 theGlauberFac[i] = 0.0; >> 76 theCoulombFac[i] = 0.0; >> 77 theA[i] = 1; >> 78 } >> 79 fNucleon = 0; >> 80 fGlauber = 0; >> 81 fHadron = 0; >> 82 fSAID = 0; 81 83 >> 84 particle = p; 82 theProton= G4Proton::Proton(); 85 theProton= G4Proton::Proton(); 83 isProton = (theProton == p); << 86 isProton = false; 84 SetForAllAtomsAndEnergies(true); << 87 isInitialized = false; 85 << 86 if (0 == theA[0]) { Initialise(); } << 87 } 88 } 88 89 89 //....oooOO0OOooo........oooOO0OOooo........oo 90 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 90 91 91 G4BGGNucleonInelasticXS::~G4BGGNucleonInelasti 92 G4BGGNucleonInelasticXS::~G4BGGNucleonInelasticXS() 92 { 93 { 93 delete fHadron; 94 delete fHadron; >> 95 delete fSAID; 94 } 96 } 95 97 96 //....oooOO0OOooo........oooOO0OOooo........oo 98 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 97 99 98 G4bool G4BGGNucleonInelasticXS::IsElementAppli 100 G4bool G4BGGNucleonInelasticXS::IsElementApplicable(const G4DynamicParticle*, 99 << 101 G4int, const G4Material*) 100 { 102 { 101 return true; 103 return true; 102 } 104 } 103 105 104 //....oooOO0OOooo........oooOO0OOooo........oo 106 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 105 107 106 G4bool G4BGGNucleonInelasticXS::IsIsoApplicabl 108 G4bool G4BGGNucleonInelasticXS::IsIsoApplicable(const G4DynamicParticle*, 107 << 109 G4int Z, G4int, 108 << 110 const G4Element*, 109 << 111 const G4Material*) 110 { 112 { 111 return (1 == Z); 113 return (1 == Z); 112 } 114 } 113 115 114 //....oooOO0OOooo........oooOO0OOooo........oo 116 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 115 117 116 G4double 118 G4double 117 G4BGGNucleonInelasticXS::GetElementCrossSectio 119 G4BGGNucleonInelasticXS::GetElementCrossSection(const G4DynamicParticle* dp, 118 << 120 G4int ZZ, const G4Material*) 119 { 121 { >> 122 // this method should be called only for Z > 1 >> 123 120 G4double cross = 0.0; 124 G4double cross = 0.0; 121 G4double ekin = dp->GetKineticEnergy(); 125 G4double ekin = dp->GetKineticEnergy(); 122 G4int Z = std::min(ZZ, 92); << 126 G4int Z = ZZ; 123 if (1 == Z) { << 127 if(1 == Z) { 124 cross = 1.0115*GetIsoCrossSection(dp,1,1); 128 cross = 1.0115*GetIsoCrossSection(dp,1,1); >> 129 } else if(2 == Z) { >> 130 if(ekin > fGlauberEnergy) { >> 131 cross = theGlauberFac[Z]*fGlauber->GetInelasticGlauberGribov(dp, Z, theA[Z]); >> 132 } else { >> 133 cross = fNucleon->GetElementCrossSection(dp, Z); >> 134 } >> 135 125 } else { 136 } else { 126 if (ekin <= fLowEnergy) { << 137 if(Z > 92) { Z = 92; } 127 cross = CoulombFactor(ekin, Z); << 138 128 cross *= (isProton) ? theCoulombFacP[Z] << 139 if(ekin <= fLowEnergy) { 129 } else if (ekin > fGlauberEnergy) { << 140 cross = theCoulombFac[Z]*CoulombFactor(ekin, Z); 130 cross = fGlauber->GetInelasticGlauberGri << 141 } else if(ekin > fGlauberEnergy) { 131 cross *= (isProton) ? theGlauberFacP[Z] << 142 cross = theGlauberFac[Z]*fGlauber->GetInelasticGlauberGribov(dp, Z, theA[Z]); 132 } else { 143 } else { 133 cross = fNucleon->GetElementCrossSection 144 cross = fNucleon->GetElementCrossSection(dp, Z); 134 } 145 } 135 } 146 } 136 147 137 #ifdef G4VERBOSE << 148 if(verboseLevel > 1) { 138 if (verboseLevel > 1) { << 139 G4cout << "G4BGGNucleonInelasticXS::GetCro 149 G4cout << "G4BGGNucleonInelasticXS::GetCrossSection for " 140 << dp->GetDefinition()->GetParticle << 150 << dp->GetDefinition()->GetParticleName() 141 << " Ekin(GeV)= " << dp->GetKineti << 151 << " Ekin(GeV)= " << dp->GetKineticEnergy()/CLHEP::GeV 142 << " in nucleus Z= " << Z << " A= << 152 << " in nucleus Z= " << Z << " A= " << theA[Z] 143 << " XS(b)= " << cross/barn << 153 << " XS(b)= " << cross/barn 144 << G4endl; << 154 << G4endl; 145 } 155 } 146 #endif << 156 //AR-18Dec2013 if(cross <= fLowestXSection) { cross = 0.0; } 147 return cross; 157 return cross; 148 } 158 } 149 159 150 //....oooOO0OOooo........oooOO0OOooo........oo 160 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 151 161 152 G4double 162 G4double 153 G4BGGNucleonInelasticXS::GetIsoCrossSection(co 163 G4BGGNucleonInelasticXS::GetIsoCrossSection(const G4DynamicParticle* dp, 154 G4int, << 164 G4int Z, G4int A, 155 const G << 165 const G4Isotope*, 156 const G << 166 const G4Element*, 157 const G << 167 const G4Material*) 158 { 168 { 159 // this method should be called only for Z = 169 // this method should be called only for Z = 1 160 fHadron->HadronNucleonXscNS(dp->GetDefinitio << 161 dp->GetKineticEnerg << 162 G4double cross = A*fHadron->GetInelasticHadr << 163 170 164 #ifdef G4VERBOSE << 171 G4double cross = 0.0; >> 172 G4double ekin = dp->GetKineticEnergy(); >> 173 >> 174 if(ekin <= fSAIDHighEnergyLimit) { >> 175 cross = fSAID->GetInelasticIsotopeCrossSection(particle, ekin, 1, 1); >> 176 } else if(ekin < fHighEnergy) { >> 177 fHadron->GetHadronNucleonXscNS(dp, theProton); >> 178 cross = (theCoulombFac[0]/ekin + 1)*fHadron->GetInelasticHadronNucleonXsc(); >> 179 } else { >> 180 fHadron->GetHadronNucleonXscPDG(dp, theProton); >> 181 cross = (theCoulombFac[1]/ekin + 1)*fHadron->GetInelasticHadronNucleonXsc(); >> 182 } >> 183 cross *= A; >> 184 165 if(verboseLevel > 1) { 185 if(verboseLevel > 1) { 166 G4cout << "G4BGGNucleonInelasticXS::GetIso << 186 G4cout << "G4BGGNucleonInelasticXS::GetCrossSection for " 167 << dp->GetDefinition()->GetParticle << 187 << dp->GetDefinition()->GetParticleName() 168 << " Ekin(GeV)= " << dp->GetKineti << 188 << " Ekin(GeV)= " << dp->GetKineticEnergy()/CLHEP::GeV 169 << " in nucleus Z=1 A=" << A << 189 << " in nucleus Z= " << Z << " A= " << A 170 << " XS(b)= " << cross/barn << 190 << " XS(b)= " << cross/barn 171 << G4endl; << 191 << G4endl; 172 } 192 } 173 #endif << 193 //AR-18Dec2013 if(cross <= fLowestXSection) { cross = 0.0; } 174 return cross; 194 return cross; 175 } 195 } 176 196 177 //....oooOO0OOooo........oooOO0OOooo........oo 197 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 178 198 179 void G4BGGNucleonInelasticXS::BuildPhysicsTabl 199 void G4BGGNucleonInelasticXS::BuildPhysicsTable(const G4ParticleDefinition& p) 180 { 200 { 181 if(&p == theProton || &p == G4Neutron::Neutr 201 if(&p == theProton || &p == G4Neutron::Neutron()) { 182 isProton = (theProton == &p); << 202 particle = &p; 183 } else { 203 } else { 184 G4ExceptionDescription ed; << 204 G4cout << "### G4BGGNucleonInelasticXS WARNING: is not applicable to " 185 ed << "This BGG cross section is applicabl << 205 << p.GetParticleName() 186 << p.GetParticleName() << G4endl; << 206 << G4endl; 187 G4Exception("G4BGGNucleonInelasticXS::Buil << 207 throw G4HadronicException(__FILE__, __LINE__, 188 FatalException, ed); << 208 "G4BGGNucleonElasticXS::BuildPhysicsTable is used for wrong particle"); 189 return; 209 return; 190 } 210 } 191 } << 192 211 193 //....oooOO0OOooo........oooOO0OOooo........oo << 212 if(isInitialized) { return; } >> 213 isInitialized = true; 194 214 195 void G4BGGNucleonInelasticXS::Initialise() << 215 fNucleon = (G4NucleonNuclearCrossSection*)G4CrossSectionDataSetRegistry::Instance()->GetCrossSectionDataSet(G4NucleonNuclearCrossSection::Default_Name()); 196 { << 216 fGlauber = (G4GlauberGribovCrossSection*)G4CrossSectionDataSetRegistry::Instance()->GetCrossSectionDataSet(G4GlauberGribovCrossSection::Default_Name()); 197 theA[0] = theA[1] = 1; << 217 >> 218 fHadron = new G4HadronNucleonXsc(); >> 219 fSAID = new G4ComponentSAIDTotalXS(); >> 220 >> 221 fNucleon->BuildPhysicsTable(*particle); >> 222 fGlauber->BuildPhysicsTable(*particle); >> 223 >> 224 if(particle == theProton) { >> 225 isProton = true; >> 226 fSAIDHighEnergyLimit = 2*GeV; >> 227 fHighEnergy = 2*GeV; >> 228 } >> 229 >> 230 G4ParticleDefinition* part = const_cast<G4ParticleDefinition*>(particle); 198 G4ThreeVector mom(0.0,0.0,1.0); 231 G4ThreeVector mom(0.0,0.0,1.0); 199 G4DynamicParticle dp(theProton, mom, fGlaube << 232 G4DynamicParticle dp(part, mom, fGlauberEnergy); 200 233 201 G4NistManager* nist = G4NistManager::Instanc 234 G4NistManager* nist = G4NistManager::Instance(); >> 235 G4int A; >> 236 202 G4double csup, csdn; 237 G4double csup, csdn; 203 238 204 for (G4int iz=2; iz<93; ++iz) { << 239 if(verboseLevel > 0) { >> 240 G4cout << "### G4BGGNucleonInelasticXS::Initialise for " >> 241 << particle->GetParticleName() << G4endl; >> 242 } >> 243 for(G4int iz=2; iz<93; iz++) { 205 244 206 G4int A = G4lrint(nist->GetAtomicMassAmu(i << 245 A = G4lrint(nist->GetAtomicMassAmu(iz)); 207 theA[iz] = A; 246 theA[iz] = A; 208 247 209 csup = fGlauber->GetInelasticGlauberGribov 248 csup = fGlauber->GetInelasticGlauberGribov(&dp, iz, A); 210 csdn = fNucleon->GetElementCrossSection(&d 249 csdn = fNucleon->GetElementCrossSection(&dp, iz); 211 theGlauberFacP[iz] = csdn/csup; << 212 } << 213 << 214 dp.SetDefinition(G4Neutron::Neutron()); << 215 for (G4int iz=2; iz<93; ++iz) { << 216 csup = fGlauber->GetInelasticGlauberGribov << 217 csdn = fNucleon->GetElementCrossSection(&d << 218 theGlauberFacN[iz] = csdn/csup; << 219 250 220 if(verboseLevel > 1) { << 251 theGlauberFac[iz] = csdn/csup; 221 G4cout << "G4BGGNucleonInelasticXS::Init << 252 if(verboseLevel > 0) { 222 << " GFactorP= " << theGlauberFacP[iz] << 253 G4cout << "Z= " << iz << " A= " << A 223 << " GFactorN= " << theGlauberFacN[iz] << 254 << " GlauberFactor= " << theGlauberFac[iz] << G4endl; 224 } 255 } 225 } 256 } >> 257 //const G4Material* mat = 0; 226 258 227 theCoulombFacP[1] = theCoulombFacN[1] = 1.0; << 259 dp.SetKineticEnergy(fSAIDHighEnergyLimit); 228 dp.SetDefinition(theProton); << 260 fHadron->GetHadronNucleonXscNS(&dp, theProton); >> 261 theCoulombFac[0] = fSAIDHighEnergyLimit* >> 262 (fSAID->GetInelasticIsotopeCrossSection(particle,fSAIDHighEnergyLimit,1,1) >> 263 /fHadron->GetInelasticHadronNucleonXsc() - 1); >> 264 >> 265 //G4cout << "Z=1 E(GeV)= " << fSAIDHighEnergyLimit/GeV >> 266 // << " xsNS(b)= " << fHadron->GetInelasticHadronNucleonXsc()/barn; >> 267 fHadron->GetHadronNucleonXscPDG(&dp, theProton); >> 268 //G4cout << " xsPDG(b)= " << fHadron->GetInelasticHadronNucleonXsc()/barn; >> 269 //G4cout << " xsSAID(b)= " << fSAID->GetInelasticIsotopeCrossSection(particle,fSAIDHighEnergyLimit,1,1)/barn << G4endl; >> 270 >> 271 dp.SetKineticEnergy(fHighEnergy); >> 272 fHadron->GetHadronNucleonXscPDG(&dp, theProton); >> 273 G4double x = fHadron->GetInelasticHadronNucleonXsc(); >> 274 >> 275 //G4cout << "Z=1 E(GeV)= " << fHighEnergy/GeV >> 276 // << " xsPDG(b)= " << fHadron->GetInelasticHadronNucleonXsc()/barn; >> 277 >> 278 fHadron->GetHadronNucleonXscNS(&dp, theProton); >> 279 theCoulombFac[1] = fHighEnergy*((theCoulombFac[0]/fHighEnergy + 1) >> 280 *fHadron->GetInelasticHadronNucleonXsc()/x - 1); >> 281 >> 282 fHadron->GetHadronNucleonXscNS(&dp, theProton); >> 283 //G4cout <<" xsNS(b)= "<<fHadron->GetInelasticHadronNucleonXsc()/barn<<G4endl; >> 284 >> 285 if(verboseLevel > 0) { >> 286 G4cout << "Z=1 A=1" << " CoulombFactor[0]= " << theCoulombFac[0] >> 287 << " CoulombFactor[1]= " << theCoulombFac[1] << G4endl; >> 288 } >> 289 theCoulombFac[2] = 1.0; >> 290 229 dp.SetKineticEnergy(fLowEnergy); 291 dp.SetKineticEnergy(fLowEnergy); 230 for (G4int iz=2; iz<93; ++iz) { << 292 for(G4int iz=3; iz<93; iz++) { 231 theCoulombFacP[iz] = fNucleon->GetElementC << 293 theCoulombFac[iz] = 232 /CoulombFactor(fLowEnergy, iz); << 294 fNucleon->GetElementCrossSection(&dp, iz)/CoulombFactor(fLowEnergy, iz); 233 } << 295 234 dp.SetDefinition(G4Neutron::Neutron()); << 296 if(verboseLevel > 0) { 235 for (G4int iz=2; iz<93; ++iz) { << 297 G4cout << "Z= " << iz << " A= " << theA[iz] 236 theCoulombFacN[iz] = fNucleon->GetElementC << 298 << " CoulombFactor= " << theCoulombFac[iz] << G4endl; 237 /CoulombFactor(fLowEnergy, iz); << 238 << 239 if (verboseLevel > 1) { << 240 G4cout << "G4BGGNucleonInelasticXS::Init << 241 << " CFactorP= " << theCoulombFacP[iz] << 242 << " CFactorN= " << theCoulombFacN[iz] << 243 } 299 } 244 } 300 } 245 } 301 } 246 302 247 //....oooOO0OOooo........oooOO0OOooo........oo 303 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 248 304 249 G4double G4BGGNucleonInelasticXS::CoulombFacto 305 G4double G4BGGNucleonInelasticXS::CoulombFactor(G4double kinEnergy, G4int Z) 250 { 306 { 251 G4double res = 0.0; << 307 G4double res= 0.0; 252 << 253 if(kinEnergy <= 0.0) { return res; } 308 if(kinEnergy <= 0.0) { return res; } 254 << 309 else if (Z <= 1) { return kinEnergy*kinEnergy; } >> 310 255 G4double elog = G4Log(kinEnergy/GeV)/llog10; 311 G4double elog = G4Log(kinEnergy/GeV)/llog10; 256 G4double aa = theA[Z]; 312 G4double aa = theA[Z]; 257 << 258 if(isProton) { << 259 313 260 res = G4NuclearRadii::CoulombFactor(Z, aa, << 314 // from G4ProtonInelasticCrossSection >> 315 if(isProton) { >> 316 >> 317 G4double ff1 = 5.6 - 0.016*aa; // slope of the drop at medium energies. >> 318 G4double ff2 = 1.37 + 1.37/aa; // start of the slope. >> 319 G4double ff3 = 0.8 + 18./aa - 0.002*aa; // stephight >> 320 res = 1.0 + ff3*(1.0 - (1.0/(1+G4Exp(-ff1*(elog + ff2))))); >> 321 >> 322 ff1 = 8. - 8./aa - 0.008*aa; // slope of the rise >> 323 ff2 = 2.34 - 5.4/aa - 0.0028*aa; // start of the rise >> 324 res /= (1.0 + G4Exp(-ff1*(elog + ff2))); 261 325 262 // from G4ProtonInelasticCrossSection << 263 if(res > 0.0) { << 264 G4double ff1 = 5.6 - 0.016*aa; // slope << 265 G4double ff2 = 1.37 + 1.37/aa; // start << 266 G4double ff3 = 0.8 + 18./aa - 0.002*aa; << 267 res *= (1.0 + ff3*(1.0 - (1.0/(1+G4Exp(- << 268 ff1 = 8. - 8./aa - 0.008*aa; // slope << 269 ff2 = 2.34 - 5.4/aa - 0.0028*aa; // star << 270 res /= (1.0 + G4Exp(-ff1*(elog + ff2))); << 271 } << 272 } else { 326 } else { >> 327 273 // from G4NeutronInelasticCrossSection 328 // from G4NeutronInelasticCrossSection 274 G4double p3 = 0.6 + 13./aa - 0.0005*aa; 329 G4double p3 = 0.6 + 13./aa - 0.0005*aa; 275 G4double p4 = 7.2449 - 0.018242*aa; 330 G4double p4 = 7.2449 - 0.018242*aa; 276 G4double p5 = 1.36 + 1.8/aa + 0.0005*aa; 331 G4double p5 = 1.36 + 1.8/aa + 0.0005*aa; 277 G4double p6 = 1. + 200./aa + 0.02*aa; 332 G4double p6 = 1. + 200./aa + 0.02*aa; 278 G4double p7 = 3.0 - (aa-70.)*(aa-200.)/110 333 G4double p7 = 3.0 - (aa-70.)*(aa-200.)/11000.; 279 334 280 G4double firstexp = G4Exp(-p4*(elog + p5) 335 G4double firstexp = G4Exp(-p4*(elog + p5)); 281 G4double secondexp = G4Exp(-p6*(elog + p7) 336 G4double secondexp = G4Exp(-p6*(elog + p7)); 282 337 283 res = (1. + p3*firstexp/(1. + firstexp))/( << 338 res = (1.+p3*firstexp/(1. + firstexp))/(1. + secondexp); >> 339 284 } 340 } 285 return res; 341 return res; 286 } 342 } 287 343 288 //....oooOO0OOooo........oooOO0OOooo........oo 344 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 289 345 290 void G4BGGNucleonInelasticXS::CrossSectionDesc 346 void G4BGGNucleonInelasticXS::CrossSectionDescription(std::ostream& outFile) const 291 { 347 { 292 outFile << "The Barashenkov-Glauber-Gribov c 348 outFile << "The Barashenkov-Glauber-Gribov cross section calculates inelastic\n" 293 << "scattering of protons and neutro 349 << "scattering of protons and neutrons from nuclei using the\n" 294 << "Barashenkov parameterization bel 350 << "Barashenkov parameterization below 91 GeV and the Glauber-Gribov\n" 295 << "parameterization above 91 GeV. 351 << "parameterization above 91 GeV. It uses the G4HadronNucleonXsc\n" 296 << "cross section component for hydr 352 << "cross section component for hydrogen targets, and the\n" 297 << "G4ComponentGGHadronNucleusXsc co << 353 << "G4GlauberGribovCrossSection component for other targets.\n"; 298 } 354 } 299 355 300 //....oooOO0OOooo........oooOO0OOooo........oo 356 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 301 357