Geant4-11
Public Member Functions | Private Member Functions | Private Attributes
G4hNuclearStoppingModel Class Reference

#include <G4hNuclearStoppingModel.hh>

Inheritance diagram for G4hNuclearStoppingModel:
G4VLowEnergyModel

Public Member Functions

 G4hNuclearStoppingModel (const G4String &name)
 
 G4hNuclearStoppingModel (G4hNuclearStoppingModel &)=delete
 
G4double HighEnergyLimit (const G4ParticleDefinition *) const override
 
G4double HighEnergyLimit (const G4ParticleDefinition *, const G4Material *) const override
 
G4bool IsInCharge (const G4DynamicParticle *, const G4Material *) const override
 
G4bool IsInCharge (const G4ParticleDefinition *, const G4Material *) const override
 
G4double LowEnergyLimit (const G4ParticleDefinition *) const override
 
G4double LowEnergyLimit (const G4ParticleDefinition *, const G4Material *) const override
 
G4hNuclearStoppingModel & operator= (const G4hNuclearStoppingModel &right)=delete
 
void SetNuclearStoppingFluctuationsOff ()
 
void SetNuclearStoppingFluctuationsOn ()
 
G4double TheValue (const G4DynamicParticle *particle, const G4Material *material) override
 
G4double TheValue (const G4ParticleDefinition *aParticle, const G4Material *material, G4double kineticEnergy) override
 
 ~G4hNuclearStoppingModel ()
 

Private Member Functions

void InitializeMe ()
 
G4double StoppingPower (const G4Material *material, G4double kineticEnergy, G4double z1, G4double m1) const
 

Private Attributes

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

Detailed Description

Definition at line 60 of file G4hNuclearStoppingModel.hh.

Constructor & Destructor Documentation

◆ G4hNuclearStoppingModel() [1/2]

G4hNuclearStoppingModel::G4hNuclearStoppingModel ( const G4String &  name)
explicit

Definition at line 69 of file G4hNuclearStoppingModel.cc.

70 :G4VLowEnergyModel(name), modelName(name)
71{
72 InitializeMe() ;
73}
G4VLowEnergyModel::G4VLowEnergyModel
G4VLowEnergyModel(const G4String &name)
Definition: G4VLowEnergyModel.cc:56
G4InuclParticleNames::name
const char * name(G4int ptype)
Definition: G4InuclParticleNames.hh:76

References InitializeMe().

◆ ~G4hNuclearStoppingModel()

G4hNuclearStoppingModel::~G4hNuclearStoppingModel ( )

Definition at line 111 of file G4hNuclearStoppingModel.cc.

112{
113 delete nStopingPowerTable;
114}
G4hNuclearStoppingModel::nStopingPowerTable
G4VhNuclearStoppingPower * nStopingPowerTable
Definition: G4hNuclearStoppingModel.hh:115

References nStopingPowerTable.

◆ G4hNuclearStoppingModel() [2/2]

G4hNuclearStoppingModel::G4hNuclearStoppingModel ( G4hNuclearStoppingModel &  )
delete

Member Function Documentation

◆ HighEnergyLimit() [1/2]

G4double G4hNuclearStoppingModel::HighEnergyLimit ( const G4ParticleDefinition *  ) const
inlineoverridevirtual

Implements G4VLowEnergyModel.

Definition at line 83 of file G4hNuclearStoppingModel.hh.

84 {return highEnergyLimit;};

References highEnergyLimit.

◆ HighEnergyLimit() [2/2]

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

Implements G4VLowEnergyModel.

Definition at line 75 of file G4hNuclearStoppingModel.hh.

77 {return highEnergyLimit;};

References highEnergyLimit.

◆ InitializeMe()

void G4hNuclearStoppingModel::InitializeMe ( )
private

Definition at line 77 of file G4hNuclearStoppingModel.cc.

78{
79 // Constants
80 highEnergyLimit = 100.*MeV ;
81 lowEnergyLimit = 1.*eV ;
82 factorPDG2AMU = 1.007276/proton_mass_c2 ;
83 theZieglerFactor= eV*cm2*1.0e-15 ;
84
85 // Registration of parametrisation models of nuclear energy losses
86 G4String blank = G4String(" ") ;
87 G4String ir49 = G4String("ICRU_R49") ;
88 G4String zi85 = G4String("Ziegler1985") ;
89 if(ir49 == modelName || blank == modelName) {
90 nStopingPowerTable = new G4hICRU49Nuclear();
91
92 } else if(zi85 == modelName) {
93 nStopingPowerTable = new G4hZiegler1985Nuclear();
94
95 } else {
96 G4cout <<
97 "G4hLowEnergyIonisation warning: There is no table with the modelName <"
98 << modelName << ">"
99 << " for nuclear stopping, <ICRU_R49> is applied "
100 << G4endl;
101 nStopingPowerTable = new G4hICRU49Nuclear();
102 }
103
104 // Default is nuclear stopping fluctuations On
105 // nStopingPowerTable->SetNuclearStoppingFluctuationsOn();
106 nStopingPowerTable->SetNuclearStoppingFluctuationsOff();
107}
cm2
static constexpr double cm2
Definition: G4SIunits.hh:100
eV
static constexpr double eV
Definition: G4SIunits.hh:201
MeV
static constexpr double MeV
Definition: G4SIunits.hh:200
G4endl
#define G4endl
Definition: G4ios.hh:57
G4cout
G4GLOB_DLL std::ostream G4cout
source.hepunit.proton_mass_c2
float proton_mass_c2
Definition: hepunit.py:274

References cm2, eV, factorPDG2AMU, G4cout, G4endl, highEnergyLimit, lowEnergyLimit, MeV, modelName, nStopingPowerTable, source.hepunit::proton_mass_c2, G4VhNuclearStoppingPower::SetNuclearStoppingFluctuationsOff(), and theZieglerFactor.

Referenced by G4hNuclearStoppingModel().

◆ IsInCharge() [1/2]

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

Implements G4VLowEnergyModel.

Definition at line 89 of file G4hNuclearStoppingModel.hh.

91 {return true;};

◆ IsInCharge() [2/2]

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

Implements G4VLowEnergyModel.

Definition at line 93 of file G4hNuclearStoppingModel.hh.

95 {return true;};

◆ LowEnergyLimit() [1/2]

G4double G4hNuclearStoppingModel::LowEnergyLimit ( const G4ParticleDefinition *  ) const
inlineoverridevirtual

Implements G4VLowEnergyModel.

Definition at line 86 of file G4hNuclearStoppingModel.hh.

87 {return lowEnergyLimit;};

References lowEnergyLimit.

◆ LowEnergyLimit() [2/2]

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

Implements G4VLowEnergyModel.

Definition at line 79 of file G4hNuclearStoppingModel.hh.

81 {return lowEnergyLimit;};

References lowEnergyLimit.

◆ operator=()

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

◆ SetNuclearStoppingFluctuationsOff()

void G4hNuclearStoppingModel::SetNuclearStoppingFluctuationsOff ( )
inline

Definition at line 100 of file G4hNuclearStoppingModel.hh.

101 {nStopingPowerTable->SetNuclearStoppingFluctuationsOff();};

References nStopingPowerTable, and G4VhNuclearStoppingPower::SetNuclearStoppingFluctuationsOff().

◆ SetNuclearStoppingFluctuationsOn()

void G4hNuclearStoppingModel::SetNuclearStoppingFluctuationsOn ( )
inline

Definition at line 97 of file G4hNuclearStoppingModel.hh.

98 {nStopingPowerTable->SetNuclearStoppingFluctuationsOn();};

References nStopingPowerTable, and G4VhNuclearStoppingPower::SetNuclearStoppingFluctuationsOn().

◆ StoppingPower()

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

Definition at line 151 of file G4hNuclearStoppingModel.cc.

155{
156 // Target nucleus
157 G4int NumberOfElements = material->GetNumberOfElements() ;
158 if(0 == NumberOfElements) return 0.0 ;
159
160 const G4ElementVector* theElementVector =
161 material->GetElementVector() ;
162 const G4double* theAtomicNumDensityVector =
163 material->GetAtomicNumDensityVector() ;
164
165 // loop for the elements in the material
166
167 G4double nloss = 0.0;
168
169 for (G4int iel=0; iel<NumberOfElements; iel++) {
170 const G4Element* element = (*theElementVector)[iel] ;
171 G4double z2 = element->GetZ();
172 G4double m2Local = element->GetA()*mole/g ;
173 nloss += (nStopingPowerTable->
174 NuclearStoppingPower(kineticEnergy, z1, z2, m1, m2Local))
175 * theAtomicNumDensityVector[iel] ;
176 }
177
178 return nloss;
179}
G4ElementVector
std::vector< const G4Element * > G4ElementVector
Definition: G4ElementVector.hh:42
mole
static constexpr double mole
Definition: G4SIunits.hh:279
g
static constexpr double g
Definition: G4SIunits.hh:168
G4double
double G4double
Definition: G4Types.hh:83
G4int
int G4int
Definition: G4Types.hh:85
G4Element::GetZ
G4double GetZ() const
Definition: G4Element.hh:131
G4Element::GetA
G4double GetA() const
Definition: G4Element.hh:139
eplot.material
string material
Definition: eplot.py:19

References g, G4Element::GetA(), G4Element::GetZ(), eplot::material, mole, and nStopingPowerTable.

Referenced by TheValue().

◆ TheValue() [1/2]

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

Implements G4VLowEnergyModel.

Definition at line 118 of file G4hNuclearStoppingModel.cc.

121{
122 // Projectile nucleus
123 G4double energy = particle->GetKineticEnergy() ;
124 G4double z1 = std::abs((particle->GetCharge())/eplus) ;
125 G4double m1 = (particle->GetMass())*factorPDG2AMU ;
126
127 G4double nloss = StoppingPower(material, energy, z1, m1) * theZieglerFactor;
128
129 return nloss;
130}
eplus
static constexpr double eplus
Definition: G4SIunits.hh:184
G4DynamicParticle::GetMass
G4double GetMass() const
G4DynamicParticle::GetCharge
G4double GetCharge() const
G4DynamicParticle::GetKineticEnergy
G4double GetKineticEnergy() const
G4hNuclearStoppingModel::StoppingPower
G4double StoppingPower(const G4Material *material, G4double kineticEnergy, G4double z1, G4double m1) const
Definition: G4hNuclearStoppingModel.cc:151
G4INCL::KinematicsUtils::energy
G4double energy(const ThreeVector &p, const G4double m)
Definition: G4INCLKinematicsUtils.cc:158

References G4INCL::KinematicsUtils::energy(), eplus, factorPDG2AMU, G4DynamicParticle::GetCharge(), G4DynamicParticle::GetKineticEnergy(), G4DynamicParticle::GetMass(), eplot::material, StoppingPower(), and theZieglerFactor.

◆ TheValue() [2/2]

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

Implements G4VLowEnergyModel.

Definition at line 134 of file G4hNuclearStoppingModel.cc.

138{
139 // Projectile nucleus
140 G4double z1 = std::abs((aParticle->GetPDGCharge())/eplus) ;
141 G4double m1 = (aParticle->GetPDGMass())*factorPDG2AMU ;
142
143 G4double nloss = StoppingPower(material, kineticEnergy, z1, m1)
144 * theZieglerFactor;
145
146 return nloss;
147}
G4ParticleDefinition::GetPDGCharge
G4double GetPDGCharge() const
Definition: G4ParticleDefinition.hh:113

References eplus, factorPDG2AMU, G4ParticleDefinition::GetPDGCharge(), G4ParticleDefinition::GetPDGMass(), eplot::material, StoppingPower(), and theZieglerFactor.

Field Documentation

◆ factorPDG2AMU

G4double G4hNuclearStoppingModel::factorPDG2AMU
private

Definition at line 118 of file G4hNuclearStoppingModel.hh.

Referenced by InitializeMe(), and TheValue().

◆ highEnergyLimit

G4double G4hNuclearStoppingModel::highEnergyLimit
private

Definition at line 125 of file G4hNuclearStoppingModel.hh.

Referenced by HighEnergyLimit(), and InitializeMe().

◆ lowEnergyLimit

G4double G4hNuclearStoppingModel::lowEnergyLimit
private

Definition at line 124 of file G4hNuclearStoppingModel.hh.

Referenced by InitializeMe(), and LowEnergyLimit().

◆ modelName

G4String G4hNuclearStoppingModel::modelName
private

Definition at line 116 of file G4hNuclearStoppingModel.hh.

Referenced by InitializeMe().

◆ nStopingPowerTable

G4VhNuclearStoppingPower* G4hNuclearStoppingModel::nStopingPowerTable
private

Definition at line 115 of file G4hNuclearStoppingModel.hh.

Referenced by InitializeMe(), SetNuclearStoppingFluctuationsOff(), SetNuclearStoppingFluctuationsOn(), StoppingPower(), and ~G4hNuclearStoppingModel().

◆ theZieglerFactor

G4double G4hNuclearStoppingModel::theZieglerFactor
private

Definition at line 121 of file G4hNuclearStoppingModel.hh.

Referenced by InitializeMe(), and TheValue().

The documentation for this class was generated from the following files:
Generated by doxygen 1.9.3