113 G4LorentzVector w0(0.,0.,0.,0.), w1(0.,0.,0.,0.), w2(0.,0.,0.,0.);
115 G4double e0, p0, e2,
p, gg, m01, m02, m12;
126 p0 = sqrt(e0*e0 - massPi0*massPi0);
161 gg = (px*px + py*py + pz*pz)*
166 G4cout <<
"G4eeTo3PiModel::SampleSecondaries WARNING matrix element g= "
167 << gg <<
" > " <<
gmax <<
" (majoranta)" <<
G4endl;
173 w0.rotateUz(direction);
174 w1.rotateUz(direction);
184 newp->push_back(dp0);
185 newp->push_back(dp1);
186 newp->push_back(dp2);
static G4PionMinus * PionMinus()
G4eeTo3PiModel(G4eeCrossSections *, G4double, G4double)
static constexpr double MeV
G4double CrossSection3pi(G4double)
Float_t x1[n_points_granero]
G4eeCrossSections * cross
static G4PionPlus * PionPlus()
G4ThreeVector G4RandomDirection()
G4double GetPDGMass() const
virtual G4double ComputeCrossSection(G4double) const override
HepLorentzVector & rotateUz(const Hep3Vector &)
virtual G4double PeakEnergy() const override
void set(double x, double y, double z, double t)
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, G4double, const G4ThreeVector &) override
G4double HighEnergy() const
G4GLOB_DLL std::ostream G4cout
static G4PionZero * PionZero()
Hep3Vector boostVector() const
CLHEP::HepLorentzVector G4LorentzVector
virtual ~G4eeTo3PiModel()
static constexpr double GeV
HepLorentzVector & boost(double, double, double)
std::complex< G4double > DpRho(G4double e)