G4IonisParamElm Class Reference

#include <G4IonisParamElm.hh>

Public Member Functions
	G4IonisParamElm (G4double Z)
virtual	~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 (G4IonisParamElm &)
G4IonisParamElm &	operator= (const G4IonisParamElm &)
G4int	operator== (const G4IonisParamElm &) const
G4int	operator!= (const G4IonisParamElm &) const
	G4IonisParamElm (__void__ &)

---

Detailed Description

Definition at line 51 of file G4IonisParamElm.hh.

---

Constructor & Destructor Documentation

G4IonisParamElm::G4IonisParamElm

(

G4double

Z

)

Definition at line 52 of file G4IonisParamElm.cc.

References FatalException, G4Exception(), G4Pow::GetInstance(), G4NistManager::GetMeanIonisationEnergy(), G4NistManager::Instance(), G4Pow::logZ(), and G4Pow::Z13().

{
  G4int Z = G4int(AtomNumber + 0.5);
  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 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 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

(

)

[virtual]

Definition at line 150 of file G4IonisParamElm.cc.

{
  if (fShellCorrectionVector) { delete [] fShellCorrectionVector; }
}

G4IonisParamElm::G4IonisParamElm

(

G4IonisParamElm &

)

Definition at line 157 of file G4IonisParamElm.cc.

{
  fShellCorrectionVector = 0;
  *this = right;
}

G4IonisParamElm::G4IonisParamElm

(

__void__ &

)

Definition at line 141 of file G4IonisParamElm.cc.

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

---

Member Function Documentation

G4double G4IonisParamElm::GetAlow

(

)

const [inline]

Definition at line 70 of file G4IonisParamElm.hh.

{return fAlow;}

G4double G4IonisParamElm::GetBlow

(

)

const [inline]

Definition at line 71 of file G4IonisParamElm.hh.

{return fBlow;}

G4double G4IonisParamElm::GetClow

(

)

const [inline]

Definition at line 72 of file G4IonisParamElm.hh.

{return fClow;}

G4double G4IonisParamElm::GetFermiVelocity

(

)

const [inline]

Definition at line 78 of file G4IonisParamElm.hh.

{return fVFermi;};

G4double G4IonisParamElm::GetLFactor

(

)

const [inline]

Definition at line 79 of file G4IonisParamElm.hh.

{return fLFactor;};

G4double G4IonisParamElm::GetlogZ3

(

)

const [inline]

Definition at line 63 of file G4IonisParamElm.hh.

Referenced by G4PairProductionRelModel::SampleSecondaries(), G4LivermorePolarizedGammaConversionModel::SampleSecondaries(), G4LivermoreNuclearGammaConversionModel::SampleSecondaries(), G4LivermoreGammaConversionModelRC::SampleSecondaries(), G4LivermoreGammaConversionModel::SampleSecondaries(), G4eBremsstrahlungModel::SampleSecondaries(), and G4BetheHeitlerModel::SampleSecondaries().

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

G4double G4IonisParamElm::GetMeanExcitationEnergy

(

)

const [inline]

Definition at line 75 of file G4IonisParamElm.hh.

{return fMeanExcitationEnergy;}

G4double* G4IonisParamElm::GetShellCorrectionVector

(

)

const [inline]

Definition at line 81 of file G4IonisParamElm.hh.

{return fShellCorrectionVector;}

G4double G4IonisParamElm::GetTau0

(

)

const [inline]

Definition at line 65 of file G4IonisParamElm.hh.

{return fTau0;};

G4double G4IonisParamElm::GetTaul

(

)

const [inline]

Definition at line 68 of file G4IonisParamElm.hh.

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

G4double G4IonisParamElm::GetZ

(

)

const [inline]

Definition at line 60 of file G4IonisParamElm.hh.

{return fZ;}       // atomic number Z

G4double G4IonisParamElm::GetZ3

(

)

const [inline]

Definition at line 61 of file G4IonisParamElm.hh.

Referenced by G4PairProductionRelModel::SampleSecondaries(), G4LivermorePolarizedGammaConversionModel::SampleSecondaries(), G4LivermoreNuclearGammaConversionModel::SampleSecondaries(), G4LivermoreGammaConversionModelRC::SampleSecondaries(), G4LivermoreGammaConversionModel::SampleSecondaries(), G4eBremsstrahlungModel::SampleSecondaries(), and G4BetheHeitlerModel::SampleSecondaries().

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

G4double G4IonisParamElm::GetZZ3

(

)

const [inline]

Definition at line 62 of file G4IonisParamElm.hh.

Referenced by G4eBremsstrahlungModel::SampleSecondaries().

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

G4int G4IonisParamElm::operator!=

(

const G4IonisParamElm &

)

const

Definition at line 200 of file G4IonisParamElm.cc.

{
  return (this != (G4IonisParamElm *) &right);
}

G4IonisParamElm & G4IonisParamElm::operator=

(

const G4IonisParamElm &

)

Definition at line 165 of file G4IonisParamElm.cc.

References fAlow, fBetheBlochLow, fBlow, fClow, fLFactor, flogZ3, fMeanExcitationEnergy, fShellCorrectionVector, fTau0, fTaul, fVFermi, fZ, fZ3, and fZZ3.

{
  if (this != &right)
    {
      fZ                     = right.fZ;
      fZ3                    = right.fZ3;
      fZZ3                   = right.fZZ3;
      flogZ3                 = right.flogZ3;
      fTau0                  = right.fTau0;
      fTaul                  = right.fTaul;
      fBetheBlochLow         = right.fBetheBlochLow;
      fAlow                  = right.fAlow;
      fBlow                  = right.fBlow;
      fClow                  = right.fClow;
      fMeanExcitationEnergy  = right.fMeanExcitationEnergy;
      fVFermi                = right.fVFermi;
      fLFactor               = right.fLFactor;
      if (fShellCorrectionVector) { delete [] fShellCorrectionVector; } 
      fShellCorrectionVector = new G4double[3];            
      fShellCorrectionVector[0] = right.fShellCorrectionVector[0];
      fShellCorrectionVector[1] = right.fShellCorrectionVector[1];
      fShellCorrectionVector[2] = right.fShellCorrectionVector[2];      
    } 
  return *this;
}

G4int G4IonisParamElm::operator==

(

const G4IonisParamElm &

)

const

Definition at line 193 of file G4IonisParamElm.cc.

{
  return (this == (G4IonisParamElm *) &right);
}

---

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

• G4IonisParamElm.hh
• G4IonisParamElm.cc

---