G4FluoData Class Reference

#include <G4FluoData.hh>

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

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;
}

~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;
    }
}

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());
 
  
  const char* path = G4FindDataDir("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;
            G4int n = (G4int)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;
}

Referenced by G4AtomicTransitionManager::Initialise().

NumberOfTransitions()

std::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;
}

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

NumberOfVacancies()

std::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;
}

Referenced by G4AtomicTransitionManager::Initialise().

operator=()

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

delete

PrintData()

void G4FluoData::PrintData

(

)

Definition at line 335 of file G4FluoData.cc.

{
  for (G4int i = 0; i <numberOfVacancies; ++i)
    {
      G4cout << "---- TransitionData for the vacancy nb "
         << i
         << " ----- "
         << G4endl;
      
      for (G4int k = 0; k<(G4int)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;
    }
}

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 = (G4int)dataSet.size();
      if (initIndex >= 0 && initIndex < nData)
    {
          n =  dataSet[initIndex];
    }
    }
  return n;
}

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 = (G4int)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 = dataSet[initIndex+1];   
       }
   }
 return n;
}

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 = (G4int)dataSet.size();
      if (initIndex >= 0 && initIndex < nData)
    {
          n =  dataSet[initIndex];
    }
    }
  return n;
}

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;
}

Referenced by G4AtomicTransitionManager::Initialise().

---

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

• G4FluoData.hh
• G4FluoData.cc