G4AugerData Class Reference

#include <G4AugerData.hh>

Public Member Functions
	G4AugerData ()
	~G4AugerData ()
size_t	NumberOfVacancies (G4int Z) const
G4int	VacancyId (G4int Z, G4int vacancyIndex) const
size_t	NumberOfTransitions (G4int Z, G4int vacancyIndex) const
G4int	StartShellId (G4int Z, G4int initialVacancyIndex, G4int transitionShellIndex) const
G4double	StartShellEnergy (G4int Z, G4int vacancyIndex, G4int transitionId, G4int augerIndex) const
G4double	StartShellProb (G4int Z, G4int vacancyIndex, G4int transitionId, G4int augerIndex) const
size_t	NumberOfAuger (G4int Z, G4int initIndex, G4int vacancyId) const
size_t	AugerShellId (G4int Z, G4int vacancyIndex, G4int transId, G4int augerIndex) const
std::vector< G4AugerTransition >	LoadData (G4int Z)
void	BuildAugerTransitionTable ()
void	PrintData (G4int Z)
G4AugerTransition *	GetAugerTransition (G4int Z, G4int vacancyShellIndex)
std::vector< G4AugerTransition > *	GetAugerTransitions (G4int Z)

Detailed Description

Definition at line 55 of file G4AugerData.hh.

Constructor & Destructor Documentation

G4AugerData()

G4AugerData::G4AugerData

(

)

Definition at line 50 of file G4AugerData.cc.

{
 
  G4int n = 0;
  G4int pos = 0;
 
    for (pos = 0 ; pos < 100; pos++) 
      {
    numberOfVacancies.push_back(n);
      }
 
  BuildAugerTransitionTable(); 
 
 
}

~G4AugerData()

G4AugerData::~G4AugerData

(

)

Definition at line 66 of file G4AugerData.cc.

{ 
  /*
  std::map<G4int,std::vector<G4AugerTransition>,std::less<G4int> >::iterator pos;
 
  for (pos = augerTransitionTable.begin(); pos != augerTransitionTable.end(); pos++)
    {
      std::vector<G4AugerTransition> dataSet = (*pos).second;
      delete dataSet;
    }
  for (pos = energyMap.begin(); pos != energyMap.end(); pos++)
    {
      std::map<G4Int,G4DataVector*,std::less<G4int>>* dataMap = (*pos).second;
      for (pos2 = newIdProbabilityMap.begin(); pos2 != idMap.end(); pos2++)
    {
      G4DataVector* dataSet = (*pos2).second;
      delete dataSet;
    }
    }
  for (pos = probabilityMap.begin(); pos != probabilityMap.end(); pos++)
    {
      std::map<G4Int,G4DataVector*,std::less<G4int>>* dataMap = (*pos).second;
      for (pos2 = newIdProbabilityMap.begin(); pos2 != idMap.end(); pos2++)
    {
      G4DataVector* dataSet = (*pos2).second;
      delete dataSet;
    }
    }
  for (pos2 = newIdMap.begin(); pos2 != idMap.end(); pos2++)
    {
      G4DataVector* dataSet = (*pos2).second;
      delete dataSet;
    }
  for (pos2 = newIdEnergyMap.begin(); pos2 != idMap.end(); pos2++)
    {
      G4DataVector* dataSet = (*pos2).second;
      delete dataSet;
    }
  for (pos2 = newIdProbabilityMap.begin(); pos2 != idMap.end(); pos2++)
    {
      G4DataVector* dataSet = (*pos2).second;
      delete dataSet;
    }
  */
 
}

Member Function Documentation

AugerShellId()

size_t G4AugerData::AugerShellId	(	G4int	Z,
		G4int	vacancyIndex,
		G4int	transId,
		G4int	augerIndex
	)		const

Definition at line 186 of file G4AugerData.cc.

{
  size_t n = 0;  
  if (vacancyIndex<0 || vacancyIndex>=numberOfVacancies[Z])
    {
      G4Exception("G4AugerData::VacancyId()","de0002", FatalErrorInArgument,"");
      return 0;
    }
  else {
    trans_Table::const_iterator element = augerTransitionTable.find(Z);
    if (element == augerTransitionTable.end()) {
      G4Exception("G4AugerData::VacancyId()","de0004", FatalErrorInArgument,  "Check element");
      return 0;
    }
    std::vector<G4AugerTransition> dataSet = (*element).second;
    n = dataSet[vacancyIndex].AugerOriginatingShellId(augerIndex,transId);
  }
  return n;
}

Referenced by PrintData().

BuildAugerTransitionTable()

void G4AugerData::BuildAugerTransitionTable

(

)

Definition at line 508 of file G4AugerData.cc.

{
 
  //  trans_Table::iterator pos = augerTransitionTable.begin();
 
  const G4MaterialTable* materialTable = G4Material::GetMaterialTable();
 
  G4int nMaterials = G4Material::GetNumberOfMaterials();
 
  G4DataVector activeZ;
  activeZ.clear();
  
  for (G4int mLocal=0; mLocal<nMaterials; mLocal++) {
 
    const G4Material* material= (*materialTable)[mLocal];        
    const G4ElementVector* elementVector = material->GetElementVector();
    const size_t nElements = material->GetNumberOfElements();
      
    for (size_t iEl=0; iEl<nElements; iEl++) {
      G4Element* element = (*elementVector)[iEl];
      G4double Z = element->GetZ();
      if (!(activeZ.contains(Z))) {
    activeZ.push_back(Z);
      }
    }
  }
 
 
  for (G4int element = 6; element < 100; element++)
    { 
      //      if(nMaterials == 0 || activeZ.contains(element)) {
      augerTransitionTable.insert(trans_Table::value_type(element,LoadData(element)));
      //    G4cout << "G4AugerData for Element no. " << element << " are loaded" << G4endl;
      // G4cout << "G4AugerData for Element no. " << element << " are loaded" << G4endl;
      //G4cout << "AugerTransitionTable complete"<< G4endl;
    }
}

Referenced by G4AugerData().

GetAugerTransition()

G4AugerTransition * G4AugerData::GetAugerTransition	(	G4int	Z,
		G4int	vacancyShellIndex
	)

Definition at line 577 of file G4AugerData.cc.

    {
      std::vector<G4AugerTransition>* dataSet = &augerTransitionTable[Z];
      std::vector<G4AugerTransition>::iterator vectorIndex = dataSet->begin() + vacancyShellIndex;
 
      G4AugerTransition* augerTransition = &(*vectorIndex);
      return augerTransition;
    }

Referenced by G4AtomicTransitionManager::ReachableAugerShell().

GetAugerTransitions()

std::vector< G4AugerTransition > * G4AugerData::GetAugerTransitions

(

G4int

Z

)

Definition at line 586 of file G4AugerData.cc.

  {
    std::vector<G4AugerTransition>* dataSet = &augerTransitionTable[Z];
    return dataSet;
  }

LoadData()

std::vector< G4AugerTransition > G4AugerData::LoadData

(

G4int

Z

)

Definition at line 276 of file G4AugerData.cc.

{ 
  // Build the complete string identifying the file with the data set
 
    std::ostringstream ost;
    if(Z != 0){
      ost << "au-tr-pr-"<< Z << ".dat";
    }
    else{
      ost << "au-tr-pr-"<<".dat"; 
    }
    G4String name(ost.str());
  
    char* path = getenv("G4LEDATA");
    if (!path)
      { 
    G4String excep = "G4AugerData::LoadData";
    G4Exception(excep,"em0006", FatalException,"" );
    std::vector<G4AugerTransition> a;
    return a;
      }
  
    G4String pathString(path);
    G4String dirFile = pathString + "/auger/" + name;
    std::ifstream file(dirFile);
    std::filebuf* lsdp = file.rdbuf();
  
    if (! (lsdp->is_open()) )
      {
    G4String excep = "G4AugerData::LoadData";
    G4String msg = "Missing" + dirFile;
    G4Exception(excep,"em0003", FatalException, msg);
      }
 
 
    G4double a = 0;
    G4int k = 1;
    G4int sLocal = 0;
  
    G4int vacId = 0;
    std::vector<G4int>* initIds = new std::vector<G4int>;
    std::vector<G4int>* newIds = new std::vector<G4int>;
    G4DataVector* transEnergies = new G4DataVector;
    G4DataVector* transProbabilities = new G4DataVector;
    std::vector<G4AugerTransition> augerTransitionVector;
    std::map<G4int,std::vector<G4int>,std::less<G4int> >* newIdMap = 
      new std::map<G4int,std::vector<G4int>,std::less<G4int> >;
    std::map<G4int,G4DataVector,std::less<G4int> >* newEnergyMap =
      new std::map<G4int,G4DataVector,std::less<G4int> >;
    std::map<G4int,G4DataVector,std::less<G4int> >* newProbabilityMap = 
      new std::map<G4int,G4DataVector,std::less<G4int> >;
 
  
    do {
      file >> a;
 
 
      G4int nColumns = 4;
 
      if (a == -1)
    {
 
 
 
      if (sLocal == 0)
        {
          // End of a shell data set
        
        
        
          std::vector<G4int>::iterator vectorIndex = initIds->begin();
 
          vacId = *vectorIndex;
        
          //initIds->erase(vectorIndex);
        
 
 
          std::vector<G4int> identifiers;
          for (vectorIndex = initIds->begin()+1 ; vectorIndex != initIds->end(); ++vectorIndex){
 
        identifiers.push_back(*vectorIndex);
          }
 
          vectorIndex = (initIds->end())-1;
 
          G4int augerShellId = *(vectorIndex);
        
        
          (*newIdMap)[augerShellId] = *newIds;
          (*newEnergyMap)[augerShellId] = *transEnergies;
          (*newProbabilityMap)[augerShellId] = *transProbabilities;
 
          augerTransitionVector.push_back(G4AugerTransition(vacId, identifiers, newIdMap, newEnergyMap, newProbabilityMap));
 
          // Now deleting all the variables I used, and creating new ones for the next shell
 
          delete newIdMap;
          delete newEnergyMap;
          delete newProbabilityMap;
        
          G4int n = initIds->size();        
          nInitShells.push_back(n);
          numberOfVacancies[Z]++;
          delete initIds;
          delete newIds;
          delete transEnergies;     
          delete transProbabilities;
          initIds = new std::vector<G4int>;
          newIds = new std::vector<G4int>;
          transEnergies = new G4DataVector;
          transProbabilities = new G4DataVector;
          newIdMap = new std::map<G4int,std::vector<G4int>,std::less<G4int> >;
          newEnergyMap = new std::map<G4int,G4DataVector,std::less<G4int> >;
          newProbabilityMap = new std::map<G4int,G4DataVector,std::less<G4int> >;   
        
 
 
        }      
      sLocal++;
      if (sLocal == nColumns)
        {
          sLocal = 0;
        }
    }
      // moved to the end in order to avoid leaks
      /*
      else if (a == -2)
    {
      // End of file; delete the empty vectors created 
      //when encountering the last -1 -1 row
      delete initIds;
      delete newIds;
      delete transEnergies;
      delete transProbabilities;
      delete newIdMap ;
      delete newEnergyMap;
      delete newProbabilityMap; 
      }*/ 
      else
    {
    
      if (k%nColumns == 3){
        // 3rd column is the transition probabilities
        transProbabilities->push_back(a);
 
        k++;}
      else if(k%nColumns == 2){  
        // 2nd column is new auger vacancy
 
        // 2nd column is new auger vacancy
 
        G4int l = (G4int)a;
        newIds->push_back(l);
 
 
        k++;
      }
      else if (k%nColumns == 1)
        {
          // 1st column is shell id
        
          if(initIds->size() == 0) {
          
        // 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
          
        initIds->push_back((G4int)a);
        // first line of initIds is the original shell of the vacancy
        file >> a;
        file >> a;
        file >> a;
        k = k+3;
          }
          else {
 
        //        std::vector<G4int>::iterator vectorIndex = (initIds->end())-1;
        if((G4int)a != initIds->back()){
 
 
          if((initIds->size()) == 1) { 
            initIds->push_back((G4int)a);
          }  
          else {
 
 
            G4int augerShellId = 0;
            augerShellId = initIds->back();
          
            (*newIdMap)[augerShellId] = *newIds;
            (*newEnergyMap)[augerShellId] = *transEnergies;
            (*newProbabilityMap)[augerShellId] = *transProbabilities;
            delete newIds;
            delete transEnergies;
            delete transProbabilities;
            newIds = new std::vector<G4int>;
            transEnergies = new G4DataVector;
            transProbabilities = new G4DataVector;
            initIds->push_back((G4int)a);
          }
        }
          }
        
          k++;    
       
        }
      else if (k%nColumns == 0)
 
        {//fourth column is transition energies
          G4double e = a * MeV; 
  
          transEnergies->push_back(e);
          k=1;
 
        }
    }
    }
 
 
    while (a != -2); // end of file
    file.close();
    delete initIds;
    delete newIds;
    delete transEnergies;
    delete transProbabilities;
    delete newIdMap ;
    delete newEnergyMap;
    delete newProbabilityMap;
    return augerTransitionVector;
 
}

Referenced by BuildAugerTransitionTable().

NumberOfAuger()

size_t G4AugerData::NumberOfAuger	(	G4int	Z,
		G4int	initIndex,
		G4int	vacancyId
	)		const

Definition at line 165 of file G4AugerData.cc.

{
  size_t n = 0;
  if (initIndex<0 || initIndex>=numberOfVacancies[Z])
    {
      G4Exception("G4AugerData::VacancyId()","de0002", FatalErrorInArgument,"");
      return 0;
    }
  else {
    trans_Table::const_iterator element = augerTransitionTable.find(Z);
    if (element == augerTransitionTable.end()) {
      G4Exception("G4AugerData::VacancyId()","de0004", FatalErrorInArgument,  "Check element");
      return 0;
    }
    std::vector<G4AugerTransition> dataSet = (*element).second;
    const std::vector<G4int>* temp =  dataSet[initIndex].AugerOriginatingShellIds(vacancyId);
    n = temp->size();
  }
  return n;
}

Referenced by PrintData().

NumberOfTransitions()

size_t G4AugerData::NumberOfTransitions	(	G4int	Z,
		G4int	vacancyIndex
	)		const

Definition at line 142 of file G4AugerData.cc.

{
  G4int n = 0;
  if (vacancyIndex<0 || vacancyIndex>=numberOfVacancies[Z])
    {
      G4Exception("G4AugerData::VacancyId()","de0002", JustWarning, "Energy deposited locally");
      return 0;
    }
  else {
    trans_Table::const_iterator element = augerTransitionTable.find(Z);
    if (element == augerTransitionTable.end()) 
      {
    G4Exception("G4AugerData::VacancyId()","de0004", FatalErrorInArgument,  "Check element");
    return 0;
      }
    std::vector<G4AugerTransition> dataSet = (*element).second;
    n = (G4int)dataSet[vacancyIndex].TransitionOriginatingShellIds()->size();
  }
 return  n;
}

Referenced by PrintData().

NumberOfVacancies()

size_t G4AugerData::NumberOfVacancies

(

G4int

Z

)

const

Definition at line 113 of file G4AugerData.cc.

{
  return numberOfVacancies[Z];
}

Referenced by G4AtomicTransitionManager::NumberOfReachableAugerShells().

PrintData()

void G4AugerData::PrintData

(

G4int

Z

)

Definition at line 546 of file G4AugerData.cc.

{
  
  for (G4int i = 0; i < numberOfVacancies[Z]; i++)
    {
      G4cout << "---- TransitionData for the vacancy nb "
         <<i
         <<" of the atomic number elemnt " 
         << Z
         <<"----- "
         <<G4endl;
      
      for (size_t k = 0; k<=NumberOfTransitions(Z,i); k++)
    { 
      G4int id = StartShellId(Z,i,k);
      
      for (size_t a = 0; a <= NumberOfAuger(Z,i,id); a++) {
        
        G4double e = StartShellEnergy(Z,i,id,a)/MeV;
        G4double p = StartShellProb(Z,i,id,a);
        G4int augerId = AugerShellId(Z, i, id, a);
        G4cout << k <<") Shell id: " << id <<G4endl;
        G4cout << "    Auger Originatig Shell Id :"<< augerId <<G4endl;
        G4cout << " - Transition energy = " << e << " MeV "<<G4endl;
        G4cout   << " - Transition probability = " << p <<G4endl;
      }
    }
      G4cout << "-------------------------------------------------" 
         << G4endl;
    }
}

StartShellEnergy()

G4double G4AugerData::StartShellEnergy	(	G4int	Z,
		G4int	vacancyIndex,
		G4int	transitionId,
		G4int	augerIndex
	)		const

Definition at line 229 of file G4AugerData.cc.

{
  G4double energy = 0;
  
  if (vacancyIndex<0 || vacancyIndex>=numberOfVacancies[Z])
    {
      G4Exception("G4AugerData::VacancyId()","de0002", FatalErrorInArgument,"");
      return 0;
    }
  else {
    trans_Table::const_iterator element = augerTransitionTable.find(Z);
    if (element == augerTransitionTable.end()) {
      G4Exception("G4AugerData::VacancyId()","de0004", FatalErrorInArgument,  "Check element");
      return 0;
    }
    std::vector<G4AugerTransition> dataSet = (*element).second;
    energy = dataSet[vacancyIndex].AugerTransitionEnergy(augerIndex,transitionId);
      
  }
  return energy;
}

Referenced by PrintData().

StartShellId()

G4int G4AugerData::StartShellId	(	G4int	Z,
		G4int	initialVacancyIndex,
		G4int	transitionShellIndex
	)		const

Definition at line 206 of file G4AugerData.cc.

{
  G4int n = 0; 
 
   if (vacancyIndex<0 || vacancyIndex>=numberOfVacancies[Z]) 
     {
       G4Exception("G4AugerData::VacancyId()","de0002", FatalErrorInArgument,"");
       return 0;
     }
  else {
    trans_Table::const_iterator element = augerTransitionTable.find(Z);
    if (element == augerTransitionTable.end()) {
      G4Exception("G4AugerData::VacancyId()","de0004", FatalErrorInArgument,  "Check element");
      return 0; 
    }
    std::vector<G4AugerTransition> dataSet = (*element).second;
    n = dataSet[vacancyIndex].TransitionOriginatingShellId(transitionShellIndex);
  }
   
 
 return n;
}

Referenced by PrintData().

StartShellProb()

G4double G4AugerData::StartShellProb	(	G4int	Z,
		G4int	vacancyIndex,
		G4int	transitionId,
		G4int	augerIndex
	)		const

Definition at line 252 of file G4AugerData.cc.

{
  G4double prob = 0;
    
    if (vacancyIndex<0 || vacancyIndex>=numberOfVacancies[Z]) 
      {
    G4Exception("G4AugerData::VacancyId()","de0002", FatalErrorInArgument,"");
    return 0;
      }
  else {
    trans_Table::const_iterator element = augerTransitionTable.find(Z);
    if (element == augerTransitionTable.end()) {
      G4Exception("G4AugerData::VacancyId()","de0004", FatalErrorInArgument,  "Check element");
      return 0;
    }
    std::vector<G4AugerTransition> dataSet = (*element).second;
    prob = dataSet[vacancyIndex].AugerTransitionProbability(augerIndex, transitionId);
 
 
 
  }
     return prob;
}

Referenced by PrintData().

VacancyId()

G4int G4AugerData::VacancyId	(	G4int	Z,
		G4int	vacancyIndex
	)		const

Definition at line 118 of file G4AugerData.cc.

{
 
  G4int n = 0;
  if (vacancyIndex<0 || vacancyIndex>=numberOfVacancies[Z])
    {
      G4Exception("G4AugerData::VacancyId()","de0002", FatalErrorInArgument,"");
    }
  else {
    trans_Table::const_iterator element = augerTransitionTable.find(Z);
    if (element == augerTransitionTable.end()) {
      G4Exception("G4AugerData::VacancyId()","de0004", FatalErrorInArgument,  "Check element");
      return 0;
    }
    std::vector<G4AugerTransition> dataSet = (*element).second;
    n = (G4int) dataSet[vacancyIndex].FinalShellId();
  }
  
  return n;
}

---

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

• G4AugerData.hh
• G4AugerData.cc