| 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 // 27 //
28 // Based on G4AtomicTransition.cc by 28 // Based on G4AtomicTransition.cc by
29 // Elena Guardincerri (Elena.Guardincerri@ge.i 29 // Elena Guardincerri (Elena.Guardincerri@ge.infn.it)
30 // 30 //
31 // Author: Alfonso Mantero (Alfonso.Mantero@ge 31 // Author: Alfonso Mantero (Alfonso.Mantero@ge.infn.it)
32 // 32 //
33 // History: 33 // History:
34 // ----------- 34 // -----------
35 // 4 Mar 2002: first implementation 35 // 4 Mar 2002: first implementation
36 // 36 //
37 // ------------------------------------------- 37 // -------------------------------------------------------------------
38 38
39 #include "G4AugerTransition.hh" 39 #include "G4AugerTransition.hh"
40 40
41 //....oooOO0OOooo........oooOO0OOooo........oo <<
42 // the final shell in wich the electron goes i 41 // 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 42 // (i.e. originating shell id, electron energy and transition probability)
>> 43
44 G4AugerTransition::G4AugerTransition(G4int fin 44 G4AugerTransition::G4AugerTransition(G4int finalShell, std::vector<G4int> transIds,
45 const std::map<G4int,std::vector< 45 const std::map<G4int,std::vector<G4int>,std::less<G4int> >* idMap,
46 const std::map<G4int,G4DataVector 46 const std::map<G4int,G4DataVector,std::less<G4int> >* energyMap,
47 const std::map<G4int,G4DataVector 47 const std::map<G4int,G4DataVector,std::less<G4int> >* probabilityMap)
48 { 48 {
49 finalShellId = finalShell; 49 finalShellId = finalShell;
50 augerOriginatingShellIdsMap = *idMap; 50 augerOriginatingShellIdsMap = *idMap;
51 augerTransitionEnergiesMap = *energyMap; 51 augerTransitionEnergiesMap = *energyMap;
52 augerTransitionProbabilitiesMap = *probabili 52 augerTransitionProbabilitiesMap = *probabilityMap;
53 transitionOriginatingShellIds = std::move(tr << 53 transitionOriginatingShellIds = transIds;
54 } << 54
55 55
56 //....oooOO0OOooo........oooOO0OOooo........oo << 56 }
57 57
58 G4AugerTransition::~G4AugerTransition() 58 G4AugerTransition::~G4AugerTransition()
59 {;} << 59 {
>> 60
>> 61 }
60 62
61 //....oooOO0OOooo........oooOO0OOooo........oo <<
62 // Returns the ids of the shells from wich an 63 // 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 64 // from wich the transition electron cames from.
>> 65
64 const std::vector<G4int>* G4AugerTransition::A 66 const std::vector<G4int>* G4AugerTransition::AugerOriginatingShellIds(G4int startShellId) const
65 { 67 {
66 auto shellId = augerOriginatingShellIdsMap.f << 68 std::map<G4int,std::vector<G4int>,std::less<G4int> >::const_iterator shellId = augerOriginatingShellIdsMap.find(startShellId);
67 69
68 const std::vector<G4int>* dataSet = &(*shell 70 const std::vector<G4int>* dataSet = &(*shellId).second;
69 if (dataSet->empty()) << 71 //const std::vector<G4int>* dataOut = 0;
70 G4cout << "Error: no auger Id found"<< G4e << 72
>> 73 if (dataSet->size() == 0) {G4cout << "Error: no auger Id found"<< G4endl;}
>> 74 else {
>> 75
>> 76 // dataOut = &dataSet;
>> 77
>> 78 }
>> 79
71 return dataSet; 80 return dataSet;
72 } 81 }
73 82
74 //....oooOO0OOooo........oooOO0OOooo........oo <<
75 // Returns the ids of the shells from wich an 83 // Returns the ids of the shells from wich an electron cuuld fill the vacancy in finalShellId
>> 84
76 const std::vector<G4int>* G4AugerTransition::T 85 const std::vector<G4int>* G4AugerTransition::TransitionOriginatingShellIds() const
77 { 86 {
>> 87
78 const std::vector<G4int>* dataSet = &transit 88 const std::vector<G4int>* dataSet = &transitionOriginatingShellIds;
79 return dataSet; 89 return dataSet;
80 } 90 }
81 91
82 //....oooOO0OOooo........oooOO0OOooo........oo << 92 // Returns the energiess of the possible auger electrons, given th shell
83 // Returns the energies of the possible auger << 93 // from wich the transition electron cames from.
84 // from which the transition electron came fro << 94
85 const G4DataVector* G4AugerTransition::AugerTr 95 const G4DataVector* G4AugerTransition::AugerTransitionEnergies(G4int startShellId) const
86 { 96 {
87 auto shellId = augerTransitionEnergiesMap.fi << 97 std::map<G4int,G4DataVector,std::less<G4int> >::const_iterator shellId = augerTransitionEnergiesMap.find(startShellId);
88 98
89 if (shellId == augerTransitionEnergiesMap.en 99 if (shellId == augerTransitionEnergiesMap.end() )
90 { 100 {
91 G4Exception("G4AugerTransition::AugerTra << 101 G4Exception("G4AugerTransition::AugerTransitionEnergies()","de0002",JustWarning,"corresponding map element not found, energy deposited locally");
92 "corresponding map element not found, en <<
93 return 0; 102 return 0;
94 } 103 }
>> 104
95 const G4DataVector* dataSet = &(*shellId).se 105 const G4DataVector* dataSet = &(*shellId).second;
>> 106
>> 107
96 return dataSet; 108 return dataSet;
97 } 109 }
98 110
99 //....oooOO0OOooo........oooOO0OOooo........oo <<
100 // Returns the emission probabilities of the a 111 // Returns the emission probabilities of the auger electrons, given the shell
101 // from wich the transition electron cames fro 112 // from wich the transition electron cames from.
>> 113
102 const G4DataVector* G4AugerTransition::AugerTr 114 const G4DataVector* G4AugerTransition::AugerTransitionProbabilities(G4int startShellId) const
103 { 115 {
104 auto shellId = augerTransitionProbabilitiesM << 116
>> 117 //debugging
>> 118 //if (startShellId == 1){G4cout <<"OI!!!"<< G4endl;}
>> 119
>> 120 std::map<G4int,G4DataVector,std::less<G4int> >::const_iterator shellId = augerTransitionProbabilitiesMap.find(startShellId);
>> 121
105 if (shellId == augerTransitionProbabilitiesM 122 if (shellId == augerTransitionProbabilitiesMap.end() )
106 { 123 {
107 124
108 G4Exception("G4AugerTransition::AugerTra << 125 G4Exception("G4AugerTransition::AugerTransitionProbabilities()","de0002",JustWarning,"corresponding map element not found, energy deposited locally");
109 JustWarning,"corresponding map element n <<
110 return 0; 126 return 0;
111 } 127 }
112 128
113 const G4DataVector* dataSet = &(*shellId).se 129 const G4DataVector* dataSet = &(*shellId).second;
>> 130 // debugging purpose:
>> 131 /* G4cout << "id: " << shellId->first << G4endl;
>> 132 G4cout << "size:" << dataSet->size() << G4endl;
>> 133 for (G4int i = 0; i < dataSet->size(); i++){
>> 134 G4cout << (dataSet[0])[i] << G4endl;
>> 135 }*/
114 return dataSet; 136 return dataSet;
115 } 137 }
116 138
117 //....oooOO0OOooo........oooOO0OOooo........oo <<
118 G4int G4AugerTransition::FinalShellId() const 139 G4int G4AugerTransition::FinalShellId() const
119 { 140 {
120 return finalShellId; 141 return finalShellId;
121 } 142 }
122 143
123 //....oooOO0OOooo........oooOO0OOooo........oo <<
124 // Returns the id of the shell from wich come 144 // Returns the id of the shell from wich come the auger electron , given the shell
125 // from wich the transition electron cames fro 145 // from wich the transition electron cames from and the index number.
>> 146
126 G4int G4AugerTransition::AugerOriginatingShell 147 G4int G4AugerTransition::AugerOriginatingShellId(G4int index, G4int startShellId) const
127 { 148 {
128 const std::vector<G4int>* ids = AugerOrigina 149 const std::vector<G4int>* ids = AugerOriginatingShellIds(startShellId);
129 // G4int i = 150 // G4int i =
130 std::vector<G4int>::const_iterator pos = ids 151 std::vector<G4int>::const_iterator pos = ids->begin();
131 G4int n = 0; 152 G4int n = 0;
132 n = *(pos+index); 153 n = *(pos+index);
133 return n; 154 return n;
134 } 155 }
135 156
136 //....oooOO0OOooo........oooOO0OOooo........oo <<
137 // Returns the energy of the auger electron, g 157 // Returns the energy of the auger electron, given the shell
138 // from wich the transition electron cames fro 158 // from wich the transition electron cames from and the index number.
>> 159
139 G4double G4AugerTransition::AugerTransitionEne 160 G4double G4AugerTransition::AugerTransitionEnergy(G4int index, G4int startShellId) const
140 { 161 {
141 const G4DataVector* energies = AugerTransiti 162 const G4DataVector* energies = AugerTransitionEnergies(startShellId);
142 G4double energy = 0; << 163 G4double energy = 0;
143 if (index < (G4int) energies->size()) { 164 if (index < (G4int) energies->size()) {
144 G4DataVector::const_iterator pos = energie 165 G4DataVector::const_iterator pos = energies->begin();
145 energy = *(pos+index); 166 energy = *(pos+index);
146 } 167 }
147 return energy; 168 return energy;
148 } 169 }
149 170
150 //....oooOO0OOooo........oooOO0OOooo........oo <<
151 // Returns the probability of the auger emissi 171 // Returns the probability of the auger emission, given the shell
152 // from wich the transition electron cames fro 172 // from wich the transition electron cames from and the index number.
>> 173
153 G4double G4AugerTransition::AugerTransitionPro 174 G4double G4AugerTransition::AugerTransitionProbability(G4int index, G4int startShellId) const
154 { 175 {
155 176
156 const G4DataVector *probabilities = AugerTra 177 const G4DataVector *probabilities = AugerTransitionProbabilities(startShellId);
157 G4DataVector::const_iterator pos = probabili 178 G4DataVector::const_iterator pos = probabilities->begin();
158 179
159 G4double probability = 0; 180 G4double probability = 0;
160 probability = *(pos+index); 181 probability = *(pos+index);
161 182
162 return probability; << 183 return probability;
>> 184
163 } 185 }
164 186
165 //....oooOO0OOooo........oooOO0OOooo........oo <<
166 G4int G4AugerTransition::TransitionOriginating 187 G4int G4AugerTransition::TransitionOriginatingShellId(G4int index) const
167 { 188 {
168 return transitionOriginatingShellIds[index] 189 return transitionOriginatingShellIds[index];
169 } 190 }
170 191
171 192
172 193
173 194
174 195
175 196
176 197
177 198
178 199
179 200
180 201
181 202
182 203
183 204
184 205
185 206
186 207
187 208
188 209
189 210
190 211
191 212