G4hNuclearStoppingModel Class Reference

#include <G4hNuclearStoppingModel.hh>

Inheritance diagram for G4hNuclearStoppingModel:

Collaboration diagram for G4hNuclearStoppingModel:

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

Private Member Functions
	G4hNuclearStoppingModel (G4hNuclearStoppingModel &)
G4hNuclearStoppingModel &	operator= (const G4hNuclearStoppingModel &right)
void	InitializeMe ()
G4double	StoppingPower (const G4Material *material, G4double kineticEnergy, G4double z1, G4double m1) const

Private Attributes
G4VhNuclearStoppingPower *	nStopingPowerTable
G4double	factorPDG2AMU
G4double	theZieglerFactor
G4String	modelName
G4double	lowEnergyLimit
G4double	highEnergyLimit

Detailed Description

Definition at line 60 of file G4hNuclearStoppingModel.hh.

Constructor & Destructor Documentation

G4hNuclearStoppingModel() [1/2]

G4hNuclearStoppingModel::G4hNuclearStoppingModel

(

const G4String &

name

)

Definition at line 69 of file G4hNuclearStoppingModel.cc.

  :G4VLowEnergyModel(name), modelName(name)
{
  InitializeMe() ;
}

Here is the call graph for this function:

Here is the caller graph for this function:

~G4hNuclearStoppingModel()

G4hNuclearStoppingModel::~G4hNuclearStoppingModel

(

)

Definition at line 111 of file G4hNuclearStoppingModel.cc.

{
  delete nStopingPowerTable;
}

G4hNuclearStoppingModel() [2/2]

G4hNuclearStoppingModel::G4hNuclearStoppingModel ( G4hNuclearStoppingModel &  )

private

Member Function Documentation

HighEnergyLimit() [1/2]

G4double G4hNuclearStoppingModel::HighEnergyLimit ( const G4ParticleDefinition *  , const G4Material *    ) const

inlinevirtual

Implements G4VLowEnergyModel.

Definition at line 76 of file G4hNuclearStoppingModel.hh.

  {return highEnergyLimit;};

HighEnergyLimit() [2/2]

G4double G4hNuclearStoppingModel::HighEnergyLimit ( const G4ParticleDefinition *  ) const

inlinevirtual

Implements G4VLowEnergyModel.

Definition at line 84 of file G4hNuclearStoppingModel.hh.

  {return highEnergyLimit;};

InitializeMe()

void G4hNuclearStoppingModel::InitializeMe ( )

private

Definition at line 77 of file G4hNuclearStoppingModel.cc.

{
  // Constants
  highEnergyLimit = 100.*MeV ;
  lowEnergyLimit  = 1.*eV ;
  factorPDG2AMU   = 1.007276/proton_mass_c2 ;
  theZieglerFactor= eV*cm2*1.0e-15 ; 

  // Registration of parametrisation models of nuclear energy losses
  G4String blank  = G4String(" ") ;
  G4String ir49   = G4String("ICRU_R49") ;
  G4String zi85   = G4String("Ziegler1985") ;
  if(ir49 == modelName || blank == modelName) {
      nStopingPowerTable = new G4hICRU49Nuclear();

  } else if(zi85 == modelName) {
      nStopingPowerTable = new G4hZiegler1985Nuclear();
        
  } else {
    G4cout << 
    "G4hLowEnergyIonisation warning: There is no table with the modelName <" 
       << modelName << ">" 
       << " for nuclear stopping, <ICRU_R49> is applied " 
       << G4endl; 
    nStopingPowerTable = new G4hICRU49Nuclear();
  }

  // Default is nuclear stopping fluctuations On
  //  nStopingPowerTable->SetNuclearStoppingFluctuationsOn();  
  nStopingPowerTable->SetNuclearStoppingFluctuationsOff();  
}

Here is the call graph for this function:

Here is the caller graph for this function:

IsInCharge() [1/2]

G4bool G4hNuclearStoppingModel::IsInCharge ( const G4DynamicParticle *  , const G4Material *    ) const

inlinevirtual

Implements G4VLowEnergyModel.

Definition at line 90 of file G4hNuclearStoppingModel.hh.

  {return true;};

IsInCharge() [2/2]

G4bool G4hNuclearStoppingModel::IsInCharge ( const G4ParticleDefinition *  , const G4Material *    ) const

inlinevirtual

Implements G4VLowEnergyModel.

Definition at line 94 of file G4hNuclearStoppingModel.hh.

  {return true;};

LowEnergyLimit() [1/2]

G4double G4hNuclearStoppingModel::LowEnergyLimit ( const G4ParticleDefinition *  , const G4Material *    ) const

inlinevirtual

Implements G4VLowEnergyModel.

Definition at line 80 of file G4hNuclearStoppingModel.hh.

  {return lowEnergyLimit;};

LowEnergyLimit() [2/2]

G4double G4hNuclearStoppingModel::LowEnergyLimit ( const G4ParticleDefinition *  ) const

inlinevirtual

Implements G4VLowEnergyModel.

Definition at line 87 of file G4hNuclearStoppingModel.hh.

  {return lowEnergyLimit;};

operator=()

G4hNuclearStoppingModel& G4hNuclearStoppingModel::operator= ( const G4hNuclearStoppingModel &  right )

private

Here is the caller graph for this function:

SetNuclearStoppingFluctuationsOff()

void G4hNuclearStoppingModel::SetNuclearStoppingFluctuationsOff ( )

inline

Definition at line 101 of file G4hNuclearStoppingModel.hh.

  {nStopingPowerTable->SetNuclearStoppingFluctuationsOff();}; 

Here is the call graph for this function:

SetNuclearStoppingFluctuationsOn()

void G4hNuclearStoppingModel::SetNuclearStoppingFluctuationsOn ( )

inline

Definition at line 98 of file G4hNuclearStoppingModel.hh.

  {nStopingPowerTable->SetNuclearStoppingFluctuationsOn();}; 

Here is the call graph for this function:

StoppingPower()

G4double G4hNuclearStoppingModel::StoppingPower ( const G4Material *  material, G4double  kineticEnergy, G4double  z1, G4double  m1  ) const

private

Definition at line 151 of file G4hNuclearStoppingModel.cc.

{
  // Target nucleus
  G4int NumberOfElements = material->GetNumberOfElements() ;
  if(0 == NumberOfElements) return 0.0 ;

  const G4ElementVector* theElementVector = 
                                 material->GetElementVector() ;
  const G4double* theAtomicNumDensityVector = 
                                 material->GetAtomicNumDensityVector() ;

  //  loop for the elements in the material

  G4double nloss = 0.0;

  for (G4int iel=0; iel<NumberOfElements; iel++) {
    const G4Element* element = (*theElementVector)[iel] ;
    G4double z2 = element->GetZ();
    G4double m2Local = element->GetA()*mole/g ;
    nloss += (nStopingPowerTable->
              NuclearStoppingPower(kineticEnergy, z1, z2, m1, m2Local))
           * theAtomicNumDensityVector[iel] ;    
  }

  return nloss;
}

Here is the call graph for this function:

Here is the caller graph for this function:

TheValue() [1/2]

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

virtual

Implements G4VLowEnergyModel.

Definition at line 118 of file G4hNuclearStoppingModel.cc.

{
  // Projectile nucleus
  G4double energy = particle->GetKineticEnergy() ;
  G4double z1 = std::abs((particle->GetCharge())/eplus) ;
  G4double m1 = (particle->GetMass())*factorPDG2AMU ;

  G4double nloss = StoppingPower(material, energy, z1, m1) * theZieglerFactor; 

  return nloss;
}

Here is the call graph for this function:

TheValue() [2/2]

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

virtual

Implements G4VLowEnergyModel.

Definition at line 134 of file G4hNuclearStoppingModel.cc.

{
  // Projectile nucleus
  G4double z1 = std::abs((aParticle->GetPDGCharge())/eplus) ;
  G4double m1 = (aParticle->GetPDGMass())*factorPDG2AMU ;

  G4double nloss = StoppingPower(material, kineticEnergy, z1, m1)
                 * theZieglerFactor; 

  return nloss;
}

Here is the call graph for this function:

Member Data Documentation

factorPDG2AMU

G4double G4hNuclearStoppingModel::factorPDG2AMU

private

Definition at line 123 of file G4hNuclearStoppingModel.hh.

highEnergyLimit

G4double G4hNuclearStoppingModel::highEnergyLimit

private

Definition at line 131 of file G4hNuclearStoppingModel.hh.

lowEnergyLimit

G4double G4hNuclearStoppingModel::lowEnergyLimit

private

Definition at line 130 of file G4hNuclearStoppingModel.hh.

modelName

G4String G4hNuclearStoppingModel::modelName

private

Definition at line 129 of file G4hNuclearStoppingModel.hh.

nStopingPowerTable

G4VhNuclearStoppingPower* G4hNuclearStoppingModel::nStopingPowerTable

private

Definition at line 121 of file G4hNuclearStoppingModel.hh.

theZieglerFactor

G4double G4hNuclearStoppingModel::theZieglerFactor

private

Definition at line 126 of file G4hNuclearStoppingModel.hh.

---

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

• G4hNuclearStoppingModel.hh
• G4hNuclearStoppingModel.cc