125 G4double xmin = 2.0*std::log(2.0*leptonmass);
126 G4double xmax = 2.0*std::log(parentmass);
133 w1 = (1.0-4.0*leptonmass*leptonmass/t);
135 w2 = ( 1.0 + 2.0*leptonmass*leptonmass/t );
136 w3 = ( 1.0 - t/parentmass/parentmass );
138 ww = w3 * w2 * std::sqrt(w1);
148 G4double sintheta = std::sqrt((1.0 - costheta)*(1.0 + costheta));
150 G4ThreeVector gdirection(sintheta*std::cos(phi),sintheta*std::sin(phi),costheta);
157 G4double beta = Pgamma/(parentmass-Pgamma);
162 G4double Elepton = std::sqrt(Plepton*Plepton + leptonmass*leptonmass );
164 sintheta = std::sqrt((1.0 - costheta)*(1.0 + costheta));
166 G4ThreeVector ldirection(sintheta*std::cos(phi),sintheta*std::sin(phi),costheta);
170 ldirection, Elepton-leptonmass );
173 -1.0*ldirection, Elepton-leptonmass );
176 p4.
boost( -1.0*gdirection.x()*beta, -1.0*gdirection.y()*beta, -1.0*gdirection.z()*beta);
179 p4.
boost( -1.0*gdirection.x()*beta, -1.0*gdirection.y()*beta, -1.0*gdirection.z()*beta);
184 delete parentparticle;
191 G4cout <<
"G4DalitzDecayChannel::DecayIt ";
192 G4cout <<
" create decay products in rest frame " <<
G4endl;
G4int PushProducts(G4DynamicParticle *aParticle)
G4ParticleDefinition * G4MT_parent
G4ParticleDefinition ** G4MT_daughters
G4GLOB_DLL std::ostream G4cout
HepLorentzVector & boost(double, double, double)
G4LorentzVector Get4Momentum() const
void Set4Momentum(const G4LorentzVector &momentum)
G4int GetVerboseLevel() const
G4double GetPDGMass() const
static G4double Pmx(G4double e, G4double p1, G4double p2)