28 #ifndef G4NeutronHPWattSpectrum_h
29 #define G4NeutronHPWattSpectrum_h 1
32 #include <CLHEP/Units/SystemOfUnits.h>
47 expm1 = std::exp(-1.);
53 inline void Init(std::istream & aDataFile)
72 G4double result = std::exp(-energy/a)*std::sinh(std::sqrt(b*energy));
G4double GetY(G4double x)
G4NeutronHPVector theFractionalProb
~G4NeutronHPWattSpectrum()
void Init(std::istream &aDataFile)
void Init(std::istream &aDataFile, G4int total, G4double ux=1., G4double uy=1.)
G4double Watt(G4double anEnergy, G4double a, G4double b)
G4NeutronHPVector theBpar
G4double energy(const ThreeVector &p, const G4double m)
G4NeutronHPWattSpectrum()
G4double Sample(G4double anEnergy)
G4NeutronHPVector theApar
G4double GetFractionalProbability(G4double anEnergy)