G4FluoData Class Reference

#include <G4FluoData.hh>

Public Member Functions
	G4FluoData (const G4FluoData &)=delete
	G4FluoData (const G4String &dir)
void	LoadData (G4int Z)
size_t	NumberOfTransitions (G4int vacancyIndex) const
size_t	NumberOfVacancies () const
G4FluoData &	operator= (const G4FluoData &right)=delete
void	PrintData ()
G4double	StartShellEnergy (G4int initIndex, G4int vacancyIndex) const
G4int	StartShellId (G4int initIndex, G4int vacancyIndex) const
G4double	StartShellProb (G4int initIndex, G4int vacancyIndex) const
G4int	VacancyId (G4int vacancyIndex) const
	Given the index of the vacancy returns its identity. More...
	~G4FluoData ()

Private Attributes
std::map< G4int, G4DataVector *, std::less< G4int > >	energyMap
G4String	fluoDirectory
std::map< G4int, std::vector< G4FluoTransition * >, std::less< G4int > >	fluoTransitionTable
std::map< G4int, G4DataVector *, std::less< G4int > >	idMap
std::vector< G4int >	nInitShells
G4int	numberOfVacancies = 0
std::map< G4int, G4DataVector *, std::less< G4int > >	probabilityMap

Detailed Description

Definition at line 51 of file G4FluoData.hh.

Constructor & Destructor Documentation

G4FluoData() [1/2]

G4FluoData::G4FluoData ( const G4String &  dir )

explicit

Definition at line 44 of file G4FluoData.cc.

{
  fluoDirectory = dir;
}

References fluoDirectory.

~G4FluoData()

G4FluoData::~G4FluoData

(

)

Definition at line 51 of file G4FluoData.cc.

{ 
  for (auto& pos : idMap)
    {
      G4DataVector* dataSet = pos.second;
      delete dataSet;
    }
  
  for (auto& pos : energyMap)
    {
      G4DataVector* dataSet = pos.second;
      delete dataSet;
    }
  
  for (auto& pos: probabilityMap)
    {
      G4DataVector* dataSet = pos.second;
      delete dataSet;
    }
}

References energyMap, idMap, pos, and probabilityMap.

G4FluoData() [2/2]

G4FluoData::G4FluoData ( const G4FluoData &  )

delete

Member Function Documentation

LoadData()

void G4FluoData::LoadData

(

G4int

Z

)

Definition at line 204 of file G4FluoData.cc.

{ 
  // Build the complete string identifying the file with the data set
  std::ostringstream ost;
  if(Z != 0){
    ost << "/fl-tr-pr-"<< Z << ".dat";
  }
  else{
    ost << "/fl-tr-pr-"<<".dat"; 
  }
  G4String name(ost.str());
 
  
  char* path = std::getenv("G4LEDATA");
  if (!path)
    { 
      G4String excep("G4FluoData::LoadData()");
      G4Exception(excep,"em0006",FatalException,"Please set G4LEDATA");
      return;
    }
  
  G4String pathString(path);
  
  G4String dirFile = pathString + fluoDirectory + name;
   
  //G4cout << "G4FluoData:: LoadData() name: " << dirFile << G4endl;
   
   
  std::ifstream file(dirFile);
  std::filebuf* lsdp = file.rdbuf();
  
  if (! (lsdp->is_open()) )
    {
      G4String excep = "G4FluoData::LoadData()";
      G4String msg = "data file: " + dirFile + " not found";
      G4Exception(excep, "em0003",FatalException, msg );
      return;
    }
  
  G4double a = 0;
  G4int k = 1;
  G4int sLocal = 0;
  
  G4int vacIndex = 0;
  G4DataVector* initIds = new G4DataVector;
  G4DataVector* transEnergies = new G4DataVector;
  G4DataVector* transProbabilities = new G4DataVector;
  
  do {
    file >> a;
    G4int nColumns = 3;
    if (a == -1)
      {
    if (sLocal == 0)
      {
        // End of a shell data set
        idMap[vacIndex] = initIds;
            energyMap[vacIndex] = transEnergies;
        probabilityMap[vacIndex] = transProbabilities;
        //      G4double size=transProbabilities->size();
            G4int n = initIds->size();
        
        nInitShells.push_back(n);
        numberOfVacancies++;
        // Start of new shell data set
        initIds = new G4DataVector;
            transEnergies = new G4DataVector;
        transProbabilities = new G4DataVector;
            vacIndex++; 
      }      
    sLocal++;
    if (sLocal == nColumns)
      {
        sLocal = 0;
      }
      }
    // moved to the end in order to avoid possible leak
    /*    else if (a == -2)
      {
    // End of file; delete the empty vectors created 
    //when encountering the last -1 -1 row
    delete initIds;
    delete transEnergies;
    delete transProbabilities;
    }*/ 
    else
      { 
    if(k%nColumns == 2)
      {  
        // 2nd column is transition  probabilities
 
       if (a != -1) transProbabilities->push_back(a);
        
        k++;
      }
    else if (k%nColumns == 1)
      {
        // 1st column is shell id
        // if this is the first data of the shell, all the colums are equal 
        // to the shell Id; so we skip the next colums ang go to the next row
        if(initIds->size() == 0) {
          if (a != -1) initIds->push_back((G4int)a);
          file >> a;
          file >> a;
          k=k+2;
        } 
        else{ 
          if (a != -1) initIds->push_back(a);
        }
        k++;    
      }
    else if (k%nColumns == 0)
 
      {//third column is transition energies
 
        if (a != -1) 
          {G4double e = a * MeV;
        transEnergies->push_back(e);}
       
        k=1;
      }
      }
  } 
  while (a != -2); // end of file
  file.close();    
  delete initIds;
  delete transEnergies;
  delete transProbabilities;
}

References energyMap, FatalException, geant4_check_module_cycles::file, fluoDirectory, G4Exception(), idMap, MeV, CLHEP::detail::n, G4InuclParticleNames::name(), nInitShells, numberOfVacancies, probabilityMap, and Z.

Referenced by G4AtomicTransitionManager::Initialise().

NumberOfTransitions()

size_t G4FluoData::NumberOfTransitions

(

G4int

vacancyIndex

)

const

Given the index of a vacancy returns the number of shells starting from wich an electrons can fill the vacancy

Definition at line 103 of file G4FluoData.cc.

{
  G4int n = 0;
  if (vacancyIndex<0 || vacancyIndex>=numberOfVacancies)
    {
      G4Exception("G4FluoData::NumberOfTransitions()","de0002",JustWarning,
          "vacancyIndex outside boundaries, energy deposited locally");
      return 0;
    }
  else
    {
      n = nInitShells[vacancyIndex]-1;
      //-1 is necessary because the elements of the vector nInitShells
      //include also the vacancy shell:
      // -1 subtracts this last one
  }
 return n;
}

References G4Exception(), JustWarning, CLHEP::detail::n, nInitShells, and numberOfVacancies.

Referenced by G4AtomicTransitionManager::Initialise(), and PrintData().

NumberOfVacancies()

size_t G4FluoData::NumberOfVacancies

(

)

const

The method returns the number of shells in wich a vacancy can be filled by a radiative transition

Definition at line 74 of file G4FluoData.cc.

{
  return numberOfVacancies;
}

References numberOfVacancies.

operator=()

G4FluoData & G4FluoData::operator= ( const G4FluoData &  right )

delete

PrintData()

void G4FluoData::PrintData

(

void

)

Definition at line 336 of file G4FluoData.cc.

{
  for (G4int i = 0; i <numberOfVacancies; i++)
    {
      G4cout << "---- TransitionData for the vacancy nb "
         <<i
         <<" ----- "
         <<G4endl;
      
      for (size_t k = 0; k<NumberOfTransitions(i); k++)
    { 
      G4int id = StartShellId(k,i);
    // let's start from 1 because the first (index = 0) element of the vector
    // is the id of the initial vacancy
      G4double e = StartShellEnergy(k,i) /MeV;
      G4double p = StartShellProb(k,i); 
      G4cout << k <<") Shell id: " << id <<G4endl;
      G4cout << " - Transition energy = " << e << " MeV "<<G4endl;
      G4cout   << " - Transition probability = " << p <<G4endl;
      
    }
      G4cout << "-------------------------------------------------" 
         << G4endl;
    }
}

References G4cout, G4endl, MeV, NumberOfTransitions(), numberOfVacancies, StartShellEnergy(), StartShellId(), and StartShellProb().

StartShellEnergy()

G4double G4FluoData::StartShellEnergy	(	G4int	initIndex,
		G4int	vacancyIndex
	)		const

Given the indexes of the starting and final shells for the transition, returns the transition energy

Definition at line 152 of file G4FluoData.cc.

{
  G4double n = -1;
  
  if (vacancyIndex<0 || vacancyIndex>=numberOfVacancies)
    {
      G4Exception("G4FluoData::StartShellEnergy()","de0002",FatalErrorInArgument,
          "vacancyIndex outside boundaries");}
  else
    {     
      auto pos = energyMap.find(vacancyIndex);
     
      G4DataVector dataSet = *((*pos).second);
     
      G4int nData = dataSet.size();
      if (initIndex >= 0 && initIndex < nData)
    {
      n =  dataSet[initIndex];
    }
    }
  return n;
}

References energyMap, FatalErrorInArgument, G4Exception(), CLHEP::detail::n, numberOfVacancies, and pos.

Referenced by G4AtomicTransitionManager::Initialise(), and PrintData().

StartShellId()

G4int G4FluoData::StartShellId	(	G4int	initIndex,
		G4int	vacancyIndex
	)		const

Given the indexes of the starting and final shells for the transition, returns the identity of the starting one

Definition at line 124 of file G4FluoData.cc.

{
 G4int n = -1;
 
 if (vacancyIndex<0 || vacancyIndex>=numberOfVacancies)
   {
     G4Exception("G4FluoData::StartShellId()","de0002",FatalErrorInArgument,
         "vacancyIndex outside boundaries");
   }
 else
   {
     auto pos = idMap.find(vacancyIndex);
     
     G4DataVector dataSet = *((*pos).second);
   
     G4int nData = dataSet.size();
     // The first Element of idMap's dataSets is the original shell of 
     // the vacancy, so we must start from the first element of dataSet
     if (initIndex >= 0 && initIndex < nData)
       {
     n =  (G4int) dataSet[initIndex+1];   
       }
   }
 return n;
}

References FatalErrorInArgument, G4Exception(), idMap, CLHEP::detail::n, numberOfVacancies, and pos.

Referenced by G4AtomicTransitionManager::Initialise(), and PrintData().

StartShellProb()

G4double G4FluoData::StartShellProb	(	G4int	initIndex,
		G4int	vacancyIndex
	)		const

Given the indexes of the starting and final shells for the transition, returns the probability of this transition

Definition at line 177 of file G4FluoData.cc.

{
  G4double n = -1;
 
  if (vacancyIndex<0 || vacancyIndex>=numberOfVacancies)
    {
      G4Exception("G4FluoData::StartShellEnergy()","de0002",JustWarning,
          "vacancyIndex outside boundaries, energy deposited locally");
      return 0;
    }
  else
    {
      auto pos = probabilityMap.find(vacancyIndex);
     
      G4DataVector dataSet = *((*pos).second);
     
      G4int nData = dataSet.size();
      if (initIndex >= 0 && initIndex < nData)
    {
      n =  dataSet[initIndex];
    }
    }
  return n;
}

References G4Exception(), JustWarning, CLHEP::detail::n, numberOfVacancies, pos, and probabilityMap.

Referenced by G4AtomicTransitionManager::Initialise(), and PrintData().

VacancyId()

G4int G4FluoData::VacancyId

(

G4int

vacancyIndex

)

const

Given the index of the vacancy returns its identity.

Definition at line 81 of file G4FluoData.cc.

{
  G4int n = -1;
  if (vacancyIndex<0 || vacancyIndex>=numberOfVacancies)
    {
      G4Exception("G4FluoData::vacancyId()","de0002",FatalErrorInArgument,
          "vacancyIndex outside boundaries");
    }
  else
    {
      auto pos = idMap.find(vacancyIndex);
      if (pos!= idMap.end())
    { G4DataVector dataSet = (*(*pos).second);
    n = (G4int) dataSet[0];
    
    }
    }
  return n;
}

References FatalErrorInArgument, G4Exception(), idMap, CLHEP::detail::n, numberOfVacancies, and pos.

Referenced by G4AtomicTransitionManager::Initialise().

Field Documentation

energyMap

std::map<G4int,G4DataVector*,std::less<G4int> > G4FluoData::energyMap

private

Definition at line 91 of file G4FluoData.hh.

Referenced by LoadData(), StartShellEnergy(), and ~G4FluoData().

fluoDirectory

G4String G4FluoData::fluoDirectory

private

Definition at line 95 of file G4FluoData.hh.

Referenced by G4FluoData(), and LoadData().

fluoTransitionTable

std::map<G4int,std::vector<G4FluoTransition*>,std::less<G4int> > G4FluoData::fluoTransitionTable

private

Definition at line 94 of file G4FluoData.hh.

idMap

std::map<G4int,G4DataVector*,std::less<G4int> > G4FluoData::idMap

private

Definition at line 90 of file G4FluoData.hh.

Referenced by LoadData(), StartShellId(), VacancyId(), and ~G4FluoData().

nInitShells

std::vector<G4int> G4FluoData::nInitShells

private

Definition at line 93 of file G4FluoData.hh.

Referenced by LoadData(), and NumberOfTransitions().

numberOfVacancies

G4int G4FluoData::numberOfVacancies = 0

private

Definition at line 97 of file G4FluoData.hh.

Referenced by LoadData(), NumberOfTransitions(), NumberOfVacancies(), PrintData(), StartShellEnergy(), StartShellId(), StartShellProb(), and VacancyId().

probabilityMap

std::map<G4int,G4DataVector*,std::less<G4int> > G4FluoData::probabilityMap

private

Definition at line 92 of file G4FluoData.hh.

Referenced by LoadData(), StartShellProb(), and ~G4FluoData().

---

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

• source/processes/electromagnetic/lowenergy/include/G4FluoData.hh
• source/processes/electromagnetic/lowenergy/src/G4FluoData.cc