G4hIonisation Class Reference

#include <G4hIonisation.hh>

Inheritance diagram for G4hIonisation:

Public Member Functions
	G4hIonisation (const G4String &name="hIoni")
virtual	~G4hIonisation ()
virtual G4bool	IsApplicable (const G4ParticleDefinition &p)
virtual G4double	MinPrimaryEnergy (const G4ParticleDefinition *p, const G4Material *, G4double cut)
virtual void	PrintInfo ()
Protected Member Functions
virtual void	InitialiseEnergyLossProcess (const G4ParticleDefinition *, const G4ParticleDefinition *)

---

Detailed Description

Definition at line 85 of file G4hIonisation.hh.

---

Constructor & Destructor Documentation

G4hIonisation::G4hIonisation

(

const G4String &

name = "hIoni"

)

Definition at line 109 of file G4hIonisation.cc.

References G4Electron::Electron(), fIonisation, G4VProcess::SetProcessSubType(), G4VEnergyLossProcess::SetSecondaryParticle(), and G4VEnergyLossProcess::SetStepFunction().

  : G4VEnergyLossProcess(name),
    isInitialised(false)
{
  SetStepFunction(0.2, 0.1*mm);
  SetProcessSubType(fIonisation);
  SetSecondaryParticle(G4Electron::Electron());
  mass = 0.0;
  ratio = 0.0;
  eth = 2*MeV;
}

G4hIonisation::~G4hIonisation

(

)

[virtual]

Definition at line 123 of file G4hIonisation.cc.

{}

---

Member Function Documentation

void G4hIonisation::InitialiseEnergyLossProcess	(	const G4ParticleDefinition *	,
		const G4ParticleDefinition *
	)			[protected, virtual]

Implements G4VEnergyLossProcess.

Definition at line 147 of file G4hIonisation.cc.

References G4VEnergyLossProcess::AddEmModel(), G4AntiProton::AntiProton(), G4VEnergyLossProcess::EmModel(), G4VEnergyLossProcess::FluctModel(), G4ParticleDefinition::GetParticleName(), G4ParticleDefinition::GetPDGCharge(), G4ParticleDefinition::GetPDGMass(), G4ParticleDefinition::GetPDGSpin(), G4KaonMinus::KaonMinus(), G4KaonPlus::KaonPlus(), G4VEnergyLossProcess::MaxKinEnergy(), G4VEnergyLossProcess::MinKinEnergy(), G4Proton::Proton(), G4VEnergyLossProcess::SetBaseParticle(), G4VEnergyLossProcess::SetEmModel(), G4VEnergyLossProcess::SetFluctModel(), G4VEmModel::SetHighEnergyLimit(), and G4VEmModel::SetLowEnergyLimit().

{
  if(!isInitialised) {

    const G4ParticleDefinition* theBaseParticle = 0;
    G4String pname = part->GetParticleName();
    G4double q = part->GetPDGCharge();

    // standard base particles
    if(part == bpart || pname == "proton" ||
       pname == "anti_proton" || 
       pname == "pi+" || pname == "pi-" || 
       pname == "kaon+" || pname == "kaon-") 
      { 
        theBaseParticle = 0;
      }
    // select base particle 
    else if(bpart == 0) {

      if(part->GetPDGSpin() == 0.0) {
        if(q > 0.0) { theBaseParticle = G4KaonPlus::KaonPlus(); }
        else { theBaseParticle = G4KaonMinus::KaonMinus(); }
      } else {
        if(q > 0.0) { theBaseParticle = G4Proton::Proton(); } 
        else { theBaseParticle = G4AntiProton::AntiProton(); }
      }

      // base particle defined by interface
    } else { 
      theBaseParticle = bpart;
    }
    SetBaseParticle(theBaseParticle);

    mass  = part->GetPDGMass();
    ratio = electron_mass_c2/mass;
    eth   = 2.0*MeV*mass/proton_mass_c2;

    if (!EmModel(1)) { 
      if(q > 0.0) { SetEmModel(new G4BraggModel(),1); }
      else { SetEmModel(new G4ICRU73QOModel(),1); }
    }
    EmModel(1)->SetLowEnergyLimit(MinKinEnergy());

    // model limit defined for protons
    //eth = (EmModel(1)->HighEnergyLimit())*mass/proton_mass_c2;
    EmModel(1)->SetHighEnergyLimit(eth);
    AddEmModel(1, EmModel(1), new G4IonFluctuations());

    if (!FluctModel()) { SetFluctModel(new G4UniversalFluctuation()); }

    if (!EmModel(2)) { SetEmModel(new G4BetheBlochModel(),2); }
    EmModel(2)->SetLowEnergyLimit(eth);
    EmModel(2)->SetHighEnergyLimit(MaxKinEnergy());
    AddEmModel(2, EmModel(2), FluctModel());  

    isInitialised = true;
  }
}

G4bool G4hIonisation::IsApplicable

(

const G4ParticleDefinition &

p

)

[virtual]

Implements G4VEnergyLossProcess.

Definition at line 128 of file G4hIonisation.cc.

References G4ParticleDefinition::GetPDGCharge(), G4ParticleDefinition::GetPDGMass(), and G4ParticleDefinition::IsShortLived().

{
  return (p.GetPDGCharge() != 0.0 && p.GetPDGMass() > 10.0*MeV &&
         !p.IsShortLived());
}

G4double G4hIonisation::MinPrimaryEnergy	(	const G4ParticleDefinition *	p,
		const G4Material *	,
		G4double	cut
	)			[virtual]

Reimplemented from G4VEnergyLossProcess.

Definition at line 136 of file G4hIonisation.cc.

References G4InuclParticleNames::gam.

{
  G4double x = 0.5*cut/electron_mass_c2;
  G4double gam = x*ratio + std::sqrt((1. + x)*(1. + x*ratio*ratio));
  return mass*(gam - 1.0);
}

void G4hIonisation::PrintInfo

(

)

[virtual]

Implements G4VEnergyLossProcess.

Definition at line 210 of file G4hIonisation.cc.

{}

---

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

• G4hIonisation.hh
• G4hIonisation.cc

---