G4eeToHadronsMultiModel Class Reference

#include <G4eeToHadronsMultiModel.hh>

Inheritance diagram for G4eeToHadronsMultiModel:

Public Member Functions
	G4eeToHadronsMultiModel (G4int ver=0, const G4String &nam="eeToHadrons")
virtual	~G4eeToHadronsMultiModel ()
virtual void	Initialise (const G4ParticleDefinition *, const G4DataVector &)
virtual G4double	CrossSectionPerVolume (const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy, G4double maxEnergy)
virtual G4double	ComputeCrossSectionPerAtom (const G4ParticleDefinition *, G4double kineticEnergy, G4double Z, G4double A, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
virtual void	SampleSecondaries (std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin=0.0, G4double maxEnergy=DBL_MAX)
virtual void	PrintInfo ()
void	SetCrossSecFactor (G4double fac)
G4double	ComputeCrossSectionPerElectron (const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)

---

Detailed Description

Definition at line 62 of file G4eeToHadronsMultiModel.hh.

---

Constructor & Destructor Documentation

G4eeToHadronsMultiModel::G4eeToHadronsMultiModel	(	G4int	ver = 0,
		const G4String &	nam = "eeToHadrons"
	)

Definition at line 65 of file G4eeToHadronsMultiModel.cc.

References DBL_MAX.

  : G4VEmModel(mname),
    csFactor(1.0),
    nModels(0),
    verbose(ver),
    isInitialised(false)
{
  thKineticEnergy  = DBL_MAX;
  maxKineticEnergy = 1.2*GeV;
  fParticleChange  = 0;
  cross = 0;
}

G4eeToHadronsMultiModel::~G4eeToHadronsMultiModel

(

)

[virtual]

Definition at line 80 of file G4eeToHadronsMultiModel.cc.

References CLHEP::detail::n.

{
  G4int n = models.size();
  if(n>0) {
    for(G4int i=0; i<n; i++) {
      delete models[i];
    }
  }
  delete cross;
}

---

Member Function Documentation

G4double G4eeToHadronsMultiModel::ComputeCrossSectionPerAtom	(	const G4ParticleDefinition *	,
		G4double	kineticEnergy,
		G4double	Z,
		G4double	A,
		G4double	cutEnergy = 0.0,
		G4double	maxEnergy = DBL_MAX
	)			[virtual]

Reimplemented from G4VEmModel.

Definition at line 168 of file G4eeToHadronsMultiModel.cc.

References ComputeCrossSectionPerElectron().

{
  return Z*ComputeCrossSectionPerElectron(p, kineticEnergy);
}

G4double G4eeToHadronsMultiModel::ComputeCrossSectionPerElectron	(	const G4ParticleDefinition *	,
		G4double	kineticEnergy,
		G4double	cutEnergy = 0.0,
		G4double	maxEnergy = DBL_MAX
	)			[inline]

Definition at line 132 of file G4eeToHadronsMultiModel.hh.

Referenced by ComputeCrossSectionPerAtom(), and CrossSectionPerVolume().

{
  G4double res = 0.0;
  if (kineticEnergy > thKineticEnergy) {
    for(G4int i=0; i<nModels; i++) {
      if(kineticEnergy >= ekinMin[i] && kineticEnergy <= ekinMax[i])
        res += (models[i])->ComputeCrossSectionPerElectron(0,kineticEnergy);
      cumSum[i] = res;
    }
  }
  return res*csFactor;
}

G4double G4eeToHadronsMultiModel::CrossSectionPerVolume	(	const G4Material *	,
		const G4ParticleDefinition *	,
		G4double	kineticEnergy,
		G4double	cutEnergy,
		G4double	maxEnergy
	)			[virtual]

Reimplemented from G4VEmModel.

Definition at line 156 of file G4eeToHadronsMultiModel.cc.

References ComputeCrossSectionPerElectron(), and G4Material::GetElectronDensity().

{
  return mat->GetElectronDensity()*
    ComputeCrossSectionPerElectron(p, kineticEnergy);
}

void G4eeToHadronsMultiModel::Initialise	(	const G4ParticleDefinition *	,
		const G4DataVector &
	)			[virtual]

Implements G4VEmModel.

Definition at line 93 of file G4eeToHadronsMultiModel.cc.

References G4VEmModel::GetParticleChangeForGamma(), G4Vee2hadrons::SetHighEnergy(), and G4Vee2hadrons::SetLowEnergy().

{
  if(!isInitialised) {
    isInitialised = true;

    cross = new G4eeCrossSections();

    G4eeToTwoPiModel* m2pi = new G4eeToTwoPiModel(cross);
    m2pi->SetHighEnergy(maxKineticEnergy);
    AddEEModel(m2pi);

    G4eeTo3PiModel* m3pi1 = new G4eeTo3PiModel(cross);
    m3pi1->SetHighEnergy(0.95*GeV);
    AddEEModel(m3pi1);

    G4eeTo3PiModel* m3pi2 = new G4eeTo3PiModel(cross);
    m3pi2->SetLowEnergy(0.95*GeV);
    m3pi2->SetHighEnergy(maxKineticEnergy);
    AddEEModel(m3pi2);

    G4ee2KChargedModel* m2kc = new G4ee2KChargedModel(cross);
    m2kc->SetHighEnergy(maxKineticEnergy);
    AddEEModel(m2kc);

    G4ee2KNeutralModel* m2kn = new G4ee2KNeutralModel(cross);
    m2kn->SetHighEnergy(maxKineticEnergy);
    AddEEModel(m2kn);

    G4eeToPGammaModel* mpg1 = new G4eeToPGammaModel(cross,"pi0");
    mpg1->SetLowEnergy(0.7*GeV);
    mpg1->SetHighEnergy(maxKineticEnergy);
    AddEEModel(mpg1);

    G4eeToPGammaModel* mpg2 = new G4eeToPGammaModel(cross,"eta");
    mpg2->SetLowEnergy(0.7*GeV);
    mpg2->SetHighEnergy(maxKineticEnergy);
    AddEEModel(mpg2);

    nModels = models.size();

    fParticleChange = GetParticleChangeForGamma();
  }
}

void G4eeToHadronsMultiModel::PrintInfo

(

)

[virtual]

Definition at line 200 of file G4eeToHadronsMultiModel.cc.

References G4cout, and G4endl.

Referenced by G4eeToHadrons::PrintInfo().

{
  if(verbose > 0) {
    G4double e1 = 0.5*thKineticEnergy*thKineticEnergy/electron_mass_c2 
      - 2.0*electron_mass_c2; 
    G4double e2 = 0.5*maxKineticEnergy*maxKineticEnergy/electron_mass_c2 
      - 2.0*electron_mass_c2; 
    G4cout << "      e+ annihilation into hadrons active from "
           << e1/GeV << " GeV to " << e2/GeV << " GeV"
           << G4endl;
  }
}

void G4eeToHadronsMultiModel::SampleSecondaries	(	std::vector< G4DynamicParticle * > *	,
		const G4MaterialCutsCouple *	,
		const G4DynamicParticle *	,
		G4double	tmin = 0.0,
		G4double	maxEnergy = DBL_MAX
	)			[virtual]

Implements G4VEmModel.

Definition at line 180 of file G4eeToHadronsMultiModel.cc.

References fStopAndKill, G4UniformRand, G4DynamicParticle::GetKineticEnergy(), and G4VParticleChange::ProposeTrackStatus().

{
  G4double kinEnergy = dp->GetKineticEnergy();
  if (kinEnergy > thKineticEnergy) {
    G4double q = cumSum[nModels-1]*G4UniformRand();
    for(G4int i=0; i<nModels; i++) {
      if(q <= cumSum[i]) {
        (models[i])->SampleSecondaries(newp, couple,dp);
        if(newp->size() > 0) fParticleChange->ProposeTrackStatus(fStopAndKill);
        break;
      }
    }
  }
}

void G4eeToHadronsMultiModel::SetCrossSecFactor

(

G4double

fac

)

Definition at line 215 of file G4eeToHadronsMultiModel.cc.

References G4cout, and G4endl.

Referenced by G4eeToHadrons::InitialiseProcess(), and G4eeToHadrons::SetCrossSecFactor().

{
  if(fac > 1.0) {
    csFactor = fac;
    if(verbose > 0)
      G4cout << "### G4eeToHadronsMultiModel: The cross section for G4eeToHadronsMultiModel "
             << " is increased by the Factor= " << csFactor << G4endl;
  }
}

---

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

• G4eeToHadronsMultiModel.hh
• G4eeToHadronsMultiModel.cc

---