G4eeToPGammaModel Class Reference

#include <G4eeToPGammaModel.hh>

Inheritance diagram for G4eeToPGammaModel:

Collaboration diagram for G4eeToPGammaModel:

Public Member Functions
	G4eeToPGammaModel (G4eeCrossSections *, const G4String &, G4double, G4double)
virtual	~G4eeToPGammaModel ()
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 60 of file G4eeToPGammaModel.hh.

Constructor & Destructor Documentation

G4eeToPGammaModel::G4eeToPGammaModel ( G4eeCrossSections *  cr, const G4String &  npart, G4double  maxkinEnergy, G4double  binWidth  )

explicit

Definition at line 65 of file G4eeToPGammaModel.cc.

: G4Vee2hadrons(cr,
        npart=="pi0" ? 782.62*MeV:1019.46*MeV,  
        maxkinEnergy,
        binWidth)
{
  G4cout << "####G4eeToPGammaModel & particle:" << npart 
     << "####" << G4endl; 

  pi0 = G4PionZero::PionZero();
  if(npart == "pi0") {
    massR = 782.62*MeV;
    particle = pi0;
  } else {
    massR = 1019.46*MeV;
    particle = G4Eta::Eta();
  }
  massP = particle->GetPDGMass();

}

Here is the call graph for this function:

G4eeToPGammaModel::~G4eeToPGammaModel ( )

virtual

Definition at line 91 of file G4eeToPGammaModel.cc.

{}

Member Function Documentation

G4double G4eeToPGammaModel::ComputeCrossSection ( G4double  e ) const

overridevirtual

Implements G4Vee2hadrons.

Definition at line 103 of file G4eeToPGammaModel.cc.

{
  G4double xs;
  if(particle == pi0) xs = cross->CrossSectionPi0G(e);
  else                xs = cross->CrossSectionEtaG(e);
  return xs;
}

Here is the call graph for this function:

G4double G4eeToPGammaModel::PeakEnergy ( ) const

overridevirtual

Implements G4Vee2hadrons.

Definition at line 96 of file G4eeToPGammaModel.cc.

{
  return massR;
}

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

overridevirtual

Implements G4Vee2hadrons.

Definition at line 113 of file G4eeToPGammaModel.cc.

{
  G4double egam = 0.5*e*(1.0 - massP*massP/(massR*massR));
  G4double tkin = e - egam - massP;
  if(tkin < 0.0) tkin = 0.0;
  G4double cost;
  do {
    cost = 2.0*G4UniformRand() - 1.0;
    // Loop checking, 07-Aug-2015, Vladimir Ivanchenko
  } while( 2.0*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* p1 = 
     new G4DynamicParticle(particle,dir,tkin);
  G4DynamicParticle* p2 = 
     new G4DynamicParticle(G4Gamma::Gamma(),-dir,egam);
  newp->push_back(p1);
  newp->push_back(p2);
}

Here is the call graph for this function:

---

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

• geant4.10.03.p01/source/processes/electromagnetic/highenergy/include/G4eeToPGammaModel.hh
• geant4.10.03.p01/source/processes/electromagnetic/highenergy/src/G4eeToPGammaModel.cc