G4eeToTwoPiModel Class Reference

#include <G4eeToTwoPiModel.hh>

Inheritance diagram for G4eeToTwoPiModel:

Public Member Functions
	G4eeToTwoPiModel (G4eeCrossSections *)
virtual	~G4eeToTwoPiModel ()
virtual G4double	ThresholdEnergy () const
virtual G4double	PeakEnergy () const
virtual G4double	ComputeCrossSection (G4double) const
virtual G4PhysicsVector *	PhysicsVector (G4double, G4double) const
virtual void	SampleSecondaries (std::vector< G4DynamicParticle * > *, G4double, const G4ThreeVector &)
Public Member Functions inherited from G4Vee2hadrons
	G4Vee2hadrons ()
virtual	~G4Vee2hadrons ()
void	SetLowEnergy (G4double val)
G4double	LowEnergy () const
void	SetHighEnergy (G4double val)
G4double	HighEnergy () const

Detailed Description

Definition at line 59 of file G4eeToTwoPiModel.hh.

Constructor & Destructor Documentation

G4eeToTwoPiModel::G4eeToTwoPiModel

(

G4eeCrossSections *

cr

)

Definition at line 64 of file G4eeToTwoPiModel.cc.

References G4ParticleDefinition::GetPDGMass(), python.hepunit::MeV, and G4PionPlus::PionPlus().

                                                       :
  cross(cr)
{
  massPi = G4PionPlus::PionPlus()->GetPDGMass();
  massRho = 775.5*MeV;
}

G4eeToTwoPiModel::~G4eeToTwoPiModel ( )

virtual

Definition at line 73 of file G4eeToTwoPiModel.cc.

{}

Member Function Documentation

G4double G4eeToTwoPiModel::ComputeCrossSection ( G4double  e ) const

virtual

Implements G4Vee2hadrons.

Definition at line 92 of file G4eeToTwoPiModel.cc.

References G4eeCrossSections::CrossSection2pi(), G4Vee2hadrons::HighEnergy(), and G4INCL::Math::min().

{
  G4double ee = std::min(HighEnergy(),e);
  return cross->CrossSection2pi(ee);
}

G4double G4eeToTwoPiModel::PeakEnergy ( ) const

virtual

Implements G4Vee2hadrons.

Definition at line 85 of file G4eeToTwoPiModel.cc.

{
  return massRho;
}

G4PhysicsVector * G4eeToTwoPiModel::PhysicsVector ( G4double  emin, G4double  emax  ) const

virtual

Implements G4Vee2hadrons.

Definition at line 100 of file G4eeToTwoPiModel.cc.

References G4INCL::Math::max(), python.hepunit::MeV, G4PhysicsVector::SetSpline(), and test::v.

{
  G4double tmin = std::max(emin, 2.0*massPi);
  G4double tmax = std::max(tmin, emax);
  G4int nbins = (G4int)((tmax - tmin)/(5.*MeV));
  G4PhysicsVector* v = new G4PhysicsLinearVector(emin,emax,nbins);
  v->SetSpline(true);
  return v;
}

void G4eeToTwoPiModel::SampleSecondaries ( std::vector< G4DynamicParticle * > *  newp, G4double  e, const G4ThreeVector &  direction  )

virtual

Implements G4Vee2hadrons.

Definition at line 113 of file G4eeToTwoPiModel.cc.

References G4UniformRand, G4PionMinus::PionMinus(), G4PionPlus::PionPlus(), G4InuclParticleNames::pip, CLHEP::Hep3Vector::rotateUz(), and python.hepunit::twopi.

{

  G4double tkin = 0.5*e - massPi;
  if(tkin < 0.0) tkin = 0.0;
  G4double cost;
  do {
    cost = 2.0*G4UniformRand() - 1.0;
  } while( G4UniformRand() > 1.0 - cost*cost );

  G4double sint = sqrt(1.0 - cost*cost);
  G4double phi  = twopi * G4UniformRand();

  G4ThreeVector dir(sint*cos(phi),sint*sin(phi), cost);
  dir.rotateUz(direction);

  // create G4DynamicParticle objects
  G4DynamicParticle* pip = 
     new G4DynamicParticle(G4PionPlus::PionPlus(),dir,tkin);
  G4DynamicParticle* pin = 
     new G4DynamicParticle(G4PionMinus::PionMinus(),-dir,tkin);
  newp->push_back(pip);
  newp->push_back(pin);
}

G4double G4eeToTwoPiModel::ThresholdEnergy ( ) const

virtual

Implements G4Vee2hadrons.

Definition at line 78 of file G4eeToTwoPiModel.cc.

{
  return 2.0*massPi;
}

---

The documentation for this class was generated from the following files:

• G4eeToTwoPiModel.hh
• G4eeToTwoPiModel.cc