39 #ifndef G4FermiPhaseSpaceDecay_hh
40 #define G4FermiPhaseSpaceDecay_hh 1
58 inline std::vector<G4LorentzVector*>*
60 const std::vector<G4double>& fragment_masses)
const;
68 std::vector<G4LorentzVector*> *
86 if (P>0.0) { res = std::sqrt(P); }
91 inline std::vector<G4LorentzVector*>*
93 const std::vector<G4double>& fragment_masses)
const
109 chi = rndmEngine->
flat();
112 }
while ( Fmax*rndmEngine->
flat() > F);
~G4FermiPhaseSpaceDecay()
static const G4double * P2[nN]
std::vector< G4LorentzVector * > * Decay(G4double parent_mass, const std::vector< G4double > &fragment_masses) const
void DumpProblem(G4double E, G4double P1, G4double P2, G4double P) const
G4double powN(G4double x, G4int n) const
G4double BetaKopylov(G4int, CLHEP::HepRandomEngine *) const
const G4FermiPhaseSpaceDecay & operator=(const G4FermiPhaseSpaceDecay &)
G4bool operator!=(const G4FermiPhaseSpaceDecay &)
G4double PtwoBody(G4double E, G4double P1, G4double P2) const
static const G4double * P1[nN]
std::vector< G4LorentzVector * > * KopylovNBodyDecay(G4double, const std::vector< G4double > &) const
G4bool operator==(const G4FermiPhaseSpaceDecay &)