G4IonisParamMat Class Reference

#include <G4IonisParamMat.hh>

Public Member Functions
	G4IonisParamMat (const G4Material *)
	~G4IonisParamMat ()
G4double	GetMeanExcitationEnergy () const
void	SetMeanExcitationEnergy (G4double value)
G4double	FindMeanExcitationEnergy (const G4Material *) const
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
void	SetDensityEffectParameters (G4double cd, G4double md, G4double ad, G4double x0, G4double x1, G4double d0)
void	SetDensityEffectParameters (const G4Material *bmat)
void	ComputeDensityEffectOnFly (G4bool)
G4double	DensityCorrection (G4double x)
G4DensityEffectCalculator *	GetDensityEffectCalculator ()
G4double	GetDensityCorrection (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 ()

Detailed Description

Definition at line 60 of file G4IonisParamMat.hh.

Constructor & Destructor Documentation

G4IonisParamMat() [1/2]

G4IonisParamMat::G4IonisParamMat

(

const G4Material *

material

)

Definition at line 59 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 == nullptr) { fDensityData = new G4DensityEffectData(); }
  fDensityEffectCalc = nullptr;
 
  // compute parameters
  ComputeMeanParameters();
  ComputeDensityEffectParameters();
  ComputeFluctModel();
  ComputeIonParameters();
}

~G4IonisParamMat()

G4IonisParamMat::~G4IonisParamMat

(

)

Definition at line 121 of file G4IonisParamMat.cc.

{
  delete fDensityEffectCalc;
  delete [] fShellCorrectionVector; 
  delete fDensityData;
  fDensityData = nullptr;
  fShellCorrectionVector = nullptr;
  fDensityEffectCalc = nullptr;
}

G4IonisParamMat() [2/2]

G4IonisParamMat::G4IonisParamMat

(

__void__ &

)

Definition at line 85 of file G4IonisParamMat.cc.

  : fMaterial(nullptr), fShellCorrectionVector(nullptr)
{
  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);
 
  fDensityEffectCalc = nullptr;
  if(fDensityData == nullptr) { fDensityData = new G4DensityEffectData(); }
}

Member Function Documentation

ComputeDensityEffectOnFly()

void G4IonisParamMat::ComputeDensityEffectOnFly

(

G4bool

val

)

Definition at line 525 of file G4IonisParamMat.cc.

{
  if(val) {
    if(!fDensityEffectCalc) { 
      G4int n = 0;
      for(size_t i=0; i<fMaterial->GetNumberOfElements(); ++i) {
    const G4int Z = fMaterial->GetElement(i)->GetZasInt();
        n += G4AtomicShells::GetNumberOfShells(Z);
      }
      // The last level is the conduction level.  If this is *not* a conductor,
      // make a dummy conductor level with zero electrons in it.
      fDensityEffectCalc = new G4DensityEffectCalculator(fMaterial, n);
    }
  } else {
    delete fDensityEffectCalc;
    fDensityEffectCalc = nullptr;
  }
}

Referenced by G4Material::ComputeDensityEffectOnFly().

DensityCorrection()

G4double G4IonisParamMat::DensityCorrection

(

G4double

x

)

Definition at line 185 of file G4IonisParamMat.cc.

{
  return (!fDensityEffectCalc) ? GetDensityCorrection(x) 
    : fDensityEffectCalc->ComputeDensityCorrection(x);
}

Referenced by G4MuBetheBlochModel::ComputeDEDXPerVolume(), G4BetheBlochModel::ComputeDEDXPerVolume(), G4LindhardSorensenIonModel::ComputeDEDXPerVolume(), and G4MollerBhabhaModel::ComputeDEDXPerVolume().

FindMeanExcitationEnergy()

G4double G4IonisParamMat::FindMeanExcitationEnergy

(

const G4Material *

mat

)

const

Definition at line 546 of file G4IonisParamMat.cc.

{
  G4double res = 0.0;
  // data from density effect data
  if(fDensityData) {
    G4int idx = fDensityData->GetIndex(mat->GetName());
    if(idx >= 0) {
      res = fDensityData->GetMeanIonisationPotential(idx);
    }
  }
 
  // The data on mean excitation energy for compaunds
  // from "Stopping Powers for Electrons and Positrons"
  // ICRU Report N#37, 1984  (energy in eV)
  // this value overwrites Density effect data
  G4String chFormula = mat->GetChemicalFormula();
  if(chFormula != "") { 
 
    static const size_t numberOfMolecula = 54; 
    static const G4String name[numberOfMolecula] = {
      // 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[numberOfMolecula] = {
 
      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 
    };
 
    for(size_t i=0; i<numberOfMolecula; i++) {
      if(chFormula == name[i]) {
    res = meanExcitation[i]*eV;
    break;
      }
    }
  }
  return res;
}

GetAdensity()

G4double G4IonisParamMat::GetAdensity ( ) const

inline

Definition at line 91 of file G4IonisParamMat.hh.

{return fAdensity;};

Referenced by G4EmCorrections::DensityCorrection(), and SetDensityEffectParameters().

GetAdjustmentFactor()

G4double G4IonisParamMat::GetAdjustmentFactor ( ) const

inline

Definition at line 85 of file G4IonisParamMat.hh.

{return fAdjustmentFactor;};

GetBirksConstant()

G4double G4IonisParamMat::GetBirksConstant ( ) const

inline

Definition at line 152 of file G4IonisParamMat.hh.

{return fBirks;};

Referenced by G4EmSaturation::DumpBirksCoefficients(), and G4EmSaturation::VisibleEnergyDeposition().

GetCdensity()

G4double G4IonisParamMat::GetCdensity ( ) const

inline

Definition at line 87 of file G4IonisParamMat.hh.

{return fCdensity;};

Referenced by G4EmCorrections::DensityCorrection(), and SetDensityEffectParameters().

GetD0density()

G4double G4IonisParamMat::GetD0density ( ) const

inline

Definition at line 97 of file G4IonisParamMat.hh.

{return fD0density;};

Referenced by SetDensityEffectParameters().

GetDensityCorrection()

G4double G4IonisParamMat::GetDensityCorrection ( G4double  x )

inline

Definition at line 240 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;
}

Referenced by G4DensityEffectCalculator::ComputeDensityCorrection(), and DensityCorrection().

GetDensityEffectCalculator()

G4DensityEffectCalculator * G4IonisParamMat::GetDensityEffectCalculator ( )

inline

Definition at line 112 of file G4IonisParamMat.hh.

  { return fDensityEffectCalc; }

GetDensityEffectData()

G4DensityEffectData * G4IonisParamMat::GetDensityEffectData ( )

static

Definition at line 178 of file G4IonisParamMat.cc.

{
  return fDensityData;
}

Referenced by G4NistMessenger::SetNewValue().

GetEnergy0fluct()

G4double G4IonisParamMat::GetEnergy0fluct ( ) const

inline

Definition at line 134 of file G4IonisParamMat.hh.

{return fEnergy0fluct;};

Referenced by G4UniversalFluctuation::SampleFluctuations().

GetEnergy1fluct()

G4double G4IonisParamMat::GetEnergy1fluct ( ) const

inline

Definition at line 126 of file G4IonisParamMat.hh.

{return fEnergy1fluct;};

Referenced by G4UniversalFluctuation::SampleFluctuations().

GetEnergy2fluct()

G4double G4IonisParamMat::GetEnergy2fluct ( ) const

inline

Definition at line 130 of file G4IonisParamMat.hh.

{return fEnergy2fluct;};

Referenced by G4UniversalFluctuation::SampleFluctuations().

GetF1fluct()

G4double G4IonisParamMat::GetF1fluct ( ) const

inline

Definition at line 122 of file G4IonisParamMat.hh.

{return fF1fluct;};

Referenced by G4UniversalFluctuation::SampleFluctuations().

GetF2fluct()

G4double G4IonisParamMat::GetF2fluct ( ) const

inline

Definition at line 124 of file G4IonisParamMat.hh.

{return fF2fluct;};

Referenced by G4UniversalFluctuation::SampleFluctuations().

GetFermiEnergy()

G4double G4IonisParamMat::GetFermiEnergy ( ) const

inline

Definition at line 142 of file G4IonisParamMat.hh.

{return fFermiEnergy;};

Referenced by G4ionEffectiveCharge::EffectiveCharge().

GetInvA23()

G4double G4IonisParamMat::GetInvA23 ( ) const

inline

Definition at line 146 of file G4IonisParamMat.hh.

{return fInvA23;};

Referenced by G4CoulombScattering::MinPrimaryEnergy(), G4WentzelOKandVIxSection::SetupKinematic(), and G4WentzelVIRelXSection::SetupKinematic().

GetLFactor()

G4double G4IonisParamMat::GetLFactor ( ) const

inline

Definition at line 144 of file G4IonisParamMat.hh.

{return fLfactor;};

GetLogEnergy1fluct()

G4double G4IonisParamMat::GetLogEnergy1fluct ( ) const

inline

Definition at line 128 of file G4IonisParamMat.hh.

{return fLogEnergy1fluct;};

Referenced by G4UniversalFluctuation::SampleFluctuations().

GetLogEnergy2fluct()

G4double G4IonisParamMat::GetLogEnergy2fluct ( ) const

inline

Definition at line 132 of file G4IonisParamMat.hh.

{return fLogEnergy2fluct;};

Referenced by G4UniversalFluctuation::SampleFluctuations().

GetLogMeanExcEnergy()

G4double G4IonisParamMat::GetLogMeanExcEnergy ( ) const

inline

Definition at line 75 of file G4IonisParamMat.hh.

{return fLogMeanExcEnergy;};

Referenced by G4UniversalFluctuation::SampleFluctuations().

GetMdensity()

G4double G4IonisParamMat::GetMdensity ( ) const

inline

Definition at line 89 of file G4IonisParamMat.hh.

{return fMdensity;};

Referenced by G4EmCorrections::DensityCorrection(), and SetDensityEffectParameters().

GetMeanEnergyPerIonPair()

G4double G4IonisParamMat::GetMeanEnergyPerIonPair ( ) const

inline

Definition at line 158 of file G4IonisParamMat.hh.

{return fMeanEnergyPerIon;};

Referenced by G4ElectronIonPair::DumpMeanEnergyPerIonPair(), G4ElectronIonPair::MeanNumberOfIonsAlongStep(), and G4eeToTwoGammaModel::SampleSecondaries().

GetMeanExcitationEnergy()

G4double G4IonisParamMat::GetMeanExcitationEnergy ( ) const

inline

Definition at line 69 of file G4IonisParamMat.hh.

{return fMeanExcitationEnergy;};

Referenced by G4EmCorrections::Bethe(), G4hImpactIonisation::BuildPhysicsTable(), G4MuBetheBlochModel::ComputeDEDXPerVolume(), G4BetheBlochModel::ComputeDEDXPerVolume(), G4LindhardSorensenIonModel::ComputeDEDXPerVolume(), G4MollerBhabhaModel::ComputeDEDXPerVolume(), G4tgbGeometryDumper::DumpMaterial(), G4DensityEffectCalculator::G4DensityEffectCalculator(), G4GDMLWriteMaterials::MaterialWrite(), G4IonParametrisedLossModel::MinEnergyCut(), G4MuBetheBlochModel::MinEnergyCut(), G4mplIonisationWithDeltaModel::MinEnergyCut(), G4AtimaEnergyLossModel::MinEnergyCut(), G4BetheBlochModel::MinEnergyCut(), G4BraggIonModel::MinEnergyCut(), G4LindhardSorensenIonModel::MinEnergyCut(), G4hImpactIonisation::PrintInfoDefinition(), and G4UniversalFluctuation::SampleFluctuations().

GetPlasmaEnergy()

G4double G4IonisParamMat::GetPlasmaEnergy ( ) const

inline

Definition at line 83 of file G4IonisParamMat.hh.

{return fPlasmaEnergy;};

Referenced by G4DensityEffectCalculator::G4DensityEffectCalculator().

GetRateionexcfluct()

G4double G4IonisParamMat::GetRateionexcfluct ( ) const

inline

Definition at line 136 of file G4IonisParamMat.hh.

{return fRateionexcfluct;};

GetShellCorrectionVector()

G4double * G4IonisParamMat::GetShellCorrectionVector ( ) const

inline

Definition at line 77 of file G4IonisParamMat.hh.

{return fShellCorrectionVector;};

Referenced by G4EmCorrections::ShellCorrectionSTD().

GetTaul()

G4double G4IonisParamMat::GetTaul ( ) const

inline

Definition at line 79 of file G4IonisParamMat.hh.

{return fTaul;};

Referenced by G4hBetheBlochModel::LowEnergyLimit(), and G4EmCorrections::ShellCorrectionSTD().

GetX0density()

G4double G4IonisParamMat::GetX0density ( ) const

inline

Definition at line 93 of file G4IonisParamMat.hh.

{return fX0density;};

Referenced by G4EmCorrections::DensityCorrection(), and SetDensityEffectParameters().

GetX1density()

G4double G4IonisParamMat::GetX1density ( ) const

inline

Definition at line 95 of file G4IonisParamMat.hh.

{return fX1density;};

Referenced by G4EmCorrections::DensityCorrection(), and SetDensityEffectParameters().

GetZeffective()

G4double G4IonisParamMat::GetZeffective ( ) const

inline

Definition at line 140 of file G4IonisParamMat.hh.

{return fZeff;};

Referenced by G4MollerBhabhaModel::ComputeDEDXPerVolume(), G4AtimaEnergyLossModel::ComputeDEDXPerVolume(), G4GoudsmitSaundersonMscModel::ComputeTruePathLengthLimit(), G4UrbanAdjointMscModel::ComputeTruePathLengthLimit(), G4AtimaFluctuations::Dispersion(), G4IonFluctuations::Dispersion(), G4ionEffectiveCharge::EffectiveCharge(), G4GoudsmitSaundersonTable::InitSCPCorrection(), and G4eDPWAElasticDCS::InitSCPCorrection().

SetBirksConstant()

void G4IonisParamMat::SetBirksConstant ( G4double  value )

inline

Definition at line 150 of file G4IonisParamMat.hh.

{fBirks = value;}; 

SetDensityEffectParameters() [1/2]

void G4IonisParamMat::SetDensityEffectParameters

(

const G4Material *

bmat

)

Definition at line 501 of file G4IonisParamMat.cc.

{
#ifdef G4MULTITHREADED
  G4MUTEXLOCK(&ionisMutex);
#endif
  const G4IonisParamMat* ipm = bmat->GetIonisation();
  fCdensity  = ipm->GetCdensity(); 
  fMdensity  = ipm->GetMdensity();
  fAdensity  = ipm->GetAdensity();
  fX0density = ipm->GetX0density();
  fX1density = ipm->GetX1density();
  fD0density = ipm->GetD0density();
 
  G4double corr = G4Log(bmat->GetDensity()/fMaterial->GetDensity());
  fCdensity  += corr;
  fX0density += corr/twoln10;
  fX1density += corr/twoln10;
#ifdef G4MULTITHREADED
  G4MUTEXUNLOCK(&ionisMutex);
#endif
}

SetDensityEffectParameters() [2/2]

void G4IonisParamMat::SetDensityEffectParameters	(	G4double	cd,
		G4double	md,
		G4double	ad,
		G4double	x0,
		G4double	x1,
		G4double	d0
	)

Definition at line 480 of file G4IonisParamMat.cc.

{
  // no check on consistence of user parameters  
#ifdef G4MULTITHREADED
  G4MUTEXLOCK(&ionisMutex);
#endif
  fCdensity  = cd; 
  fMdensity  = md;
  fAdensity  = ad;
  fX0density = x0;
  fX1density = x1;
  fD0density = d0;
#ifdef G4MULTITHREADED
  G4MUTEXUNLOCK(&ionisMutex);
#endif
}

SetMeanEnergyPerIonPair()

void G4IonisParamMat::SetMeanEnergyPerIonPair ( G4double  value )

inline

Definition at line 156 of file G4IonisParamMat.hh.

{fMeanEnergyPerIon = value;}; 

Referenced by G4ElectronIonPair::FindG4MeanEnergyPerIonPair().

SetMeanExcitationEnergy()

void G4IonisParamMat::SetMeanExcitationEnergy

(

G4double

value

)

Definition at line 453 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();
}

Referenced by G4tgbMaterialMgr::FindOrBuildG4Material(), and G4GDMLReadMaterials::MaterialRead().

---

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

• G4IonisParamMat.hh
• G4IonisParamMat.cc