G4eeToPGammaModel Class Reference

#include <G4eeToPGammaModel.hh>

Inheritance diagram for G4eeToPGammaModel:

Public Member Functions
	G4eeToPGammaModel (G4eeCrossSections *, const G4String &)
virtual	~G4eeToPGammaModel ()
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 ()
virtual G4double	ThresholdEnergy () const =0
virtual G4double	PeakEnergy () const =0
virtual G4double	ComputeCrossSection (G4double) const =0
virtual G4PhysicsVector *	PhysicsVector (G4double, G4double) const =0
virtual void	SampleSecondaries (std::vector< G4DynamicParticle * > *, G4double, const G4ThreeVector &)=0
void	SetLowEnergy (G4double val)
G4double	LowEnergy () const
void	SetHighEnergy (G4double val)
G4double	HighEnergy () const

Detailed Description

Definition at line 60 of file G4eeToPGammaModel.hh.

Constructor & Destructor Documentation

G4eeToPGammaModel()

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

Definition at line 65 of file G4eeToPGammaModel.cc.

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

~G4eeToPGammaModel()

G4eeToPGammaModel::~G4eeToPGammaModel ( )

virtual

Definition at line 81 of file G4eeToPGammaModel.cc.

{}

Member Function Documentation

ComputeCrossSection()

G4double G4eeToPGammaModel::ComputeCrossSection ( G4double  e ) const

virtual

Implements G4Vee2hadrons.

Definition at line 100 of file G4eeToPGammaModel.cc.

{
  G4double ee = std::min(HighEnergy(),e);
  G4double xs;
  if(particle == pi0) xs = cross->CrossSectionPi0G(ee);
  else                xs = cross->CrossSectionEtaG(ee);
  return xs;
}

PeakEnergy()

G4double G4eeToPGammaModel::PeakEnergy ( ) const

virtual

Implements G4Vee2hadrons.

Definition at line 93 of file G4eeToPGammaModel.cc.

{
  return massR;
}

PhysicsVector()

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

virtual

Implements G4Vee2hadrons.

Definition at line 111 of file G4eeToPGammaModel.cc.

{
  G4double tmin = std::max(emin, ThresholdEnergy());
  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;
}

SampleSecondaries()

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

virtual

Implements G4Vee2hadrons.

Definition at line 124 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;
  } 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);
}

ThresholdEnergy()

G4double G4eeToPGammaModel::ThresholdEnergy ( ) const

virtual

Implements G4Vee2hadrons.

Definition at line 86 of file G4eeToPGammaModel.cc.

{
  return LowEnergy();
}

Referenced by PhysicsVector().

---

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

• G4eeToPGammaModel.hh
• G4eeToPGammaModel.cc