| 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 // INCL++ intra-nuclear cascade model 26 // INCL++ intra-nuclear cascade model
27 // Alain Boudard, CEA-Saclay, France 27 // Alain Boudard, CEA-Saclay, France
28 // Joseph Cugnon, University of Liege, Belgium 28 // Joseph Cugnon, University of Liege, Belgium
29 // Jean-Christophe David, CEA-Saclay, France 29 // Jean-Christophe David, CEA-Saclay, France
30 // Pekka Kaitaniemi, CEA-Saclay, France, and H 30 // Pekka Kaitaniemi, CEA-Saclay, France, and Helsinki Institute of Physics, Finland
31 // Sylvie Leray, CEA-Saclay, France 31 // Sylvie Leray, CEA-Saclay, France
32 // Davide Mancusi, CEA-Saclay, France 32 // Davide Mancusi, CEA-Saclay, France
33 // 33 //
34 #define INCLXX_IN_GEANT4_MODE 1 34 #define INCLXX_IN_GEANT4_MODE 1
35 35
36 #include "globals.hh" 36 #include "globals.hh"
37 37
38 /** \file G4INCLNaturalIsotopicDistributions.c 38 /** \file G4INCLNaturalIsotopicDistributions.cc
39 * \brief Classes that stores isotopic abundan 39 * \brief Classes that stores isotopic abundances
40 * 40 *
41 * \date 21st October 2012 41 * \date 21st October 2012
42 * \author Davide Mancusi 42 * \author Davide Mancusi
43 */ 43 */
44 44
45 #include "G4INCLNaturalIsotopicDistributions.h 45 #include "G4INCLNaturalIsotopicDistributions.hh"
46 #include "G4INCLRandom.hh" 46 #include "G4INCLRandom.hh"
47 #include "G4INCLLogger.hh" 47 #include "G4INCLLogger.hh"
48 // #include <cassert> 48 // #include <cassert>
49 #include <utility> 49 #include <utility>
50 #include <iostream> 50 #include <iostream>
51 51
52 namespace G4INCL { 52 namespace G4INCL {
53 53
54 Isotope::Isotope(const G4int A, const G4doub 54 Isotope::Isotope(const G4int A, const G4double abundance) :
55 theA(A), 55 theA(A),
56 theAbundance(abundance) 56 theAbundance(abundance)
57 {} 57 {}
58 58
59 IsotopicDistribution::IsotopicDistribution(I 59 IsotopicDistribution::IsotopicDistribution(IsotopeVector const &aVector) :
60 theIsotopes(aVector) 60 theIsotopes(aVector)
61 { 61 {
62 G4double previousAbundance = 0.; 62 G4double previousAbundance = 0.;
63 // Cumulate the abundances 63 // Cumulate the abundances
64 for(IsotopeIter i=theIsotopes.begin(), e=t 64 for(IsotopeIter i=theIsotopes.begin(), e=theIsotopes.end(); i!=e; ++i) {
65 i->theAbundance += previousAbundance; 65 i->theAbundance += previousAbundance;
66 previousAbundance = i->theAbundance; 66 previousAbundance = i->theAbundance;
67 } 67 }
68 // Normalize the abundances to 1 68 // Normalize the abundances to 1
69 const G4double normalisation = 1./theIsoto 69 const G4double normalisation = 1./theIsotopes.back().theAbundance;
70 for(IsotopeIter i=theIsotopes.begin(), e=t 70 for(IsotopeIter i=theIsotopes.begin(), e=theIsotopes.end(); i!=e; ++i)
71 i->theAbundance *= normalisation; 71 i->theAbundance *= normalisation;
72 } 72 }
73 73
74 G4int IsotopicDistribution::drawRandomIsotop 74 G4int IsotopicDistribution::drawRandomIsotope() const {
75 const G4double r = Random::shoot(); 75 const G4double r = Random::shoot();
76 for(unsigned int i=0; i<theIsotopes.size() 76 for(unsigned int i=0; i<theIsotopes.size()-1; ++i) {
77 if(r<=theIsotopes.at(i).theAbundance) 77 if(r<=theIsotopes.at(i).theAbundance)
78 return theIsotopes.at(i).theA; 78 return theIsotopes.at(i).theA;
79 } 79 }
80 return theIsotopes.back().theA; 80 return theIsotopes.back().theA;
81 } 81 }
82 82
83 IsotopeVector const &IsotopicDistribution::g 83 IsotopeVector const &IsotopicDistribution::getIsotopes() const {
84 return theIsotopes; 84 return theIsotopes;
85 } 85 }
86 86
87 IsotopicDistribution const &NaturalIsotopicD 87 IsotopicDistribution const &NaturalIsotopicDistributions::getIsotopicDistribution(G4int Z) const {
88 std::map<G4int, IsotopicDistribution>::con 88 std::map<G4int, IsotopicDistribution>::const_iterator i = theDistributions.find(Z);
89 if(i!=theDistributions.end()) 89 if(i!=theDistributions.end())
90 return i->second; 90 return i->second;
91 else { 91 else {
92 INCL_FATAL("Requested natural isotopic d 92 INCL_FATAL("Requested natural isotopic distribution for synthetic element Z = " << Z << '\n');
93 return theDistributions.begin()->second; 93 return theDistributions.begin()->second;
94 } 94 }
95 } 95 }
96 96
97 G4int NaturalIsotopicDistributions::drawRand 97 G4int NaturalIsotopicDistributions::drawRandomIsotope(G4int Z) const {
98 return getIsotopicDistribution(Z).drawRand 98 return getIsotopicDistribution(Z).drawRandomIsotope();
99 } 99 }
100 100
101 namespace { 101 namespace {
102 std::pair<G4int, Isotope> theRawDistributi 102 std::pair<G4int, Isotope> theRawDistributions[] = {
103 std::pair<G4int, Isotope>(1, Isotope(1, 103 std::pair<G4int, Isotope>(1, Isotope(1, 99.985)),
104 std::pair<G4int, Isotope>(1, Isotope(2, 104 std::pair<G4int, Isotope>(1, Isotope(2, 0.015)),
105 std::pair<G4int, Isotope>(2, Isotope(3, 105 std::pair<G4int, Isotope>(2, Isotope(3, 0.000137)),
106 std::pair<G4int, Isotope>(2, Isotope(4, 106 std::pair<G4int, Isotope>(2, Isotope(4, 99.999863)),
107 std::pair<G4int, Isotope>(3, Isotope(6, 107 std::pair<G4int, Isotope>(3, Isotope(6, 7.5)),
108 std::pair<G4int, Isotope>(3, Isotope(7, 108 std::pair<G4int, Isotope>(3, Isotope(7, 92.5)),
109 std::pair<G4int, Isotope>(4, Isotope(9, 109 std::pair<G4int, Isotope>(4, Isotope(9, 100.0)),
110 std::pair<G4int, Isotope>(5, Isotope(10, 110 std::pair<G4int, Isotope>(5, Isotope(10, 19.9)),
111 std::pair<G4int, Isotope>(5, Isotope(11, 111 std::pair<G4int, Isotope>(5, Isotope(11, 80.1)),
112 std::pair<G4int, Isotope>(6, Isotope(12, 112 std::pair<G4int, Isotope>(6, Isotope(12, 98.90)),
113 std::pair<G4int, Isotope>(6, Isotope(13, 113 std::pair<G4int, Isotope>(6, Isotope(13, 1.10)),
114 std::pair<G4int, Isotope>(7, Isotope(14, 114 std::pair<G4int, Isotope>(7, Isotope(14, 99.634)),
115 std::pair<G4int, Isotope>(7, Isotope(15, 115 std::pair<G4int, Isotope>(7, Isotope(15, 0.366)),
116 std::pair<G4int, Isotope>(8, Isotope(16, 116 std::pair<G4int, Isotope>(8, Isotope(16, 99.762)),
117 std::pair<G4int, Isotope>(8, Isotope(17, 117 std::pair<G4int, Isotope>(8, Isotope(17, 0.038)),
118 std::pair<G4int, Isotope>(8, Isotope(18, 118 std::pair<G4int, Isotope>(8, Isotope(18, 0.200)),
119 std::pair<G4int, Isotope>(9, Isotope(19, 119 std::pair<G4int, Isotope>(9, Isotope(19, 100.0)),
120 std::pair<G4int, Isotope>(10, Isotope(20 120 std::pair<G4int, Isotope>(10, Isotope(20, 90.48)),
121 std::pair<G4int, Isotope>(10, Isotope(21 121 std::pair<G4int, Isotope>(10, Isotope(21, 0.27)),
122 std::pair<G4int, Isotope>(10, Isotope(22 122 std::pair<G4int, Isotope>(10, Isotope(22, 9.25)),
123 std::pair<G4int, Isotope>(11, Isotope(23 123 std::pair<G4int, Isotope>(11, Isotope(23, 100.0)),
124 std::pair<G4int, Isotope>(12, Isotope(24 124 std::pair<G4int, Isotope>(12, Isotope(24, 78.99)),
125 std::pair<G4int, Isotope>(12, Isotope(25 125 std::pair<G4int, Isotope>(12, Isotope(25, 10.00)),
126 std::pair<G4int, Isotope>(12, Isotope(26 126 std::pair<G4int, Isotope>(12, Isotope(26, 11.01)),
127 std::pair<G4int, Isotope>(13, Isotope(27 127 std::pair<G4int, Isotope>(13, Isotope(27, 100.0)),
128 std::pair<G4int, Isotope>(14, Isotope(28 128 std::pair<G4int, Isotope>(14, Isotope(28, 92.23)),
129 std::pair<G4int, Isotope>(14, Isotope(29 129 std::pair<G4int, Isotope>(14, Isotope(29, 4.67)),
130 std::pair<G4int, Isotope>(14, Isotope(30 130 std::pair<G4int, Isotope>(14, Isotope(30, 3.10)),
131 std::pair<G4int, Isotope>(15, Isotope(31 131 std::pair<G4int, Isotope>(15, Isotope(31, 100.0)),
132 std::pair<G4int, Isotope>(16, Isotope(32 132 std::pair<G4int, Isotope>(16, Isotope(32, 95.02)),
133 std::pair<G4int, Isotope>(16, Isotope(33 133 std::pair<G4int, Isotope>(16, Isotope(33, 0.75)),
134 std::pair<G4int, Isotope>(16, Isotope(34 134 std::pair<G4int, Isotope>(16, Isotope(34, 4.21)),
135 std::pair<G4int, Isotope>(16, Isotope(36 135 std::pair<G4int, Isotope>(16, Isotope(36, 0.02)),
136 std::pair<G4int, Isotope>(17, Isotope(35 136 std::pair<G4int, Isotope>(17, Isotope(35, 75.77)),
137 std::pair<G4int, Isotope>(17, Isotope(37 137 std::pair<G4int, Isotope>(17, Isotope(37, 24.23)),
138 std::pair<G4int, Isotope>(18, Isotope(36 138 std::pair<G4int, Isotope>(18, Isotope(36, 0.337)),
139 std::pair<G4int, Isotope>(18, Isotope(38 139 std::pair<G4int, Isotope>(18, Isotope(38, 0.063)),
140 std::pair<G4int, Isotope>(18, Isotope(40 140 std::pair<G4int, Isotope>(18, Isotope(40, 99.600)),
141 std::pair<G4int, Isotope>(19, Isotope(39 141 std::pair<G4int, Isotope>(19, Isotope(39, 93.2581)),
142 std::pair<G4int, Isotope>(19, Isotope(40 142 std::pair<G4int, Isotope>(19, Isotope(40, 0.0117)),
143 std::pair<G4int, Isotope>(19, Isotope(41 143 std::pair<G4int, Isotope>(19, Isotope(41, 6.7302)),
144 std::pair<G4int, Isotope>(20, Isotope(40 144 std::pair<G4int, Isotope>(20, Isotope(40, 96.941)),
145 std::pair<G4int, Isotope>(20, Isotope(42 145 std::pair<G4int, Isotope>(20, Isotope(42, 0.647)),
146 std::pair<G4int, Isotope>(20, Isotope(43 146 std::pair<G4int, Isotope>(20, Isotope(43, 0.135)),
147 std::pair<G4int, Isotope>(20, Isotope(44 147 std::pair<G4int, Isotope>(20, Isotope(44, 2.086)),
148 std::pair<G4int, Isotope>(20, Isotope(46 148 std::pair<G4int, Isotope>(20, Isotope(46, 0.004)),
149 std::pair<G4int, Isotope>(20, Isotope(48 149 std::pair<G4int, Isotope>(20, Isotope(48, 0.187)),
150 std::pair<G4int, Isotope>(21, Isotope(45 150 std::pair<G4int, Isotope>(21, Isotope(45, 100.0)),
151 std::pair<G4int, Isotope>(22, Isotope(46 151 std::pair<G4int, Isotope>(22, Isotope(46, 8.0)),
152 std::pair<G4int, Isotope>(22, Isotope(47 152 std::pair<G4int, Isotope>(22, Isotope(47, 7.3)),
153 std::pair<G4int, Isotope>(22, Isotope(48 153 std::pair<G4int, Isotope>(22, Isotope(48, 73.8)),
154 std::pair<G4int, Isotope>(22, Isotope(49 154 std::pair<G4int, Isotope>(22, Isotope(49, 5.5)),
155 std::pair<G4int, Isotope>(22, Isotope(50 155 std::pair<G4int, Isotope>(22, Isotope(50, 5.4)),
156 std::pair<G4int, Isotope>(23, Isotope(50 156 std::pair<G4int, Isotope>(23, Isotope(50, 0.250)),
157 std::pair<G4int, Isotope>(23, Isotope(51 157 std::pair<G4int, Isotope>(23, Isotope(51, 99.750)),
158 std::pair<G4int, Isotope>(24, Isotope(50 158 std::pair<G4int, Isotope>(24, Isotope(50, 4.345)),
159 std::pair<G4int, Isotope>(24, Isotope(52 159 std::pair<G4int, Isotope>(24, Isotope(52, 83.789)),
160 std::pair<G4int, Isotope>(24, Isotope(53 160 std::pair<G4int, Isotope>(24, Isotope(53, 9.501)),
161 std::pair<G4int, Isotope>(24, Isotope(54 161 std::pair<G4int, Isotope>(24, Isotope(54, 2.365)),
162 std::pair<G4int, Isotope>(25, Isotope(55 162 std::pair<G4int, Isotope>(25, Isotope(55, 100.0)),
163 std::pair<G4int, Isotope>(26, Isotope(54 163 std::pair<G4int, Isotope>(26, Isotope(54, 5.8)),
164 std::pair<G4int, Isotope>(26, Isotope(56 164 std::pair<G4int, Isotope>(26, Isotope(56, 91.72)),
165 std::pair<G4int, Isotope>(26, Isotope(57 165 std::pair<G4int, Isotope>(26, Isotope(57, 2.2)),
166 std::pair<G4int, Isotope>(26, Isotope(58 166 std::pair<G4int, Isotope>(26, Isotope(58, 0.28)),
167 std::pair<G4int, Isotope>(27, Isotope(59 167 std::pair<G4int, Isotope>(27, Isotope(59, 100.0)),
168 std::pair<G4int, Isotope>(28, Isotope(58 168 std::pair<G4int, Isotope>(28, Isotope(58, 68.077)),
169 std::pair<G4int, Isotope>(28, Isotope(60 169 std::pair<G4int, Isotope>(28, Isotope(60, 26.223)),
170 std::pair<G4int, Isotope>(28, Isotope(61 170 std::pair<G4int, Isotope>(28, Isotope(61, 1.140)),
171 std::pair<G4int, Isotope>(28, Isotope(62 171 std::pair<G4int, Isotope>(28, Isotope(62, 3.634)),
172 std::pair<G4int, Isotope>(28, Isotope(64 172 std::pair<G4int, Isotope>(28, Isotope(64, 0.926)),
173 std::pair<G4int, Isotope>(29, Isotope(63 173 std::pair<G4int, Isotope>(29, Isotope(63, 69.17)),
174 std::pair<G4int, Isotope>(29, Isotope(65 174 std::pair<G4int, Isotope>(29, Isotope(65, 30.83)),
175 std::pair<G4int, Isotope>(30, Isotope(64 175 std::pair<G4int, Isotope>(30, Isotope(64, 48.6)),
176 std::pair<G4int, Isotope>(30, Isotope(66 176 std::pair<G4int, Isotope>(30, Isotope(66, 27.9)),
177 std::pair<G4int, Isotope>(30, Isotope(67 177 std::pair<G4int, Isotope>(30, Isotope(67, 4.1)),
178 std::pair<G4int, Isotope>(30, Isotope(68 178 std::pair<G4int, Isotope>(30, Isotope(68, 18.8)),
179 std::pair<G4int, Isotope>(30, Isotope(70 179 std::pair<G4int, Isotope>(30, Isotope(70, 0.6)),
180 std::pair<G4int, Isotope>(31, Isotope(69 180 std::pair<G4int, Isotope>(31, Isotope(69, 60.108)),
181 std::pair<G4int, Isotope>(31, Isotope(71 181 std::pair<G4int, Isotope>(31, Isotope(71, 39.892)),
182 std::pair<G4int, Isotope>(32, Isotope(70 182 std::pair<G4int, Isotope>(32, Isotope(70, 21.23)),
183 std::pair<G4int, Isotope>(32, Isotope(72 183 std::pair<G4int, Isotope>(32, Isotope(72, 27.66)),
184 std::pair<G4int, Isotope>(32, Isotope(73 184 std::pair<G4int, Isotope>(32, Isotope(73, 7.73)),
185 std::pair<G4int, Isotope>(32, Isotope(74 185 std::pair<G4int, Isotope>(32, Isotope(74, 35.94)),
186 std::pair<G4int, Isotope>(32, Isotope(76 186 std::pair<G4int, Isotope>(32, Isotope(76, 7.44)),
187 std::pair<G4int, Isotope>(33, Isotope(75 187 std::pair<G4int, Isotope>(33, Isotope(75, 100.0)),
188 std::pair<G4int, Isotope>(34, Isotope(74 188 std::pair<G4int, Isotope>(34, Isotope(74, 0.89)),
189 std::pair<G4int, Isotope>(34, Isotope(76 189 std::pair<G4int, Isotope>(34, Isotope(76, 9.36)),
190 std::pair<G4int, Isotope>(34, Isotope(77 190 std::pair<G4int, Isotope>(34, Isotope(77, 7.63)),
191 std::pair<G4int, Isotope>(34, Isotope(78 191 std::pair<G4int, Isotope>(34, Isotope(78, 23.78)),
192 std::pair<G4int, Isotope>(34, Isotope(80 192 std::pair<G4int, Isotope>(34, Isotope(80, 49.61)),
193 std::pair<G4int, Isotope>(34, Isotope(82 193 std::pair<G4int, Isotope>(34, Isotope(82, 8.73)),
194 std::pair<G4int, Isotope>(35, Isotope(79 194 std::pair<G4int, Isotope>(35, Isotope(79, 50.69)),
195 std::pair<G4int, Isotope>(35, Isotope(81 195 std::pair<G4int, Isotope>(35, Isotope(81, 49.31)),
196 std::pair<G4int, Isotope>(36, Isotope(78 196 std::pair<G4int, Isotope>(36, Isotope(78, 0.35)),
197 std::pair<G4int, Isotope>(36, Isotope(80 197 std::pair<G4int, Isotope>(36, Isotope(80, 2.25)),
198 std::pair<G4int, Isotope>(36, Isotope(82 198 std::pair<G4int, Isotope>(36, Isotope(82, 11.6)),
199 std::pair<G4int, Isotope>(36, Isotope(83 199 std::pair<G4int, Isotope>(36, Isotope(83, 11.5)),
200 std::pair<G4int, Isotope>(36, Isotope(84 200 std::pair<G4int, Isotope>(36, Isotope(84, 57.0)),
201 std::pair<G4int, Isotope>(36, Isotope(86 201 std::pair<G4int, Isotope>(36, Isotope(86, 17.3)),
202 std::pair<G4int, Isotope>(37, Isotope(85 202 std::pair<G4int, Isotope>(37, Isotope(85, 72.165)),
203 std::pair<G4int, Isotope>(37, Isotope(87 203 std::pair<G4int, Isotope>(37, Isotope(87, 27.835)),
204 std::pair<G4int, Isotope>(38, Isotope(84 204 std::pair<G4int, Isotope>(38, Isotope(84, 0.56)),
205 std::pair<G4int, Isotope>(38, Isotope(86 205 std::pair<G4int, Isotope>(38, Isotope(86, 9.86)),
206 std::pair<G4int, Isotope>(38, Isotope(87 206 std::pair<G4int, Isotope>(38, Isotope(87, 7.00)),
207 std::pair<G4int, Isotope>(38, Isotope(88 207 std::pair<G4int, Isotope>(38, Isotope(88, 82.58)),
208 std::pair<G4int, Isotope>(39, Isotope(89 208 std::pair<G4int, Isotope>(39, Isotope(89, 100.0)),
209 std::pair<G4int, Isotope>(40, Isotope(90 209 std::pair<G4int, Isotope>(40, Isotope(90, 51.45)),
210 std::pair<G4int, Isotope>(40, Isotope(91 210 std::pair<G4int, Isotope>(40, Isotope(91, 11.22)),
211 std::pair<G4int, Isotope>(40, Isotope(92 211 std::pair<G4int, Isotope>(40, Isotope(92, 17.15)),
212 std::pair<G4int, Isotope>(40, Isotope(94 212 std::pair<G4int, Isotope>(40, Isotope(94, 17.38)),
213 std::pair<G4int, Isotope>(40, Isotope(96 213 std::pair<G4int, Isotope>(40, Isotope(96, 2.80)),
214 std::pair<G4int, Isotope>(41, Isotope(93 214 std::pair<G4int, Isotope>(41, Isotope(93, 100.0)),
215 std::pair<G4int, Isotope>(42, Isotope(92 215 std::pair<G4int, Isotope>(42, Isotope(92, 14.84)),
216 std::pair<G4int, Isotope>(42, Isotope(94 216 std::pair<G4int, Isotope>(42, Isotope(94, 9.25)),
217 std::pair<G4int, Isotope>(42, Isotope(95 217 std::pair<G4int, Isotope>(42, Isotope(95, 15.92)),
218 std::pair<G4int, Isotope>(42, Isotope(96 218 std::pair<G4int, Isotope>(42, Isotope(96, 16.68)),
219 std::pair<G4int, Isotope>(42, Isotope(97 219 std::pair<G4int, Isotope>(42, Isotope(97, 9.55)),
220 std::pair<G4int, Isotope>(42, Isotope(98 220 std::pair<G4int, Isotope>(42, Isotope(98, 24.13)),
221 std::pair<G4int, Isotope>(42, Isotope(10 221 std::pair<G4int, Isotope>(42, Isotope(100, 9.63)),
222 std::pair<G4int, Isotope>(44, Isotope(96 222 std::pair<G4int, Isotope>(44, Isotope(96, 5.52)),
223 std::pair<G4int, Isotope>(44, Isotope(98 223 std::pair<G4int, Isotope>(44, Isotope(98, 1.88)),
224 std::pair<G4int, Isotope>(44, Isotope(99 224 std::pair<G4int, Isotope>(44, Isotope(99, 12.7)),
225 std::pair<G4int, Isotope>(44, Isotope(10 225 std::pair<G4int, Isotope>(44, Isotope(100, 12.6)),
226 std::pair<G4int, Isotope>(44, Isotope(10 226 std::pair<G4int, Isotope>(44, Isotope(101, 17.0)),
227 std::pair<G4int, Isotope>(44, Isotope(10 227 std::pair<G4int, Isotope>(44, Isotope(102, 31.6)),
228 std::pair<G4int, Isotope>(44, Isotope(10 228 std::pair<G4int, Isotope>(44, Isotope(104, 18.7)),
229 std::pair<G4int, Isotope>(45, Isotope(10 229 std::pair<G4int, Isotope>(45, Isotope(103, 100.0)),
230 std::pair<G4int, Isotope>(46, Isotope(10 230 std::pair<G4int, Isotope>(46, Isotope(102, 1.02)),
231 std::pair<G4int, Isotope>(46, Isotope(10 231 std::pair<G4int, Isotope>(46, Isotope(104, 11.14)),
232 std::pair<G4int, Isotope>(46, Isotope(10 232 std::pair<G4int, Isotope>(46, Isotope(105, 22.33)),
233 std::pair<G4int, Isotope>(46, Isotope(10 233 std::pair<G4int, Isotope>(46, Isotope(106, 27.33)),
234 std::pair<G4int, Isotope>(46, Isotope(10 234 std::pair<G4int, Isotope>(46, Isotope(108, 26.46)),
235 std::pair<G4int, Isotope>(46, Isotope(11 235 std::pair<G4int, Isotope>(46, Isotope(110, 11.72)),
236 std::pair<G4int, Isotope>(47, Isotope(10 236 std::pair<G4int, Isotope>(47, Isotope(107, 51.839)),
237 std::pair<G4int, Isotope>(47, Isotope(10 237 std::pair<G4int, Isotope>(47, Isotope(109, 48.161)),
238 std::pair<G4int, Isotope>(48, Isotope(10 238 std::pair<G4int, Isotope>(48, Isotope(106, 1.25)),
239 std::pair<G4int, Isotope>(48, Isotope(10 239 std::pair<G4int, Isotope>(48, Isotope(108, 0.89)),
240 std::pair<G4int, Isotope>(48, Isotope(11 240 std::pair<G4int, Isotope>(48, Isotope(110, 12.49)),
241 std::pair<G4int, Isotope>(48, Isotope(11 241 std::pair<G4int, Isotope>(48, Isotope(111, 12.80)),
242 std::pair<G4int, Isotope>(48, Isotope(11 242 std::pair<G4int, Isotope>(48, Isotope(112, 24.13)),
243 std::pair<G4int, Isotope>(48, Isotope(11 243 std::pair<G4int, Isotope>(48, Isotope(113, 12.22)),
244 std::pair<G4int, Isotope>(48, Isotope(11 244 std::pair<G4int, Isotope>(48, Isotope(114, 28.73)),
245 std::pair<G4int, Isotope>(48, Isotope(11 245 std::pair<G4int, Isotope>(48, Isotope(116, 7.49)),
246 std::pair<G4int, Isotope>(49, Isotope(11 246 std::pair<G4int, Isotope>(49, Isotope(113, 4.3)),
247 std::pair<G4int, Isotope>(49, Isotope(11 247 std::pair<G4int, Isotope>(49, Isotope(115, 95.7)),
248 std::pair<G4int, Isotope>(50, Isotope(11 248 std::pair<G4int, Isotope>(50, Isotope(112, 0.97)),
249 std::pair<G4int, Isotope>(50, Isotope(11 249 std::pair<G4int, Isotope>(50, Isotope(114, 0.65)),
250 std::pair<G4int, Isotope>(50, Isotope(11 250 std::pair<G4int, Isotope>(50, Isotope(115, 0.34)),
251 std::pair<G4int, Isotope>(50, Isotope(11 251 std::pair<G4int, Isotope>(50, Isotope(116, 14.53)),
252 std::pair<G4int, Isotope>(50, Isotope(11 252 std::pair<G4int, Isotope>(50, Isotope(117, 7.68)),
253 std::pair<G4int, Isotope>(50, Isotope(11 253 std::pair<G4int, Isotope>(50, Isotope(118, 24.23)),
254 std::pair<G4int, Isotope>(50, Isotope(11 254 std::pair<G4int, Isotope>(50, Isotope(119, 8.59)),
255 std::pair<G4int, Isotope>(50, Isotope(12 255 std::pair<G4int, Isotope>(50, Isotope(120, 32.59)),
256 std::pair<G4int, Isotope>(50, Isotope(12 256 std::pair<G4int, Isotope>(50, Isotope(122, 4.63)),
257 std::pair<G4int, Isotope>(50, Isotope(12 257 std::pair<G4int, Isotope>(50, Isotope(124, 5.79)),
258 std::pair<G4int, Isotope>(51, Isotope(12 258 std::pair<G4int, Isotope>(51, Isotope(121, 57.36)),
259 std::pair<G4int, Isotope>(51, Isotope(12 259 std::pair<G4int, Isotope>(51, Isotope(123, 42.64)),
260 std::pair<G4int, Isotope>(52, Isotope(12 260 std::pair<G4int, Isotope>(52, Isotope(120, 0.096)),
261 std::pair<G4int, Isotope>(52, Isotope(12 261 std::pair<G4int, Isotope>(52, Isotope(122, 2.603)),
262 std::pair<G4int, Isotope>(52, Isotope(12 262 std::pair<G4int, Isotope>(52, Isotope(123, 0.908)),
263 std::pair<G4int, Isotope>(52, Isotope(12 263 std::pair<G4int, Isotope>(52, Isotope(124, 4.816)),
264 std::pair<G4int, Isotope>(52, Isotope(12 264 std::pair<G4int, Isotope>(52, Isotope(125, 7.139)),
265 std::pair<G4int, Isotope>(52, Isotope(12 265 std::pair<G4int, Isotope>(52, Isotope(126, 18.95)),
266 std::pair<G4int, Isotope>(52, Isotope(12 266 std::pair<G4int, Isotope>(52, Isotope(128, 31.69)),
267 std::pair<G4int, Isotope>(52, Isotope(13 267 std::pair<G4int, Isotope>(52, Isotope(130, 33.80)),
268 std::pair<G4int, Isotope>(53, Isotope(12 268 std::pair<G4int, Isotope>(53, Isotope(127, 100.0)),
269 std::pair<G4int, Isotope>(54, Isotope(12 269 std::pair<G4int, Isotope>(54, Isotope(124, 0.10)),
270 std::pair<G4int, Isotope>(54, Isotope(12 270 std::pair<G4int, Isotope>(54, Isotope(126, 0.09)),
271 std::pair<G4int, Isotope>(54, Isotope(12 271 std::pair<G4int, Isotope>(54, Isotope(128, 1.91)),
272 std::pair<G4int, Isotope>(54, Isotope(12 272 std::pair<G4int, Isotope>(54, Isotope(129, 26.4)),
273 std::pair<G4int, Isotope>(54, Isotope(13 273 std::pair<G4int, Isotope>(54, Isotope(130, 4.1)),
274 std::pair<G4int, Isotope>(54, Isotope(13 274 std::pair<G4int, Isotope>(54, Isotope(131, 21.2)),
275 std::pair<G4int, Isotope>(54, Isotope(13 275 std::pair<G4int, Isotope>(54, Isotope(132, 26.9)),
276 std::pair<G4int, Isotope>(54, Isotope(13 276 std::pair<G4int, Isotope>(54, Isotope(134, 10.4)),
277 std::pair<G4int, Isotope>(54, Isotope(13 277 std::pair<G4int, Isotope>(54, Isotope(136, 8.9)),
278 std::pair<G4int, Isotope>(55, Isotope(13 278 std::pair<G4int, Isotope>(55, Isotope(133, 100.0)),
279 std::pair<G4int, Isotope>(56, Isotope(13 279 std::pair<G4int, Isotope>(56, Isotope(130, 0.106)),
280 std::pair<G4int, Isotope>(56, Isotope(13 280 std::pair<G4int, Isotope>(56, Isotope(132, 0.101)),
281 std::pair<G4int, Isotope>(56, Isotope(13 281 std::pair<G4int, Isotope>(56, Isotope(134, 2.417)),
282 std::pair<G4int, Isotope>(56, Isotope(13 282 std::pair<G4int, Isotope>(56, Isotope(135, 6.592)),
283 std::pair<G4int, Isotope>(56, Isotope(13 283 std::pair<G4int, Isotope>(56, Isotope(136, 7.854)),
284 std::pair<G4int, Isotope>(56, Isotope(13 284 std::pair<G4int, Isotope>(56, Isotope(137, 11.23)),
285 std::pair<G4int, Isotope>(56, Isotope(13 285 std::pair<G4int, Isotope>(56, Isotope(138, 71.70)),
286 std::pair<G4int, Isotope>(57, Isotope(13 286 std::pair<G4int, Isotope>(57, Isotope(138, 0.0902)),
287 std::pair<G4int, Isotope>(57, Isotope(13 287 std::pair<G4int, Isotope>(57, Isotope(139, 99.9098)),
288 std::pair<G4int, Isotope>(58, Isotope(13 288 std::pair<G4int, Isotope>(58, Isotope(136, 0.19)),
289 std::pair<G4int, Isotope>(58, Isotope(13 289 std::pair<G4int, Isotope>(58, Isotope(138, 0.25)),
290 std::pair<G4int, Isotope>(58, Isotope(14 290 std::pair<G4int, Isotope>(58, Isotope(140, 88.48)),
291 std::pair<G4int, Isotope>(58, Isotope(14 291 std::pair<G4int, Isotope>(58, Isotope(142, 11.08)),
292 std::pair<G4int, Isotope>(59, Isotope(14 292 std::pair<G4int, Isotope>(59, Isotope(141, 100.0)),
293 std::pair<G4int, Isotope>(60, Isotope(14 293 std::pair<G4int, Isotope>(60, Isotope(142, 27.13)),
294 std::pair<G4int, Isotope>(60, Isotope(14 294 std::pair<G4int, Isotope>(60, Isotope(143, 12.18)),
295 std::pair<G4int, Isotope>(60, Isotope(14 295 std::pair<G4int, Isotope>(60, Isotope(144, 23.80)),
296 std::pair<G4int, Isotope>(60, Isotope(14 296 std::pair<G4int, Isotope>(60, Isotope(145, 8.30)),
297 std::pair<G4int, Isotope>(60, Isotope(14 297 std::pair<G4int, Isotope>(60, Isotope(146, 17.19)),
298 std::pair<G4int, Isotope>(60, Isotope(14 298 std::pair<G4int, Isotope>(60, Isotope(148, 5.76)),
299 std::pair<G4int, Isotope>(60, Isotope(15 299 std::pair<G4int, Isotope>(60, Isotope(150, 5.64)),
300 std::pair<G4int, Isotope>(62, Isotope(14 300 std::pair<G4int, Isotope>(62, Isotope(144, 3.1)),
301 std::pair<G4int, Isotope>(62, Isotope(14 301 std::pair<G4int, Isotope>(62, Isotope(147, 15.0)),
302 std::pair<G4int, Isotope>(62, Isotope(14 302 std::pair<G4int, Isotope>(62, Isotope(148, 11.3)),
303 std::pair<G4int, Isotope>(62, Isotope(14 303 std::pair<G4int, Isotope>(62, Isotope(149, 13.8)),
304 std::pair<G4int, Isotope>(62, Isotope(15 304 std::pair<G4int, Isotope>(62, Isotope(150, 7.4)),
305 std::pair<G4int, Isotope>(62, Isotope(15 305 std::pair<G4int, Isotope>(62, Isotope(152, 26.7)),
306 std::pair<G4int, Isotope>(62, Isotope(15 306 std::pair<G4int, Isotope>(62, Isotope(154, 22.7)),
307 std::pair<G4int, Isotope>(63, Isotope(15 307 std::pair<G4int, Isotope>(63, Isotope(151, 47.8)),
308 std::pair<G4int, Isotope>(63, Isotope(15 308 std::pair<G4int, Isotope>(63, Isotope(153, 52.2)),
309 std::pair<G4int, Isotope>(64, Isotope(15 309 std::pair<G4int, Isotope>(64, Isotope(152, 0.20)),
310 std::pair<G4int, Isotope>(64, Isotope(15 310 std::pair<G4int, Isotope>(64, Isotope(154, 2.18)),
311 std::pair<G4int, Isotope>(64, Isotope(15 311 std::pair<G4int, Isotope>(64, Isotope(155, 14.80)),
312 std::pair<G4int, Isotope>(64, Isotope(15 312 std::pair<G4int, Isotope>(64, Isotope(156, 20.47)),
313 std::pair<G4int, Isotope>(64, Isotope(15 313 std::pair<G4int, Isotope>(64, Isotope(157, 15.65)),
314 std::pair<G4int, Isotope>(64, Isotope(15 314 std::pair<G4int, Isotope>(64, Isotope(158, 24.84)),
315 std::pair<G4int, Isotope>(64, Isotope(16 315 std::pair<G4int, Isotope>(64, Isotope(160, 21.86)),
316 std::pair<G4int, Isotope>(65, Isotope(15 316 std::pair<G4int, Isotope>(65, Isotope(159, 100.0)),
317 std::pair<G4int, Isotope>(66, Isotope(15 317 std::pair<G4int, Isotope>(66, Isotope(156, 0.06)),
318 std::pair<G4int, Isotope>(66, Isotope(15 318 std::pair<G4int, Isotope>(66, Isotope(158, 0.10)),
319 std::pair<G4int, Isotope>(66, Isotope(16 319 std::pair<G4int, Isotope>(66, Isotope(160, 2.34)),
320 std::pair<G4int, Isotope>(66, Isotope(16 320 std::pair<G4int, Isotope>(66, Isotope(161, 18.9)),
321 std::pair<G4int, Isotope>(66, Isotope(16 321 std::pair<G4int, Isotope>(66, Isotope(162, 25.5)),
322 std::pair<G4int, Isotope>(66, Isotope(16 322 std::pair<G4int, Isotope>(66, Isotope(163, 24.9)),
323 std::pair<G4int, Isotope>(66, Isotope(16 323 std::pair<G4int, Isotope>(66, Isotope(164, 28.2)),
324 std::pair<G4int, Isotope>(67, Isotope(16 324 std::pair<G4int, Isotope>(67, Isotope(165, 100.0)),
325 std::pair<G4int, Isotope>(68, Isotope(16 325 std::pair<G4int, Isotope>(68, Isotope(162, 0.14)),
326 std::pair<G4int, Isotope>(68, Isotope(16 326 std::pair<G4int, Isotope>(68, Isotope(164, 1.61)),
327 std::pair<G4int, Isotope>(68, Isotope(16 327 std::pair<G4int, Isotope>(68, Isotope(166, 33.6)),
328 std::pair<G4int, Isotope>(68, Isotope(16 328 std::pair<G4int, Isotope>(68, Isotope(167, 22.95)),
329 std::pair<G4int, Isotope>(68, Isotope(16 329 std::pair<G4int, Isotope>(68, Isotope(168, 26.8)),
330 std::pair<G4int, Isotope>(68, Isotope(17 330 std::pair<G4int, Isotope>(68, Isotope(170, 14.9)),
331 std::pair<G4int, Isotope>(69, Isotope(16 331 std::pair<G4int, Isotope>(69, Isotope(169, 100.0)),
332 std::pair<G4int, Isotope>(70, Isotope(16 332 std::pair<G4int, Isotope>(70, Isotope(168, 0.13)),
333 std::pair<G4int, Isotope>(70, Isotope(17 333 std::pair<G4int, Isotope>(70, Isotope(170, 3.05)),
334 std::pair<G4int, Isotope>(70, Isotope(17 334 std::pair<G4int, Isotope>(70, Isotope(171, 14.3)),
335 std::pair<G4int, Isotope>(70, Isotope(17 335 std::pair<G4int, Isotope>(70, Isotope(172, 21.9)),
336 std::pair<G4int, Isotope>(70, Isotope(17 336 std::pair<G4int, Isotope>(70, Isotope(173, 16.12)),
337 std::pair<G4int, Isotope>(70, Isotope(17 337 std::pair<G4int, Isotope>(70, Isotope(174, 31.8)),
338 std::pair<G4int, Isotope>(70, Isotope(17 338 std::pair<G4int, Isotope>(70, Isotope(176, 12.7)),
339 std::pair<G4int, Isotope>(71, Isotope(17 339 std::pair<G4int, Isotope>(71, Isotope(175, 97.41)),
340 std::pair<G4int, Isotope>(71, Isotope(17 340 std::pair<G4int, Isotope>(71, Isotope(176, 2.59)),
341 std::pair<G4int, Isotope>(72, Isotope(17 341 std::pair<G4int, Isotope>(72, Isotope(174, 0.162)),
342 std::pair<G4int, Isotope>(72, Isotope(17 342 std::pair<G4int, Isotope>(72, Isotope(176, 5.206)),
343 std::pair<G4int, Isotope>(72, Isotope(17 343 std::pair<G4int, Isotope>(72, Isotope(177, 18.606)),
344 std::pair<G4int, Isotope>(72, Isotope(17 344 std::pair<G4int, Isotope>(72, Isotope(178, 27.297)),
345 std::pair<G4int, Isotope>(72, Isotope(17 345 std::pair<G4int, Isotope>(72, Isotope(179, 13.629)),
346 std::pair<G4int, Isotope>(72, Isotope(18 346 std::pair<G4int, Isotope>(72, Isotope(180, 35.100)),
347 std::pair<G4int, Isotope>(73, Isotope(18 347 std::pair<G4int, Isotope>(73, Isotope(180, 0.012)),
348 std::pair<G4int, Isotope>(73, Isotope(18 348 std::pair<G4int, Isotope>(73, Isotope(181, 99.988)),
349 std::pair<G4int, Isotope>(74, Isotope(18 349 std::pair<G4int, Isotope>(74, Isotope(180, 0.13)),
350 std::pair<G4int, Isotope>(74, Isotope(18 350 std::pair<G4int, Isotope>(74, Isotope(182, 26.3)),
351 std::pair<G4int, Isotope>(74, Isotope(18 351 std::pair<G4int, Isotope>(74, Isotope(183, 14.3)),
352 std::pair<G4int, Isotope>(74, Isotope(18 352 std::pair<G4int, Isotope>(74, Isotope(184, 30.67)),
353 std::pair<G4int, Isotope>(74, Isotope(18 353 std::pair<G4int, Isotope>(74, Isotope(186, 28.6)),
354 std::pair<G4int, Isotope>(75, Isotope(18 354 std::pair<G4int, Isotope>(75, Isotope(185, 37.40)),
355 std::pair<G4int, Isotope>(75, Isotope(18 355 std::pair<G4int, Isotope>(75, Isotope(187, 62.60)),
356 std::pair<G4int, Isotope>(76, Isotope(18 356 std::pair<G4int, Isotope>(76, Isotope(184, 0.02)),
357 std::pair<G4int, Isotope>(76, Isotope(18 357 std::pair<G4int, Isotope>(76, Isotope(186, 1.58)),
358 std::pair<G4int, Isotope>(76, Isotope(18 358 std::pair<G4int, Isotope>(76, Isotope(187, 1.6)),
359 std::pair<G4int, Isotope>(76, Isotope(18 359 std::pair<G4int, Isotope>(76, Isotope(188, 13.3)),
360 std::pair<G4int, Isotope>(76, Isotope(18 360 std::pair<G4int, Isotope>(76, Isotope(189, 16.1)),
361 std::pair<G4int, Isotope>(76, Isotope(19 361 std::pair<G4int, Isotope>(76, Isotope(190, 26.4)),
362 std::pair<G4int, Isotope>(76, Isotope(19 362 std::pair<G4int, Isotope>(76, Isotope(192, 41.0)),
363 std::pair<G4int, Isotope>(77, Isotope(19 363 std::pair<G4int, Isotope>(77, Isotope(191, 37.3)),
364 std::pair<G4int, Isotope>(77, Isotope(19 364 std::pair<G4int, Isotope>(77, Isotope(193, 62.7)),
365 std::pair<G4int, Isotope>(78, Isotope(19 365 std::pair<G4int, Isotope>(78, Isotope(190, 0.01)),
366 std::pair<G4int, Isotope>(78, Isotope(19 366 std::pair<G4int, Isotope>(78, Isotope(192, 0.79)),
367 std::pair<G4int, Isotope>(78, Isotope(19 367 std::pair<G4int, Isotope>(78, Isotope(194, 32.9)),
368 std::pair<G4int, Isotope>(78, Isotope(19 368 std::pair<G4int, Isotope>(78, Isotope(195, 33.8)),
369 std::pair<G4int, Isotope>(78, Isotope(19 369 std::pair<G4int, Isotope>(78, Isotope(196, 25.3)),
370 std::pair<G4int, Isotope>(78, Isotope(19 370 std::pair<G4int, Isotope>(78, Isotope(198, 7.2)),
371 std::pair<G4int, Isotope>(79, Isotope(19 371 std::pair<G4int, Isotope>(79, Isotope(197, 100.0)),
372 std::pair<G4int, Isotope>(80, Isotope(19 372 std::pair<G4int, Isotope>(80, Isotope(196, 0.15)),
373 std::pair<G4int, Isotope>(80, Isotope(19 373 std::pair<G4int, Isotope>(80, Isotope(198, 9.97)),
374 std::pair<G4int, Isotope>(80, Isotope(19 374 std::pair<G4int, Isotope>(80, Isotope(199, 16.87)),
375 std::pair<G4int, Isotope>(80, Isotope(20 375 std::pair<G4int, Isotope>(80, Isotope(200, 23.10)),
376 std::pair<G4int, Isotope>(80, Isotope(20 376 std::pair<G4int, Isotope>(80, Isotope(201, 13.18)),
377 std::pair<G4int, Isotope>(80, Isotope(20 377 std::pair<G4int, Isotope>(80, Isotope(202, 29.86)),
378 std::pair<G4int, Isotope>(80, Isotope(20 378 std::pair<G4int, Isotope>(80, Isotope(204, 6.87)),
379 std::pair<G4int, Isotope>(81, Isotope(20 379 std::pair<G4int, Isotope>(81, Isotope(203, 29.524)),
380 std::pair<G4int, Isotope>(81, Isotope(20 380 std::pair<G4int, Isotope>(81, Isotope(205, 70.476)),
381 std::pair<G4int, Isotope>(82, Isotope(20 381 std::pair<G4int, Isotope>(82, Isotope(204, 1.4)),
382 std::pair<G4int, Isotope>(82, Isotope(20 382 std::pair<G4int, Isotope>(82, Isotope(206, 24.1)),
383 std::pair<G4int, Isotope>(82, Isotope(20 383 std::pair<G4int, Isotope>(82, Isotope(207, 22.1)),
384 std::pair<G4int, Isotope>(82, Isotope(20 384 std::pair<G4int, Isotope>(82, Isotope(208, 52.4)),
385 std::pair<G4int, Isotope>(83, Isotope(20 385 std::pair<G4int, Isotope>(83, Isotope(209, 100.0)),
386 std::pair<G4int, Isotope>(90, Isotope(23 386 std::pair<G4int, Isotope>(90, Isotope(232, 100.0)),
387 std::pair<G4int, Isotope>(92, Isotope(23 387 std::pair<G4int, Isotope>(92, Isotope(234, 0.0055)),
388 std::pair<G4int, Isotope>(92, Isotope(23 388 std::pair<G4int, Isotope>(92, Isotope(235, 0.7200)),
389 std::pair<G4int, Isotope>(92, Isotope(23 389 std::pair<G4int, Isotope>(92, Isotope(238, 99.2745))
390 }; 390 };
391 391
392 // Cool hack to get the size of an array i 392 // Cool hack to get the size of an array in C++
393 template<typename T, ::std::size_t N> ::st 393 template<typename T, ::std::size_t N> ::std::size_t sizeOfArray(const T(&)[ N ] ) {
394 return N; 394 return N;
395 } 395 }
396 } 396 }
397 397
398 NaturalIsotopicDistributions::NaturalIsotopi 398 NaturalIsotopicDistributions::NaturalIsotopicDistributions() {
399 G4int oldZ = -1; 399 G4int oldZ = -1;
400 IsotopeVector aVector; 400 IsotopeVector aVector;
401 for(unsigned int i=0; i<sizeOfArray(theRaw 401 for(unsigned int i=0; i<sizeOfArray(theRawDistributions); ++i) {
402 std::pair<G4int, Isotope> const &aPair = 402 std::pair<G4int, Isotope> const &aPair = theRawDistributions[i];
403 if(aPair.first == oldZ) { 403 if(aPair.first == oldZ) {
404 aVector.push_back(aPair.second); 404 aVector.push_back(aPair.second);
405 } else { 405 } else {
406 if(oldZ!=-1) 406 if(oldZ!=-1)
407 theDistributions.insert(std::pair<G4 407 theDistributions.insert(std::pair<G4int, IsotopicDistribution>(oldZ, IsotopicDistribution(aVector)));
408 oldZ = aPair.first; 408 oldZ = aPair.first;
409 aVector.clear(); 409 aVector.clear();
410 aVector.push_back(aPair.second); 410 aVector.push_back(aPair.second);
411 } 411 }
412 } 412 }
413 // last element 413 // last element
414 theDistributions.insert(std::pair<G4int, I 414 theDistributions.insert(std::pair<G4int, IsotopicDistribution>(oldZ, IsotopicDistribution(aVector)));
415 } 415 }
416 416
417 } 417 }
418 418
419 419