66 std::map<G4int,std::vector<G4AtomicShell*>,std::less<G4int> >::iterator
pos;
70 std::vector<G4AtomicShell*>vec = (*pos).second;
71 G4int vecSize = vec.size();
73 for (
G4int i=0; i< vecSize; ++i){
79 std::map<G4int,std::vector<G4FluoTransition*>,std::less<G4int> >::iterator ppos;
82 std::vector<G4FluoTransition*>vec = (*ppos).second;
84 G4int vecSize=vec.size();
86 for (
G4int i=0; i< vecSize; i++){
96 std::map<G4int,std::vector<G4AtomicShell*>,std::less<G4int> >::const_iterator
pos;
102 std::vector<G4AtomicShell*> v = (*pos).second;
103 if (shellIndex < v.size()) {
return v[shellIndex]; }
107 size_t lastShell = v.size();
109 ed <<
"No de-excitation for Z= " << Z
110 <<
" shellIndex= " << shellIndex
111 <<
">= numberOfShells= " << lastShell;
112 G4Exception(
"G4AtomicTransitionManager::Shell()",
"de0001",
114 if (lastShell > 0) {
return v[lastShell - 1]; }
120 ed <<
"No de-excitation for Z= " << Z
121 <<
" shellIndex= " << shellIndex
122 <<
". AtomicShell not found - check if data are uploaded";
123 G4Exception(
"G4AtomicTransitionManager::Shell()",
"de0001",
136 std::map<G4int,std::vector<G4FluoTransition*>,std::less<G4int> >::const_iterator
pos;
140 std::vector<G4FluoTransition*> v = (*pos).second;
141 if (shellIndex < v.size()) {
return(v[shellIndex]); }
145 ed <<
"No fluo transition for Z= " << Z
146 <<
" shellIndex= " << shellIndex;
147 G4Exception(
"G4AtomicTransitionManager::ReachebleShell()",
"de0002",
154 ed <<
"No transition table for Z= " << Z
155 <<
" shellIndex= " << shellIndex;
156 G4Exception(
"G4AtomicTransitionManager::ReachableShell()",
"de0001",
164 G4int vacancyShellIndex)
const
171 std::map<G4int,std::vector<G4AtomicShell*>,std::less<G4int> >::const_iterator
pos;
177 res = ((*pos).second).size();
181 ed <<
"No deexcitation for Z= " << Z;
182 G4Exception(
"G4AtomicTransitionManager::NumberOfShells()",
"de0001",
193 std::map<G4int,std::vector<G4FluoTransition*>,std::less<G4int> >::const_iterator
pos;
198 res = ((*pos).second).size();
203 ed <<
"No deexcitation for Z= " << Z
204 <<
", so energy deposited locally";
205 G4Exception(
"G4AtomicTransitionManager::NumberOfReachebleShells()",
221 G4int Z,
size_t shellIndex)
const
223 std::map<G4int,std::vector<G4FluoTransition*>,std::less<G4int> >::const_iterator
pos;
230 std::vector<G4FluoTransition*> v = (*pos).second;
232 if (shellIndex < v.size())
237 for (
size_t j=0; j<transProb.size(); ++j)
239 totalRadTransProb += transProb[j];
245 ed <<
"Zero transition probability for Z=" << Z
246 <<
" shellIndex= " << shellIndex;
248 "G4AtomicTransitionManager::TotalRadiativeTransitionProbability()",
255 ed <<
"No deexcitation for Z=" << Z
256 <<
" shellIndex= " << shellIndex;
258 "G4AtomicTransitionManager::TotalRadiativeTransitionProbability()",
259 "de0001",
FatalException,ed,
"Cannot compute transition probability");
261 return totalRadTransProb;
265 G4int Z,
size_t shellIndex)
const
268 if(prob > 1.0 || prob < 0.0) {
270 ed <<
"Total probability mismatch Z= " << Z
271 <<
" shellIndex= " << shellIndex
272 <<
" prob= " << prob;
274 "G4AtomicTransitionManager::TotalNonRadiativeTransitionProbability()",
275 "de0003",
FatalException,ed,
"Cannot compute non-radiative probability");
295 "/fluor_Bearden":
"/fluor");
299 shellManager->
LoadData(fluoDirectory+
"/binding");
308 std::vector<G4AtomicShell*> vectorOfShells;
309 size_t shellIndex = 0;
313 for (shellIndex = 0; shellIndex<numberOfShells; ++shellIndex)
320 vectorOfShells.push_back(shell);
332 std::vector<G4FluoTransition*> vectorOfTransitions;
335 size_t numberOfVacancies = fluoManager-> NumberOfVacancies();
337 for(
size_t vacancyIndex = 0; vacancyIndex<numberOfVacancies;
340 std::vector<G4int> vectorOfIds;
345 size_t numberOfTransitions =
347 for (
size_t origShellIndex = 0; origShellIndex < numberOfTransitions;
350 G4int originatingShellId =
352 vectorOfIds.push_back(originatingShellId);
356 vectorOfEnergies.push_back(transitionEnergy);
359 vectorOfProbabilities.push_back(transitionProbability);
363 vectorOfEnergies,vectorOfProbabilities);
364 vectorOfTransitions.push_back(transition);
G4int NumberOfShells(G4int Z) const
G4double StartShellProb(G4int initIndex, G4int vacancyIndex) const
std::ostringstream G4ExceptionDescription
const G4AugerTransition * ReachableAugerShell(G4int Z, G4int shellIndex) const
G4int NumberOfReachableShells(G4int Z) const
std::map< G4int, std::vector< G4FluoTransition * >, std::less< G4int > > transitionTable
std::map< G4int, std::vector< G4AtomicShell * >, std::less< G4int > > shellTable
G4AtomicTransitionManager()
const G4DataVector & TransitionProbabilities() const
#define G4MUTEX_INITIALIZER
G4int ShellId(G4int Z, G4int shellIndex) const
void LoadData(const G4String &fileName)
G4double TotalNonRadiativeTransitionProbability(G4int Z, size_t shellIndex) const
G4double BindingEnergy(G4int Z, G4int shellIndex) const
G4int VacancyId(G4int vacancyIndex) const
G4int NumberOfReachableAugerShells(G4int Z) const
G4int StartShellId(G4int initIndex, G4int vacancyIndex) const
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
G4double TotalRadiativeTransitionProbability(G4int Z, size_t shellIndex) const
G4double StartShellEnergy(G4int initIndex, G4int vacancyIndex) const
const G4FluoTransition * ReachableShell(G4int Z, size_t shellIndex) const
size_t NumberOfShells(G4int Z) const
static G4EmParameters * Instance()
size_t NumberOfTransitions(G4int vacancyIndex) const
~G4AtomicTransitionManager()
static G4AtomicTransitionManager * instance
G4AugerTransition * GetAugerTransition(G4int Z, G4int vacancyShellIndex)
size_t NumberOfVacancies(G4int Z) const
G4double bindingEnergy(G4int A, G4int Z)
static G4AtomicTransitionManager * Instance()
G4AtomicShell * Shell(G4int Z, size_t shellIndex) const
G4bool isInitialized()
Check if the generator is initialized.
static const G4double pos
G4bool BeardenFluoDir() const