69 fGlauberEnergy = 91.*
GeV;
72 fSAIDHighEnergyLimit = 1.3*
GeV;
74 for (
G4int i = 0; i < 93; ++i) {
75 theGlauberFac[i] = 0.0;
76 theCoulombFac[i] = 0.0;
87 isInitialized =
false;
130 if(ekin > fGlauberEnergy) {
137 if(Z > 92) { Z = 92; }
139 if(ekin <= fLowEnergy) {
140 cross = theCoulombFac[Z]*CoulombFactor(ekin, Z);
141 }
else if(ekin > fGlauberEnergy) {
149 G4cout <<
"G4BGGNucleonInelasticXS::GetCrossSection for "
152 <<
" in nucleus Z= " << Z <<
" A= " << theA[Z]
153 <<
" XS(b)= " << cross/
barn
174 if(ekin <= fSAIDHighEnergyLimit) {
176 }
else if(ekin < fHighEnergy) {
186 G4cout <<
"G4BGGNucleonInelasticXS::GetCrossSection for "
189 <<
" in nucleus Z= " << Z <<
" A= " << A
190 <<
" XS(b)= " << cross/
barn
204 G4cout <<
"### G4BGGNucleonInelasticXS WARNING: is not applicable to "
208 "G4BGGNucleonElasticXS::BuildPhysicsTable is used for wrong particle");
212 if(isInitialized) {
return; }
213 isInitialized =
true;
224 if(particle == theProton) {
226 fSAIDHighEnergyLimit = 2*
GeV;
240 G4cout <<
"### G4BGGNucleonInelasticXS::Initialise for "
251 theGlauberFac[
iz] = csdn/csup;
253 G4cout <<
"Z= " <<
iz <<
" A= " << A
254 <<
" GlauberFactor= " << theGlauberFac[
iz] <<
G4endl;
261 theCoulombFac[0] = fSAIDHighEnergyLimit*
279 theCoulombFac[1] = fHighEnergy*((theCoulombFac[0]/fHighEnergy + 1)
286 G4cout <<
"Z=1 A=1" <<
" CoulombFactor[0]= " << theCoulombFac[0]
287 <<
" CoulombFactor[1]= " << theCoulombFac[1] <<
G4endl;
289 theCoulombFac[2] = 1.0;
298 <<
" CoulombFactor= " << theCoulombFac[
iz] <<
G4endl;
308 if(kinEnergy <= 0.0) {
return res; }
309 else if (Z <= 1) {
return kinEnergy*kinEnergy; }
319 G4double ff3 = 0.8 + 18./aa - 0.002*aa;
320 res = 1.0 + ff3*(1.0 - (1.0/(1+
G4Exp(-ff1*(elog + ff2)))));
322 ff1 = 8. - 8./aa - 0.008*aa;
323 ff2 = 2.34 - 5.4/aa - 0.0028*aa;
324 res /= (1.0 +
G4Exp(-ff1*(elog + ff2)));
329 G4double p3 = 0.6 + 13./aa - 0.0005*aa;
331 G4double p5 = 1.36 + 1.8/aa + 0.0005*aa;
332 G4double p6 = 1. + 200./aa + 0.02*aa;
333 G4double p7 = 3.0 - (aa-70.)*(aa-200.)/11000.;
338 res = (1.+p3*firstexp/(1. + firstexp))/(1. + secondexp);
348 outFile <<
"The Barashenkov-Glauber-Gribov cross section calculates inelastic\n"
349 <<
"scattering of protons and neutrons from nuclei using the\n"
350 <<
"Barashenkov parameterization below 91 GeV and the Glauber-Gribov\n"
351 <<
"parameterization above 91 GeV. It uses the G4HadronNucleonXsc\n"
352 <<
"cross section component for hydrogen targets, and the\n"
353 <<
"G4GlauberGribovCrossSection component for other targets.\n";
G4VCrossSectionDataSet * GetCrossSectionDataSet(const G4String &name, G4bool warning=true)
virtual G4double GetInelasticIsotopeCrossSection(const G4ParticleDefinition *, G4double kinEnergy, G4int, G4int)
G4double GetKineticEnergy() const
virtual void BuildPhysicsTable(const G4ParticleDefinition &)
virtual G4double GetElementCrossSection(const G4DynamicParticle *aParticle, G4int Z, const G4Material *mat=0)
G4ParticleDefinition * GetDefinition() const
static G4NistManager * Instance()
const G4String & GetParticleName() const
G4double GetHadronNucleonXscNS(const G4DynamicParticle *, const G4ParticleDefinition *)
G4BGGNucleonInelasticXS(const G4ParticleDefinition *)
G4GLOB_DLL std::ostream G4cout
virtual void CrossSectionDescription(std::ostream &) const
static const char * Default_Name()
static G4CrossSectionDataSetRegistry * Instance()
static G4Proton * Proton()
static G4Neutron * Neutron()
void SetKineticEnergy(G4double aEnergy)
virtual G4double GetElementCrossSection(const G4DynamicParticle *, G4int Z, const G4Material *mat=0)
G4double G4Log(G4double x)
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
G4double GetInelasticGlauberGribov(const G4DynamicParticle *, G4int Z, G4int A)
virtual G4bool IsElementApplicable(const G4DynamicParticle *, G4int Z, const G4Material *mat=0)
virtual G4double GetIsoCrossSection(const G4DynamicParticle *, G4int Z, G4int A, const G4Isotope *iso=0, const G4Element *elm=0, const G4Material *mat=0)
G4double GetAtomicMassAmu(const G4String &symb) const
virtual void BuildPhysicsTable(const G4ParticleDefinition &)
G4double GetHadronNucleonXscPDG(const G4DynamicParticle *, const G4ParticleDefinition *)
virtual ~G4BGGNucleonInelasticXS()
static const char * Default_Name()
virtual G4bool IsIsoApplicable(const G4DynamicParticle *, G4int Z, G4int A, const G4Element *elm=0, const G4Material *mat=0)
G4double GetInelasticHadronNucleonXsc()