source: trunk/source/processes/electromagnetic/lowenergy/src/G4hShellCrossSectionExpData.cc@ 1197

Last change on this file since 1197 was 1196, checked in by garnier, 16 years ago

update CVS release candidate geant4.9.3.01
File size: 12.4 KB
Line 
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//
27//
28// Author: Simona Saliceti (simona.saliceti@ge.infn.it)
29//
30// History:
31// -----------
32// 22 Apr 2004 First committed to cvs
33//
34// -------------------------------------------------------------------
35// $Id: G4hShellCrossSectionExpData.cc,v 1.4 2009/06/10 13:32:36 mantero Exp $
36// GEANT4 tag $Name: geant4-09-03-cand-01 $
37
38#include "G4hShellCrossSectionExpData.hh"
39#include <fstream>
40
41G4hShellCrossSectionExpData::G4hShellCrossSectionExpData()
42{
43
44
45 FillVectorValues();
46
47 FillParameterMap();
48
49}
50
51G4hShellCrossSectionExpData::~G4hShellCrossSectionExpData()
52{
53 std::map< G4int,std::vector<G4double>*,std::less<G4int> >::iterator pos;
54 for (pos = parameterMap.begin(); pos != parameterMap.end(); ++pos)
55 {
56 std::vector<G4double>* dataSet = (*pos).second;
57 delete dataSet;
58 }
59}
60
61inline void G4hShellCrossSectionExpData::InitializeVector(std::vector<G4double> &vect, G4double value1, G4double value2, G4double value3) const
62{
63 vect.push_back(value1);
64 vect.push_back(value2);
65 vect.push_back(value3);
66}
67
68void G4hShellCrossSectionExpData::FillVectorValues()
69{
70 InitializeVector(parameter6C, 25965819., 0.1974268719, -0.9746342);
71 InitializeVector(parameter7N, 14254204., 0.09584816998, -0.8908825);
72 InitializeVector(parameter8O, 14257000., 0.9582500E-01, -0.8909900);
73 InitializeVector(parameter9F, 5628098., 0.2221402E-01, -0.7857742);
74 InitializeVector(parameter10Ne, 3866500., 0.9298651E-02, -0.7569925);
75 InitializeVector(parameter11Na, 2816614., 0.3161790E-02, -0.7429561);
76 InitializeVector(parameter12Mg, 2113709., 0.9267687E-03, -0.7342144);
77 InitializeVector(parameter13Al, 1647127., 0.2344573E-03, -0.7305032);
78 InitializeVector(parameter14Si, 1267490., 0.5787061E-04, -0.7184879);
79 InitializeVector(parameter15P, 1035792., 0.1101495E-04, -0.7199858);
80 InitializeVector(parameter16S, 825871.6, 0.2326387E-05, -0.7095044);
81 InitializeVector(parameter17Cl, 13.41653, -14.66877, -0.7017217);
82 InitializeVector(parameter18Ar, 555038.1, 0.6543986E-07, -0.6968985);
83 InitializeVector(parameter19K, 469063.8, 0.8876799E-08, -0.6940991);
84 InitializeVector(parameter20Ca, 405918.4, 0.9881247E-09, -0.6953619);
85 InitializeVector(parameter21Sc, 351611.0, 0.1054022E-09, -0.6950487);
86 InitializeVector(parameter22Ti, 313223.2, 0.8541349E-11, -0.6984913);
87 InitializeVector(parameter23V, 257752.3, 0.1092866E-11, -0.6841001);
88 InitializeVector(parameter24Cr, 230064.8, 0.7874358E-13, -0.6862135);
89 InitializeVector(parameter25Mn, 151740.8, 0.7988999E-13, -0.6279536);
90 InitializeVector(parameter26Fe, 12.21083, -35.54018, -0.7140945);
91 InitializeVector(parameter27Co, 12.14600, -38.73178, -0.7229655);
92 InitializeVector(parameter28Ni, 12.10570, -42.25900, -0.7354265);
93 InitializeVector(parameter29Cu, 11.87911, -44.37304, -0.7089269);
94 InitializeVector(parameter30Zn, 11.86411, -48.35208, -0.7252861);
95 InitializeVector(parameter31Ga, 11.66042, -50.73798, -0.7029260);
96 InitializeVector(parameter32Ge, 11.66824, -55.26208, -0.7218891);
97 InitializeVector(parameter33As, 11.47804, -57.94218, -0.7016547);
98 InitializeVector(parameter34Se, 11.30017, -60.70304, -0.6836499);
99 InitializeVector(parameter35Br, 11.32571, -65.77174, -0.7042627);
100 InitializeVector(parameter36Kr, 10.44025, -61.74093, -0.5383769);
101 InitializeVector(parameter37Rb, 10.27499, -64.44025, -0.5236434);
102 InitializeVector(parameter38Sr, 10.31786, -69.83844, -0.5475903);
103 InitializeVector(parameter39Y, 10.14157, -72.55550, -0.5302060);
104 InitializeVector(parameter40Zr, 10.03310, -76.52999, -0.5243382);
105 InitializeVector(parameter41Nb, 10.03804, -81.38335, -0.5404782);
106 InitializeVector(parameter42Mo, 9.872166, -84.24110, -0.5245091);
107 InitializeVector(parameter43Tc, 9.941685, -90.72401, -0.5503877);
108 InitializeVector(parameter44Ru, 9.771618, -93.50847, -0.5334178);
109 InitializeVector(parameter45Rh, 9.819192, -99.99145, -0.5543751);
110 InitializeVector(parameter46Pd, 9.680034, -103.5450, -0.5387755);
111 InitializeVector(parameter47Ag, 9.703977, -109.8104, -0.5550457);
112 InitializeVector(parameter48Cd, 9.548004, -113.0565, -0.5399088);
113 InitializeVector(parameter49In, 9.378392, -115.8466, -0.5226109);
114 InitializeVector(parameter50Sn, 9.406608, -122.6137, -0.5386317);
115 InitializeVector(parameter51Sb, 9.254029, -125.8533, -0.5238908);
116 InitializeVector(parameter52Te, 9.311451, -133.6049, -0.5440237);
117 InitializeVector(parameter53I, 9.158569, -136.7913, -0.5291739);
118 InitializeVector(parameter54Xe, 9.011420, -140.1908, -0.5150210);
119 InitializeVector(parameter55Cs, 9.112158, -149.6160, -0.5411919);
120 InitializeVector(parameter56Ba, 7.696734, -124.4350, -0.3181673);
121 InitializeVector(parameter57La, 7.728227, -130.9067, -0.3348193);
122 InitializeVector(parameter58Ce, 7.595518, -134.1356, -0.3233585);
123 InitializeVector(parameter59Pr, 7.451400, -136.9318, -0.3102762);
124 InitializeVector(parameter60Nd, 7.493138, -144.0968, -0.3273972);
125 InitializeVector(parameter61Pm, 7.358854, -147.1354, -0.3155824);
126 InitializeVector(parameter62Sm, 7.231769, -150.3722, -0.3048225);
127 InitializeVector(parameter63Eu, 7.060385, -152.3517, -0.2873322);
128 InitializeVector(parameter64Gd, 7.179881, -161.9544, -0.3162093);
129 InitializeVector(parameter65Tb, 7.048202, -164.9670, -0.3046019);
130 InitializeVector(parameter66Dy, 1.962871, -4.620275, 0.5748519E-02);
131 InitializeVector(parameter67Ho, 1.992797, -4.863559, 0.5772170E-02);
132 InitializeVector(parameter68Er, 2.013621, -5.068997, 0.5749294E-02);
133 InitializeVector(parameter69Tm, 2.026489, -5.243042, 0.5671808E-02);
134 InitializeVector(parameter70Yb, 1.946252, -5.047087, 0.5011610E-02);
135 InitializeVector(parameter71Lu, 1.968246, -5.253083, 0.5007004E-02);
136 InitializeVector(parameter72Hf, 1.980888, -5.421401, 0.4949681E-02);
137 InitializeVector(parameter73Ta, 1.925043, -5.310529, 0.4504530E-02);
138 InitializeVector(parameter74W, 1.939689, -5.484359, 0.4464923E-02);
139 InitializeVector(parameter75Re, 1.951746, -5.645718, 0.4415765E-02);
140 InitializeVector(parameter76Os, 2.249802, -6.974045, 0.5965646E-02);
141 InitializeVector(parameter77Ir, 2.257475, -7.121614, 0.5873831E-02);
142 InitializeVector(parameter78Pt, 2.262231, -7.255240, 0.5769194E-02);
143 InitializeVector(parameter79Au, 2.188981, -7.080360, 0.5161316E-02);
144 InitializeVector(parameter80Hg, 2.199681, -7.233214, 0.5111361E-02);
145 InitializeVector(parameter81Tl, 2.211119, -7.389515, 0.5064001E-02);
146 InitializeVector(parameter82Pb, 2.204573, -7.468263, 0.4934075E-02);
147 InitializeVector(parameter83Bi, 2.148083, -7.348053, 0.4494417E-02);
148 InitializeVector(parameter84Po, 2.154512, -7.479838, 0.4433237E-02);
149 InitializeVector(parameter85At, 2.164621, -7.626572, 0.4388465E-02);
150 InitializeVector(parameter86Rn, 2.288465, -8.212229, 0.5807703E-02);
151 InitializeVector(parameter87Fr, 2.294291, -8.336506, 0.5753172E-02);
152 InitializeVector(parameter88Ra, 2.318430, -8.525764, 0.5849121E-02);
153 InitializeVector(parameter89Ac, 2.306637, -8.578354, 0.5672658E-02);
154 InitializeVector(parameter90Th, 2.254484, -8.479531, 0.5163647E-02);
155 InitializeVector(parameter91Pa, 2.246310, -8.540535, 0.5044748E-02);
156 InitializeVector(parameter92U, 2.276106, -8.748476, 0.5186345E-02);
157}
158
159void G4hShellCrossSectionExpData::FillParameterMap()
160{
161 parameterMap [6] = &parameter6C ;
162 parameterMap [7] = &parameter7N ;
163 parameterMap [8] = &parameter8O ;
164 parameterMap [9] = &parameter9F ;
165 parameterMap [10] = &parameter10Ne;
166 parameterMap [11] = &parameter11Na;
167 parameterMap [12] = &parameter12Mg;
168 parameterMap [13] = &parameter13Al;
169 parameterMap [14] = &parameter14Si;
170 parameterMap [15] = &parameter15P ;
171 parameterMap [16] = &parameter16S ;
172 parameterMap [17] = &parameter17Cl;
173 parameterMap [18] = &parameter18Ar;
174 parameterMap [19] = &parameter19K ;
175 parameterMap [20] = &parameter20Ca;
176 parameterMap [21] = &parameter21Sc;
177 parameterMap [22] = &parameter22Ti;
178 parameterMap [23] = &parameter23V ;
179 parameterMap [24] = &parameter24Cr;
180 parameterMap [25] = &parameter25Mn;
181 parameterMap [26] = &parameter26Fe;
182 parameterMap [27] = &parameter27Co;
183 parameterMap [28] = &parameter28Ni;
184 parameterMap [29] = &parameter29Cu;
185 parameterMap [30] = &parameter30Zn;
186 parameterMap [31] = &parameter31Ga;
187 parameterMap [32] = &parameter32Ge;
188 parameterMap [33] = &parameter33As;
189 parameterMap [34] = &parameter34Se;
190 parameterMap [35] = &parameter35Br;
191 parameterMap [36] = &parameter36Kr;
192 parameterMap [37] = &parameter37Rb;
193 parameterMap [38] = &parameter38Sr;
194 parameterMap [39] = &parameter39Y ;
195 parameterMap [40] = &parameter40Zr;
196 parameterMap [41] = &parameter41Nb;
197 parameterMap [42] = &parameter42Mo;
198 parameterMap [43] = &parameter43Tc;
199 parameterMap [44] = &parameter44Ru;
200 parameterMap [45] = &parameter45Rh;
201 parameterMap [46] = &parameter46Pd;
202 parameterMap [47] = &parameter47Ag;
203 parameterMap [48] = &parameter48Cd;
204 parameterMap [49] = &parameter49In;
205 parameterMap [50] = &parameter50Sn;
206 parameterMap [51] = &parameter51Sb;
207 parameterMap [52] = &parameter52Te;
208 parameterMap [53] = &parameter53I ;
209 parameterMap [54] = &parameter54Xe;
210 parameterMap [55] = &parameter55Cs;
211 parameterMap [56] = &parameter56Ba;
212 parameterMap [57] = &parameter57La;
213 parameterMap [58] = &parameter58Ce;
214 parameterMap [59] = &parameter59Pr;
215 parameterMap [60] = &parameter60Nd;
216 parameterMap [61] = &parameter61Pm;
217 parameterMap [62] = &parameter62Sm;
218 parameterMap [63] = &parameter63Eu;
219 parameterMap [64] = &parameter64Gd;
220 parameterMap [65] = &parameter65Tb;
221 parameterMap [66] = &parameter66Dy;
222 parameterMap [67] = &parameter67Ho;
223 parameterMap [68] = &parameter68Er;
224 parameterMap [69] = &parameter69Tm;
225 parameterMap [70] = &parameter70Yb;
226 parameterMap [71] = &parameter71Lu;
227 parameterMap [72] = &parameter72Hf;
228 parameterMap [73] = &parameter73Ta;
229 parameterMap [74] = &parameter74W ;
230 parameterMap [75] = &parameter75Re;
231 parameterMap [76] = &parameter76Os;
232 parameterMap [77] = &parameter77Ir;
233 parameterMap [78] = &parameter78Pt;
234 parameterMap [79] = &parameter79Au;
235 parameterMap [80] = &parameter80Hg;
236 parameterMap [81] = &parameter81Tl;
237 parameterMap [82] = &parameter82Pb;
238 parameterMap [83] = &parameter83Bi;
239 parameterMap [84] = &parameter84Po;
240 parameterMap [85] = &parameter85At;
241 parameterMap [86] = &parameter86Rn;
242 parameterMap [87] = &parameter87Fr;
243 parameterMap [88] = &parameter88Ra;
244 parameterMap [89] = &parameter89Ac;
245 parameterMap [90] = &parameter90Th;
246 parameterMap [91] = &parameter91Pa;
247 parameterMap [92] = &parameter92U ;
248}
249
250std::vector<G4double>* G4hShellCrossSectionExpData::GetParam(G4int Z)
251{
252 return parameterMap[Z];
253}
Note: See TracBrowser for help on using the repository browser.