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 82 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 244 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 217 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 flimitElastic[1]; }

Referenced by G4MicroElecElasticModel_new::Initialise().

GetElasticModelLowLimit()

G4double G4MicroElecMaterialStructure::GetElasticModelLowLimit ( )

inline

Definition at line 100 of file G4MicroElecMaterialStructure.hh.

{return flimitElastic[0];}

Referenced by G4MicroElecElasticModel_new::Initialise().

GetEnergyGap()

G4double G4MicroElecMaterialStructure::GetEnergyGap ( )

inline

Definition at line 94 of file G4MicroElecMaterialStructure.hh.

{ return energyGap; }

Referenced by G4MicroElecCapture::PostStepDoIt(), and G4MicroElecInelasticModel_new::SampleSecondaries().

GetInelasticModelHighLimit()

G4double G4MicroElecMaterialStructure::GetInelasticModelHighLimit

(

G4int

pdg

)

Definition at line 298 of file G4MicroElecMaterialStructure.cc.

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

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

GetInelasticModelLowLimit()

G4double G4MicroElecMaterialStructure::GetInelasticModelLowLimit

(

G4int

pdg

)

Definition at line 282 of file G4MicroElecMaterialStructure.cc.

{
  G4double res = 0.0;
  if(pdg == 11)
  {
    res = flimitInelastic[0];
  }
  else if(pdg == 2212)
  {
    res = flimitInelastic[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; }

Referenced by G4MicroElecCapture::PostStepDoIt().

GetLimitEnergy()

G4double G4MicroElecMaterialStructure::GetLimitEnergy

(

G4int

level

)

Definition at line 273 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(), G4MicroElecElasticModel_new::Initialise(), and G4MicroElecSurface::Initialise().

GetZ()

G4double G4MicroElecMaterialStructure::GetZ

(

G4int

Shell

)

Definition at line 224 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 314 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 97 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")
        {
          flimitInelastic[0] = data;
        }
        if(nameParameter == "DielectricModelHighEnergyLimit_e")
        {
          flimitInelastic[1] = data;
        }
        if(nameParameter == "DielectricModelLowEnergyLimit_p")
        {
          flimitInelastic[2] = data;
        }
        if(nameParameter == "DielectricModelHighEnergyLimit_p")
        {
          flimitInelastic[3] = data;
        }
 
        if(nameParameter == "ElasticModelLowEnergyLimit")
        {
          flimitElastic[0] = data;
        }
        if(nameParameter == "ElasticModelHighEnergyLimit")
        {
          flimitElastic[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