28 #ifndef G4NeutronHPMadlandNixSpectrum_h
29 #define G4NeutronHPMadlandNixSpectrum_h 1
33 #include <CLHEP/Units/PhysicalConstants.h>
52 expm1 = std::exp(-1.);
58 inline void Init(std::istream & aDataFile)
62 theAvarageKineticPerNucleonForLightFragments*=
CLHEP::eV;
64 theAvarageKineticPerNucleonForHeavyFragments*=
CLHEP::eV;
93 result = 1- (0.3480242*t - 0.0958798*t*t + 0.7478556*t*t*t)*std::exp(-aValue);
102 result = 0.5*
Gamma05(aValue) - std::sqrt(aValue)*std::exp(-aValue);
109 result = 1.5*
Gamma15(aValue) - std::pow(aValue,1.5)*std::exp(aValue);
119 G4double result =-gamma - std::log(aValue);
130 term = -term*aValue*(count-1)/(count*count);
132 if(std::fabs(term)/std::fabs(result)<precision)
break;
G4double theAvarageKineticPerNucleonForHeavyFragments
G4double GetY(G4double x)
void Init(std::istream &aDataFile)
G4NeutronHPVector theFractionalProb
void Init(std::istream &aDataFile, G4int total, G4double ux=1., G4double uy=1.)
G4NeutronHPMadlandNixSpectrum()
G4double Sample(G4double anEnergy)
G4double Madland(G4double aSecEnergy, G4double tm)
G4double GetFractionalProbability(G4double anEnergy)
G4double theAvarageKineticPerNucleonForLightFragments
G4double FissionIntegral(G4double tm, G4double anEnergy)
G4double GIntegral(G4double tm, G4double anEnergy, G4double aMean)
~G4NeutronHPMadlandNixSpectrum()
G4double Gamma15(G4double aValue)
G4double Gamma05(G4double aValue)
G4double E1(G4double aValue)
G4double Gamma25(G4double aValue)
G4NeutronHPVector theMaxTemp