G4MicroElecMaterialStructure Class Reference

#include <G4MicroElecMaterialStructure.hh>

Public Member Functions
	G4MicroElecMaterialStructure (const G4String &matName="")
virtual	~G4MicroElecMaterialStructure ()=default
void	ReadMaterialFile ()
G4double	Energy (G4int level)
G4int	NumberOfLevels ()
G4double	GetZ (G4int Shell)
G4double	ConvertUnit (const G4String &unitName)
G4double	GetEnergyGap ()
G4double	GetInitialEnergy ()
G4int	GetEADL_Enumerator (G4int shell)
G4double	GetWorkFunction ()
G4String	GetMaterialName ()
G4double	GetLimitEnergy (G4int level)
G4double	GetElasticModelLowLimit ()
G4double	GetElasticModelHighLimit ()
G4double	GetInelasticModelLowLimit (G4int pdg)
G4double	GetInelasticModelHighLimit (G4int pdg)
G4bool	IsShellWeaklyBound (G4int level)

Detailed Description

Definition at line 83 of file G4MicroElecMaterialStructure.hh.

Constructor & Destructor Documentation

G4MicroElecMaterialStructure()

G4MicroElecMaterialStructure::G4MicroElecMaterialStructure

(

const G4String &

matName = ""

)

Definition at line 81 of file G4MicroElecMaterialStructure.cc.

{
  materialName = matName;
  if (matName == "Vacuum" || matName == "uum") {
    workFunction = 0;
    initialEnergy = 0;
  }
  else {
    ReadMaterialFile();
  }
  nLevels = (G4int)energyConstant.size();
}

~G4MicroElecMaterialStructure()

virtual G4MicroElecMaterialStructure::~G4MicroElecMaterialStructure ( )

virtualdefault

Member Function Documentation

ConvertUnit()

G4double G4MicroElecMaterialStructure::ConvertUnit

(

const G4String &

unitName

)

Definition at line 243 of file G4MicroElecMaterialStructure.cc.

{
  G4double unitValue = 0;
  if(unitName == "meV")
  {
    unitValue = 1e-3 * CLHEP::eV;
  }
  else if(unitName == "eV")
  {
    unitValue = CLHEP::eV;
  }
  else if(unitName == "keV")
  {
    unitValue = CLHEP::keV;
  }
  else if(unitName == "MeV")
  {
    unitValue = CLHEP::MeV;
  }
  else if(unitName == "noUnit")
  {
    unitValue = 1;
  }
 
  return unitValue;
}

Referenced by ReadMaterialFile().

Energy()

G4double G4MicroElecMaterialStructure::Energy

(

G4int

level

)

Definition at line 216 of file G4MicroElecMaterialStructure.cc.

{
  return (level >= 0 && level < nLevels) ? energyConstant[level] : 0.0;
}

Referenced by G4MicroElecInelasticModel_new::CrossSectionPerVolume(), and G4MicroElecInelasticModel_new::SampleSecondaries().

GetEADL_Enumerator()

G4int G4MicroElecMaterialStructure::GetEADL_Enumerator ( G4int  shell )

inline

Definition at line 96 of file G4MicroElecMaterialStructure.hh.

{ return EADL_Enumerator[shell]; };

Referenced by G4MicroElecInelasticModel_new::SampleSecondaries().

GetElasticModelHighLimit()

G4double G4MicroElecMaterialStructure::GetElasticModelHighLimit ( )

inline

Definition at line 101 of file G4MicroElecMaterialStructure.hh.

{ return limitElastic[1]; }

Referenced by G4MicroElecElasticModel_new::Initialise().

GetElasticModelLowLimit()

G4double G4MicroElecMaterialStructure::GetElasticModelLowLimit ( )

inline

Definition at line 100 of file G4MicroElecMaterialStructure.hh.

{return limitElastic[0];}

Referenced by G4MicroElecElasticModel_new::Initialise().

GetEnergyGap()

G4double G4MicroElecMaterialStructure::GetEnergyGap ( )

inline

Definition at line 94 of file G4MicroElecMaterialStructure.hh.

{ return energyGap; }

Referenced by G4MicroElecInelasticModel_new::SampleSecondaries().

GetInelasticModelHighLimit()

G4double G4MicroElecMaterialStructure::GetInelasticModelHighLimit

(

G4int

pdg

)

Definition at line 297 of file G4MicroElecMaterialStructure.cc.

{
  G4double res = 0.0;
  if(pdg == 11)
  {
    res = limitInelastic[1];
  }
  else if(pdg == 2212)
  {
    res = limitInelastic[3];
  }
  return res;
}

Referenced by G4MicroElecInelasticModel_new::CrossSectionPerVolume(), and G4MicroElecInelasticModel_new::SampleSecondaries().

GetInelasticModelLowLimit()

G4double G4MicroElecMaterialStructure::GetInelasticModelLowLimit

(

G4int

pdg

)

Definition at line 281 of file G4MicroElecMaterialStructure.cc.

{
  G4double res = 0.0;
  if(pdg == 11)
  {
    res = limitInelastic[0];
  }
  else if(pdg == 2212)
  {
    res = limitInelastic[2];
  }
  return res;
}

Referenced by G4MicroElecInelasticModel_new::CrossSectionPerVolume(), and G4MicroElecInelasticModel_new::SampleSecondaries().

GetInitialEnergy()

G4double G4MicroElecMaterialStructure::GetInitialEnergy ( )

inline

Definition at line 95 of file G4MicroElecMaterialStructure.hh.

{ return initialEnergy; }

GetLimitEnergy()

G4double G4MicroElecMaterialStructure::GetLimitEnergy

(

G4int

level

)

Definition at line 272 of file G4MicroElecMaterialStructure.cc.

{
  G4double E = LimitEnergy[level];
  if (IsShellWeaklyBound(level)) { E = energyGap+ initialEnergy; }
  return E;
}

Referenced by G4MicroElecInelasticModel_new::DifferentialCrossSection(), and G4MicroElecInelasticModel_new::SampleSecondaries().

GetMaterialName()

G4String G4MicroElecMaterialStructure::GetMaterialName ( )

inline

Definition at line 98 of file G4MicroElecMaterialStructure.hh.

{ return materialName; };

GetWorkFunction()

G4double G4MicroElecMaterialStructure::GetWorkFunction ( )

inline

Definition at line 97 of file G4MicroElecMaterialStructure.hh.

{ return workFunction; };

Referenced by G4MicroElecSurface::BuildPhysicsTable(), and G4MicroElecElasticModel_new::Initialise().

GetZ()

G4double G4MicroElecMaterialStructure::GetZ

(

G4int

Shell

)

Definition at line 223 of file G4MicroElecMaterialStructure.cc.

{
  if (Shell >= 0 && Shell < nLevels) {
    if(!isCompound)
    {
      return Z;
    }
    else
    {
      return compoundShellZ[Shell];
    }
  }
  else
  {
    return 0;
  }
}

Referenced by G4MicroElecInelasticModel_new::SampleSecondaries().

IsShellWeaklyBound()

G4bool G4MicroElecMaterialStructure::IsShellWeaklyBound

(

G4int

level

)

Definition at line 313 of file G4MicroElecMaterialStructure.cc.

{
  return isShellWeaklyBoundVector[level];
}

Referenced by GetLimitEnergy(), and G4MicroElecInelasticModel_new::SampleSecondaries().

NumberOfLevels()

G4int G4MicroElecMaterialStructure::NumberOfLevels ( )

inline

Definition at line 91 of file G4MicroElecMaterialStructure.hh.

{ return nLevels; }

Referenced by G4MicroElecInelasticModel_new::CrossSectionPerVolume(), and G4MicroElecInelasticModel_new::Initialise().

ReadMaterialFile()

void G4MicroElecMaterialStructure::ReadMaterialFile

(

)

Definition at line 96 of file G4MicroElecMaterialStructure.cc.

{
  const char* path = G4FindDataDir("G4LEDATA");
 
  if (materialName[0] == 'G' && materialName[1] == '4') {
    //in the case the NIST database is used
    materialName.erase(0, 1);
    materialName.erase(0, 1);
    materialName.erase(0, 1);
  }
  
  std::ostringstream fileName;
  fileName << path << "/microelec/Structure/Data_" + materialName + ".dat";
  std::ifstream fichier(fileName.str().c_str());
  
  int varLength = 0;
  G4String nameParameter;
  
  G4String unitName;    
  G4double unitValue;
  G4double data;
  G4String filler;  
  G4String type;
  
  if (fichier)
    {
      fichier >> filler >> type;
      materialName = filler;
      if (type == "Compound") {isCompound = true; Z = 0; }
      else { isCompound = false; Z = std::stoi(type); }
      while(!fichier.eof()) {
    
    getline(fichier, filler);
    std::stringstream line(filler);
    
    if (filler[0] == '#' || filler.empty()) {continue;}
    
    line >> varLength;
    line >> nameParameter;
    line >> unitName;
    unitValue = ConvertUnit(unitName);
    
    for (int i = 0; i < varLength; i++)
      {
        line >> data;   data = data*unitValue;
 
      if(nameParameter == "WorkFunction")
      {
        workFunction = data;
      }
      if(nameParameter == "EnergyGap")
      {
        energyGap = data;
      }
 
      if(nameParameter == "EnergyPeak")
      {
        energyConstant.push_back(data);
      }
      if(nameParameter == "EnergyLimit")
      {
        LimitEnergy.push_back(data);
      }
      if(nameParameter == "EADL")
      {
        EADL_Enumerator.push_back(data);
      }
 
      if (nameParameter == "WeaklyBoundShell")
          {if (data == 0) { isShellWeaklyBoundVector.push_back(false); }
        else {isShellWeaklyBoundVector.push_back(true);}}
 
        if(nameParameter == "WeaklyBoundInitialEnergy")
        {
          initialEnergy = data;
        }
 
        if(nameParameter == "ShellAtomicNumber")
        {
          compoundShellZ.push_back(data);
        }
 
        if(nameParameter == "DielectricModelLowEnergyLimit_e")
        {
          limitInelastic[0] = data;
        }
        if(nameParameter == "DielectricModelHighEnergyLimit_e")
        {
          limitInelastic[1] = data;
        }
        if(nameParameter == "DielectricModelLowEnergyLimit_p")
        {
          limitInelastic[2] = data;
        }
        if(nameParameter == "DielectricModelHighEnergyLimit_p")
        {
          limitInelastic[3] = data;
        }
 
        if(nameParameter == "ElasticModelLowEnergyLimit")
        {
          limitElastic[0] = data;
        }
        if(nameParameter == "ElasticModelHighEnergyLimit")
        {
          limitElastic[1] = data;
        }
    }
      }
      fichier.close();  // on ferme le fichier
    }
  else {
    G4String str = "file ";
    str += fileName.str() + " not found!";
    G4Exception("G4MicroElecMaterialStructure::ReadMaterialFile", "em0002", FatalException, str);
  }
}

Referenced by G4MicroElecMaterialStructure().

---

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

• G4MicroElecMaterialStructure.hh
• G4MicroElecMaterialStructure.cc