48 using namespace CLHEP;
107 if( pName ==
"neutron" &&
119 G4double Tkin, rosxsc, xsc, delta, err=1.e-5;
133 delta = 0.5*std::abs( (rosxsc-xsc) )/(rosxsc+xsc);
194 znq2 = 1. + 2.*
fee*x/
fM;
202 result /= ( x +
fAm )*znq2*znq2*znf4;
204 result *= ( 1 -
x )/( 1 + q2/4./
fM2 ) + 2.*
x;
214 1.52681, 1.54903, 1.57123, 1.59341, 1.61556, 1.63769, 1.6598, 1.68189, 1.70396,
215 1.72601, 1.74805, 1.77007, 1.79208, 1.81407, 1.83605, 1.85801, 1.87997, 1.90192,
216 1.92385, 1.94578, 1.96771, 1.98962, 2.01154, 2.03345, 2.05535, 2.07725, 2.09915,
217 2.12105, 2.14295, 2.16485, 2.18675, 2.20865, 2.23055, 2.25244, 2.27433, 2.29621,
218 2.31807, 2.33992, 2.36173, 2.38351, 2.40524, 2.42691, 2.4485, 2.47, 2.49138,
219 2.51262, 2.53369, 2.55457, 2.57524, 2.59565, 2.61577, 2.63559, 2.65505, 2.67414,
220 2.69281, 2.71104, 2.72881, 2.74607, 2.76282, 2.77903, 2.79467, 2.80974, 2.82422,
221 2.83811, 2.85139, 2.86408, 2.87616, 2.88764, 2.89854, 2.90885, 2.91859, 2.92777,
222 2.93641, 2.94453, 2.95213, 2.95924, 2.96588, 2.97207, 2.97782, 2.98316, 2.98811,
223 2.99268, 2.9969, 3.00078, 3.00435, 3.00761, 3.01059, 3.01331, 3.01578, 3.01801,
224 3.02003, 3.02185, 3.02347, 3.02491, 3.02619, 3.02732, 3.0283, 3.02915, 3.02988,
225 3.03049, 3.03099, 3.03139, 3.03169, 3.03191, 3.03203, 3.03208, 3.03205, 3.03195,
226 3.03177, 3.03152, 3.0312, 3.03081, 3.03034, 3.0298, 3.02919, 3.02849, 3.02771,
227 3.02684, 3.02588, 3.02482, 3.02365, 3.02237, 3.02097, 3.01943, 3.01775, 3.0159,
228 3.01389, 3.01169, 3.00929, 3.00666, 3.00379, 3.00065, 2.99722, 2.99347, 2.98936,
229 2.98487, 2.97996, 2.97459, 2.9687, 2.96226, 2.9552, 2.94748, 2.93903, 2.92977,
230 2.91965, 2.90858, 2.89649, 2.88329, 2.86889, 2.85321, 2.83615, 2.81764, 2.7976,
231 2.77594, 2.7526, 2.72754, 2.70071, 2.67209, 2.64171, 2.60957, 2.57575, 2.54031,
232 2.50336, 2.46504, 2.42548, 2.38484, 2.34328, 2.30099, 2.2581, 2.21478, 2.17115,
233 2.12735, 2.08345, 2.03954, 1.99569, 1.95191, 1.90825, 1.86471, 1.82129, 1.77799,
234 1.7348, 1.69171, 1.64869, 1.60575, 1.56286, 1.52, 1.47718, 1.43437, 1.39157,
235 1.34877, 1.30596, 1.26314, 1.22031, 1.17746, 1.13459, 1.0917, 1.04879, 1.00585,
236 0.962892, 0.919908 };
CLHEP::Hep3Vector G4ThreeVector
static constexpr double MeV
G4double GetLowEdgeEnergy(size_t binNumber) const
G4double CalculateAm(G4double momentum)
static constexpr double hbarc
const G4String & GetParticleName() const
G4double Value(G4double theEnergy, size_t &lastidx) const
G4double XscIntegrand(G4double x)
virtual G4double GetElementCrossSection(const G4DynamicParticle *, G4int Z, const G4Material *)
static constexpr double neutron_mass_c2
G4ParticleDefinition * GetDefinition() const
~G4NeutronElectronElXsc()
static const G4double fXscArray[200]
static constexpr double electron_mass_c2
G4double GetRosenbluthXsc(const G4DynamicParticle *, G4int Z, const G4Material *)
G4Material * FindOrBuildMaterial(const G4String &name, G4bool isotopes=true, G4bool warning=false)
G4PhysicsLogVector * fEnergyXscVector
G4double G4ParticleHPJENDLHEData::G4double result
static constexpr double microbarn
virtual G4bool IsElementApplicable(const G4DynamicParticle *, G4int Z, const G4Material *)
static G4Neutron * Neutron()
G4double GetKineticEnergy() const
G4GLOB_DLL std::ostream G4cout
static constexpr double pi
G4double GetTotalEnergy() const
static constexpr double fine_structure_const
static constexpr double GeV
void PutValue(size_t index, G4double theValue)
static G4NistManager * Instance()