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