37 #ifndef G4GeneralPhaseSpaceDecay_h
38 #define G4GeneralPhaseSpaceDecay_h 1
53 G4int theNumberOfDaughters,
55 const G4String& theDaughterName2 =
"",
56 const G4String& theDaughterName3 =
"");
61 G4int theNumberOfDaughters,
63 const G4String& theDaughterName2 =
"",
64 const G4String& theDaughterName3 =
"");
69 G4int theNumberOfDaughters,
78 G4int theNumberOfDaughters,
126 G4HadronicException(__FILE__, __LINE__,
"G4GeneralPhaseSpaceDecay::Pmx energy in cms > mass1+mass2");
128 G4double ppp = (e+p1+p2)*(e+p1-p2)*(e-p1+p2)*(e-p1-p2)/(4.0*e*e);
129 if (ppp>0)
return std::sqrt(ppp);
G4DecayProducts * ManyBodyDecayIt()
G4double GetParentMass() const
G4GeneralPhaseSpaceDecay(G4int Verbose=1)
G4DecayProducts * ThreeBodyDecayIt()
const G4double * theDaughterMasses
virtual ~G4GeneralPhaseSpaceDecay()
G4DecayProducts * TwoBodyDecayIt()
virtual G4DecayProducts * DecayIt(G4double mass=0.0)
static G4double Pmx(G4double e, G4double p1, G4double p2)
G4DecayProducts * OneBodyDecayIt()
void SetParentMass(const G4double aParentMass)