62 instanceOfWorker =
false;
64 instanceOfWorker =
true;
66 element_cache =
nullptr;
67 material_cache =
nullptr;
75 if ( theCrossSections !=
nullptr && instanceOfWorker !=
true ) {
77 delete theCrossSections;
78 theCrossSections =
nullptr;
102 if ( dp->
GetKineticEnergy() == ke_cache && element == element_cache && material == material_cache )
return xs_cache;
105 element_cache = element;
106 material_cache = material;
126 throw G4HadronicException(__FILE__, __LINE__,
"Attempt to use NeutronHP data for particles other than neutrons!!!");
136 if ( theCrossSections ==
nullptr )
144 for( std::size_t i=0; i<numberOfElements; ++i )
157 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;
182 for ( std::size_t i = 0 ; i < numberOfElements ; ++i )
187 for ( ie = 0 ; ie < 130 ; ++ie )
190 G4bool outOfRange =
false;
192 if ( eKinetic < 20*MeV )
194 G4cout << eKinetic/eV <<
" " << (*((*theCrossSections)(i))).GetValue(eKinetic, outOfRange)/barn <<
G4endl;
216 if ( eKinetic < aT * k_Boltzmann )
223 return ( (*((*theCrossSections)(index))).GetValue(eKinetic, outOfRange) )* factor;
247 G4int size =
G4int(std::max(10., aT/60*kelvin));
251 while(counter == 0 || std::abs(buffer-result/std::max(1,counter)) > 0.03*buffer)
253 if(counter) buffer = result/counter;
258 boosted.
Lorentz(theNeutron, aThermalNuc);
260 aXsection = (*((*theCrossSections)(index))).GetValue(theEkin, outOfRange);
263 aXsection *= (targetVelocity-neutronVelocity).mag()/neutronVMag;
292 outFile <<
"High Precision cross data based on Evaluated Nuclear Data Files (ENDF) for elastic reaction of neutrons below 20MeV\n" ;
std::vector< G4Element * > G4ElementTable
G4GLOB_DLL std::ostream G4cout
G4ParticleDefinition * GetDefinition() const
G4double GetKineticEnergy() const
G4ThreeVector GetMomentum() const
static G4ElementTable * GetElementTable()
static size_t GetNumberOfElements()
static G4HadronicParameters * Instance()
G4double GetTemperature() const
static G4Neutron * Neutron()
static G4double GetNuclearMass(const G4double A, const G4double Z)
G4ReactionProduct GetThermalNucleus(G4double aMass, G4double temp=-1) const
G4double GetPDGMass() const
static G4ParticleHPData * Instance(G4ParticleDefinition *projectile)
G4PhysicsVector * MakePhysicsVector(G4Element *thE, G4ParticleHPFissionData *theP)
G4ParticleHPElasticData()
G4bool IsIsoApplicable(const G4DynamicParticle *, G4int, G4int, const G4Element *, const G4Material *)
~G4ParticleHPElasticData()
void SetVerboseLevel(G4int)
void BuildPhysicsTable(const G4ParticleDefinition &)
G4double GetIsoCrossSection(const G4DynamicParticle *, G4int, G4int, const G4Isotope *, const G4Element *, const G4Material *)
G4int GetVerboseLevel() const
void DumpPhysicsTable(const G4ParticleDefinition &)
virtual void CrossSectionDescription(std::ostream &) const
G4double GetCrossSection(const G4DynamicParticle *, const G4Element *, G4double aT)
G4bool GetNeglectDoppler()
G4PhysicsTable * GetElasticCrossSections()
void SetVerboseLevel(G4int i)
void RegisterElasticCrossSections(G4PhysicsTable *val)
static G4ParticleHPManager * GetInstance()
void push_back(G4PhysicsVector *)
static G4Pow * GetInstance()
G4double powA(G4double A, G4double y) const
void SetMomentum(const G4double x, const G4double y, const G4double z)
G4double GetKineticEnergy() const
G4ThreeVector GetMomentum() const
void Lorentz(const G4ReactionProduct &p1, const G4ReactionProduct &p2)
void SetKineticEnergy(const G4double en)
G4double GetMaxKinEnergy() const
void SetMaxKinEnergy(G4double value)
void SetMinKinEnergy(G4double value)
G4double GetMinKinEnergy() const