| 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: G4RDAugerTransition.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 "G4RDAugerTransition.hh" 40 #include "G4RDAugerTransition.hh"
40 41
41 // 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
42 // (i.e. originating shell id, electron energy 43 // (i.e. originating shell id, electron energy and transition probability)
43 44
44 G4RDAugerTransition::G4RDAugerTransition(G4int 45 G4RDAugerTransition::G4RDAugerTransition(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 = transIds; 54 transitionOriginatingShellIds = transIds;
54 55
55 56
56 } 57 }
57 58
58 G4RDAugerTransition::~G4RDAugerTransition() 59 G4RDAugerTransition::~G4RDAugerTransition()
59 { 60 {
60 61
61 } 62 }
62 63
63 // 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
64 // from wich the transition electron cames fro 65 // from wich the transition electron cames from.
65 66
66 const std::vector<G4int>* G4RDAugerTransition: 67 const std::vector<G4int>* G4RDAugerTransition::AugerOriginatingShellIds(G4int startShellId) const
67 { 68 {
68 std::map<G4int,std::vector<G4int>,std::less< 69 std::map<G4int,std::vector<G4int>,std::less<G4int> >::const_iterator shellId = augerOriginatingShellIdsMap.find(startShellId);
69 70
70 const std::vector<G4int>* dataSet = &(*shell 71 const std::vector<G4int>* dataSet = &(*shellId).second;
71 //const std::vector<G4int>* dataOut = 0; 72 //const std::vector<G4int>* dataOut = 0;
72 73
73 if (dataSet->size() == 0) {G4cout << "Error: 74 if (dataSet->size() == 0) {G4cout << "Error: no auger Id found"<< G4endl;}
74 else { 75 else {
75 76
76 // dataOut = &dataSet; 77 // dataOut = &dataSet;
77 78
78 } 79 }
79 80
80 return dataSet; 81 return dataSet;
81 } 82 }
82 83
83 // 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
84 85
85 const std::vector<G4int>* G4RDAugerTransition: 86 const std::vector<G4int>* G4RDAugerTransition::TransitionOriginatingShellIds() const
86 { 87 {
87 88
88 const std::vector<G4int>* dataSet = &transit 89 const std::vector<G4int>* dataSet = &transitionOriginatingShellIds;
89 return dataSet; 90 return dataSet;
90 } 91 }
91 92
92 // Returns the energiess of the possible auger 93 // Returns the energiess of the possible auger electrons, given th shell
93 // from wich the transition electron cames fro 94 // from wich the transition electron cames from.
94 95
95 const G4DataVector* G4RDAugerTransition::Auger 96 const G4DataVector* G4RDAugerTransition::AugerTransitionEnergies(G4int startShellId) const
96 { 97 {
97 std::map<G4int,G4DataVector,std::less<G4int> 98 std::map<G4int,G4DataVector,std::less<G4int> >::const_iterator shellId = augerTransitionEnergiesMap.find(startShellId);
98 const G4DataVector* dataSet = &(*shellId).se 99 const G4DataVector* dataSet = &(*shellId).second;
99 100
100 101
101 return dataSet; 102 return dataSet;
102 } 103 }
103 104
104 // Returns the emission probabilities of the a 105 // Returns the emission probabilities of the auger electrons, given th shell
105 // from wich the transition electron cames fro 106 // from wich the transition electron cames from.
106 107
107 const G4DataVector* G4RDAugerTransition::Auger 108 const G4DataVector* G4RDAugerTransition::AugerTransitionProbabilities(G4int startShellId) const
108 { 109 {
109 std::map<G4int,G4DataVector,std::less<G4int> 110 std::map<G4int,G4DataVector,std::less<G4int> >::const_iterator shellId = augerTransitionProbabilitiesMap.find(startShellId);
110 const G4DataVector* dataSet = &(*shellId).se 111 const G4DataVector* dataSet = &(*shellId).second;
111 return dataSet; 112 return dataSet;
112 } 113 }
113 114
114 G4int G4RDAugerTransition::FinalShellId() cons 115 G4int G4RDAugerTransition::FinalShellId() const
115 { 116 {
116 return finalShellId; 117 return finalShellId;
117 } 118 }
118 119
119 // 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
120 // from wich the transition electron cames fro 121 // from wich the transition electron cames from and the index number.
121 122
122 G4int G4RDAugerTransition::AugerOriginatingShe 123 G4int G4RDAugerTransition::AugerOriginatingShellId(G4int index, G4int startShellId) const
123 { 124 {
124 const std::vector<G4int>* ids = AugerOrigina 125 const std::vector<G4int>* ids = AugerOriginatingShellIds(startShellId);
125 // G4int i = 126 // G4int i =
126 std::vector<G4int>::const_iterator pos = ids 127 std::vector<G4int>::const_iterator pos = ids->begin();
127 G4int n = 0; 128 G4int n = 0;
128 n = *(pos+index); 129 n = *(pos+index);
129 return n; 130 return n;
130 } 131 }
131 132
132 // Returns the energy of the auger electron, g 133 // Returns the energy of the auger electron, given the shell
133 // from wich the transition electron cames fro 134 // from wich the transition electron cames from and the index number.
134 135
135 G4double G4RDAugerTransition::AugerTransitionE 136 G4double G4RDAugerTransition::AugerTransitionEnergy(G4int index, G4int startShellId) const
136 { 137 {
137 const G4DataVector* energies = AugerTransiti 138 const G4DataVector* energies = AugerTransitionEnergies(startShellId);
138 G4double energy = 0; 139 G4double energy = 0;
139 if (index < (G4int) energies->size()) { 140 if (index < (G4int) energies->size()) {
140 G4DataVector::const_iterator pos = energie 141 G4DataVector::const_iterator pos = energies->begin();
141 energy = *(pos+index); 142 energy = *(pos+index);
142 } 143 }
143 return energy; 144 return energy;
144 } 145 }
145 146
146 // Returns the probability of the auger emissi 147 // Returns the probability of the auger emission, given the shell
147 // from wich the transition electron cames fro 148 // from wich the transition electron cames from and the index number.
148 149
149 G4double G4RDAugerTransition::AugerTransitionP 150 G4double G4RDAugerTransition::AugerTransitionProbability(G4int index, G4int startShellId) const
150 { 151 {
151 152
152 const G4DataVector *probabilities = AugerTra 153 const G4DataVector *probabilities = AugerTransitionProbabilities(startShellId);
153 G4DataVector::const_iterator pos = probabili 154 G4DataVector::const_iterator pos = probabilities->begin();
154 155
155 G4double probability = 0; 156 G4double probability = 0;
156 probability = *(pos+index); 157 probability = *(pos+index);
157 158
158 return probability; 159 return probability;
159 160
160 } 161 }
161 162
162 G4int G4RDAugerTransition::TransitionOriginati 163 G4int G4RDAugerTransition::TransitionOriginatingShellId(G4int index) const
163 { 164 {
164 return transitionOriginatingShellIds[index] 165 return transitionOriginatingShellIds[index];
165 } 166 }
166 167
167 168
168 169
169 170
170 171
171 172
172 173
173 174
174 175
175 176
176 177
177 178
178 179
179 180
180 181
181 182
182 183
183 184
184 185
185 186
186 187
187 188