29 #ifndef G4ParticleHPMadlandNixSpectrum_h
30 #define G4ParticleHPMadlandNixSpectrum_h 1
64 inline void Init(std::istream & aDataFile)
68 theAvarageKineticPerNucleonForLightFragments*=
CLHEP::eV;
70 theAvarageKineticPerNucleonForHeavyFragments*=
CLHEP::eV;
99 result = 1- (0.3480242*t - 0.0958798*t*t + 0.7478556*t*t*
t)*
G4Exp(-aValue);
108 result = 0.5*
Gamma05(aValue) - std::sqrt(aValue)*
G4Exp(-aValue);
136 term = -term*aValue*(count-1)/(count*count);
138 if(std::fabs(term)/std::fabs(result)<precision)
break;