G4AtomicTransitionManager Class Reference

#include <G4AtomicTransitionManager.hh>

Public Member Functions
void	Initialise ()
G4AtomicShell *	Shell (G4int Z, size_t shellIndex) const
const G4FluoTransition *	ReachableShell (G4int Z, size_t shellIndex) const
const G4AugerTransition *	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) const
G4double	TotalNonRadiativeTransitionProbability (G4int Z, size_t shellIndex) const
void	SetVerboseLevel (G4int vl)
G4int	GetVerboseLevel ()

Static Public Member Functions
static G4AtomicTransitionManager *	Instance ()

Detailed Description

Definition at line 59 of file G4AtomicTransitionManager.hh.

Member Function Documentation

GetVerboseLevel()

G4int G4AtomicTransitionManager::GetVerboseLevel ( )

inline

Definition at line 115 of file G4AtomicTransitionManager.hh.

{return verboseLevel;};

Initialise()

void G4AtomicTransitionManager::Initialise

(

)

Definition at line 285 of file G4AtomicTransitionManager.cc.

{
  G4AutoLock l(&AtomicTransitionManagerMutex);
 
  //G4cout << "!!! G4AtomicTransitionManager::Initialise " << isInitialized 
  // << G4endl;
  if(isInitialized) { return; }
  isInitialized = true;
 
  // Selection of fluorescence files
  G4String fluoDirectory = (G4EmParameters::Instance()->BeardenFluoDir()?
                "/fluor_Bearden":"/fluor");
 
  // infTableLimit is initialized to 6 because EADL lacks data for Z<=5
  G4ShellData* shellManager = new G4ShellData;
  shellManager->LoadData(fluoDirectory+"/binding");
 
  // initialization of the data for auger effect
  augerData = new G4AugerData;
  
  // Fills shellTable with the data from EADL, identities and binding 
  // energies of shells
  for (G4int Z = zMin; Z<= zMax; ++Z) 
    {
      std::vector<G4AtomicShell*> vectorOfShells;  
      size_t shellIndex = 0; 
 
      size_t numberOfShells = shellManager->NumberOfShells(Z);
      // G4cout << "For Z= " << Z << "  " << numberOfShells << " shells" << G4endl;
      for (shellIndex = 0; shellIndex<numberOfShells; ++shellIndex) 
    { 
      G4int shellId = shellManager->ShellId(Z,shellIndex);
      G4double bindingEnergy = shellManager->BindingEnergy(Z,shellIndex);
     
      G4AtomicShell * shell = new G4AtomicShell(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)
    {  
      G4FluoData* fluoManager = new G4FluoData(fluoDirectory);
      std::vector<G4FluoTransition*> 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);
        }
      G4FluoTransition* transition = 
        new G4FluoTransition (finalShell,vectorOfIds,
                  vectorOfEnergies,vectorOfProbabilities);
      vectorOfTransitions.push_back(transition); 
    }
      transitionTable[Znum] = vectorOfTransitions;
      delete fluoManager;
    }
  delete shellManager;
}

Referenced by G4LivermoreIonisationModel::Initialise(), and G4UAtomicDeexcitation::InitialiseForNewRun().

Instance()

G4AtomicTransitionManager * G4AtomicTransitionManager::Instance ( )

static

Definition at line 44 of file G4AtomicTransitionManager.cc.

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

Referenced by G4eIonisationSpectrum::AverageEnergy(), G4ecpssrBaseKxsModel::CalculateCrossSection(), G4ecpssrBaseLixsModel::CalculateL1CrossSection(), G4ecpssrBaseLixsModel::CalculateL2CrossSection(), G4ecpssrBaseLixsModel::CalculateL3CrossSection(), G4ecpssrBaseLixsModel::CalculateVelocity(), G4LivermoreIonisationCrossSection::G4LivermoreIonisationCrossSection(), G4LivermoreIonisationModel::G4LivermoreIonisationModel(), G4OrlicLiXsModel::G4OrlicLiXsModel(), G4PenelopeComptonModel::G4PenelopeComptonModel(), G4PenelopeIonisationCrossSection::G4PenelopeIonisationCrossSection(), G4PenelopePhotoElectricModel::G4PenelopePhotoElectricModel(), G4UAtomicDeexcitation::G4UAtomicDeexcitation(), G4hImpactIonisation::PostStepDoIt(), G4eIonisationSpectrum::Probability(), G4eIonisationSpectrum::SampleEnergy(), G4PenelopeIonisationModel::SampleSecondaries(), and G4PenelopePhotoElectricModel::SampleSecondaries().

NumberOfReachableAugerShells()

G4int G4AtomicTransitionManager::NumberOfReachableAugerShells

(

G4int

Z

)

const

Definition at line 216 of file G4AtomicTransitionManager.cc.

{
  return augerData->NumberOfVacancies(Z);
}

NumberOfReachableShells()

G4int G4AtomicTransitionManager::NumberOfReachableShells

(

G4int

Z

)

const

Definition at line 192 of file G4AtomicTransitionManager.cc.

{
  std::map<G4int,std::vector<G4FluoTransition*>,std::less<G4int> >::const_iterator pos;
  pos = transitionTable.find(Z);
  G4int res = 0;
  if (pos!= transitionTable.end())
    {
      res = ((*pos).second).size();
    }
  else
    {
      G4ExceptionDescription ed;
      ed << "No deexcitation for Z= " << Z
     << ", so energy deposited locally";
      G4Exception("G4AtomicTransitionManager::NumberOfReachebleShells()",
          "de0001",FatalException,ed,"");
    } 
  return res;
}

NumberOfShells()

G4int G4AtomicTransitionManager::NumberOfShells

(

G4int

Z

)

const

Definition at line 170 of file G4AtomicTransitionManager.cc.

{
  std::map<G4int,std::vector<G4AtomicShell*>,std::less<G4int> >::const_iterator pos;
  pos = shellTable.find(Z);
 
  G4int res = 0;
  if (pos != shellTable.end()){
 
    res = ((*pos).second).size();
 
  } else {
    G4ExceptionDescription ed;
    ed << "No deexcitation for Z= " << Z;
    G4Exception("G4AtomicTransitionManager::NumberOfShells()","de0001",
        FatalException, ed, "");
  } 
  return res;
}

Referenced by G4LivermoreIonisationModel::ComputeDEDXPerVolume(), G4PenelopeIonisationCrossSection::CrossSection(), G4LivermoreIonisationCrossSection::CrossSection(), G4PenelopeIonisationCrossSection::GetCrossSection(), G4LivermoreIonisationCrossSection::GetCrossSection(), and G4PenelopePhotoElectricModel::SampleSecondaries().

ReachableAugerShell()

const G4AugerTransition * G4AtomicTransitionManager::ReachableAugerShell	(	G4int	Z,
		G4int	shellIndex
	)		const

Definition at line 164 of file G4AtomicTransitionManager.cc.

{
  return augerData->GetAugerTransition(Z,vacancyShellIndex);
}

ReachableShell()

const G4FluoTransition * G4AtomicTransitionManager::ReachableShell	(	G4int	Z,
		size_t	shellIndex
	)		const

Definition at line 135 of file G4AtomicTransitionManager.cc.

{
  std::map<G4int,std::vector<G4FluoTransition*>,std::less<G4int> >::const_iterator pos;
  pos = transitionTable.find(Z);
  if (pos!= transitionTable.end())
    {
      std::vector<G4FluoTransition*> v = (*pos).second;      
      if (shellIndex < v.size()) { return(v[shellIndex]); }
 
      else {
    G4ExceptionDescription ed;
    ed << "No fluo transition for Z= " << Z 
       << "  shellIndex= " << shellIndex;
    G4Exception("G4AtomicTransitionManager::ReachebleShell()","de0002", 
            FatalException,ed,"");
      }
    } 
  else 
    {
      G4ExceptionDescription ed;
      ed << "No transition table for Z= " << Z 
     << "  shellIndex= " << shellIndex;
      G4Exception("G4AtomicTransitionManager::ReachableShell()","de0001",
          FatalException,ed,"");
    } 
  return 0;
}

SetVerboseLevel()

void G4AtomicTransitionManager::SetVerboseLevel ( G4int  vl )

inline

Definition at line 114 of file G4AtomicTransitionManager.hh.

{verboseLevel = vl;};

Shell()

G4AtomicShell * G4AtomicTransitionManager::Shell	(	G4int	Z,
		size_t	shellIndex
	)		const

Definition at line 94 of file G4AtomicTransitionManager.cc.

{ 
  std::map<G4int,std::vector<G4AtomicShell*>,std::less<G4int> >::const_iterator pos;
  
  pos = shellTable.find(Z);
  
  if (pos!= shellTable.end())
    {
      std::vector<G4AtomicShell*> v = (*pos).second;
      if (shellIndex < v.size()) { return v[shellIndex]; }
 
      else 
    {
      size_t lastShell = v.size();
      G4ExceptionDescription ed;
      ed << "No de-excitation for Z= " << Z 
         << "  shellIndex= " << shellIndex
         << ">=  numberOfShells= " << lastShell;
      if (verboseLevel > 0) 
              G4Exception("G4AtomicTransitionManager::Shell()","de0001",
        JustWarning,ed," AtomicShell not found");
      if (lastShell > 0) { return v[lastShell - 1]; }
    }
    }
  else
    {
      G4ExceptionDescription ed;
      ed << "No de-excitation for Z= " << Z 
     << "  shellIndex= " << shellIndex
     << ". AtomicShell not found - check if data are uploaded";
      G4Exception("G4AtomicTransitionManager::Shell()","de0001",
          FatalException,ed,"");
    } 
  return 0;
}

Referenced by G4ecpssrBaseKxsModel::CalculateCrossSection(), G4OrlicLiXsModel::CalculateL1CrossSection(), G4ecpssrBaseLixsModel::CalculateL1CrossSection(), G4OrlicLiXsModel::CalculateL2CrossSection(), G4ecpssrBaseLixsModel::CalculateL2CrossSection(), G4OrlicLiXsModel::CalculateL3CrossSection(), G4ecpssrBaseLixsModel::CalculateL3CrossSection(), G4ecpssrBaseLixsModel::CalculateVelocity(), G4UAtomicDeexcitation::GetAtomicShell(), G4hImpactIonisation::PostStepDoIt(), G4LivermoreIonisationModel::SampleSecondaries(), G4PenelopeComptonModel::SampleSecondaries(), G4PenelopeIonisationModel::SampleSecondaries(), and G4PenelopePhotoElectricModel::SampleSecondaries().

TotalNonRadiativeTransitionProbability()

G4double G4AtomicTransitionManager::TotalNonRadiativeTransitionProbability	(	G4int	Z,
		size_t	shellIndex
	)		const

Definition at line 265 of file G4AtomicTransitionManager.cc.

{
  G4double prob = 1.0 - TotalRadiativeTransitionProbability(Z, shellIndex);
  if(prob > 1.0 || prob < 0.0) {
    G4ExceptionDescription ed;
    ed << "Total probability mismatch Z= " << Z 
       << "  shellIndex= " << shellIndex
       << "  prob= " << prob;
    G4Exception( 
    "G4AtomicTransitionManager::TotalNonRadiativeTransitionProbability()",
    "de0003",FatalException,ed,"Cannot compute non-radiative probability");
    return 0.0;
  }
  return prob;    
}

TotalRadiativeTransitionProbability()

G4double G4AtomicTransitionManager::TotalRadiativeTransitionProbability	(	G4int	Z,
		size_t	shellIndex
	)		const

Definition at line 221 of file G4AtomicTransitionManager.cc.

{
  std::map<G4int,std::vector<G4FluoTransition*>,std::less<G4int> >::const_iterator pos;
 
  pos = transitionTable.find(Z);
  G4double totalRadTransProb = 0.0;
 
  if (pos!= transitionTable.end())
    {
      std::vector<G4FluoTransition*> v = (*pos).second;
      
      if (shellIndex < v.size())
    {
      G4FluoTransition* transition = v[shellIndex];
      G4DataVector transProb = transition->TransitionProbabilities();
    
      for (size_t j=0; j<transProb.size(); ++j) // AM -- corrected, it was 1
        {
          totalRadTransProb += transProb[j];
        }
    }
      else 
    {
      G4ExceptionDescription ed;
      ed << "Zero transition probability for Z=" << Z 
         << "  shellIndex= " << shellIndex;
      G4Exception(
      "G4AtomicTransitionManager::TotalRadiativeTransitionProbability()",
      "de0002",FatalException,"Incorrect de-excitation");
    }
    }
  else
    {
      G4ExceptionDescription ed;
      ed << "No deexcitation for Z=" << Z 
     << "  shellIndex= " << shellIndex;
      G4Exception(
      "G4AtomicTransitionManager::TotalRadiativeTransitionProbability()",
      "de0001",FatalException,ed,"Cannot compute transition probability");
    } 
  return totalRadTransProb;
}

Referenced by TotalNonRadiativeTransitionProbability().

---

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

• G4AtomicTransitionManager.hh
• G4AtomicTransitionManager.cc