29 #ifndef G4NeutronHPNBodyPhaseSpace_h
30 #define G4NeutronHPNBodyPhaseSpace_h 1
56 void Init(std::ifstream & aDataFile)
58 aDataFile >> theTotalMass >> theTotalCount;
69 result = std::sqrt(anEnergy)*std::pow(eMax-anEnergy, 3.*n/2.-4.);
76 if(theTotalCount==3) result = 4./CLHEP::pi/std::pow(GetEmax(anEnergy, mass),2);
77 if(theTotalCount==4) result = 105./32./std::pow(GetEmax(anEnergy, mass), 3.5);
78 if(theTotalCount==5) result = 256./14./CLHEP::pi/std::pow(GetEmax(anEnergy, mass), 5.);
88 result = availableEnergy*(theTotalMass-mass)/theTotalMass;
92 G4double MeanEnergyOfThisInteraction() {
return -1; }