Geant4-11
G4AtomicTransitionManager.cc
Go to the documentation of this file.
1//
2// ********************************************************************
3// * License and Disclaimer *
4// * *
5// * The Geant4 software is copyright of the Copyright Holders of *
6// * the Geant4 Collaboration. It is provided under the terms and *
7// * conditions of the Geant4 Software License, included in the file *
8// * LICENSE and available at http://cern.ch/geant4/license . These *
9// * include a list of copyright holders. *
10// * *
11// * Neither the authors of this software system, nor their employing *
12// * institutes,nor the agencies providing financial support for this *
13// * work make any representation or warranty, express or implied, *
14// * regarding this software system or assume any liability for its *
15// * use. Please see the license in the file LICENSE and URL above *
16// * for the full disclaimer and the limitation of liability. *
17// * *
18// * This code implementation is the result of the scientific and *
19// * technical work of the GEANT4 collaboration. *
20// * By using, copying, modifying or distributing the software (or *
21// * any work based on the software) you agree to acknowledge its *
22// * use in resulting scientific publications, and indicate your *
23// * acceptance of all terms of the Geant4 Software license. *
24// ********************************************************************
25//
26//
27//
28// Authors: Elena Guardincerri (Elena.Guardincerri@ge.infn.it)
29// Alfonso Mantero (Alfonso.Mantero@ge.infn.it)
30//
31// History:
32// -----------
33// 16 Sep 2001 E. Guardincerri First Committed to cvs
34//
35// -------------------------------------------------------------------
36
38#include "G4EmParameters.hh"
39#include "G4FluoData.hh"
40#include "G4AugerData.hh"
41#include "G4AutoLock.hh"
43
44//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
46
47//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
49{
50 if (instance == nullptr) {
52 }
53 return instance;
54}
55
56//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
58 : augerData(nullptr),
59 verboseLevel(0),
60 isInitialized(false)
61{}
62
63//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
65{
66 delete augerData;
67
68 for (auto& pos : shellTable){
69 std::vector<G4AtomicShell*>vec = pos.second;
70 G4int vecSize = vec.size();
71 for (G4int i=0; i< vecSize; ++i){
72 G4AtomicShell* shell = vec[i];
73 delete shell;
74 }
75 }
76
77 for (auto& ppos : transitionTable)
78 {
79 std::vector<G4FluoTransition*>vec = ppos.second;
80 G4int vecSize=vec.size();
81
82 for (G4int i=0; i< vecSize; ++i){
83 G4FluoTransition* transition = vec[i];
84 delete transition;
85 }
86 }
87}
88//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
90G4AtomicTransitionManager::Shell(G4int Z, size_t shellIndex) const
91{
92 auto pos = shellTable.find(Z);
93
94 if (pos!= shellTable.end())
95 {
96 std::vector<G4AtomicShell*> v = (*pos).second;
97 if (shellIndex < v.size()) { return v[shellIndex]; }
98
99 else
100 {
101 size_t lastShell = v.size();
103 ed << "No de-excitation for Z= " << Z
104 << " shellIndex= " << shellIndex
105 << ">= numberOfShells= " << lastShell;
106 if (verboseLevel > 0)
107 G4Exception("G4AtomicTransitionManager::Shell()","de0001",
108 JustWarning,ed," AtomicShell not found");
109 if (lastShell > 0) { return v[lastShell - 1]; }
110 }
111 }
112 else
113 {
115 ed << "No de-excitation for Z= " << Z
116 << " shellIndex= " << shellIndex
117 << ". AtomicShell not found - check if data are uploaded";
118 G4Exception("G4AtomicTransitionManager::Shell()","de0001",
119 FatalException,ed,"");
120 }
121 return 0;
122}
123
124//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
125
126// This function gives, upon Z and the Index of the initial shell where
127// the vacancy is, the radiative transition that can happen (originating
128// shell, energy, probability)
129const G4FluoTransition*
131{
132 auto pos = transitionTable.find(Z);
133 if (pos!= transitionTable.end())
134 {
135 std::vector<G4FluoTransition*> v = (*pos).second;
136 if (shellIndex < v.size()) { return(v[shellIndex]); }
137
138 else {
140 ed << "No fluo transition for Z= " << Z
141 << " shellIndex= " << shellIndex;
142 G4Exception("G4AtomicTransitionManager::ReachebleShell()","de0002",
143 FatalException,ed,"");
144 }
145 }
146 else
147 {
149 ed << "No transition table for Z= " << Z
150 << " shellIndex= " << shellIndex;
151 G4Exception("G4AtomicTransitionManager::ReachableShell()","de0001",
152 FatalException,ed,"");
153 }
154 return 0;
155}
156
157//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
158const G4AugerTransition*
160 G4int vacancyShellIndex) const
161{
162 return augerData->GetAugerTransition(Z,vacancyShellIndex);
163}
164
165//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
167{
168 auto pos = shellTable.find(Z);
169
170 G4int res = 0;
171 if (pos != shellTable.end()){
172
173 res = ((*pos).second).size();
174
175 } else {
177 ed << "No deexcitation for Z= " << Z;
178 G4Exception("G4AtomicTransitionManager::NumberOfShells()","de0001",
179 FatalException, ed, "");
180 }
181 return res;
182}
183
184//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
185// This function returns the number of possible radiative transitions for
186// the atom with atomic number Z i.e. the number of shell in wich a vacancy
187// can be filled with a radiative transition
189{
190 auto pos = transitionTable.find(Z);
191 G4int res = 0;
192 if (pos!= transitionTable.end())
193 {
194 res = ((*pos).second).size();
195 }
196 else
197 {
199 ed << "No deexcitation for Z= " << Z
200 << ", so energy deposited locally";
201 G4Exception("G4AtomicTransitionManager::NumberOfReachebleShells()",
202 "de0001",FatalException,ed,"");
203 }
204 return res;
205}
206
207//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
208// This function returns the number of possible NON-radiative transitions
209// for the atom with atomic number Z i.e. the number of shell in wich a
210// vacancy can be filled with a NON-radiative transition
212{
214}
215
216//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
218 G4int Z, size_t shellIndex) const
219{
220 auto pos = transitionTable.find(Z);
221 G4double totalRadTransProb = 0.0;
222
223 if (pos!= transitionTable.end())
224 {
225 std::vector<G4FluoTransition*> v = (*pos).second;
226
227 if (shellIndex < v.size())
228 {
229 G4FluoTransition* transition = v[shellIndex];
230 G4DataVector transProb = transition->TransitionProbabilities();
231
232 for (size_t j=0; j<transProb.size(); ++j)
233 {
234 totalRadTransProb += transProb[j];
235 }
236 }
237 else
238 {
240 ed << "Zero transition probability for Z=" << Z
241 << " shellIndex= " << shellIndex;
243 "G4AtomicTransitionManager::TotalRadiativeTransitionProbability()",
244 "de0002",FatalException,"Incorrect de-excitation");
245 }
246 }
247 else
248 {
250 ed << "No deexcitation for Z=" << Z
251 << " shellIndex= " << shellIndex;
253 "G4AtomicTransitionManager::TotalRadiativeTransitionProbability()",
254 "de0001",FatalException,ed,"Cannot compute transition probability");
255 }
256 return totalRadTransProb;
257}
258
259//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
261 G4int Z, size_t shellIndex) const
262{
263 G4double prob = 1.0 - TotalRadiativeTransitionProbability(Z, shellIndex);
264 if(prob > 1.0 || prob < 0.0) {
266 ed << "Total probability mismatch Z= " << Z
267 << " shellIndex= " << shellIndex
268 << " prob= " << prob;
270 "G4AtomicTransitionManager::TotalNonRadiativeTransitionProbability()",
271 "de0003",FatalException,ed,"Cannot compute non-radiative probability");
272 return 0.0;
273 }
274 return prob;
275}
276
277//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
279{
281
282 if(isInitialized) { return; }
283 isInitialized = true;
284
285 // Selection of fluorescence files
286
287 G4String fluoDirectory = "/fluor";
288
289 fluoDirectory = (G4EmParameters::Instance()->BeardenFluoDir()?
290 "/fluor_Bearden":"/fluor");
291
292 // infTableLimit is initialized to 6 because EADL lacks data for Z<=5
293 G4ShellData* shellManager = new G4ShellData;
294 shellManager->LoadData(fluoDirectory+"/binding");
295
296 // initialization of the data for auger effect
298
300 {
301 fluoDirectory = "/fluor_ANSTO";
302 G4cout << "The ANSTO fluorescence radiation yields substitute the EADL data libraries" << G4endl;
303 }
304 // Fills shellTable with the data from EADL, identities and binding
305 // energies of shells
306 for (G4int Z = zMin; Z<= zMax; ++Z)
307 {
308
309 if ((Z > 92) && (G4EmParameters::Instance()->ANSTOFluoDir())) fluoDirectory ="/fluor";
310
311 std::vector<G4AtomicShell*> vectorOfShells;
312 size_t shellIndex = 0;
313
314 size_t numberOfShells = shellManager->NumberOfShells(Z);
315 // G4cout << fluoDirectory <<" For Z= " << Z << " " << numberOfShells << " shells" << G4endl;
316 for (shellIndex = 0; shellIndex<numberOfShells; ++shellIndex)
317 {
318 G4int shellId = shellManager->ShellId(Z,shellIndex);
319 G4double bindingEnergy = shellManager->BindingEnergy(Z,shellIndex);
320
321 G4AtomicShell * shell = new G4AtomicShell(shellId,bindingEnergy);
322
323 vectorOfShells.push_back(shell);
324 }
325 shellTable[Z] = vectorOfShells;
326 }
327
328 // Fills transitionTable with the data from EADL, identities, transition
329 // energies and transition probabilities
330 if (G4EmParameters::Instance()->ANSTOFluoDir()) fluoDirectory = "/fluor_ANSTO";
331
332
333 for (G4int Znum= infTableLimit; Znum<=supTableLimit; ++Znum)
334 {
335 if ((Znum > 92) && (G4EmParameters::Instance()->ANSTOFluoDir())) fluoDirectory ="/fluor";
336 G4FluoData* fluoManager = new G4FluoData(fluoDirectory);
337 std::vector<G4FluoTransition*> vectorOfTransitions;
338 fluoManager->LoadData(Znum);
339
340 size_t numberOfVacancies = fluoManager-> NumberOfVacancies();
341
342 for(size_t vacancyIndex = 0; vacancyIndex<numberOfVacancies;
343 ++vacancyIndex)
344 {
345 std::vector<G4int> vectorOfIds;
346 G4DataVector vectorOfEnergies;
347 G4DataVector vectorOfProbabilities;
348
349 G4int finalShell = fluoManager->VacancyId(vacancyIndex);
350 size_t numberOfTransitions =
351 fluoManager->NumberOfTransitions(vacancyIndex);
352 for (size_t origShellIndex = 0; origShellIndex < numberOfTransitions;
353 ++origShellIndex)
354 {
355 G4int originatingShellId =
356 fluoManager->StartShellId(origShellIndex,vacancyIndex);
357 vectorOfIds.push_back(originatingShellId);
358
359 G4double transitionEnergy =
360 fluoManager->StartShellEnergy(origShellIndex,vacancyIndex);
361 vectorOfEnergies.push_back(transitionEnergy);
362 G4double transitionProbability =
363 fluoManager->StartShellProb(origShellIndex,vacancyIndex);
364 vectorOfProbabilities.push_back(transitionProbability);
365 }
366 G4FluoTransition* transition =
367 new G4FluoTransition (finalShell,vectorOfIds,
368 vectorOfEnergies,vectorOfProbabilities);
369 vectorOfTransitions.push_back(transition);
370 }
371 transitionTable[Znum] = vectorOfTransitions;
372 delete fluoManager;
373 }
374 delete shellManager;
375}
376
377
378
379
380
static const G4double pos
@ JustWarning
@ FatalException
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
Definition: G4Exception.cc:35
std::ostringstream G4ExceptionDescription
Definition: G4Exception.hh:40
#define G4MUTEX_INITIALIZER
Definition: G4Threading.hh:85
std::mutex G4Mutex
Definition: G4Threading.hh:81
double G4double
Definition: G4Types.hh:83
int G4int
Definition: G4Types.hh:85
const G4int Z[17]
#define G4endl
Definition: G4ios.hh:57
G4GLOB_DLL std::ostream G4cout
G4int NumberOfReachableShells(G4int Z) const
G4double TotalRadiativeTransitionProbability(G4int Z, size_t shellIndex) const
const G4AugerTransition * ReachableAugerShell(G4int Z, G4int shellIndex) const
G4AtomicShell * Shell(G4int Z, size_t shellIndex) const
std::map< G4int, std::vector< G4AtomicShell * >, std::less< G4int > > shellTable
static G4AtomicTransitionManager * instance
const G4FluoTransition * ReachableShell(G4int Z, size_t shellIndex) const
void Initialise()
needs to be called once from other code before start of run
G4double TotalNonRadiativeTransitionProbability(G4int Z, size_t shellIndex) const
static G4AtomicTransitionManager * Instance()
G4int NumberOfReachableAugerShells(G4int Z) const
std::map< G4int, std::vector< G4FluoTransition * >, std::less< G4int > > transitionTable
G4AugerTransition * GetAugerTransition(G4int Z, G4int vacancyShellIndex)
Definition: G4AugerData.cc:450
size_t NumberOfVacancies(G4int Z) const
Definition: G4AugerData.cc:66
static G4EmParameters * Instance()
G4bool BeardenFluoDir() const
G4bool ANSTOFluoDir() const
G4double StartShellEnergy(G4int initIndex, G4int vacancyIndex) const
Definition: G4FluoData.cc:152
size_t NumberOfTransitions(G4int vacancyIndex) const
Definition: G4FluoData.cc:103
G4int VacancyId(G4int vacancyIndex) const
Given the index of the vacancy returns its identity.
Definition: G4FluoData.cc:81
void LoadData(G4int Z)
Definition: G4FluoData.cc:204
G4int StartShellId(G4int initIndex, G4int vacancyIndex) const
Definition: G4FluoData.cc:124
G4double StartShellProb(G4int initIndex, G4int vacancyIndex) const
Definition: G4FluoData.cc:177
const G4DataVector & TransitionProbabilities() const
Return the probabilities of the transitions.
G4int ShellId(G4int Z, G4int shellIndex) const
Definition: G4ShellData.cc:118
G4double BindingEnergy(G4int Z, G4int shellIndex) const
Definition: G4ShellData.cc:161
size_t NumberOfShells(G4int Z) const
Definition: G4ShellData.cc:81
void LoadData(const G4String &fileName)
Definition: G4ShellData.cc:228
G4bool isInitialized()
G4double bindingEnergy(G4int A, G4int Z)