42 minEnergy(19.9*
MeV), maxEnergy(19.9*
GeV)
51 outFile <<
"G4NeutronInelasticCrossSection calculates the inelastic neutron\n"
52 <<
"scattering cross section for nuclei using the Wellisch-Laidlaw\n"
53 <<
"parameterization between 19.9 MeV and 19.9 GeV. Above 19.9 GeV\n"
54 <<
"the cross section is assumed to be constant.\n";
79 if(anEnergy <
keV) {
return cross_section; }
85 G4int nOfNeutrons = A - Z;
88 G4double p2 = 1. + 10./atomicNumber - 0.0006*atomicNumber;
89 G4double p3 = 0.6+ 13./atomicNumber - 0.0005*atomicNumber;
90 G4double p4 = 7.2449 - 0.018242*atomicNumber;
91 G4double p5 = 1.64 - 1.8/atomicNumber - 0.0005*atomicNumber;
92 G4double p6 = 1. + 200./atomicNumber + 0.02*atomicNumber;
93 G4double p7 = (atomicNumber-70.)*(atomicNumber-200.)/11000.;
101 G4double corr = 1. + p3*(1.-1./first);
108 if(xsec < 0.0) { xsec = 0.0; }
static G4Pow * GetInstance()
virtual void CrossSectionDescription(std::ostream &) const
G4double GetKineticEnergy() const
G4double GetCrossSection(G4double kineticEnergy, G4int Z, G4int A)
virtual G4double GetElementCrossSection(const G4DynamicParticle *, G4int Z, const G4Material *)
static G4NistManager * Instance()
~G4NeutronInelasticCrossSection()
G4double logZ(G4int Z) const
G4double Z13(G4int Z) const
double A(double temperature)
G4NeutronInelasticCrossSection()
virtual G4bool IsElementApplicable(const G4DynamicParticle *aPart, G4int Z, const G4Material *)
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
static const double millibarn