| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 // 26 //
>> 27 // $Id: G4AugerTransition.cc,v 1.2 ????
27 // 28 //
28 // Based on G4AtomicTransition.cc by 29 // Based on G4AtomicTransition.cc by
29 // Elena Guardincerri (Elena.Guardincerri@ge.i 30 // Elena Guardincerri (Elena.Guardincerri@ge.infn.it)
30 // 31 //
31 // Author: Alfonso Mantero (Alfonso.Mantero@ge 32 // Author: Alfonso Mantero (Alfonso.Mantero@ge.infn.it)
32 // 33 //
33 // History: 34 // History:
34 // ----------- 35 // -----------
35 // 4 Mar 2002: first implementation 36 // 4 Mar 2002: first implementation
36 // 37 //
37 // ------------------------------------------- 38 // -------------------------------------------------------------------
38 39
39 #include "G4AugerTransition.hh" 40 #include "G4AugerTransition.hh"
40 41
41 //....oooOO0OOooo........oooOO0OOooo........oo <<
42 // the final shell in wich the electron goes i 42 // the final shell in wich the electron goes is needed, to know the data for the auger electron emitted
43 // (i.e. originating shell id, electron energy 43 // (i.e. originating shell id, electron energy and transition probability)
>> 44
44 G4AugerTransition::G4AugerTransition(G4int fin 45 G4AugerTransition::G4AugerTransition(G4int finalShell, std::vector<G4int> transIds,
45 const std::map<G4int,std::vector< 46 const std::map<G4int,std::vector<G4int>,std::less<G4int> >* idMap,
46 const std::map<G4int,G4DataVector 47 const std::map<G4int,G4DataVector,std::less<G4int> >* energyMap,
47 const std::map<G4int,G4DataVector 48 const std::map<G4int,G4DataVector,std::less<G4int> >* probabilityMap)
48 { 49 {
49 finalShellId = finalShell; 50 finalShellId = finalShell;
50 augerOriginatingShellIdsMap = *idMap; 51 augerOriginatingShellIdsMap = *idMap;
51 augerTransitionEnergiesMap = *energyMap; 52 augerTransitionEnergiesMap = *energyMap;
52 augerTransitionProbabilitiesMap = *probabili 53 augerTransitionProbabilitiesMap = *probabilityMap;
53 transitionOriginatingShellIds = std::move(tr << 54 transitionOriginatingShellIds = transIds;
54 } << 55
55 56
56 //....oooOO0OOooo........oooOO0OOooo........oo << 57 }
57 58
58 G4AugerTransition::~G4AugerTransition() 59 G4AugerTransition::~G4AugerTransition()
59 {;} << 60 {
>> 61
>> 62 }
60 63
61 //....oooOO0OOooo........oooOO0OOooo........oo <<
62 // Returns the ids of the shells from wich an 64 // Returns the ids of the shells from wich an auger electron culd came from, given th shell
63 // from wich the transition electron cames fro 65 // from wich the transition electron cames from.
>> 66
64 const std::vector<G4int>* G4AugerTransition::A 67 const std::vector<G4int>* G4AugerTransition::AugerOriginatingShellIds(G4int startShellId) const
65 { 68 {
66 auto shellId = augerOriginatingShellIdsMap.f << 69 std::map<G4int,std::vector<G4int>,std::less<G4int> >::const_iterator shellId = augerOriginatingShellIdsMap.find(startShellId);
67 70
68 const std::vector<G4int>* dataSet = &(*shell 71 const std::vector<G4int>* dataSet = &(*shellId).second;
69 if (dataSet->empty()) << 72 //const std::vector<G4int>* dataOut = 0;
70 G4cout << "Error: no auger Id found"<< G4e << 73
>> 74 if (dataSet->size() == 0) {G4cout << "Error: no auger Id found"<< G4endl;}
>> 75 else {
>> 76
>> 77 // dataOut = &dataSet;
>> 78
>> 79 }
>> 80
71 return dataSet; 81 return dataSet;
72 } 82 }
73 83
74 //....oooOO0OOooo........oooOO0OOooo........oo <<
75 // Returns the ids of the shells from wich an 84 // Returns the ids of the shells from wich an electron cuuld fill the vacancy in finalShellId
>> 85
76 const std::vector<G4int>* G4AugerTransition::T 86 const std::vector<G4int>* G4AugerTransition::TransitionOriginatingShellIds() const
77 { 87 {
>> 88
78 const std::vector<G4int>* dataSet = &transit 89 const std::vector<G4int>* dataSet = &transitionOriginatingShellIds;
79 return dataSet; 90 return dataSet;
80 } 91 }
81 92
82 //....oooOO0OOooo........oooOO0OOooo........oo << 93 // Returns the energiess of the possible auger electrons, given th shell
83 // Returns the energies of the possible auger << 94 // from wich the transition electron cames from.
84 // from which the transition electron came fro << 95
85 const G4DataVector* G4AugerTransition::AugerTr 96 const G4DataVector* G4AugerTransition::AugerTransitionEnergies(G4int startShellId) const
86 { 97 {
87 auto shellId = augerTransitionEnergiesMap.fi << 98 std::map<G4int,G4DataVector,std::less<G4int> >::const_iterator shellId = augerTransitionEnergiesMap.find(startShellId);
88 <<
89 if (shellId == augerTransitionEnergiesMap.en <<
90 { <<
91 G4Exception("G4AugerTransition::AugerTra <<
92 "corresponding map element not found, en <<
93 return 0; <<
94 } <<
95 const G4DataVector* dataSet = &(*shellId).se 99 const G4DataVector* dataSet = &(*shellId).second;
>> 100
>> 101
96 return dataSet; 102 return dataSet;
97 } 103 }
98 104
99 //....oooOO0OOooo........oooOO0OOooo........oo << 105 // Returns the emission probabilities of the auger electrons, given th shell
100 // Returns the emission probabilities of the a <<
101 // from wich the transition electron cames fro 106 // from wich the transition electron cames from.
>> 107
102 const G4DataVector* G4AugerTransition::AugerTr 108 const G4DataVector* G4AugerTransition::AugerTransitionProbabilities(G4int startShellId) const
103 { 109 {
104 auto shellId = augerTransitionProbabilitiesM << 110 std::map<G4int,G4DataVector,std::less<G4int> >::const_iterator shellId = augerTransitionProbabilitiesMap.find(startShellId);
105 if (shellId == augerTransitionProbabilitiesM <<
106 { <<
107 <<
108 G4Exception("G4AugerTransition::AugerTra <<
109 JustWarning,"corresponding map element n <<
110 return 0; <<
111 } <<
112 <<
113 const G4DataVector* dataSet = &(*shellId).se 111 const G4DataVector* dataSet = &(*shellId).second;
114 return dataSet; 112 return dataSet;
115 } 113 }
116 114
117 //....oooOO0OOooo........oooOO0OOooo........oo <<
118 G4int G4AugerTransition::FinalShellId() const 115 G4int G4AugerTransition::FinalShellId() const
119 { 116 {
120 return finalShellId; 117 return finalShellId;
121 } 118 }
122 119
123 //....oooOO0OOooo........oooOO0OOooo........oo <<
124 // Returns the id of the shell from wich come 120 // Returns the id of the shell from wich come the auger electron , given the shell
125 // from wich the transition electron cames fro 121 // from wich the transition electron cames from and the index number.
>> 122
126 G4int G4AugerTransition::AugerOriginatingShell 123 G4int G4AugerTransition::AugerOriginatingShellId(G4int index, G4int startShellId) const
127 { 124 {
128 const std::vector<G4int>* ids = AugerOrigina 125 const std::vector<G4int>* ids = AugerOriginatingShellIds(startShellId);
129 // G4int i = 126 // G4int i =
130 std::vector<G4int>::const_iterator pos = ids 127 std::vector<G4int>::const_iterator pos = ids->begin();
131 G4int n = 0; 128 G4int n = 0;
132 n = *(pos+index); 129 n = *(pos+index);
133 return n; 130 return n;
134 } 131 }
135 132
136 //....oooOO0OOooo........oooOO0OOooo........oo <<
137 // Returns the energy of the auger electron, g 133 // Returns the energy of the auger electron, given the shell
138 // from wich the transition electron cames fro 134 // from wich the transition electron cames from and the index number.
>> 135
139 G4double G4AugerTransition::AugerTransitionEne 136 G4double G4AugerTransition::AugerTransitionEnergy(G4int index, G4int startShellId) const
140 { 137 {
141 const G4DataVector* energies = AugerTransiti 138 const G4DataVector* energies = AugerTransitionEnergies(startShellId);
142 G4double energy = 0; << 139 G4double energy = 0;
143 if (index < (G4int) energies->size()) { 140 if (index < (G4int) energies->size()) {
144 G4DataVector::const_iterator pos = energie 141 G4DataVector::const_iterator pos = energies->begin();
145 energy = *(pos+index); 142 energy = *(pos+index);
146 } 143 }
147 return energy; 144 return energy;
148 } 145 }
149 146
150 //....oooOO0OOooo........oooOO0OOooo........oo <<
151 // Returns the probability of the auger emissi 147 // Returns the probability of the auger emission, given the shell
152 // from wich the transition electron cames fro 148 // from wich the transition electron cames from and the index number.
>> 149
153 G4double G4AugerTransition::AugerTransitionPro 150 G4double G4AugerTransition::AugerTransitionProbability(G4int index, G4int startShellId) const
154 { 151 {
155 152
156 const G4DataVector *probabilities = AugerTra 153 const G4DataVector *probabilities = AugerTransitionProbabilities(startShellId);
157 G4DataVector::const_iterator pos = probabili 154 G4DataVector::const_iterator pos = probabilities->begin();
158 155
159 G4double probability = 0; 156 G4double probability = 0;
160 probability = *(pos+index); 157 probability = *(pos+index);
161 158
162 return probability; << 159 return probability;
>> 160
163 } 161 }
164 162
165 //....oooOO0OOooo........oooOO0OOooo........oo <<
166 G4int G4AugerTransition::TransitionOriginating 163 G4int G4AugerTransition::TransitionOriginatingShellId(G4int index) const
167 { 164 {
168 return transitionOriginatingShellIds[index] 165 return transitionOriginatingShellIds[index];
169 } 166 }
170 167
171 168
172 169
173 170
174 171
175 172
176 173
177 174
178 175
179 176
180 177
181 178
182 179
183 180
184 181
185 182
186 183
187 184
188 185
189 186
190 187
191 188