121 throw G4HadronicException(__FILE__, __LINE__,
"Attempt to use NeutronHP data for particles other than neutrons!!!");
126 G4cout <<
"Find a flag of \"G4NEUTRONHP_NEGLECT_DOPPLER\"." <<
G4endl;
127 G4cout <<
"On the fly Doppler broadening will be neglect in the cross section calculation of elastic scattering of neutrons (<20MeV)." <<
G4endl;
147 for(
size_t i=0; i<numberOfElements; ++i )
160 throw G4HadronicException(__FILE__, __LINE__,
"Attempt to use NeutronHP data for particles other than neutrons!!!");
171 G4cout <<
"Elastic Cross Section of Neutron HP"<<
G4endl;
172 G4cout <<
"(Pointwise cross-section at 0 Kelvin.)" <<
G4endl;
181 for (
size_t i = 0 ; i < numberOfElements ; ++i )
188 for ( ie = 0 ; ie < 130 ; ie++ )
191 G4bool outOfRange =
false;
193 if ( eKinetic < 20*
MeV )
195 G4cout << eKinetic/
eV <<
" " << (*((*theCrossSections)(i))).GetValue(eKinetic, outOfRange)/
barn <<
G4endl;
233 return ( (*((*
theCrossSections)(index))).GetValue(eKinetic, outOfRange) )* factor;
238 theNeutron.SetKineticEnergy( eKinetic );
261 while(counter == 0 || std::abs(buffer-result/
std::max(1,counter)) > 0.03*buffer)
263 if(counter) buffer = result/counter;
268 boosted.
Lorentz(theNeutron, aThermalNuc);
270 aXsection = (*((*theCrossSections)(index))).GetValue(theEkin, outOfRange);
273 aXsection *= (targetVelocity-neutronVelocity).mag()/neutronVMag;
301 outFile <<
"High Precision cross data based on Evaluated Nuclear Data Files (ENDF) for elastic reaction of neutrons below 20MeV\n" ;
static constexpr double kelvin
G4double GetMaxKinEnergy() const
T max(const T t1, const T t2)
brief Return the largest of the two arguments
static G4ParticleHPManager * GetInstance()
G4double GetIsoCrossSection(const G4DynamicParticle *, G4int, G4int, const G4Isotope *, const G4Element *, const G4Material *)
void SetMinKinEnergy(G4double value)
static constexpr double MeV
void SetMaxKinEnergy(G4double value)
const G4Material * material_cache
void SetMomentum(const G4double x, const G4double y, const G4double z)
void SetVerboseLevel(G4int i)
G4ReactionProduct GetThermalNucleus(G4double aMass, G4double temp=-1) const
void RegisterElasticCrossSections(G4PhysicsTable *val)
G4bool IsIsoApplicable(const G4DynamicParticle *, G4int, G4int, const G4Element *, const G4Material *)
~G4ParticleHPElasticData()
G4double GetPDGMass() const
G4PhysicsTable * theCrossSections
const G4Element * element_cache
void BuildPhysicsTable(const G4ParticleDefinition &)
virtual void CrossSectionDescription(std::ostream &) const
G4double GetCrossSection(const G4DynamicParticle *, const G4Element *, G4double aT)
G4double GetKineticEnergy() const
static G4Pow * GetInstance()
G4ParticleDefinition * GetDefinition() const
G4double powA(G4double A, G4double y) const
G4ParticleHPElasticData()
std::vector< G4Element * > G4ElementTable
G4PhysicsVector * MakePhysicsVector(G4Element *thE, G4ParticleHPFissionData *theP)
static G4double GetNuclearMass(const G4double A, const G4double Z)
static constexpr double eV
G4ThreeVector GetMomentum() const
G4double GetMinKinEnergy() const
void SetVerboseLevel(G4int)
G4double G4ParticleHPJENDLHEData::G4double result
G4ThreeVector GetMomentum() const
void push_back(G4PhysicsVector *)
static G4Neutron * Neutron()
static constexpr double k_Boltzmann
static constexpr double barn
G4double GetKineticEnergy() const
G4GLOB_DLL std::ostream G4cout
static size_t GetNumberOfElements()
static const G4double eps
G4PhysicsTable * GetElasticCrossSections()
G4int GetVerboseLevel() const
void Lorentz(const G4ReactionProduct &p1, const G4ReactionProduct &p2)
void DumpPhysicsTable(const G4ParticleDefinition &)
static G4ParticleHPData * Instance(G4ParticleDefinition *projectile)
static G4ElementTable * GetElementTable()