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__ &)
G4IonisParamElm &	operator= (const G4IonisParamElm &)=delete
G4bool	operator== (const G4IonisParamElm &) const =delete
G4bool	operator!= (const G4IonisParamElm &) const =delete
	G4IonisParamElm (G4IonisParamElm &)=delete

Detailed Description

Definition at line 50 of file G4IonisParamElm.hh.

Constructor & Destructor Documentation

G4IonisParamElm() [1/3]

G4IonisParamElm::G4IonisParamElm

(

G4double

Z

)

Definition at line 51 of file G4IonisParamElm.cc.

{
  G4int Z = G4lrint(AtomNumber);
  if (Z < 1) {
    G4Exception("G4IonisParamElm::G4IonisParamElm()",  "mat501", FatalException,
                "It is not allowed to create an Element with Z<1");
  }
  G4Pow* g4pow = G4Pow::GetInstance();
 
  // some basic functions of the atomic number
  fZ     = Z;
  fZ3    = g4pow->Z13(Z);
  fZZ3   = fZ3*g4pow->Z13(Z+1);
  flogZ3 = g4pow->logZ(Z)/3.;
   
  fMeanExcitationEnergy = 
    G4NistManager::Instance()->GetMeanIonisationEnergy(Z);
 
  // compute parameters for ion transport
  // The aproximation from:
  // J.F.Ziegler, J.P. Biersack, U. Littmark
  // The Stopping and Range of Ions in Matter,
  // Vol.1, Pergamon Press, 1985
  // Fast ions or hadrons
 
  G4int iz = Z - 1;
  if(91 < iz) { iz = 91; }
 
  static const G4double vFermi[92] = {
    1.0309,  0.15976, 0.59782, 1.0781,  1.0486,  1.0,     1.058,   0.93942, 0.74562, 0.3424,
    0.45259, 0.71074, 0.90519, 0.97411, 0.97184, 0.89852, 0.70827, 0.39816, 0.36552, 0.62712,
    0.81707, 0.9943,  1.1423,  1.2381,  1.1222,  0.92705, 1.0047,  1.2,     1.0661,  0.97411,
    0.84912, 0.95,    1.0903,  1.0429,  0.49715, 0.37755, 0.35211, 0.57801, 0.77773, 1.0207,
    1.029,   1.2542,  1.122,   1.1241,  1.0882,  1.2709,  1.2542,  0.90094, 0.74093, 0.86054,
    0.93155, 1.0047,  0.55379, 0.43289, 0.32636, 0.5131,  0.695,   0.72591, 0.71202, 0.67413,
    0.71418, 0.71453, 0.5911,  0.70263, 0.68049, 0.68203, 0.68121, 0.68532, 0.68715, 0.61884,
    0.71801, 0.83048, 1.1222,  1.2381,  1.045,   1.0733,  1.0953,  1.2381,  1.2879,  0.78654,
    0.66401, 0.84912, 0.88433, 0.80746, 0.43357, 0.41923, 0.43638, 0.51464, 0.73087, 0.81065,
    1.9578,  1.0257} ;
 
  static const G4double lFactor[92] = {
    1.0,  1.0,  1.1,  1.06, 1.01, 1.03, 1.04, 0.99, 0.95, 0.9,
    0.82, 0.81, 0.83, 0.88, 1.0,  0.95, 0.97, 0.99, 0.98, 0.97,
    0.98, 0.97, 0.96, 0.93, 0.91, 0.9,  0.88, 0.9,  0.9,  0.9,
    0.9,  0.85, 0.9,  0.9,  0.91, 0.92, 0.9,  0.9,  0.9,  0.9,
    0.9,  0.88, 0.9,  0.88, 0.88, 0.9,  0.9,  0.88, 0.9,  0.9,
    0.9,  0.9,  0.96, 1.2,  0.9,  0.88, 0.88, 0.85, 0.9,  0.9,
    0.92, 0.95, 0.99, 1.03, 1.05, 1.07, 1.08, 1.1,  1.08, 1.08,
    1.08, 1.08, 1.09, 1.09, 1.1,  1.11, 1.12, 1.13, 1.14, 1.15,
    1.17, 1.2,  1.18, 1.17, 1.17, 1.16, 1.16, 1.16, 1.16, 1.16,
    1.16, 1.16} ;
 
  fVFermi  = vFermi[iz];
  fLFactor = lFactor[iz];
 
  // obsolete parameters for ionisation
  fTau0 = 0.1*fZ3*MeV/proton_mass_c2;
  fTaul = 2.*MeV/proton_mass_c2;
 
  // compute the Bethe-Bloch formula for energy = fTaul*particle mass
  G4double rate = fMeanExcitationEnergy/electron_mass_c2 ;
  G4double w = fTaul*(fTaul+2.) ;
  fBetheBlochLow = (fTaul+1.)*(fTaul+1.)*std::log(2.*w/rate)/w - 1. ;
  fBetheBlochLow = 2.*fZ*twopi_mc2_rcl2*fBetheBlochLow ; 
  
  fClow = std::sqrt(fTaul)*fBetheBlochLow;
  fAlow = 6.458040 * fClow/fTau0;
  G4double Taum = 0.035*fZ3*MeV/proton_mass_c2;
  fBlow =-3.229020*fClow/(fTau0*std::sqrt(Taum));
 
  // Shell correction parameterization
  fShellCorrectionVector = new G4double[3]; //[3]
  rate = 0.001*fMeanExcitationEnergy/eV;
  G4double rate2 = rate*rate;
    /*    
    fShellCorrectionVector[0] = ( 1.10289e5 + 5.14781e8*rate)*rate2 ;
    fShellCorrectionVector[1] = ( 7.93805e3 - 2.22565e7*rate)*rate2 ;
    fShellCorrectionVector[2] = (-9.92256e1 + 2.10823e5*rate)*rate2 ;
    */
  fShellCorrectionVector[0] = ( 0.422377   + 3.858019*rate)*rate2 ;
  fShellCorrectionVector[1] = ( 0.0304043  - 0.1667989*rate)*rate2 ;
  fShellCorrectionVector[2] = (-0.00038106 + 0.00157955*rate)*rate2 ;
}

~G4IonisParamElm()

G4IonisParamElm::~G4IonisParamElm

(

)

Definition at line 149 of file G4IonisParamElm.cc.

{ delete[] fShellCorrectionVector; }

G4IonisParamElm() [2/3]

G4IonisParamElm::G4IonisParamElm

(

__void__ &

)

Definition at line 140 of file G4IonisParamElm.cc.

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

G4IonisParamElm() [3/3]

G4IonisParamElm::G4IonisParamElm ( G4IonisParamElm &  )

delete

Member Function Documentation

GetAlow()

G4double G4IonisParamElm::GetAlow ( ) const

inline

Definition at line 69 of file G4IonisParamElm.hh.

{return fAlow;}

GetBlow()

G4double G4IonisParamElm::GetBlow ( ) const

inline

Definition at line 70 of file G4IonisParamElm.hh.

{return fBlow;}

GetClow()

G4double G4IonisParamElm::GetClow ( ) const

inline

Definition at line 71 of file G4IonisParamElm.hh.

{return fClow;}

GetFermiVelocity()

G4double G4IonisParamElm::GetFermiVelocity ( ) const

inline

Definition at line 77 of file G4IonisParamElm.hh.

{return fVFermi;};

GetLFactor()

G4double G4IonisParamElm::GetLFactor ( ) const

inline

Definition at line 78 of file G4IonisParamElm.hh.

{return fLFactor;};

GetlogZ3()

G4double G4IonisParamElm::GetlogZ3 ( ) const

inline

Definition at line 62 of file G4IonisParamElm.hh.

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

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

GetMeanExcitationEnergy()

G4double G4IonisParamElm::GetMeanExcitationEnergy ( ) const

inline

Definition at line 74 of file G4IonisParamElm.hh.

{return fMeanExcitationEnergy;}

GetShellCorrectionVector()

G4double * G4IonisParamElm::GetShellCorrectionVector ( ) const

inline

Definition at line 80 of file G4IonisParamElm.hh.

{return fShellCorrectionVector;}

GetTau0()

G4double G4IonisParamElm::GetTau0 ( ) const

inline

Definition at line 64 of file G4IonisParamElm.hh.

{return fTau0;};

GetTaul()

G4double G4IonisParamElm::GetTaul ( ) const

inline

Definition at line 67 of file G4IonisParamElm.hh.

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

GetZ()

G4double G4IonisParamElm::GetZ ( ) const

inline

Definition at line 59 of file G4IonisParamElm.hh.

{return fZ;}       // atomic number Z

GetZ3()

G4double G4IonisParamElm::GetZ3 ( ) const

inline

Definition at line 60 of file G4IonisParamElm.hh.

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

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

GetZZ3()

G4double G4IonisParamElm::GetZZ3 ( ) const

inline

Definition at line 61 of file G4IonisParamElm.hh.

{return fZZ3;}     // std::pow (Z(Z+1),1/3)

operator!=()

G4bool G4IonisParamElm::operator!= ( const G4IonisParamElm &  ) const

delete

operator=()

G4IonisParamElm & G4IonisParamElm::operator= ( const G4IonisParamElm &  )

delete

operator==()

G4bool G4IonisParamElm::operator== ( const G4IonisParamElm &  ) const

delete

---

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

• G4IonisParamElm.hh
• G4IonisParamElm.cc