53 material_cache = NULL;
62 delete theCrossSections;
84 if ( dp->
GetKineticEnergy() == ke_cache && element == element_cache && material == material_cache )
return xs_cache;
87 element_cache = element;
107 throw G4HadronicException(__FILE__, __LINE__,
"Attempt to use NeutronHP data for particles other than neutrons!!!");
112 if ( theCrossSections == NULL )
120 for(
size_t i=0; i<numberOfElements; ++i )
131 throw G4HadronicException(__FILE__, __LINE__,
"Attempt to use NeutronHP data for particles other than neutrons!!!");
142 G4cout <<
"Fission Cross Section of Neutron HP"<<
G4endl;
143 G4cout <<
"(Pointwise cross-section at 0 Kelvin.)" <<
G4endl;
152 for (
size_t i = 0 ; i < numberOfElements ; ++i )
157 if ( (*((*theCrossSections)(i))).GetVectorLength() == 0 )
159 G4cout <<
"The cross-section data of the fission of this element is not available." <<
G4endl;
166 for ( ie = 0 ; ie < 130 ; ie++ )
168 G4double eKinetic = 1.0e-5 * std::pow ( 10.0 , ie/10.0 ) *
eV;
169 G4bool outOfRange =
false;
171 if ( eKinetic < 20*
MeV )
173 G4cout << eKinetic/
eV <<
" " << (*((*theCrossSections)(i))).GetValue(eKinetic, outOfRange)/
barn <<
G4endl;
190 if(anE->
GetZ()<90)
return result;
195 if ( ( ( *theCrossSections )( index ) )->GetVectorLength() == 0 )
return result;
201 theNeutron.SetKineticEnergy( eKinetic );
222 while(counter == 0 || std::abs(buffer-result/
std::max(1,counter)) > 0.01*buffer)
224 if(counter) buffer = result/counter;
229 boosted.
Lorentz(theNeutron, aThermalNuc);
231 aXsection = (*((*theCrossSections)(
index))).GetValue(theEkin, outOfRange);
234 aXsection *= (targetVelocity-neutronVelocity).mag()/neutronVMag;
G4PhysicsVector * MakePhysicsVector(G4Element *thE, G4NeutronHPFissionData *theP)
static G4double GetNuclearMass(const G4double A, const G4double Z)
void Lorentz(const G4ReactionProduct &p1, const G4ReactionProduct &p2)
G4double GetKineticEnergy() const
G4bool IsIsoApplicable(const G4DynamicParticle *, G4int, G4int, const G4Element *, const G4Material *)
static G4NeutronHPManager * GetInstance()
void SetMomentum(const G4double x, const G4double y, const G4double z)
void push_back(G4PhysicsVector *)
G4double G4NeutronHPJENDLHEData::G4double result
G4ParticleDefinition * GetDefinition() const
G4ReactionProduct GetThermalNucleus(G4double aMass, G4double temp=-1) const
~G4NeutronHPFissionData()
G4GLOB_DLL std::ostream G4cout
G4double GetCrossSection(const G4DynamicParticle *, const G4Element *, G4double aT)
static size_t GetNumberOfElements()
void SetMinKinEnergy(G4double value)
G4int GetVerboseLevel() const
G4double GetIsoCrossSection(const G4DynamicParticle *, G4int, G4int, const G4Isotope *, const G4Element *, const G4Material *)
static G4Neutron * Neutron()
static G4NeutronHPData * Instance()
G4double GetKineticEnergy() const
void BuildPhysicsTable(const G4ParticleDefinition &)
void SetVerboseLevel(G4int i)
void SetVerboseLevel(G4int)
G4double GetPDGMass() const
void SetMaxKinEnergy(G4double value)
T max(const T t1, const T t2)
brief Return the largest of the two arguments
void DumpPhysicsTable(const G4ParticleDefinition &)
G4ThreeVector GetMomentum() const
std::vector< G4Element * > G4ElementTable
static G4ElementTable * GetElementTable()
G4double GetMinKinEnergy() const
G4ThreeVector GetMomentum() const
G4double GetMaxKinEnergy() const