G4eeTo3PiModel Class Reference

#include <G4eeTo3PiModel.hh>

Inheritance diagram for G4eeTo3PiModel:

Collaboration diagram for G4eeTo3PiModel:

Public Member Functions
	G4eeTo3PiModel (G4eeCrossSections *, G4double, G4double)
virtual	~G4eeTo3PiModel ()
virtual G4double	PeakEnergy () const override
virtual G4double	ComputeCrossSection (G4double) const override
virtual void	SampleSecondaries (std::vector< G4DynamicParticle * > *, G4double, const G4ThreeVector &) override
Public Member Functions inherited from G4Vee2hadrons
	G4Vee2hadrons (G4eeCrossSections *cr, G4double vlowEnergy, G4double vhighEnergy, G4double vdelta)
virtual	~G4Vee2hadrons ()
G4PhysicsVector *	PhysicsVector () const
G4double	LowEnergy () const
G4double	HighEnergy () const

Additional Inherited Members
Protected Attributes inherited from G4Vee2hadrons
G4eeCrossSections *	cross

Detailed Description

Definition at line 59 of file G4eeTo3PiModel.hh.

Constructor & Destructor Documentation

G4eeTo3PiModel::G4eeTo3PiModel ( G4eeCrossSections *  cr, G4double  maxkinEnergy, G4double  binWidth  )

explicit

Definition at line 66 of file G4eeTo3PiModel.cc.

:  G4Vee2hadrons(cr,
             0.41612*GeV,    //threshold
         maxkinEnergy,
         binWidth)
{
  G4cout << "####G4eeTo3PiModel####" << G4endl;

  massPi  = G4PionPlus::PionPlus()->GetPDGMass();
  massPi0 = G4PionZero::PionZero()->GetPDGMass();
  massOm  = 782.62*MeV;
  massPhi = 1019.46*MeV;
  gmax    = 3.0e-8;
}

Here is the call graph for this function:

G4eeTo3PiModel::~G4eeTo3PiModel ( )

virtual

Definition at line 85 of file G4eeTo3PiModel.cc.

{}

Member Function Documentation

G4double G4eeTo3PiModel::ComputeCrossSection ( G4double  e ) const

overridevirtual

Implements G4Vee2hadrons.

Definition at line 99 of file G4eeTo3PiModel.cc.

{
  return cross->CrossSection3pi(e);
}

Here is the call graph for this function:

G4double G4eeTo3PiModel::PeakEnergy ( ) const

overridevirtual

Implements G4Vee2hadrons.

Definition at line 90 of file G4eeTo3PiModel.cc.

{
  G4double e = massOm;
  if(HighEnergy() > massPhi) { e = massPhi; } 
  return e;
}

Here is the call graph for this function:

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

overridevirtual

Implements G4Vee2hadrons.

Definition at line 106 of file G4eeTo3PiModel.cc.

{

  G4double x0 = massPi0/e;
  G4double x1 = massPi/e;

  G4LorentzVector w0(0.,0.,0.,0.), w1(0.,0.,0.,0.), w2(0.,0.,0.,0.);
  G4ThreeVector dir0, dir1;
  G4double e0, p0, e2, p, gg, m01, m02, m12;

  // max pi0 energy
  G4double edel  = 0.5*e*(1.0 + x0*x0 - 4.0*x1*x1) - massPi0;

  const G4int nmax = 200;
  G4int nn = 0;
  do {
    ++nn;
    // pi0 sample
    e0 = edel*G4UniformRand() + massPi0;
    p0 = sqrt(e0*e0 - massPi0*massPi0);
    dir0 = G4RandomDirection();
    w0 = G4LorentzVector(p0*dir0.x(),p0*dir0.y(),p0*dir0.z(),e0);

    // pi+pi- pair
    w1 = G4LorentzVector(-p0*dir0.x(),-p0*dir0.y(),-p0*dir0.z(),e-e0);
    G4ThreeVector bst = w1.boostVector();
    e2 = 0.25*w1.m2();

    // pi+ 
    p = sqrt(e2 - massPi*massPi);
    dir1 = G4RandomDirection();
    w2 = G4LorentzVector(p*dir1.x(),p*dir1.y(),p*dir1.z(),sqrt(e2));
    // pi- 
    w1.set(-w2.px(), -w2.py(), -w2.pz(), w2.e());

    w1.boost(bst);
    w2.boost(bst);

    G4double px2 = w2.x();
    G4double py2 = w2.y();
    G4double pz2 = w2.z();

    G4double px1 = w1.x();
    G4double py1 = w1.y();
    G4double pz1 = w1.z();

    m01 = w0*w1;
    m02 = w0*w2;
    m12 = w1*w2;

    G4double px = py1*pz2 - py2*pz1;
    G4double py = pz1*px2 - pz2*px1;
    G4double pz = px1*py2 - px2*py1;

    gg = (px*px + py*py + pz*pz)*
      norm( 1.0/cross->DpRho(m01) +  1.0/cross->DpRho(m02)
        + 1.0/cross->DpRho(m12) );

    if(gg > gmax) {
      G4cout << "G4eeTo3PiModel::SampleSecondaries WARNING matrix element g= "
         << gg << " > " << gmax << " (majoranta)" << G4endl;
      gmax = gg;
    }
    // Loop checking, 07-Aug-2015, Vladimir Ivanchenko
  } while( gmax*G4UniformRand() > gg || nn < nmax);

  w0.rotateUz(direction);
  w1.rotateUz(direction);
  w2.rotateUz(direction);

  // create G4DynamicParticle objects 
  G4DynamicParticle* dp0 = 
     new G4DynamicParticle(G4PionZero::PionZero(), w0);
  G4DynamicParticle* dp1 = 
     new G4DynamicParticle(G4PionPlus::PionPlus(), w1);
  G4DynamicParticle* dp2 = 
     new G4DynamicParticle(G4PionMinus::PionMinus(), w2);
  newp->push_back(dp0);
  newp->push_back(dp1);
  newp->push_back(dp2);
}

Here is the call graph for this function:

---

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

• source/geant4.10.03.p03/source/processes/electromagnetic/highenergy/include/G4eeTo3PiModel.hh
• source/geant4.10.03.p03/source/processes/electromagnetic/highenergy/src/G4eeTo3PiModel.cc