Geant4 10.7.0
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
Public Member Functions | List of all members
G4IonisParamElm Class Reference

#include <G4IonisParamElm.hh>

Public Member Functions

 G4IonisParamElm (G4double Z)
 
 ~G4IonisParamElm ()
 
G4double GetZ () const
 
G4double GetZ3 () const
 
G4double GetZZ3 () const
 
G4double GetlogZ3 () const
 
G4double GetTau0 () const
 
G4double GetTaul () const
 
G4double GetAlow () const
 
G4double GetBlow () const
 
G4double GetClow () const
 
G4double GetMeanExcitationEnergy () const
 
G4double GetFermiVelocity () const
 
G4double GetLFactor () const
 
G4double * GetShellCorrectionVector () const
 
 G4IonisParamElm (__void__ &)
 

Detailed Description

Definition at line 50 of file G4IonisParamElm.hh.

Constructor & Destructor Documentation

◆ G4IonisParamElm() [1/2]

G4IonisParamElm::G4IonisParamElm ( G4double  Z)

Definition at line 51 of file G4IonisParamElm.cc.

52{
53 G4int Z = G4lrint(AtomNumber);
54 if (Z < 1) {
55 G4Exception("G4IonisParamElm::G4IonisParamElm()", "mat501", FatalException,
56 "It is not allowed to create an Element with Z<1");
57 }
58 G4Pow* g4pow = G4Pow::GetInstance();
59
60 // some basic functions of the atomic number
61 fZ = Z;
62 fZ3 = g4pow->Z13(Z);
63 fZZ3 = fZ3*g4pow->Z13(Z+1);
64 flogZ3 = g4pow->logZ(Z)/3.;
65
66 fMeanExcitationEnergy =
67 G4NistManager::Instance()->GetMeanIonisationEnergy(Z);
68
69 // compute parameters for ion transport
70 // The aproximation from:
71 // J.F.Ziegler, J.P. Biersack, U. Littmark
72 // The Stopping and Range of Ions in Matter,
73 // Vol.1, Pergamon Press, 1985
74 // Fast ions or hadrons
75
76 G4int iz = Z - 1;
77 if(91 < iz) { iz = 91; }
78
79 static const G4double vFermi[92] = {
80 1.0309, 0.15976, 0.59782, 1.0781, 1.0486, 1.0, 1.058, 0.93942, 0.74562, 0.3424,
81 0.45259, 0.71074, 0.90519, 0.97411, 0.97184, 0.89852, 0.70827, 0.39816, 0.36552, 0.62712,
82 0.81707, 0.9943, 1.1423, 1.2381, 1.1222, 0.92705, 1.0047, 1.2, 1.0661, 0.97411,
83 0.84912, 0.95, 1.0903, 1.0429, 0.49715, 0.37755, 0.35211, 0.57801, 0.77773, 1.0207,
84 1.029, 1.2542, 1.122, 1.1241, 1.0882, 1.2709, 1.2542, 0.90094, 0.74093, 0.86054,
85 0.93155, 1.0047, 0.55379, 0.43289, 0.32636, 0.5131, 0.695, 0.72591, 0.71202, 0.67413,
86 0.71418, 0.71453, 0.5911, 0.70263, 0.68049, 0.68203, 0.68121, 0.68532, 0.68715, 0.61884,
87 0.71801, 0.83048, 1.1222, 1.2381, 1.045, 1.0733, 1.0953, 1.2381, 1.2879, 0.78654,
88 0.66401, 0.84912, 0.88433, 0.80746, 0.43357, 0.41923, 0.43638, 0.51464, 0.73087, 0.81065,
89 1.9578, 1.0257} ;
90
91 static const G4double lFactor[92] = {
92 1.0, 1.0, 1.1, 1.06, 1.01, 1.03, 1.04, 0.99, 0.95, 0.9,
93 0.82, 0.81, 0.83, 0.88, 1.0, 0.95, 0.97, 0.99, 0.98, 0.97,
94 0.98, 0.97, 0.96, 0.93, 0.91, 0.9, 0.88, 0.9, 0.9, 0.9,
95 0.9, 0.85, 0.9, 0.9, 0.91, 0.92, 0.9, 0.9, 0.9, 0.9,
96 0.9, 0.88, 0.9, 0.88, 0.88, 0.9, 0.9, 0.88, 0.9, 0.9,
97 0.9, 0.9, 0.96, 1.2, 0.9, 0.88, 0.88, 0.85, 0.9, 0.9,
98 0.92, 0.95, 0.99, 1.03, 1.05, 1.07, 1.08, 1.1, 1.08, 1.08,
99 1.08, 1.08, 1.09, 1.09, 1.1, 1.11, 1.12, 1.13, 1.14, 1.15,
100 1.17, 1.2, 1.18, 1.17, 1.17, 1.16, 1.16, 1.16, 1.16, 1.16,
101 1.16, 1.16} ;
102
103 fVFermi = vFermi[iz];
104 fLFactor = lFactor[iz];
105
106 // obsolete parameters for ionisation
107 fTau0 = 0.1*fZ3*MeV/proton_mass_c2;
108 fTaul = 2.*MeV/proton_mass_c2;
109
110 // compute the Bethe-Bloch formula for energy = fTaul*particle mass
111 G4double rate = fMeanExcitationEnergy/electron_mass_c2 ;
112 G4double w = fTaul*(fTaul+2.) ;
113 fBetheBlochLow = (fTaul+1.)*(fTaul+1.)*std::log(2.*w/rate)/w - 1. ;
114 fBetheBlochLow = 2.*fZ*twopi_mc2_rcl2*fBetheBlochLow ;
115
116 fClow = std::sqrt(fTaul)*fBetheBlochLow;
117 fAlow = 6.458040 * fClow/fTau0;
118 G4double Taum = 0.035*fZ3*MeV/proton_mass_c2;
119 fBlow =-3.229020*fClow/(fTau0*std::sqrt(Taum));
120
121 // Shell correction parameterization
122 fShellCorrectionVector = new G4double[3]; //[3]
123 rate = 0.001*fMeanExcitationEnergy/eV;
124 G4double rate2 = rate*rate;
125 /*
126 fShellCorrectionVector[0] = ( 1.10289e5 + 5.14781e8*rate)*rate2 ;
127 fShellCorrectionVector[1] = ( 7.93805e3 - 2.22565e7*rate)*rate2 ;
128 fShellCorrectionVector[2] = (-9.92256e1 + 2.10823e5*rate)*rate2 ;
129 */
130 fShellCorrectionVector[0] = ( 0.422377 + 3.858019*rate)*rate2 ;
131 fShellCorrectionVector[1] = ( 0.0304043 - 0.1667989*rate)*rate2 ;
132 fShellCorrectionVector[2] = (-0.00038106 + 0.00157955*rate)*rate2 ;
133}
FatalException
@ FatalException
Definition: G4ExceptionSeverity.hh:61
G4Exception
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
Definition: G4Exception.cc:35
G4double
double G4double
Definition: G4Types.hh:83
G4int
int G4int
Definition: G4Types.hh:85
G4NistManager::GetMeanIonisationEnergy
G4double GetMeanIonisationEnergy(G4int Z) const
Definition: G4NistManager.hh:446
G4NistManager::Instance
static G4NistManager * Instance()
Definition: G4NistManager.cc:70
G4Pow
Definition: G4Pow.hh:49
G4Pow::GetInstance
static G4Pow * GetInstance()
Definition: G4Pow.cc:41
G4Pow::logZ
G4double logZ(G4int Z) const
Definition: G4Pow.hh:137
G4Pow::Z13
G4double Z13(G4int Z) const
Definition: G4Pow.hh:123
G4lrint
int G4lrint(double ad)
Definition: templates.hh:134

◆ ~G4IonisParamElm()

G4IonisParamElm::~G4IonisParamElm ( )

Definition at line 149 of file G4IonisParamElm.cc.

150{
151 if (fShellCorrectionVector) { delete [] fShellCorrectionVector; }
152}

◆ G4IonisParamElm() [2/2]

G4IonisParamElm::G4IonisParamElm ( __void__ &  )

Definition at line 140 of file G4IonisParamElm.cc.

141 : fShellCorrectionVector(nullptr)
142{
143 fZ=fZ3=fZZ3=flogZ3=fTau0=fTaul=fBetheBlochLow=fAlow=fBlow=fClow
144 =fMeanExcitationEnergy=fVFermi=fLFactor=0.0;
145}

Member Function Documentation

◆ GetAlow()

G4double G4IonisParamElm::GetAlow ( ) const
inline

Definition at line 69 of file G4IonisParamElm.hh.

69{return fAlow;}

◆ GetBlow()

G4double G4IonisParamElm::GetBlow ( ) const
inline

Definition at line 70 of file G4IonisParamElm.hh.

70{return fBlow;}

◆ GetClow()

G4double G4IonisParamElm::GetClow ( ) const
inline

Definition at line 71 of file G4IonisParamElm.hh.

71{return fClow;}

◆ GetFermiVelocity()

G4double G4IonisParamElm::GetFermiVelocity ( ) const
inline

Definition at line 77 of file G4IonisParamElm.hh.

77{return fVFermi;};

◆ GetLFactor()

G4double G4IonisParamElm::GetLFactor ( ) const
inline

Definition at line 78 of file G4IonisParamElm.hh.

78{return fLFactor;};

◆ GetlogZ3()

G4double G4IonisParamElm::GetlogZ3 ( ) const
inline

Definition at line 62 of file G4IonisParamElm.hh.

62{return flogZ3;} // std::log(Z)/3

Referenced by G4BetheHeitlerModel::InitialiseElementData(), G4LivermoreGammaConversionModelRC::SampleSecondaries(), G4LivermoreNuclearGammaConversionModel::SampleSecondaries(), G4LivermorePolarizedGammaConversionModel::SampleSecondaries(), and G4BetheHeitlerModel::SampleSecondaries().

◆ GetMeanExcitationEnergy()

G4double G4IonisParamElm::GetMeanExcitationEnergy ( ) const
inline

Definition at line 74 of file G4IonisParamElm.hh.

74{return fMeanExcitationEnergy;}

◆ GetShellCorrectionVector()

G4double * G4IonisParamElm::GetShellCorrectionVector ( ) const
inline

Definition at line 80 of file G4IonisParamElm.hh.

80{return fShellCorrectionVector;}

◆ GetTau0()

G4double G4IonisParamElm::GetTau0 ( ) const
inline

Definition at line 64 of file G4IonisParamElm.hh.

64{return fTau0;};

◆ GetTaul()

G4double G4IonisParamElm::GetTaul ( ) const
inline

Definition at line 67 of file G4IonisParamElm.hh.

67{return fTaul;} // 2*MeV/proton mass

◆ GetZ()

G4double G4IonisParamElm::GetZ ( ) const
inline

Definition at line 59 of file G4IonisParamElm.hh.

59{return fZ;} // atomic number Z

◆ GetZ3()

G4double G4IonisParamElm::GetZ3 ( ) const
inline

Definition at line 60 of file G4IonisParamElm.hh.

60{return fZ3;} // std::pow (Z,1/3)

Referenced by G4LivermoreGammaConversionModelRC::SampleSecondaries(), G4LivermoreNuclearGammaConversionModel::SampleSecondaries(), G4LivermorePolarizedGammaConversionModel::SampleSecondaries(), G4BetheHeitlerModel::SampleSecondaries(), and G4BetheHeitler5DModel::SampleSecondaries().

◆ GetZZ3()

G4double G4IonisParamElm::GetZZ3 ( ) const
inline

Definition at line 61 of file G4IonisParamElm.hh.

61{return fZZ3;} // std::pow (Z(Z+1),1/3)
The documentation for this class was generated from the following files: