#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() [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 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::~G4IonisParamElm ( )

virtual

Definition at line 150 of file G4IonisParamElm.cc.

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

◆ G4IonisParamElm() [2/3]

G4IonisParamElm::G4IonisParamElm ( G4IonisParamElm & right )

Definition at line 157 of file G4IonisParamElm.cc.

{
  fShellCorrectionVector = 0;
  *this = right;
}

◆ G4IonisParamElm() [3/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;
}

Member Function Documentation

◆ GetAlow()

G4double G4IonisParamElm::GetAlow ( ) const

inline

Definition at line 70 of file G4IonisParamElm.hh.

70{return fAlow;}

◆ GetBlow()

G4double G4IonisParamElm::GetBlow ( ) const

inline

Definition at line 71 of file G4IonisParamElm.hh.

71{return fBlow;}

◆ GetClow()

G4double G4IonisParamElm::GetClow ( ) const

inline

Definition at line 72 of file G4IonisParamElm.hh.

72{return fClow;}

◆ GetFermiVelocity()

G4double G4IonisParamElm::GetFermiVelocity ( ) const

inline

Definition at line 78 of file G4IonisParamElm.hh.

78{return fVFermi;};

◆ GetLFactor()

G4double G4IonisParamElm::GetLFactor ( ) const

inline

Definition at line 79 of file G4IonisParamElm.hh.

79{return fLFactor;};

◆ GetlogZ3()

G4double G4IonisParamElm::GetlogZ3 ( ) const

inline

Definition at line 63 of file G4IonisParamElm.hh.

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

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

◆ GetMeanExcitationEnergy()

G4double G4IonisParamElm::GetMeanExcitationEnergy ( ) const

inline

Definition at line 75 of file G4IonisParamElm.hh.

75{return fMeanExcitationEnergy;}

◆ GetShellCorrectionVector()

G4double * G4IonisParamElm::GetShellCorrectionVector ( ) const

inline

Definition at line 81 of file G4IonisParamElm.hh.

81{return fShellCorrectionVector;}

◆ GetTau0()

G4double G4IonisParamElm::GetTau0 ( ) const

inline

Definition at line 65 of file G4IonisParamElm.hh.

65{return fTau0;};

◆ GetTaul()

G4double G4IonisParamElm::GetTaul ( ) const

inline

Definition at line 68 of file G4IonisParamElm.hh.

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

◆ GetZ()

G4double G4IonisParamElm::GetZ ( ) const

inline

Definition at line 60 of file G4IonisParamElm.hh.

60{return fZ;} // atomic number Z

◆ GetZ3()

G4double G4IonisParamElm::GetZ3 ( ) const

inline

Definition at line 61 of file G4IonisParamElm.hh.

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

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

◆ GetZZ3()

G4double G4IonisParamElm::GetZZ3 ( ) const

inline

Definition at line 62 of file G4IonisParamElm.hh.

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

Referenced by G4eBremsstrahlungModel::SampleSecondaries().

◆ operator!=()

G4int G4IonisParamElm::operator!= ( const G4IonisParamElm & right ) const

Definition at line 200 of file G4IonisParamElm.cc.

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

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

◆ operator==()

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

Definition at line 193 of file G4IonisParamElm.cc.

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

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

Public Member Functions

Detailed Description

Constructor & Destructor Documentation

◆ G4IonisParamElm() [1/3]

◆ ~G4IonisParamElm()

◆ G4IonisParamElm() [2/3]

◆ G4IonisParamElm() [3/3]

Member Function Documentation

◆ GetAlow()

◆ GetBlow()

◆ GetClow()

◆ GetFermiVelocity()

◆ GetLFactor()

◆ GetlogZ3()

◆ GetMeanExcitationEnergy()

◆ GetShellCorrectionVector()

◆ GetTau0()

◆ GetTaul()

◆ GetZ()

◆ GetZ3()

◆ GetZZ3()

◆ operator!=()

◆ operator=()

◆ operator==()