G4ee2KNeutralModel Class Reference

#include <G4ee2KNeutralModel.hh>

Inheritance diagram for G4ee2KNeutralModel:

Public Member Functions
	G4ee2KNeutralModel (G4eeCrossSections *)
virtual	~G4ee2KNeutralModel ()
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 G4ee2KNeutralModel.hh.

Constructor & Destructor Documentation

G4ee2KNeutralModel::G4ee2KNeutralModel

(

G4eeCrossSections *

cr

)

Definition at line 64 of file G4ee2KNeutralModel.cc.

References G4ParticleDefinition::GetPDGMass(), G4KaonZeroLong::KaonZeroLong(), and python.hepunit::MeV.

                                                           :
  cross(cr)
{
  massK = G4KaonZeroLong::KaonZeroLong()->GetPDGMass();
  massPhi = 1019.46*MeV;
}

G4ee2KNeutralModel::~G4ee2KNeutralModel ( )

virtual

Definition at line 73 of file G4ee2KNeutralModel.cc.

{}

Member Function Documentation

G4double G4ee2KNeutralModel::ComputeCrossSection ( G4double  e ) const

virtual

Implements G4Vee2hadrons.

Definition at line 92 of file G4ee2KNeutralModel.cc.

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

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

G4double G4ee2KNeutralModel::PeakEnergy ( ) const

virtual

Implements G4Vee2hadrons.

Definition at line 85 of file G4ee2KNeutralModel.cc.

{
  return massPhi;
}

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

virtual

Implements G4Vee2hadrons.

Definition at line 100 of file G4ee2KNeutralModel.cc.

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

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

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

virtual

Implements G4Vee2hadrons.

Definition at line 113 of file G4ee2KNeutralModel.cc.

References G4UniformRand, G4KaonZeroLong::KaonZeroLong(), G4KaonZeroShort::KaonZeroShort(), CLHEP::Hep3Vector::rotateUz(), and python.hepunit::twopi.

{

  G4double tkin = 0.5*e - massK;
  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* p1 = 
     new G4DynamicParticle(G4KaonZeroLong::KaonZeroLong(),dir,tkin);
  G4DynamicParticle* p2 = 
     new G4DynamicParticle(G4KaonZeroShort::KaonZeroShort(),-dir,tkin);
  newp->push_back(p1);
  newp->push_back(p2);
}

G4double G4ee2KNeutralModel::ThresholdEnergy ( ) const

virtual

Implements G4Vee2hadrons.

Definition at line 78 of file G4ee2KNeutralModel.cc.

{
  return 2.0*massK;
}

---

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

• G4ee2KNeutralModel.hh
• G4ee2KNeutralModel.cc