104 if(tkin < 0.0) tkin = 0.0;
111 G4double sint = sqrt(1.0 - cost*cost);
115 dir.rotateUz(direction);
122 newp->push_back(pip);
123 newp->push_back(pin);
CLHEP::Hep3Vector G4ThreeVector
G4double CrossSection2pi(G4double)
G4GLOB_DLL std::ostream G4cout
static const double twopi
static G4PionPlus * PionPlus()
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, G4double, const G4ThreeVector &)
G4double GetPDGMass() const
static G4PionMinus * PionMinus()
virtual ~G4eeToTwoPiModel()
G4eeCrossSections * cross
G4eeToTwoPiModel(G4eeCrossSections *, G4double, G4double)
virtual G4double ComputeCrossSection(G4double) const
virtual G4double PeakEnergy() const