G4hBetheBlochModel Class Reference

#include <G4hBetheBlochModel.hh>

Inheritance diagram for G4hBetheBlochModel:

Collaboration diagram for G4hBetheBlochModel:

Public Member Functions
	G4hBetheBlochModel (const G4String &name)
	~G4hBetheBlochModel ()
G4double	TheValue (const G4DynamicParticle *particle, const G4Material *material)
G4double	TheValue (const G4ParticleDefinition *aParticle, const G4Material *material, G4double kineticEnergy)
G4double	HighEnergyLimit (const G4ParticleDefinition *aParticle, const G4Material *material) const
G4double	LowEnergyLimit (const G4ParticleDefinition *aParticle, const G4Material *material) const
G4double	HighEnergyLimit (const G4ParticleDefinition *aParticle) const
G4double	LowEnergyLimit (const G4ParticleDefinition *aParticle) const
G4bool	IsInCharge (const G4DynamicParticle *particle, const G4Material *material) const
G4bool	IsInCharge (const G4ParticleDefinition *aParticle, const G4Material *material) const
Public Member Functions inherited from G4VLowEnergyModel
	G4VLowEnergyModel (const G4String &name)
virtual	~G4VLowEnergyModel ()

Private Member Functions
G4double	BetheBlochFormula (const G4Material *material, G4double kineticEnergy, G4double particleMass) const

Private Attributes
G4double	lowEnergyLimit
G4double	highEnergyLimit
const G4double	twoln10
const G4double	bg2lim
const G4double	taulim

Detailed Description

Definition at line 59 of file G4hBetheBlochModel.hh.

Constructor & Destructor Documentation

G4hBetheBlochModel()

G4hBetheBlochModel::G4hBetheBlochModel

(

const G4String &

name

)

Definition at line 63 of file G4hBetheBlochModel.cc.

  : G4VLowEnergyModel(name), 
    lowEnergyLimit(1.*MeV),
    highEnergyLimit(100.*GeV),
    twoln10(2.*std::log(10.)),
    bg2lim(0.0169), 
    taulim(8.4146e-3)
{;}

~G4hBetheBlochModel()

G4hBetheBlochModel::~G4hBetheBlochModel

(

)

Definition at line 74 of file G4hBetheBlochModel.cc.

{;}

Member Function Documentation

BetheBlochFormula()

G4double G4hBetheBlochModel::BetheBlochFormula ( const G4Material *  material, G4double  kineticEnergy, G4double  particleMass  ) const

private

Definition at line 156 of file G4hBetheBlochModel.cc.

{
  // This member function is applied normally to proton/antiproton
  G4double ionloss ;

  G4double rateMass = electron_mass_c2/particleMass ;

  G4double taul = material->GetIonisation()->GetTaul() ;
  G4double tau  = kineticEnergy/particleMass ;    // tau is relative energy
  
  // It is not normal case for this function
  // for low energy parametrisation have to be applied
  if ( tau < taul ) tau = taul ; 
    
  // some local variables 
    
  G4double gamma,bg2,beta2,tmax,x,delta,sh ;
  G4double electronDensity = material->GetElectronDensity();
  G4double eexc  = material->GetIonisation()->GetMeanExcitationEnergy();
  G4double eexc2 = eexc*eexc ;
  G4double cden  = material->GetIonisation()->GetCdensity();
  G4double mden  = material->GetIonisation()->GetMdensity();
  G4double aden  = material->GetIonisation()->GetAdensity();
  G4double x0den = material->GetIonisation()->GetX0density();
  G4double x1den = material->GetIonisation()->GetX1density();
  G4double* shellCorrectionVector =
            material->GetIonisation()->GetShellCorrectionVector();
    
  gamma = tau + 1.0 ;
  bg2 = tau*(tau+2.0) ;
  beta2 = bg2/(gamma*gamma) ;
  tmax = 2.*electron_mass_c2*bg2/(1.+2.*gamma*rateMass+rateMass*rateMass) ;
        
  ionloss = std::log(2.0*electron_mass_c2*bg2*tmax/eexc2)-2.0*beta2 ;
    
  // density correction     
  x = std::log(bg2)/twoln10 ;
  if ( x < x0den ) {
    delta = 0.0 ;

  } else { 
    delta = twoln10*x - cden ;
    if ( x < x1den ) delta += aden*std::pow((x1den-x),mden) ;
  }
    
  // shell correction 
  sh = 0.0 ;      
  x  = 1.0 ;

  if ( bg2 > bg2lim ) {
    for (G4int k=0; k<=2; k++) {
    x *= bg2 ;
    sh += shellCorrectionVector[k]/x;
    }

  } else {
    for (G4int k=0; k<=2; k++) {
    x *= bg2lim ;
    sh += shellCorrectionVector[k]/x;
    }
    sh *= std::log(tau/taul)/std::log(taulim/taul) ;     
  }
    
  // now compute the total ionization loss
    
  ionloss -= delta + sh ;
  ionloss *= twopi_mc2_rcl2*electronDensity/beta2 ;

  if ( ionloss < 0.0) ionloss = 0.0 ;
  
  return ionloss;
}

Here is the call graph for this function:

Here is the caller graph for this function:

HighEnergyLimit() [1/2]

G4double G4hBetheBlochModel::HighEnergyLimit ( const G4ParticleDefinition *  aParticle, const G4Material *  material  ) const

virtual

Implements G4VLowEnergyModel.

Definition at line 104 of file G4hBetheBlochModel.cc.

{
  return highEnergyLimit ;
}

HighEnergyLimit() [2/2]

G4double G4hBetheBlochModel::HighEnergyLimit ( const G4ParticleDefinition *  aParticle ) const

virtual

Implements G4VLowEnergyModel.

Definition at line 124 of file G4hBetheBlochModel.cc.

{
  return highEnergyLimit ;
}

IsInCharge() [1/2]

G4bool G4hBetheBlochModel::IsInCharge ( const G4DynamicParticle *  particle, const G4Material *  material  ) const

virtual

Implements G4VLowEnergyModel.

Definition at line 140 of file G4hBetheBlochModel.cc.

{
  return true ;
}

IsInCharge() [2/2]

G4bool G4hBetheBlochModel::IsInCharge ( const G4ParticleDefinition *  aParticle, const G4Material *  material  ) const

virtual

Implements G4VLowEnergyModel.

Definition at line 148 of file G4hBetheBlochModel.cc.

{
  return true ;
}

LowEnergyLimit() [1/2]

G4double G4hBetheBlochModel::LowEnergyLimit ( const G4ParticleDefinition *  aParticle, const G4Material *  material  ) const

virtual

Implements G4VLowEnergyModel.

Definition at line 113 of file G4hBetheBlochModel.cc.

{
  G4double taul = (material->GetIonisation()->GetTaul())*
                  (aParticle->GetPDGMass()) ;
  return taul ;
}

Here is the call graph for this function:

LowEnergyLimit() [2/2]

G4double G4hBetheBlochModel::LowEnergyLimit ( const G4ParticleDefinition *  aParticle ) const

virtual

Implements G4VLowEnergyModel.

Definition at line 132 of file G4hBetheBlochModel.cc.

{
  return lowEnergyLimit ;
}

TheValue() [1/2]

G4double G4hBetheBlochModel::TheValue ( const G4DynamicParticle *  particle, const G4Material *  material  )

virtual

Implements G4VLowEnergyModel.

Definition at line 79 of file G4hBetheBlochModel.cc.

{
  G4double energy = particle->GetKineticEnergy() ;
  G4double particleMass = particle->GetMass() ;

  G4double eloss  = BetheBlochFormula(material,energy,particleMass) ;

  return eloss ;
}

Here is the call graph for this function:

TheValue() [2/2]

G4double G4hBetheBlochModel::TheValue ( const G4ParticleDefinition *  aParticle, const G4Material *  material, G4double  kineticEnergy  )

virtual

Implements G4VLowEnergyModel.

Definition at line 92 of file G4hBetheBlochModel.cc.

{
  G4double particleMass = aParticle->GetPDGMass() ;
  G4double eloss  = BetheBlochFormula(material,kineticEnergy,particleMass) ;

  return eloss ;
}

Here is the call graph for this function:

Member Data Documentation

bg2lim

const G4double G4hBetheBlochModel::bg2lim

private

Definition at line 106 of file G4hBetheBlochModel.hh.

highEnergyLimit

G4double G4hBetheBlochModel::highEnergyLimit

private

Definition at line 101 of file G4hBetheBlochModel.hh.

lowEnergyLimit

G4double G4hBetheBlochModel::lowEnergyLimit

private

Definition at line 100 of file G4hBetheBlochModel.hh.

taulim

const G4double G4hBetheBlochModel::taulim

private

Definition at line 107 of file G4hBetheBlochModel.hh.

twoln10

const G4double G4hBetheBlochModel::twoln10

private

Definition at line 105 of file G4hBetheBlochModel.hh.

---

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

• G4hBetheBlochModel.hh
• G4hBetheBlochModel.cc