118 G4cout <<
"Find a flag of \"G4NEUTRONHP_NEGLECT_DOPPLER\"." <<
G4endl;
119 G4cout <<
"On the fly Doppler broadening will be neglect in the cross section calculation of fission reaction of neutrons (<20MeV)." <<
G4endl;
124 throw G4HadronicException(__FILE__, __LINE__,
"Attempt to use NeutronHP data for particles other than neutrons!!!");
143 for(
size_t i=0; i<numberOfElements; ++i )
156 throw G4HadronicException(__FILE__, __LINE__,
"Attempt to use NeutronHP data for particles other than neutrons!!!");
167 G4cout <<
"Fission Cross Section of Neutron HP"<<
G4endl;
168 G4cout <<
"(Pointwise cross-section at 0 Kelvin.)" <<
G4endl;
177 for (
size_t i = 0 ; i < numberOfElements ; ++i )
184 G4cout <<
"The cross-section data of the fission of this element is not available." <<
G4endl;
191 for ( ie = 0 ; ie < 130 ; ie++ )
194 G4bool outOfRange =
false;
196 if ( eKinetic < 20*
MeV )
198 G4cout << eKinetic/
eV <<
" " << (*((*theCrossSections)(i))).GetValue(eKinetic, outOfRange)/
barn <<
G4endl;
215 if(anE->
GetZ()<88)
return result;
220 if ( ( ( *
theCrossSections )( index ) )->GetVectorLength() == 0 )
return result;
226 theNeutronRP.SetKineticEnergy( eKinetic );
237 return ( (*((*
theCrossSections)(index))).GetValue(eKinetic, outOfRange) )* factor;
259 while(counter == 0 ||
std::abs(buffer-result/
std::max(1,counter)) > 0.01*buffer)
261 if(counter) buffer = result/counter;
266 boosted.
Lorentz(theNeutronRP, aThermalNuc);
268 aXsection = (*((*theCrossSections)(index))).GetValue(theEkin, outOfRange);
271 aXsection *= (targetVelocity-neutronVelocity).mag()/neutronVMag;
290 outFile <<
"High Precision cross data based on Evaluated Nuclear Data Files (ENDF) for induced fission reaction of neutrons below 20MeV\n" ;