28 #ifndef G4NeutronHPSimpleEvapSpectrum_h
29 #define G4NeutronHPSimpleEvapSpectrum_h 1
32 #include <CLHEP/Units/SystemOfUnits.h>
47 expm1 = std::exp(-1.);
53 inline void Init(std::istream & aDataFile)
72 result = -theta*std::log(random);
75 while(cut>result/max);
void Init(std::istream &aDataFile)
G4double GetY(G4double x)
~G4NeutronHPSimpleEvapSpectrum()
G4double Evapo(G4double anEnergy, G4double theta)
void Init(std::istream &aDataFile, G4int total, G4double ux=1., G4double uy=1.)
T max(const T t1, const T t2)
brief Return the largest of the two arguments
G4NeutronHPSimpleEvapSpectrum()
G4double Sample(G4double anEnergy)
G4NeutronHPVector theThetaDist
G4NeutronHPVector theFractionalProb
G4double GetFractionalProbability(G4double anEnergy)