#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 python.hepunit::electron_mass_c2, python.hepunit::eV, FatalException, G4Exception(), G4Pow::GetInstance(), G4NistManager::GetMeanIonisationEnergy(), G4NistManager::Instance(), iz, G4Pow::logZ(), python.hepunit::MeV, python.hepunit::proton_mass_c2, python.hepunit::twopi_mc2_rcl2, 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 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 ( )

virtual

Definition at line 150 of file G4IonisParamElm.cc.

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

G4IonisParamElm::G4IonisParamElm ( G4IonisParamElm & right )

Definition at line 157 of file G4IonisParamElm.cc.

References right.

 {
   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.

70 {return fAlow;}

G4double G4IonisParamElm::GetBlow ( ) const

inline

Definition at line 71 of file G4IonisParamElm.hh.

71 {return fBlow;}

G4double G4IonisParamElm::GetClow ( ) const

inline

Definition at line 72 of file G4IonisParamElm.hh.

72 {return fClow;}

G4double G4IonisParamElm::GetFermiVelocity ( ) const

inline

Definition at line 78 of file G4IonisParamElm.hh.

78 {return fVFermi;};

G4double G4IonisParamElm::GetLFactor ( ) const

inline

Definition at line 79 of file G4IonisParamElm.hh.

79 {return fLFactor;};

G4double G4IonisParamElm::GetlogZ3 ( ) const

inline

Definition at line 63 of file G4IonisParamElm.hh.

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

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

G4double G4IonisParamElm::GetMeanExcitationEnergy ( ) const

inline

Definition at line 75 of file G4IonisParamElm.hh.

75 {return fMeanExcitationEnergy;}

G4double* G4IonisParamElm::GetShellCorrectionVector ( ) const

inline

Definition at line 81 of file G4IonisParamElm.hh.

81 {return fShellCorrectionVector;}

G4double G4IonisParamElm::GetTau0 ( ) const

inline

Definition at line 65 of file G4IonisParamElm.hh.

65 {return fTau0;};

G4double G4IonisParamElm::GetTaul ( ) const

inline

Definition at line 68 of file G4IonisParamElm.hh.

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

G4double G4IonisParamElm::GetZ ( ) const

inline

Definition at line 60 of file G4IonisParamElm.hh.

60 {return fZ;} // atomic number Z

G4double G4IonisParamElm::GetZ3 ( ) const

inline

Definition at line 61 of file G4IonisParamElm.hh.

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

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

G4double G4IonisParamElm::GetZZ3 ( ) const

inline

Definition at line 62 of file G4IonisParamElm.hh.

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

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

Definition at line 200 of file G4IonisParamElm.cc.

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

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

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

Member Function Documentation