#include <G4IonisParamMat.hh>

Collaboration diagram for G4IonisParamMat:

Public Member Functions
	G4IonisParamMat (G4Material *)

virtual	~G4IonisParamMat ()

G4double	GetMeanExcitationEnergy () const

void	SetMeanExcitationEnergy (G4double value)

G4double	FindMeanExcitationEnergy (const G4String &chFormula)

G4double	GetLogMeanExcEnergy () const

G4double *	GetShellCorrectionVector () const

G4double	GetTaul () const

G4double	GetPlasmaEnergy () const

G4double	GetAdjustmentFactor () const

G4double	GetCdensity () const

G4double	GetMdensity () const

G4double	GetAdensity () const

G4double	GetX0density () const

G4double	GetX1density () const

G4double	GetD0density () const

G4double	DensityCorrection (G4double x)

G4double	GetF1fluct () const

G4double	GetF2fluct () const

G4double	GetEnergy1fluct () const

G4double	GetLogEnergy1fluct () const

G4double	GetEnergy2fluct () const

G4double	GetLogEnergy2fluct () const

G4double	GetEnergy0fluct () const

G4double	GetRateionexcfluct () const

G4double	GetZeffective () const

G4double	GetFermiEnergy () const

G4double	GetLFactor () const

G4double	GetInvA23 () const

void	SetBirksConstant (G4double value)

G4double	GetBirksConstant () const

void	SetMeanEnergyPerIonPair (G4double value)

G4double	GetMeanEnergyPerIonPair () const

	G4IonisParamMat (__void__ &)

Static Public Member Functions
static G4DensityEffectData *	GetDensityEffectData ()

Private Member Functions
void	ComputeMeanParameters ()

void	ComputeDensityEffect ()

void	ComputeFluctModel ()

void	ComputeIonParameters ()

G4IonisParamMat &	operator= (const G4IonisParamMat &)

G4int	operator== (const G4IonisParamMat &) const

G4int	operator!= (const G4IonisParamMat &) const

	G4IonisParamMat (const G4IonisParamMat &)

Static Private Attributes
static G4DensityEffectData *	fDensityData = 0

Detailed Description

Definition at line 58 of file G4IonisParamMat.hh.

Constructor & Destructor Documentation

◆ G4IonisParamMat() [1/3]

G4IonisParamMat::G4IonisParamMat ( G4Material * material )

Definition at line 56 of file G4IonisParamMat.cc.

   : fMaterial(material)
 {
   fBirks = 0.;
   fMeanEnergyPerIon = 0.0;
   twoln10 = 2.*G4Pow::GetInstance()->logZ(10);
 
   // minimal set of default parameters for density effect
   fCdensity = 0.0;
   fD0density = 0.0;
   fAdjustmentFactor = 1.0;
   if(!fDensityData) { fDensityData = new G4DensityEffectData(); }
 
   // compute parameters
   ComputeMeanParameters();
   ComputeDensityEffect();
   ComputeFluctModel();
   ComputeIonParameters();
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ ~G4IonisParamMat()

G4IonisParamMat::~G4IonisParamMat ( )

virtual

Definition at line 116 of file G4IonisParamMat.cc.

 {
   if (fShellCorrectionVector) { delete [] fShellCorrectionVector; }
   if (fDensityData) { delete fDensityData; }
   fDensityData = 0;
   fShellCorrectionVector = 0;
 }

◆ G4IonisParamMat() [2/3]

G4IonisParamMat::G4IonisParamMat ( __void__ & )

Definition at line 81 of file G4IonisParamMat.cc.

   : fMaterial(0), fShellCorrectionVector(0)
 {
   fMeanExcitationEnergy = 0.0;
   fLogMeanExcEnergy = 0.0;
   fTaul = 0.0;
   fCdensity = 0.0;
   fMdensity = 0.0;
   fAdensity = 0.0;
   fX0density = 0.0;
   fX1density = 0.0;
   fD0density = 0.0;
   fPlasmaEnergy = 0.0;
   fAdjustmentFactor = 0.0;
   fF1fluct = 0.0;          
   fF2fluct = 0.0;                       
   fEnergy1fluct = 0.0;
   fLogEnergy1fluct = 0.0;
   fEnergy2fluct = 0.0;
   fLogEnergy2fluct = 0.0;
   fEnergy0fluct = 0.0;
   fRateionexcfluct = 0.0;
   fZeff = 0.0;
   fFermiEnergy = 0.0;
   fLfactor = 0.0;
   fInvA23 = 0.0;
   fBirks = 0.0;
   fMeanEnergyPerIon = 0.0;
   twoln10 = 2.*G4Pow::GetInstance()->logZ(10);
 
   if(!fDensityData) { fDensityData = new G4DensityEffectData(); }
 }

Here is the call graph for this function:

◆ G4IonisParamMat() [3/3]

G4IonisParamMat::G4IonisParamMat ( const G4IonisParamMat & )

private

Member Function Documentation

◆ ComputeDensityEffect()

void G4IonisParamMat::ComputeDensityEffect ( )

private

Definition at line 180 of file G4IonisParamMat.cc.

 {
   G4State State = fMaterial->GetState();
 
   // Check if density effect data exist in the table
   // R.M. Sternheimer, Atomic Data and Nuclear Data Tables, 30: 261 (1984)
   G4int idx = fDensityData->GetIndex(fMaterial->GetName());
   G4int nelm= fMaterial->GetNumberOfElements();
   G4int Z0  = G4lrint((*(fMaterial->GetElementVector()))[0]->GetZ());
 
   // for simple non-NIST materials 
   G4double corr = 0.0;
   if(idx < 0 && 1 == nelm) {
     idx = fDensityData->GetElementIndex(Z0, fMaterial->GetState());
 
     // Correction for base material or for non-nominal density
     // Except cases of very different density defined in user code
     if(idx >= 0) {
       const G4Material* bmat = fMaterial->GetBaseMaterial();
       if(bmat) {
     corr = G4Log(bmat->GetDensity()/fMaterial->GetDensity());
       } else {
     G4double dens = G4NistManager::Instance()->GetNominalDensity(Z0);
     if(dens <= 0.0) { idx = -1; }
     else {
       corr = G4Log(dens/fMaterial->GetDensity());
     }
       }
       // 1.0 is an arbitrary empirical limit 
       // parameterisation with very different density is not applicable
       if(std::abs(corr) > 1.0) { idx = -1; }
     }
   }
 
   //G4cout<<"DensityEffect for "<<fMaterial->GetName()<<"  "<< idx << G4endl; 
 
   if(idx >= 0) {
 
     // Take parameters for the density effect correction from
     // R.M. Sternheimer et al. Density Effect For The Ionization Loss 
     // of Charged Particles in Various Substances. 
     // Atom. Data Nucl. Data Tabl. 30 (1984) 261-271. 
 
     fCdensity = fDensityData->GetCdensity(idx); 
     //G4double Cdensity = fCdensity; 
     fMdensity = fDensityData->GetMdensity(idx);
     fAdensity = fDensityData->GetAdensity(idx);
     fX0density = fDensityData->GetX0density(idx);
     fX1density = fDensityData->GetX1density(idx);
     fD0density = fDensityData->GetDelta0density(idx);
     fPlasmaEnergy = fDensityData->GetPlasmaEnergy(idx);
     fAdjustmentFactor = fDensityData->GetAdjustmentFactor(idx);
 
     // parameter C is computed and not taken from Sternheimer tables
     //fCdensity = 1. + 2*G4Log(fMeanExcitationEnergy/fPlasmaEnergy);
     //G4cout << "IonisParamMat: " << fMaterial->GetName() 
     //     << "  Cst= " << Cdensity << " C= " << fCdensity << G4endl;
 
     // correction on nominal density
     fCdensity  += corr;
     fX0density += corr/twoln10;
     fX1density += corr/twoln10;
 
   } else {
 
     static const G4double Cd2 = 4*pi*hbarc_squared*classic_electr_radius;
     fPlasmaEnergy = std::sqrt(Cd2*fMaterial->GetTotNbOfElectPerVolume());
 
     // Compute parameters for the density effect correction in DE/Dx formula.
     // The parametrization is from R.M. Sternheimer, Phys. Rev.B,3:3681 (1971)
     G4int icase;
     
     fCdensity = 1. + 2*G4Log(fMeanExcitationEnergy/fPlasmaEnergy);
     //
     // condensed materials
     //  
     if ((State == kStateSolid)||(State == kStateLiquid)) {
 
       static const G4double E100eV  = 100.*eV; 
       static const G4double ClimiS[] = {3.681 , 5.215 };
       static const G4double X0valS[] = {1.0   , 1.5   };
       static const G4double X1valS[] = {2.0   , 3.0   };
                                 
       if(fMeanExcitationEnergy < E100eV) { icase = 0; }
       else                               { icase = 1; } 
 
       if(fCdensity < ClimiS[icase])    { fX0density = 0.2; }
       else { fX0density = 0.326*fCdensity - X0valS[icase]; }
 
       fX1density = X1valS[icase]; fMdensity = 3.0;
       
       //special: Hydrogen
       if (1 == nelm && 1 == Z0) {
          fX0density = 0.425; fX1density = 2.0; fMdensity = 5.949;
       }
     }
     //
     // gases
     //
     if (State == kStateGas) { 
 
       fMdensity = 3.;
       fX1density = 4.0;
       //static const G4double ClimiG[] = {10.,10.5,11.,11.5,12.25,13.804};
       //static const G4double X0valG[] = {1.6,1.7,1.8,1.9,2.0,2.0};
       //static const G4double X1valG[] = {4.0,4.0,4.0,4.0,4.0,5.0};
 
       if(fCdensity < 10.) {
     fX0density = 1.6; 
       } else if(fCdensity < 11.5) { 
     fX0density = 1.6 + 0.2*(fCdensity - 10.); 
       } else if(fCdensity < 12.25) { 
     fX0density = 1.9 + (fCdensity - 11.5)/7.5; 
       } else if(fCdensity < 13.804) { 
     fX0density = 2.0; 
     fX1density = 4.0 + (fCdensity - 12.25)/1.554;
       } else {
     fX0density = 0.326*fCdensity-2.5; fX1density = 5.0; 
       }
       
       //special: Hydrogen
       if (1 == nelm && 1 == Z0) {
          fX0density = 1.837; fX1density = 3.0; fMdensity = 4.754;
       }
       
       //special: Helium
       if (1 == nelm && 2 == Z0) {
          fX0density = 2.191; fX1density = 3.0; fMdensity = 3.297;
       }
     }
   }
 
   // change parameters if the gas is not in STP.
   // For the correction the density(STP) is needed. 
   // Density(STP) is calculated here : 
   
     
   if (State == kStateGas) { 
     G4double Density  = fMaterial->GetDensity();
     G4double Pressure = fMaterial->GetPressure();
     G4double Temp     = fMaterial->GetTemperature();
       
     G4double DensitySTP = Density*STP_Pressure*Temp/(Pressure*NTP_Temperature);
 
     G4double ParCorr = G4Log(Density/DensitySTP);
   
     fCdensity  -= ParCorr;
     fX0density -= ParCorr/twoln10;
     fX1density -= ParCorr/twoln10;
   }
 
   // fAdensity parameter can be fixed for not conductive materials 
   if(0.0 == fD0density) {
     G4double Xa = fCdensity/twoln10;
     fAdensity = twoln10*(Xa-fX0density)
       /std::pow((fX1density-fX0density),fMdensity);
   }
   /*  
   G4cout << "G4IonisParamMat: density effect data for <" 
          << fMaterial->GetName() 
      << "> " << G4endl;
   G4cout << "Eplasma(eV)= " << fPlasmaEnergy/eV
      << " rho= " << fAdjustmentFactor
      << " -C= " << fCdensity 
      << " x0= " << fX0density
      << " x1= " << fX1density
      << " a= " << fAdensity
      << " m= " << fMdensity
      << G4endl;
   */
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ ComputeFluctModel()

void G4IonisParamMat::ComputeFluctModel ( )

private

Definition at line 354 of file G4IonisParamMat.cc.

 {
   // compute parameters for the energy loss fluctuation model
   // needs an 'effective Z' 
   G4double Zeff = 0.;
   for (size_t i=0;i<fMaterial->GetNumberOfElements();i++) {
      Zeff += (fMaterial->GetFractionVector())[i]
              *((*(fMaterial->GetElementVector()))[i]->GetZ());
   }
   if (Zeff > 2.) { fF2fluct = 2./Zeff; }
   else           { fF2fluct = 0.; }
 
   fF1fluct         = 1. - fF2fluct;
   fEnergy2fluct    = 10.*Zeff*Zeff*eV;
   fLogEnergy2fluct = G4Log(fEnergy2fluct);
   fLogEnergy1fluct = (fLogMeanExcEnergy - fF2fluct*fLogEnergy2fluct)
                      /fF1fluct;
   fEnergy1fluct    = G4Exp(fLogEnergy1fluct);
   fEnergy0fluct    = 10.*eV;
   fRateionexcfluct = 0.4;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ ComputeIonParameters()

void G4IonisParamMat::ComputeIonParameters ( )

private

Definition at line 378 of file G4IonisParamMat.cc.

 {
   // get elements in the actual material,
   const G4ElementVector* theElementVector = fMaterial->GetElementVector() ;
   const G4double* theAtomicNumDensityVector =
                          fMaterial->GetAtomicNumDensityVector() ;
   const G4int NumberOfElements = fMaterial->GetNumberOfElements() ;
 
   //  loop for the elements in the material
   //  to find out average values Z, vF, lF
   G4double z(0.0), vF(0.0), lF(0.0), norm(0.0), a23(0.0);
 
   G4Pow* g4pow = G4Pow::GetInstance();
   if( 1 == NumberOfElements ) {
     const G4Element* element = (*theElementVector)[0];
     z = element->GetZ();
     vF= element->GetIonisation()->GetFermiVelocity();
     lF= element->GetIonisation()->GetLFactor();
     a23 = 1.0/g4pow->A23(element->GetN());
 
   } else {
     for (G4int iel=0; iel<NumberOfElements; iel++)
       {
         const G4Element* element = (*theElementVector)[iel];
         const G4double weight = theAtomicNumDensityVector[iel];
         norm += weight ;
         z    += element->GetZ() * weight;
         vF   += element->GetIonisation()->GetFermiVelocity() * weight;
         lF   += element->GetIonisation()->GetLFactor() * weight;
     a23  += weight/g4pow->A23(element->GetN());
       }
     z  /= norm;
     vF /= norm;
     lF /= norm;
     a23 /= norm;
   }  
   fZeff        = z;
   fLfactor     = lF;
   fFermiEnergy = 25.*keV*vF*vF;
   fInvA23      = a23;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ ComputeMeanParameters()

void G4IonisParamMat::ComputeMeanParameters ( )

private

Definition at line 126 of file G4IonisParamMat.cc.

 {
   // compute mean excitation energy and shell correction vector
   fTaul = (*(fMaterial->GetElementVector()))[0]->GetIonisation()->GetTaul();
 
   fMeanExcitationEnergy = 0.;
   fLogMeanExcEnergy = 0.;
 
   size_t nElements = fMaterial->GetNumberOfElements();
   const G4ElementVector* elmVector = fMaterial->GetElementVector();
   const G4double* nAtomsPerVolume = fMaterial->GetVecNbOfAtomsPerVolume();
  
   G4String ch = fMaterial->GetChemicalFormula();
 
   if(ch != "") { fMeanExcitationEnergy = FindMeanExcitationEnergy(ch); }
 
   // Chemical formula defines mean excitation energy
   if(fMeanExcitationEnergy > 0.0) {
     fLogMeanExcEnergy = G4Log(fMeanExcitationEnergy);
 
     // Compute average 
   } else {
     for (size_t i=0; i < nElements; i++) {
       const G4Element* elm = (*elmVector)[i];
       fLogMeanExcEnergy += nAtomsPerVolume[i]*elm->GetZ()
     *G4Log(elm->GetIonisation()->GetMeanExcitationEnergy());
     }
     fLogMeanExcEnergy /= fMaterial->GetTotNbOfElectPerVolume();
     fMeanExcitationEnergy = G4Exp(fLogMeanExcEnergy);
   }
 
   fShellCorrectionVector = new G4double[3]; 
 
   for (G4int j=0; j<=2; j++)
   {
     fShellCorrectionVector[j] = 0.;
 
     for (size_t k=0; k<nElements; k++) {
       fShellCorrectionVector[j] += nAtomsPerVolume[k]
     *(((*elmVector)[k])->GetIonisation()->GetShellCorrectionVector())[j];
     }
     fShellCorrectionVector[j] *= 2.0/fMaterial->GetTotNbOfElectPerVolume();
   } 
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ DensityCorrection()

G4double G4IonisParamMat::DensityCorrection ( G4double x )

inline

Definition at line 222 of file G4IonisParamMat.hh.

 {
   // x = log10(beta*gamma)  
 
   G4double y = 0.0;
   if(x < fX0density) {
     if(fD0density > 0.0) { y = fD0density*G4Exp(twoln10*(x - fX0density)); }
   } else if(x >= fX1density) { y = twoln10*x - fCdensity; }
   else {y = twoln10*x - fCdensity + fAdensity*G4Exp(G4Log(fX1density - x)*fMdensity);}
   return y;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ FindMeanExcitationEnergy()

G4double G4IonisParamMat::FindMeanExcitationEnergy ( const G4String & chFormula )

Definition at line 449 of file G4IonisParamMat.cc.

 {
   // The data on mean excitation energy for compaunds
   // from "Stopping Powers for Electrons and Positrons"
   // ICRU Report N#37, 1984  (energy in eV)
 
   size_t numberOfMolecula = 54; 
   static const G4String name[54] = {
     // gas 0 - 12
     "NH_3",       "C_4H_10",     "CO_2",      "C_2H_6",      "C_7H_16-Gas",
     // "G4_AMMONIA", "G4_BUTANE","G4_CARBON_DIOXIDE","G4_ETHANE", "G4_N-HEPTANE"
     "C_6H_14-Gas",   "CH_4",     "NO",        "N_2O",        "C_8H_18-Gas",      
     // "G4_N-HEXANE" , "G4_METHANE", "x", "G4_NITROUS_OXIDE", "G4_OCTANE"
     "C_5H_12-Gas",   "C_3H_8",   "H_2O-Gas",                        
     // "G4_N-PENTANE", "G4_PROPANE", "G4_WATER_VAPOR"
 
     // liquid 13 - 39
     "C_3H_6O",    "C_6H_5NH_2",  "C_6H_6",    "C_4H_9OH",    "CCl_4", 
     //"G4_ACETONE","G4_ANILINE","G4_BENZENE","G4_N-BUTYL_ALCOHOL","G4_CARBON_TETRACHLORIDE"  
     "C_6H_5Cl",   "CHCl_3",      "C_6H_12",   "C_6H_4Cl_2",  "C_4Cl_2H_8O",
     //"G4_CHLOROBENZENE","G4_CHLOROFORM","G4_CYCLOHEXANE","G4_1,2-DICHLOROBENZENE","G4_DICHLORODIETHYL_ETHER"
     "C_2Cl_2H_4", "(C_2H_5)_2O", "C_2H_5OH",  "C_3H_5(OH)_3","C_7H_16",
     //"G4_1,2-DICHLOROETHANE","G4_DIETHYL_ETHER","G4_ETHYL_ALCOHOL","G4_GLYCEROL","G4_N-HEPTANE"
     "C_6H_14",    "CH_3OH",      "C_6H_5NO_2","C_5H_12",     "C_3H_7OH",
     //"G4_N-HEXANE","G4_METHANOL","G4_NITROBENZENE","G4_N-PENTANE","G4_N-PROPYL_ALCOHOL",  
     "C_5H_5N",    "C_8H_8",      "C_2Cl_4",   "C_7H_8",      "C_2Cl_3H",
     //"G4_PYRIDINE","G4_POLYSTYRENE","G4_TETRACHLOROETHYLENE","G4_TOLUENE","G4_TRICHLOROETHYLENE"
     "H_2O",       "C_8H_10",                                                 
     // "G4_WATER", "G4_XYLENE"
 
     // solid 40 - 53
     "C_5H_5N_5",  "C_5H_5N_5O",  "(C_6H_11NO)-nylon",  "C_25H_52",
     // "G4_ADENINE", "G4_GUANINE", "G4_NYLON-6-6", "G4_PARAFFIN"
     "(C_2H_4)-Polyethylene",     "(C_5H_8O_2)-Polymethil_Methacrylate",
     // "G4_ETHYLENE", "G4_PLEXIGLASS"
     "(C_8H_8)-Polystyrene",      "A-150-tissue",       "Al_2O_3",  "CaF_2",
     // "G4_POLYSTYRENE", "G4_A-150_TISSUE", "G4_ALUMINUM_OXIDE", "G4_CALCIUM_FLUORIDE"
     "LiF",        "Photo_Emulsion",  "(C_2F_4)-Teflon",  "SiO_2"
     // "G4_LITHIUM_FLUORIDE", "G4_PHOTO_EMULSION", "G4_TEFLON", "G4_SILICON_DIOXIDE"
   } ;
     
   static const G4double meanExcitation[54] = {
 
     53.7,   48.3,  85.0,  45.4,  49.2,
     49.1,   41.7,  87.8,  84.9,  49.5,
     48.2,   47.1,  71.6,
 
     64.2,   66.2,  63.4,  59.9,  166.3,
     89.1,  156.0,  56.4, 106.5,  103.3, 
    111.9,   60.0,  62.9,  72.6,   54.4,  
     54.0,  67.6,   75.8,  53.6,   61.1,  
     66.2,  64.0,  159.2,  62.5,  148.1,  
     75.0,  61.8,
 
     71.4,  75.0,   63.9,  48.3,   57.4,
     74.0,  68.7,   65.1, 145.2,  166.,
     94.0, 331.0,   99.1, 139.2 
   };
 
   G4double x = fMeanExcitationEnergy;
 
   for(size_t i=0; i<numberOfMolecula; i++) {
     if(chFormula == name[i]) {
       x = meanExcitation[i]*eV;
       break;
     }
   }
   return x;
 }

Here is the caller graph for this function:

◆ GetAdensity()

G4double G4IonisParamMat::GetAdensity ( ) const

inline

Definition at line 93 of file G4IonisParamMat.hh.

93 {return fAdensity;};

G4IonisParamMat::fAdensity

G4double fAdensity

Definition: G4IonisParamMat.hh:186

Here is the caller graph for this function:

◆ GetAdjustmentFactor()

G4double G4IonisParamMat::GetAdjustmentFactor ( ) const

inline

Definition at line 87 of file G4IonisParamMat.hh.

87 {return fAdjustmentFactor;};

G4IonisParamMat::fAdjustmentFactor

G4double fAdjustmentFactor

Definition: G4IonisParamMat.hh:192

◆ GetBirksConstant()

G4double G4IonisParamMat::GetBirksConstant ( ) const

inline

Definition at line 139 of file G4IonisParamMat.hh.

139 {return fBirks;};

G4IonisParamMat::fBirks

G4double fBirks

Definition: G4IonisParamMat.hh:211

Here is the caller graph for this function:

◆ GetCdensity()

G4double G4IonisParamMat::GetCdensity ( ) const

inline

Definition at line 89 of file G4IonisParamMat.hh.

89 {return fCdensity;};

G4IonisParamMat::fCdensity

G4double fCdensity

Definition: G4IonisParamMat.hh:184

Here is the caller graph for this function:

◆ GetD0density()

G4double G4IonisParamMat::GetD0density ( ) const

inline

Definition at line 99 of file G4IonisParamMat.hh.

99 {return fD0density;};

G4IonisParamMat::fD0density

G4double fD0density

Definition: G4IonisParamMat.hh:189

Here is the call graph for this function:

◆ GetDensityEffectData()

G4DensityEffectData * G4IonisParamMat::GetDensityEffectData ( )

static

Definition at line 173 of file G4IonisParamMat.cc.

 {
   return fDensityData;
 }

Here is the caller graph for this function:

◆ GetEnergy0fluct()

G4double G4IonisParamMat::GetEnergy0fluct ( ) const

inline

Definition at line 121 of file G4IonisParamMat.hh.

121 {return fEnergy0fluct;};

G4IonisParamMat::fEnergy0fluct

G4double fEnergy0fluct

Definition: G4IonisParamMat.hh:201

Here is the caller graph for this function:

◆ GetEnergy1fluct()

G4double G4IonisParamMat::GetEnergy1fluct ( ) const

inline

Definition at line 113 of file G4IonisParamMat.hh.

113 {return fEnergy1fluct;};

G4IonisParamMat::fEnergy1fluct

G4double fEnergy1fluct

Definition: G4IonisParamMat.hh:197

Here is the caller graph for this function:

◆ GetEnergy2fluct()

G4double G4IonisParamMat::GetEnergy2fluct ( ) const

inline

Definition at line 117 of file G4IonisParamMat.hh.

117 {return fEnergy2fluct;};

G4IonisParamMat::fEnergy2fluct

G4double fEnergy2fluct

Definition: G4IonisParamMat.hh:199

Here is the caller graph for this function:

◆ GetF1fluct()

G4double G4IonisParamMat::GetF1fluct ( ) const

inline

Definition at line 109 of file G4IonisParamMat.hh.

109 {return fF1fluct;};

G4IonisParamMat::fF1fluct

G4double fF1fluct

Definition: G4IonisParamMat.hh:195

Here is the caller graph for this function:

◆ GetF2fluct()

G4double G4IonisParamMat::GetF2fluct ( ) const

inline

Definition at line 111 of file G4IonisParamMat.hh.

111 {return fF2fluct;};

G4IonisParamMat::fF2fluct

G4double fF2fluct

Definition: G4IonisParamMat.hh:196

Here is the caller graph for this function:

◆ GetFermiEnergy()

G4double G4IonisParamMat::GetFermiEnergy ( ) const

inline

Definition at line 129 of file G4IonisParamMat.hh.

129 {return fFermiEnergy;};

G4IonisParamMat::fFermiEnergy

G4double fFermiEnergy

Definition: G4IonisParamMat.hh:206

Here is the caller graph for this function:

◆ GetInvA23()

G4double G4IonisParamMat::GetInvA23 ( ) const

inline

Definition at line 133 of file G4IonisParamMat.hh.

133 {return fInvA23;};

G4IonisParamMat::fInvA23

G4double fInvA23

Definition: G4IonisParamMat.hh:208

Here is the caller graph for this function:

◆ GetLFactor()

G4double G4IonisParamMat::GetLFactor ( ) const

inline

Definition at line 131 of file G4IonisParamMat.hh.

131 {return fLfactor;};

G4IonisParamMat::fLfactor

G4double fLfactor

Definition: G4IonisParamMat.hh:207

◆ GetLogEnergy1fluct()

G4double G4IonisParamMat::GetLogEnergy1fluct ( ) const

inline

Definition at line 115 of file G4IonisParamMat.hh.

115 {return fLogEnergy1fluct;};

G4IonisParamMat::fLogEnergy1fluct

G4double fLogEnergy1fluct

Definition: G4IonisParamMat.hh:198

Here is the caller graph for this function:

◆ GetLogEnergy2fluct()

G4double G4IonisParamMat::GetLogEnergy2fluct ( ) const

inline

Definition at line 119 of file G4IonisParamMat.hh.

119 {return fLogEnergy2fluct;};

G4IonisParamMat::fLogEnergy2fluct

G4double fLogEnergy2fluct

Definition: G4IonisParamMat.hh:200

Here is the caller graph for this function:

◆ GetLogMeanExcEnergy()

G4double G4IonisParamMat::GetLogMeanExcEnergy ( ) const

inline

Definition at line 77 of file G4IonisParamMat.hh.

77 {return fLogMeanExcEnergy;};

G4IonisParamMat::fLogMeanExcEnergy

G4double fLogMeanExcEnergy

Definition: G4IonisParamMat.hh:179

Here is the caller graph for this function:

◆ GetMdensity()

G4double G4IonisParamMat::GetMdensity ( ) const

inline

Definition at line 91 of file G4IonisParamMat.hh.

91 {return fMdensity;};

G4IonisParamMat::fMdensity

G4double fMdensity

Definition: G4IonisParamMat.hh:185

Here is the caller graph for this function:

◆ GetMeanEnergyPerIonPair()

G4double G4IonisParamMat::GetMeanEnergyPerIonPair ( ) const

inline

Definition at line 145 of file G4IonisParamMat.hh.

145 {return fMeanEnergyPerIon;};

G4IonisParamMat::fMeanEnergyPerIon

G4double fMeanEnergyPerIon

Definition: G4IonisParamMat.hh:213

Here is the call graph for this function:

Here is the caller graph for this function:

◆ GetMeanExcitationEnergy()

G4double G4IonisParamMat::GetMeanExcitationEnergy ( ) const

inline

Definition at line 71 of file G4IonisParamMat.hh.

71 {return fMeanExcitationEnergy;};

G4IonisParamMat::fMeanExcitationEnergy

G4double fMeanExcitationEnergy

Definition: G4IonisParamMat.hh:178

Here is the call graph for this function:

Here is the caller graph for this function:

◆ GetPlasmaEnergy()

G4double G4IonisParamMat::GetPlasmaEnergy ( ) const

inline

Definition at line 85 of file G4IonisParamMat.hh.

85 {return fPlasmaEnergy;};

G4IonisParamMat::fPlasmaEnergy

G4double fPlasmaEnergy

Definition: G4IonisParamMat.hh:191

◆ GetRateionexcfluct()

G4double G4IonisParamMat::GetRateionexcfluct ( ) const

inline

Definition at line 123 of file G4IonisParamMat.hh.

123 {return fRateionexcfluct;};

G4IonisParamMat::fRateionexcfluct

G4double fRateionexcfluct

Definition: G4IonisParamMat.hh:202

Here is the caller graph for this function:

◆ GetShellCorrectionVector()

G4double* G4IonisParamMat::GetShellCorrectionVector ( ) const

inline

Definition at line 79 of file G4IonisParamMat.hh.

79 {return fShellCorrectionVector;};

G4IonisParamMat::fShellCorrectionVector

G4double * fShellCorrectionVector

Definition: G4IonisParamMat.hh:180

Here is the caller graph for this function:

◆ GetTaul()

G4double G4IonisParamMat::GetTaul ( ) const

inline

Definition at line 81 of file G4IonisParamMat.hh.

81 {return fTaul;};

G4IonisParamMat::fTaul

G4double fTaul

Definition: G4IonisParamMat.hh:181

Here is the caller graph for this function:

◆ GetX0density()

G4double G4IonisParamMat::GetX0density ( ) const

inline

Definition at line 95 of file G4IonisParamMat.hh.

95 {return fX0density;};

G4IonisParamMat::fX0density

G4double fX0density

Definition: G4IonisParamMat.hh:187

Here is the caller graph for this function:

◆ GetX1density()

G4double G4IonisParamMat::GetX1density ( ) const

inline

Definition at line 97 of file G4IonisParamMat.hh.

97 {return fX1density;};

G4IonisParamMat::fX1density

G4double fX1density

Definition: G4IonisParamMat.hh:188

Here is the caller graph for this function:

◆ GetZeffective()

G4double G4IonisParamMat::GetZeffective ( ) const

inline

Definition at line 127 of file G4IonisParamMat.hh.

127 {return fZeff;};

G4IonisParamMat::fZeff

G4double fZeff

Definition: G4IonisParamMat.hh:205

Here is the caller graph for this function:

◆ operator!=()

G4int G4IonisParamMat::operator!= ( const G4IonisParamMat & ) const

private

Here is the caller graph for this function:

◆ operator=()

G4IonisParamMat& G4IonisParamMat::operator= ( const G4IonisParamMat & )

private

Here is the caller graph for this function:

◆ operator==()

G4int G4IonisParamMat::operator== ( const G4IonisParamMat & ) const

private

Here is the caller graph for this function:

◆ SetBirksConstant()

void G4IonisParamMat::SetBirksConstant ( G4double value )

inline

Definition at line 137 of file G4IonisParamMat.hh.

137 {fBirks = value;};

G4IonisParamMat::fBirks

G4double fBirks

Definition: G4IonisParamMat.hh:211

Here is the caller graph for this function:

◆ SetMeanEnergyPerIonPair()

void G4IonisParamMat::SetMeanEnergyPerIonPair ( G4double value )

inline

Definition at line 143 of file G4IonisParamMat.hh.

143 {fMeanEnergyPerIon = value;};

G4IonisParamMat::fMeanEnergyPerIon

G4double fMeanEnergyPerIon

Definition: G4IonisParamMat.hh:213

Here is the caller graph for this function:

◆ SetMeanExcitationEnergy()

void G4IonisParamMat::SetMeanExcitationEnergy ( G4double value )

Definition at line 422 of file G4IonisParamMat.cc.

 {
   if(value == fMeanExcitationEnergy || value <= 0.0) { return; }
   if (G4NistManager::Instance()->GetVerbose() > 1) {
     G4cout << "G4Material: Mean excitation energy is changed for "
            << fMaterial->GetName()
            << " Iold= " << fMeanExcitationEnergy/eV
            << "eV; Inew= " << value/eV << " eV;"
            << G4endl;
   }
   
   fMeanExcitationEnergy = value;
 
   // add corrections to density effect
   G4double newlog = G4Log(value);
   G4double corr = 2*(newlog - fLogMeanExcEnergy);
   fCdensity  += corr;
   fX0density += corr/twoln10;
   fX1density += corr/twoln10;
 
   // recompute parameters of fluctuation model
   fLogMeanExcEnergy = newlog;
   ComputeFluctModel();
 }

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

◆ fAdensity

G4double G4IonisParamMat::fAdensity

private

Definition at line 186 of file G4IonisParamMat.hh.

◆ fAdjustmentFactor

G4double G4IonisParamMat::fAdjustmentFactor

private

Definition at line 192 of file G4IonisParamMat.hh.

◆ fBirks

G4double G4IonisParamMat::fBirks

private

Definition at line 211 of file G4IonisParamMat.hh.

◆ fCdensity

G4double G4IonisParamMat::fCdensity

private

Definition at line 184 of file G4IonisParamMat.hh.

◆ fD0density

G4double G4IonisParamMat::fD0density

private

Definition at line 189 of file G4IonisParamMat.hh.

◆ fDensityData

G4DensityEffectData * G4IonisParamMat::fDensityData = 0

staticprivate

Definition at line 216 of file G4IonisParamMat.hh.

◆ fEnergy0fluct

G4double G4IonisParamMat::fEnergy0fluct

private

Definition at line 201 of file G4IonisParamMat.hh.

◆ fEnergy1fluct

G4double G4IonisParamMat::fEnergy1fluct

private

Definition at line 197 of file G4IonisParamMat.hh.

◆ fEnergy2fluct

G4double G4IonisParamMat::fEnergy2fluct

private

Definition at line 199 of file G4IonisParamMat.hh.

◆ fF1fluct

G4double G4IonisParamMat::fF1fluct

private

Definition at line 195 of file G4IonisParamMat.hh.

◆ fF2fluct

G4double G4IonisParamMat::fF2fluct

private

Definition at line 196 of file G4IonisParamMat.hh.

◆ fFermiEnergy

G4double G4IonisParamMat::fFermiEnergy

private

Definition at line 206 of file G4IonisParamMat.hh.

◆ fInvA23

G4double G4IonisParamMat::fInvA23

private

Definition at line 208 of file G4IonisParamMat.hh.

◆ fLfactor

G4double G4IonisParamMat::fLfactor

private

Definition at line 207 of file G4IonisParamMat.hh.

◆ fLogEnergy1fluct

G4double G4IonisParamMat::fLogEnergy1fluct

private

Definition at line 198 of file G4IonisParamMat.hh.

◆ fLogEnergy2fluct

G4double G4IonisParamMat::fLogEnergy2fluct

private

Definition at line 200 of file G4IonisParamMat.hh.

◆ fLogMeanExcEnergy

G4double G4IonisParamMat::fLogMeanExcEnergy

private

Definition at line 179 of file G4IonisParamMat.hh.

◆ fMaterial

G4Material* G4IonisParamMat::fMaterial

private

Definition at line 175 of file G4IonisParamMat.hh.

◆ fMdensity

G4double G4IonisParamMat::fMdensity

private

Definition at line 185 of file G4IonisParamMat.hh.

◆ fMeanEnergyPerIon

G4double G4IonisParamMat::fMeanEnergyPerIon

private

Definition at line 213 of file G4IonisParamMat.hh.

◆ fMeanExcitationEnergy

G4double G4IonisParamMat::fMeanExcitationEnergy

private

Definition at line 178 of file G4IonisParamMat.hh.

◆ fPlasmaEnergy

G4double G4IonisParamMat::fPlasmaEnergy

private

Definition at line 191 of file G4IonisParamMat.hh.

◆ fRateionexcfluct

G4double G4IonisParamMat::fRateionexcfluct

private

Definition at line 202 of file G4IonisParamMat.hh.

◆ fShellCorrectionVector

G4double* G4IonisParamMat::fShellCorrectionVector

private

Definition at line 180 of file G4IonisParamMat.hh.

◆ fTaul

G4double G4IonisParamMat::fTaul

private

Definition at line 181 of file G4IonisParamMat.hh.

◆ fX0density

G4double G4IonisParamMat::fX0density

private

Definition at line 187 of file G4IonisParamMat.hh.

◆ fX1density

G4double G4IonisParamMat::fX1density

private

Definition at line 188 of file G4IonisParamMat.hh.

◆ fZeff

G4double G4IonisParamMat::fZeff

private

Definition at line 205 of file G4IonisParamMat.hh.

◆ twoln10

G4double G4IonisParamMat::twoln10

private

Definition at line 217 of file G4IonisParamMat.hh.

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

Private Attributes
G4Material *	fMaterial

G4double	fMeanExcitationEnergy

G4double	fLogMeanExcEnergy

G4double *	fShellCorrectionVector

G4double	fTaul

G4double	fCdensity

G4double	fMdensity

G4double	fAdensity

G4double	fX0density

G4double	fX1density

G4double	fD0density

G4double	fPlasmaEnergy

G4double	fAdjustmentFactor

G4double	fF1fluct

G4double	fF2fluct

G4double	fEnergy1fluct

G4double	fLogEnergy1fluct

G4double	fEnergy2fluct

G4double	fLogEnergy2fluct

G4double	fEnergy0fluct

G4double	fRateionexcfluct

G4double	fZeff

G4double	fFermiEnergy

G4double	fLfactor

G4double	fInvA23

G4double	fBirks

G4double	fMeanEnergyPerIon

G4double	twoln10

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4IonisParamMat() [1/3]

◆ ~G4IonisParamMat()

◆ G4IonisParamMat() [2/3]

◆ G4IonisParamMat() [3/3]

Member Function Documentation

◆ ComputeDensityEffect()

◆ ComputeFluctModel()

◆ ComputeIonParameters()

◆ ComputeMeanParameters()

◆ DensityCorrection()

◆ FindMeanExcitationEnergy()

◆ GetAdensity()

◆ GetAdjustmentFactor()

◆ GetBirksConstant()

◆ GetCdensity()

◆ GetD0density()

◆ GetDensityEffectData()

◆ GetEnergy0fluct()

◆ GetEnergy1fluct()

◆ GetEnergy2fluct()

◆ GetF1fluct()

◆ GetF2fluct()

◆ GetFermiEnergy()

◆ GetInvA23()

◆ GetLFactor()

◆ GetLogEnergy1fluct()

◆ GetLogEnergy2fluct()

◆ GetLogMeanExcEnergy()

◆ GetMdensity()

◆ GetMeanEnergyPerIonPair()

◆ GetMeanExcitationEnergy()

◆ GetPlasmaEnergy()

◆ GetRateionexcfluct()

◆ GetShellCorrectionVector()

◆ GetTaul()

◆ GetX0density()

◆ GetX1density()

◆ GetZeffective()

◆ operator!=()

◆ operator=()

◆ operator==()

◆ SetBirksConstant()

◆ SetMeanEnergyPerIonPair()

◆ SetMeanExcitationEnergy()

Member Data Documentation

◆ fAdensity

◆ fAdjustmentFactor

◆ fBirks

◆ fCdensity

◆ fD0density

◆ fDensityData

◆ fEnergy0fluct

◆ fEnergy1fluct

◆ fEnergy2fluct

◆ fF1fluct

◆ fF2fluct

◆ fFermiEnergy

◆ fInvA23

◆ fLfactor

◆ fLogEnergy1fluct

◆ fLogEnergy2fluct

◆ fLogMeanExcEnergy

◆ fMaterial

◆ fMdensity

◆ fMeanEnergyPerIon

◆ fMeanExcitationEnergy

◆ fPlasmaEnergy

◆ fRateionexcfluct

◆ fShellCorrectionVector

◆ fTaul

◆ fX0density

◆ fX1density

◆ fZeff