G4IonisParamElm Class Reference

#include <G4IonisParamElm.hh>

Collaboration diagram for G4IonisParamElm:

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 &	operator= (const G4IonisParamElm &)
G4int	operator== (const G4IonisParamElm &) const
G4int	operator!= (const G4IonisParamElm &) const
	G4IonisParamElm (__void__ &)

Private Member Functions
	G4IonisParamElm (G4IonisParamElm &)

Private Attributes
G4double	fZ
G4double	fZ3
G4double	fZZ3
G4double	flogZ3
G4double	fTau0
G4double	fTaul
G4double	fBetheBlochLow
G4double	fAlow
G4double	fBlow
G4double	fClow
G4double	fMeanExcitationEnergy
G4double *	fShellCorrectionVector
G4double	fVFermi
G4double	fLFactor

Detailed Description

Definition at line 51 of file G4IonisParamElm.hh.

Constructor & Destructor Documentation

G4IonisParamElm() [1/3]

G4IonisParamElm::G4IonisParamElm

(

G4double

Z

)

Definition at line 52 of file G4IonisParamElm.cc.

{
  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 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 ;
}

Here is the call graph for this function:

Here is the caller graph for this function:

~G4IonisParamElm()

G4IonisParamElm::~G4IonisParamElm ( )

virtual

Definition at line 150 of file G4IonisParamElm.cc.

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

G4IonisParamElm() [2/3]

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;
}

G4IonisParamElm() [3/3]

G4IonisParamElm::G4IonisParamElm ( G4IonisParamElm &  )

private

Member Function Documentation

GetAlow()

G4double G4IonisParamElm::GetAlow ( ) const

inline

Definition at line 70 of file G4IonisParamElm.hh.

{return fAlow;}

GetBlow()

G4double G4IonisParamElm::GetBlow ( ) const

inline

Definition at line 71 of file G4IonisParamElm.hh.

{return fBlow;}

GetClow()

G4double G4IonisParamElm::GetClow ( ) const

inline

Definition at line 72 of file G4IonisParamElm.hh.

{return fClow;}

GetFermiVelocity()

G4double G4IonisParamElm::GetFermiVelocity ( ) const

inline

Definition at line 78 of file G4IonisParamElm.hh.

{return fVFermi;};

Here is the caller graph for this function:

GetLFactor()

G4double G4IonisParamElm::GetLFactor ( ) const

inline

Definition at line 79 of file G4IonisParamElm.hh.

{return fLFactor;};

Here is the caller graph for this function:

GetlogZ3()

G4double G4IonisParamElm::GetlogZ3 ( ) const

inline

Definition at line 63 of file G4IonisParamElm.hh.

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

Here is the caller graph for this function:

GetMeanExcitationEnergy()

G4double G4IonisParamElm::GetMeanExcitationEnergy ( ) const

inline

Definition at line 75 of file G4IonisParamElm.hh.

{return fMeanExcitationEnergy;}

Here is the caller graph for this function:

GetShellCorrectionVector()

G4double* G4IonisParamElm::GetShellCorrectionVector ( ) const

inline

Definition at line 81 of file G4IonisParamElm.hh.

{return fShellCorrectionVector;}

Here is the call graph for this function:

GetTau0()

G4double G4IonisParamElm::GetTau0 ( ) const

inline

Definition at line 65 of file G4IonisParamElm.hh.

{return fTau0;};

GetTaul()

G4double G4IonisParamElm::GetTaul ( ) const

inline

Definition at line 68 of file G4IonisParamElm.hh.

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

GetZ()

G4double G4IonisParamElm::GetZ ( ) const

inline

Definition at line 60 of file G4IonisParamElm.hh.

{return fZ;}       // atomic number Z

GetZ3()

G4double G4IonisParamElm::GetZ3 ( ) const

inline

Definition at line 61 of file G4IonisParamElm.hh.

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

Here is the caller graph for this function:

GetZZ3()

G4double G4IonisParamElm::GetZZ3 ( ) const

inline

Definition at line 62 of file G4IonisParamElm.hh.

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

operator!=()

G4int G4IonisParamElm::operator!=

(

const G4IonisParamElm &

right

)

const

Definition at line 200 of file G4IonisParamElm.cc.

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

Here is the caller graph for this function:

operator=()

G4IonisParamElm & G4IonisParamElm::operator=

(

const G4IonisParamElm &

right

)

Definition at line 165 of file G4IonisParamElm.cc.

{
  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;
}

Here is the caller graph for this function:

operator==()

G4int G4IonisParamElm::operator==

(

const G4IonisParamElm &

right

)

const

Definition at line 193 of file G4IonisParamElm.cc.

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

Here is the caller graph for this function:

Member Data Documentation

fAlow

G4double G4IonisParamElm::fAlow

private

Definition at line 110 of file G4IonisParamElm.hh.

fBetheBlochLow

G4double G4IonisParamElm::fBetheBlochLow

private

Definition at line 109 of file G4IonisParamElm.hh.

fBlow

G4double G4IonisParamElm::fBlow

private

Definition at line 110 of file G4IonisParamElm.hh.

fClow

G4double G4IonisParamElm::fClow

private

Definition at line 110 of file G4IonisParamElm.hh.

fLFactor

G4double G4IonisParamElm::fLFactor

private

Definition at line 116 of file G4IonisParamElm.hh.

flogZ3

G4double G4IonisParamElm::flogZ3

private

Definition at line 103 of file G4IonisParamElm.hh.

fMeanExcitationEnergy

G4double G4IonisParamElm::fMeanExcitationEnergy

private

Definition at line 111 of file G4IonisParamElm.hh.

fShellCorrectionVector

G4double* G4IonisParamElm::fShellCorrectionVector

private

Definition at line 112 of file G4IonisParamElm.hh.

fTau0

G4double G4IonisParamElm::fTau0

private

Definition at line 107 of file G4IonisParamElm.hh.

fTaul

G4double G4IonisParamElm::fTaul

private

Definition at line 108 of file G4IonisParamElm.hh.

fVFermi

G4double G4IonisParamElm::fVFermi

private

Definition at line 115 of file G4IonisParamElm.hh.

fZ

G4double G4IonisParamElm::fZ

private

Definition at line 100 of file G4IonisParamElm.hh.

fZ3

G4double G4IonisParamElm::fZ3

private

Definition at line 101 of file G4IonisParamElm.hh.

fZZ3

G4double G4IonisParamElm::fZZ3

private

Definition at line 102 of file G4IonisParamElm.hh.

---

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

• G4IonisParamElm.hh
• G4IonisParamElm.cc