G4RDAtomicTransitionManager Class Reference

#include <G4RDAtomicTransitionManager.hh>

Public Member Functions
G4RDAtomicShell *	Shell (G4int Z, size_t shellIndex) const
const G4RDFluoTransition *	ReachableShell (G4int Z, size_t shellIndex) const
const G4RDAugerTransition *	ReachableAugerShell (G4int Z, G4int shellIndex) const
G4int	NumberOfShells (G4int Z) const
G4int	NumberOfReachableShells (G4int Z) const
G4int	NumberOfReachableAugerShells (G4int Z) const
G4double	TotalRadiativeTransitionProbability (G4int Z, size_t shellIndex)
G4double	TotalNonRadiativeTransitionProbability (G4int Z, size_t shellIndex)

Static Public Member Functions
static G4RDAtomicTransitionManager *	Instance ()

Protected Member Functions
	G4RDAtomicTransitionManager (G4int minZ=1, G4int maxZ=100, G4int limitInfTable=6, G4int limitSupTable=100)
	~G4RDAtomicTransitionManager ()

Detailed Description

Definition at line 61 of file G4RDAtomicTransitionManager.hh.

Constructor & Destructor Documentation

G4RDAtomicTransitionManager::G4RDAtomicTransitionManager ( G4int  minZ = 1, G4int  maxZ = 100, G4int  limitInfTable = 6, G4int  limitSupTable = 100  )

protected

Definition at line 41 of file G4RDAtomicTransitionManager.cc.

  :zMin(minZ), 
  zMax(maxZ),
  infTableLimit(limitInfTable),
  supTableLimit(limitSupTable)
{
  // infTableLimit is initialized to 6 because EADL lacks data for Z<=5
  G4RDShellData* shellManager = new G4RDShellData;

  // initialization of the data for auger effect
  
  augerData = new G4RDAugerData;

  shellManager->LoadData("/fluor/binding");
  
  // Fills shellTable with the data from EADL, identities and binding 
  // energies of shells
  for (G4int Z = zMin; Z<= zMax; Z++)
    {
      std::vector<G4RDAtomicShell*> vectorOfShells;  
      size_t shellIndex = 0; 

      size_t numberOfShells=shellManager->NumberOfShells(Z);
      for (shellIndex = 0; shellIndex<numberOfShells; shellIndex++) 
    { 
      G4int shellId = shellManager->ShellId(Z,shellIndex);
      G4double bindingEnergy = shellManager->BindingEnergy(Z,shellIndex);
     
      G4RDAtomicShell * shell = new G4RDAtomicShell(shellId,bindingEnergy);
    
      vectorOfShells.push_back(shell);
    }
    
      //     shellTable.insert(std::make_pair(Z, vectorOfShells));
      shellTable[Z] = vectorOfShells;
    }
  
  // Fills transitionTable with the data from EADL, identities, transition 
  // energies and transition probabilities
  for (G4int Znum= infTableLimit; Znum<=supTableLimit; Znum++)
    {  G4RDFluoData* fluoManager = new G4RDFluoData;
    std::vector<G4RDFluoTransition*> vectorOfTransitions;
    fluoManager->LoadData(Znum);
    
    size_t numberOfVacancies = fluoManager-> NumberOfVacancies();
    
    for (size_t vacancyIndex = 0; vacancyIndex<numberOfVacancies;  vacancyIndex++)
      
      {
    std::vector<G4int>  vectorOfIds;
    G4DataVector vectorOfEnergies;
    G4DataVector vectorOfProbabilities;
    
    G4int finalShell = fluoManager->VacancyId(vacancyIndex);
    size_t numberOfTransitions = fluoManager->NumberOfTransitions(vacancyIndex);
    for (size_t origShellIndex = 0; origShellIndex < numberOfTransitions;
         origShellIndex++)
        
      {
        
        G4int originatingShellId = fluoManager->StartShellId(origShellIndex,vacancyIndex);
        
        vectorOfIds.push_back(originatingShellId);
        
        G4double transitionEnergy = fluoManager->StartShellEnergy(origShellIndex,vacancyIndex);
        vectorOfEnergies.push_back(transitionEnergy);
        G4double transitionProbability = fluoManager->StartShellProb(origShellIndex,vacancyIndex);
        vectorOfProbabilities.push_back(transitionProbability);
      }
      G4RDFluoTransition * transition = new G4RDFluoTransition (finalShell,vectorOfIds,
                                    vectorOfEnergies,vectorOfProbabilities);
      vectorOfTransitions.push_back(transition); 
      }
    //      transitionTable.insert(std::make_pair(Znum, vectorOfTransitions));
    transitionTable[Znum] = vectorOfTransitions;
      
      delete fluoManager;
    }
  delete shellManager;
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4RDAtomicTransitionManager::~G4RDAtomicTransitionManager ( )

protected

Definition at line 123 of file G4RDAtomicTransitionManager.cc.

{ 

  delete augerData;

std::map<G4int,std::vector<G4RDAtomicShell*>,std::less<G4int> >::iterator pos;
 
 for (pos = shellTable.begin(); pos != shellTable.end(); pos++){
   
   std::vector< G4RDAtomicShell*>vec = (*pos).second;
   
   G4int vecSize=vec.size();
   
   for (G4int i=0; i< vecSize; i++){
     G4RDAtomicShell* shell = vec[i];
     delete shell;     
   }
   
 }
 
 std::map<G4int,std::vector<G4RDFluoTransition*>,std::less<G4int> >::iterator ppos;
 
 for (ppos = transitionTable.begin(); ppos != transitionTable.end(); ppos++){
   
   std::vector<G4RDFluoTransition*>vec = (*ppos).second;
   
   G4int vecSize=vec.size();
   
   for (G4int i=0; i< vecSize; i++){
     G4RDFluoTransition* transition = vec[i];
     delete transition;      
   }
   
 }   
 
}

Member Function Documentation

G4RDAtomicTransitionManager * G4RDAtomicTransitionManager::Instance ( void  )

static

Definition at line 163 of file G4RDAtomicTransitionManager.cc.

{
  if (instance == 0)
    {
      instance = new G4RDAtomicTransitionManager;
     
    }
  return instance;
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4int G4RDAtomicTransitionManager::NumberOfReachableAugerShells

(

G4int

Z

)

const

Definition at line 298 of file G4RDAtomicTransitionManager.cc.

{
  G4int n = augerData->NumberOfVacancies(Z);
  return n;
}

Here is the call graph for this function:

G4int G4RDAtomicTransitionManager::NumberOfReachableShells

(

G4int

Z

)

const

Definition at line 274 of file G4RDAtomicTransitionManager.cc.

{
std::map<G4int,std::vector<G4RDFluoTransition*>,std::less<G4int> >::const_iterator pos;

  pos = transitionTable.find(Z);

  if (pos!= transitionTable.end())
    {
      std::vector<G4RDFluoTransition*> v = (*pos).second;
      return v.size();
    }
  else
    {
      G4cout << "G4AtomicTransitionMagare warning: No fluorescence or Auger for Z=" << Z << G4endl;
      G4cout << "Absorbed enrgy deposited locally" << G4endl;

      //    G4Exception("G4RDAtomicTransitionManager:Z not found");
      return 0;
    } 
}

G4int G4RDAtomicTransitionManager::NumberOfShells

(

G4int

Z

)

const

Definition at line 248 of file G4RDAtomicTransitionManager.cc.

{

std::map<G4int,std::vector<G4RDAtomicShell*>,std::less<G4int> >::const_iterator pos;

  pos = shellTable.find(Z);

  if (pos!= shellTable.end()){

    std::vector<G4RDAtomicShell*> v = (*pos).second;

    return v.size();
  }

  else{
    G4cout << "G4AtomicTransitionMagare warning: No fluorescence or Auger for Z=" << Z << G4endl;
    G4cout << "Absorbed enrgy deposited locally" << G4endl;

    //    G4Exception("G4RDAtomicTransitionManager:Z not found");
    return 0;
  } 
}

const G4RDAugerTransition * G4RDAtomicTransitionManager::ReachableAugerShell	(	G4int	Z,
		G4int	shellIndex
	)		const

Definition at line 239 of file G4RDAtomicTransitionManager.cc.

{
  
  G4RDAugerTransition* augerTransition = augerData->GetAugerTransition(Z,vacancyShellIndex);
  return augerTransition;
}

Here is the call graph for this function:

const G4RDFluoTransition * G4RDAtomicTransitionManager::ReachableShell	(	G4int	Z,
		size_t	shellIndex
	)		const

Definition at line 215 of file G4RDAtomicTransitionManager.cc.

{
  std::map<G4int,std::vector<G4RDFluoTransition*>,std::less<G4int> >::const_iterator pos;
  pos = transitionTable.find(Z);
  if (pos!= transitionTable.end())
    {
      std::vector<G4RDFluoTransition*> v = (*pos).second;      
      if (shellIndex < v.size()) return(v[shellIndex]);
      else {
    G4Exception("G4RDAtomicTransitionManager::ReachableShell()",
                    "InvalidCondition", FatalException,
                    "Reachable shell not found!");
    return 0;
      }
  }
  else{
    G4cout << "G4AtomicTransitionMagare warning: No fluorescence or Auger for Z=" << Z << G4endl;
    G4cout << "Absorbed enrgy deposited locally" << G4endl;

    //    G4Exception("G4RDAtomicTransitionManager:Z not found");
    return 0;
  } 
}

Here is the call graph for this function:

G4RDAtomicShell * G4RDAtomicTransitionManager::Shell	(	G4int	Z,
		size_t	shellIndex
	)		const

Definition at line 174 of file G4RDAtomicTransitionManager.cc.

{ 
  std::map<G4int,std::vector<G4RDAtomicShell*>,std::less<G4int> >::const_iterator pos;
  
  pos = shellTable.find(Z);
  
  if (pos!= shellTable.end())
    {
      std::vector<G4RDAtomicShell*> v = (*pos).second;
      if (shellIndex<v.size())
    {
      return(v[shellIndex]);
    }
      else 
    {
      size_t lastShell = v.size();
      G4cout << "G4RDAtomicTransitionManager::Shell - Z = " 
         << Z << ", shellIndex = " << shellIndex 
         << " not found; number of shells = " << lastShell << G4endl;
      //  G4Exception("G4RDAtomicTransitionManager:shell not found");
      if (lastShell > 0)
        {
          return v[lastShell - 1];
        }
      else
        {       
          return 0;
        }
    }
    }
  else
    {
      G4Exception("G4RDAtomicTransitionManager::Shell()",
                  "InvalidSetup", FatalException, "Z not found!");
      return 0;
    } 
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4RDAtomicTransitionManager::TotalNonRadiativeTransitionProbability	(	G4int	Z,
		size_t	shellIndex
	)

Definition at line 348 of file G4RDAtomicTransitionManager.cc.

{

  std::map<G4int,std::vector<G4RDFluoTransition*>,std::less<G4int> >::iterator pos;
  
  pos = transitionTable.find(Z);
  
  if (pos!= transitionTable.end()){
    
    std::vector<G4RDFluoTransition*> v = (*pos).second;
  
    
    if (shellIndex<v.size()){

      G4RDFluoTransition* transition=v[shellIndex];
      G4DataVector transProb = transition->TransitionProbabilities();
      G4double totalRadTransProb = 0;
      
      for(size_t j = 0; j<transProb.size(); j++) // AM -- Corrected, was 1
    {
      totalRadTransProb = totalRadTransProb + transProb[j];
    }
      
      G4double totalNonRadTransProb= (1 - totalRadTransProb);
      
      return totalNonRadTransProb;    }
    
    else {
      G4Exception("G4RDAtomicTransitionManager::TotalNonRadiativeTransitionProbability()",
                  "InvalidCondition", FatalException, "Shell not found!");
      return 0;
    }
  }
  else{
    G4cout << "G4AtomicTransitionMagare warning: No fluorescence or Auger for Z=" << Z << G4endl;
    G4cout << "Absorbed enrgy deposited locally" << G4endl;

    //    G4Exception("G4RDAtomicTransitionManager:Z not found");
    return 0;
  } 
}

Here is the call graph for this function:

G4double G4RDAtomicTransitionManager::TotalRadiativeTransitionProbability	(	G4int	Z,
		size_t	shellIndex
	)

Definition at line 306 of file G4RDAtomicTransitionManager.cc.

{
std::map<G4int,std::vector<G4RDFluoTransition*>,std::less<G4int> >::iterator pos;

  pos = transitionTable.find(Z);

  if (pos!= transitionTable.end())
    {
      std::vector<G4RDFluoTransition*> v = (*pos).second;
      
    if (shellIndex < v.size())
      {
    G4RDFluoTransition* transition = v[shellIndex];
    G4DataVector transProb = transition->TransitionProbabilities();
    G4double totalRadTransProb = 0;
    
    for (size_t j = 0; j<transProb.size(); j++) // AM -- corrected, it was 1
    {
      totalRadTransProb = totalRadTransProb + transProb[j];
    }
      return totalRadTransProb;   
      
    }
    else {
      G4Exception("G4RDAtomicTransitionManager::TotalRadiativeTransitionProbability()",
                  "InvalidCondition", FatalException, "Shell not found!" );
      return 0;
      
    }
  }
  else{
    G4cout << "G4AtomicTransitionMagare warning: No fluorescence or Auger for Z=" << Z << G4endl;
    G4cout << "Absorbed enrgy deposited locally" << G4endl;

    //    G4Exception("G4RDAtomicTransitionManager:Z not found");

    return 0;
  } 
}

Here is the call graph for this function:

---

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

• source/geant4.10.03.p03/examples/advanced/eRosita/physics/include/G4RDAtomicTransitionManager.hh
• source/geant4.10.03.p03/examples/advanced/eRosita/physics/src/G4RDAtomicTransitionManager.cc