source: trunk/source/processes/electromagnetic/lowenergy/src/G4hShellCrossSectionDoubleExpData.cc@ 968

Last change on this file since 968 was 961, checked in by garnier, 17 years ago

update processes
File size: 34.4 KB
RevLine 
[819]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: G4hShellCrossSectionDoubleExpData.cc,v 1.4 2006/06/29 19:42:38 gunter Exp $
[961]36// GEANT4 tag $Name: geant4-09-02-ref-02 $
[819]37
38#include "G4hShellCrossSectionDoubleExpData.hh"
39
40G4hShellCrossSectionDoubleExpData::G4hShellCrossSectionDoubleExpData()
41{
42
43 FillVectorValuesEnergy();
44 FillVectorValuesPar1();
45 FillVectorValuesPar2();
46 FillParameterMapEnergy();
47 FillParameterMapPar1();
48 FillParameterMapPar2();
49
50}
51
52G4hShellCrossSectionDoubleExpData::~G4hShellCrossSectionDoubleExpData()
53{
54
55 std::map< G4int,std::vector<G4double>*,std::less<G4int> >::iterator pos1;
56 for (pos1 = parameterMapEnergy.begin(); pos1 != parameterMapEnergy.end(); ++pos1)
57 {
58 std::vector<G4double>* dataSet1 = (*pos1).second;
59 delete dataSet1;
60 }
61
62 std::map< G4int,std::vector<G4double>*,std::less<G4int> >::iterator pos2;
63 for (pos2 = parameterMapPar1.begin(); pos2 != parameterMapPar1.end(); ++pos2)
64 {
65 std::vector<G4double>* dataSet2 = (*pos2).second;
66 delete dataSet2;
67 }
68
69 std::map< G4int,std::vector<G4double>*,std::less<G4int> >::iterator pos3;
70 for (pos3 = parameterMapPar2.begin(); pos3 != parameterMapPar2.end(); ++pos3)
71 {
72 std::vector<G4double>* dataSet3 = (*pos3).second;
73 delete dataSet3;
74 }
75}
76
77inline void G4hShellCrossSectionDoubleExpData::InitializeVectorEnergy(std::vector<G4double> &vectEnergy, G4double value) const
78{
79 vectEnergy.push_back(value);
80}
81
82void G4hShellCrossSectionDoubleExpData::FillVectorValuesEnergy()
83{
84 InitializeVectorEnergy(energy6C, 0.35);
85 InitializeVectorEnergy(energy7N, 0.6);
86 InitializeVectorEnergy(energy8O, 0.7);
87 InitializeVectorEnergy(energy9F, 1.0);
88 InitializeVectorEnergy(energy10Ne,1.0);
89 InitializeVectorEnergy(energy11Na,1.0);
90 InitializeVectorEnergy(energy12Mg,2.0);
91 InitializeVectorEnergy(energy13Al,2.0);
92 InitializeVectorEnergy(energy14Si,3.0);
93 InitializeVectorEnergy(energy15P, 3.0);
94 InitializeVectorEnergy(energy16S, 5.0);
95 InitializeVectorEnergy(energy17Cl,5.0);
96 InitializeVectorEnergy(energy18Ar,5.0);
97 InitializeVectorEnergy(energy19K, 5.0);
98 InitializeVectorEnergy(energy20Ca,5.0);
99 InitializeVectorEnergy(energy21Sc,5.0);
100 InitializeVectorEnergy(energy22Ti,6.0);
101 InitializeVectorEnergy(energy23V, 6.0);
102 InitializeVectorEnergy(energy24Cr,6.0);
103 InitializeVectorEnergy(energy25Mn,6.0);
104 InitializeVectorEnergy(energy26Fe,3.0);
105 InitializeVectorEnergy(energy27Co,3.5);
106 InitializeVectorEnergy(energy28Ni,3.5);
107 InitializeVectorEnergy(energy29Cu,4.0);
108 InitializeVectorEnergy(energy30Zn,4.0);
109 InitializeVectorEnergy(energy31Ga,4.5);
110 InitializeVectorEnergy(energy32Ge,5.0);
111 InitializeVectorEnergy(energy33As,5.0);
112 InitializeVectorEnergy(energy34Se,5.5);
113 InitializeVectorEnergy(energy35Br,6.0);
114 InitializeVectorEnergy(energy36Kr,6.0);
115 InitializeVectorEnergy(energy37Rb,6.0);
116 InitializeVectorEnergy(energy38Sr,7.0);
117 InitializeVectorEnergy(energy39Y, 7.0);
118 InitializeVectorEnergy(energy40Zr,8.0);
119 InitializeVectorEnergy(energy41Nb,8.0);
120 InitializeVectorEnergy(energy42Mo,8.0);
121 InitializeVectorEnergy(energy43Tc,9.0);
122 InitializeVectorEnergy(energy44Ru,9.0);
123 InitializeVectorEnergy(energy45Rh,9.0);
124 InitializeVectorEnergy(energy46Pd,9.0);
125 InitializeVectorEnergy(energy47Ag,10.);
126 InitializeVectorEnergy(energy48Cd,10.);
127 InitializeVectorEnergy(energy49In,11.);
128 InitializeVectorEnergy(energy50Sn,12.);
129 InitializeVectorEnergy(energy51Sb,12.);
130 InitializeVectorEnergy(energy52Te,12.);
131 InitializeVectorEnergy(energy53I, 14.);
132 InitializeVectorEnergy(energy54Xe,14.);
133 InitializeVectorEnergy(energy55Cs,16.);
134 InitializeVectorEnergy(energy56Ba,16.);
135 InitializeVectorEnergy(energy57La,16.);
136 InitializeVectorEnergy(energy58Ce,18.);
137 InitializeVectorEnergy(energy59Pr,18.);
138 InitializeVectorEnergy(energy60Nd,20.);
139 InitializeVectorEnergy(energy61Pm,20.);
140 InitializeVectorEnergy(energy62Sm,20.);
141 InitializeVectorEnergy(energy63Eu,23.);
142 InitializeVectorEnergy(energy64Gd,23.);
143 InitializeVectorEnergy(energy65Tb,23.);
144 InitializeVectorEnergy(energy66Dy,23.);
145 InitializeVectorEnergy(energy67Ho,23.);
146 InitializeVectorEnergy(energy68Er,23.);
147 InitializeVectorEnergy(energy69Tm,26.);
148 InitializeVectorEnergy(energy70Yb,30.);
149 InitializeVectorEnergy(energy71Lu,30.);
150 InitializeVectorEnergy(energy72Hf,30.);
151 InitializeVectorEnergy(energy73Ta,26.);
152 InitializeVectorEnergy(energy74W, 30.);
153 InitializeVectorEnergy(energy75Re,30.);
154 InitializeVectorEnergy(energy76Os,30.);
155 InitializeVectorEnergy(energy77Ir,30.);
156 InitializeVectorEnergy(energy78Pt,35.);
157 InitializeVectorEnergy(energy79Au,35.);
158 InitializeVectorEnergy(energy80Hg,35.);
159 InitializeVectorEnergy(energy81Tl,35.);
160 InitializeVectorEnergy(energy82Pb,40.);
161 InitializeVectorEnergy(energy83Bi,45.);
162 InitializeVectorEnergy(energy84Po,45.);
163 InitializeVectorEnergy(energy85At,45.);
164 InitializeVectorEnergy(energy86Rn,45.);
165 InitializeVectorEnergy(energy87Fr,45.);
166 InitializeVectorEnergy(energy88Ra,45.);
167 InitializeVectorEnergy(energy89Ac,50.);
168 InitializeVectorEnergy(energy90Th,50.);
169 InitializeVectorEnergy(energy91Pa,50.);
170 InitializeVectorEnergy(energy92U, 55.);
171}
172
173inline void G4hShellCrossSectionDoubleExpData::InitializeVectorPar1(std::vector<G4double> &vect, G4double value1, G4double value2, G4double value3) const
174{
175 vect.push_back(value1);
176 vect.push_back(value2);
177 vect.push_back(value3);
178}
179
180// First parameters
181
182void G4hShellCrossSectionDoubleExpData::FillVectorValuesPar1()
183{
184 InitializeVectorPar1(parlow6C, 6.028754 , 22.76324 , 8.217967 );
185 InitializeVectorPar1(parlow7N, 5.477011 , 18.81141 , 4.858727 );
186 InitializeVectorPar1(parlow8O, 5.425050 , 16.80088 , 3.910342 );
187 InitializeVectorPar1(parlow9F, 5.017089 , 14.56733 , 2.689939 );
188 InitializeVectorPar1(parlow10Ne,4.668802 , 12.75720 , 1.897390 );
189 InitializeVectorPar1(parlow11Na,4.705349 , 11.54904 , 1.700995 );
190 InitializeVectorPar1(parlow12Mg,4.263983 , 10.00729 , 1.117699 );
191 InitializeVectorPar1(parlow13Al,4.324701 , 9.009748 , 1.013421 );
192 InitializeVectorPar1(parlow14Si,4.090029 , 7.939616 , 0.7610449 );
193 InitializeVectorPar1(parlow15P, 4.098163 , 7.072802 , 0.6808610 );
194 InitializeVectorPar1(parlow16S, 3.632073 , 6.268209 , 0.4380426 );
195 InitializeVectorPar1(parlow17Cl,3.769954 , 5.545212 , 0.4345047 );
196 InitializeVectorPar1(parlow18Ar,3.803460 , 4.852357 , 0.4018530 );
197 InitializeVectorPar1(parlow19K, 3.818340 , 4.215003 , 0.3707851 );
198 InitializeVectorPar1(parlow20Ca,3.767600 , 3.621560 , 0.3257078 );
199 InitializeVectorPar1(parlow21Sc,3.851859 , 3.047405 , 0.3256020 );
200 InitializeVectorPar1(parlow22Ti,3.739032 , 2.540971 , 0.2726244 );
201 InitializeVectorPar1(parlow23V, 3.762692 , 2.028288 , 0.2595049 );
202 InitializeVectorPar1(parlow24Cr,3.830856 , 1.517501 , 0.2585256 );
203 InitializeVectorPar1(parlow25Mn,3.889344 , 1.032337 , 0.2568166 );
204 InitializeVectorPar1(parlow26Fe,4.610710686 , 0.275389472 , 2.764554447 );
205 InitializeVectorPar1(parlow27Co,3.55117935 , 0.228178229 , 2.702047259 );
206 InitializeVectorPar1(parlow28Ni,2.334937582 , 0.232156421 , 2.738149196 );
207 InitializeVectorPar1(parlow29Cu,1.915920665 , 0.182589984 , 2.664081065 );
208 InitializeVectorPar1(parlow30Zn,1.258082776 , 0.190338228 , 2.710680405 );
209 InitializeVectorPar1(parlow31Ga,1.090692834 , 0.140413681 , 2.624294005 );
210 InitializeVectorPar1(parlow32Ge,0.868211657 , 0.115779844 , 2.587823968 );
211 InitializeVectorPar1(parlow33As,0.588895163 , 0.121427557 , 2.630013394 );
212 InitializeVectorPar1(parlow34Se,0.464947416 , 0.103777355 , 2.60655728 );
213 InitializeVectorPar1(parlow35Br,0.372685678 , 0.087649176 , 2.580702817 );
214 InitializeVectorPar1(parlow36Kr,0.261618542 , 0.071657004 , 2.628062571 );
215 InitializeVectorPar1(parlow37Rb,0.203966062 , 0.063471976 , 2.621118231 );
216 InitializeVectorPar1(parlow38Sr,0.151309868 , 0.062650983 , 2.635670955 );
217 InitializeVectorPar1(parlow39Y, 0.130831069 , 0.046546373 , 2.596438847 );
218 InitializeVectorPar1(parlow40Zr,0.113465174 , 0.034961204 , 2.557081599 );
219 InitializeVectorPar1(parlow41Nb,0.089305327 , 0.031631918 , 2.560378036 );
220 InitializeVectorPar1(parlow42Mo,0.066940214 , 0.032387014 , 2.580578224 );
221 InitializeVectorPar1(parlow43Tc,0.054005087 , 0.028625811 , 2.578181335 );
222 InitializeVectorPar1(parlow44Ru,0.047435043 , 0.021858688 , 2.544955569 );
223 InitializeVectorPar1(parlow45Rh,0.037036289 , 0.021226071 , 2.557147432 );
224 InitializeVectorPar1(parlow46Pd,0.031108855 , 0.017082016 , 2.560638611 );
225 InitializeVectorPar1(parlow47Ag,0.030039299 , 0.010687233 , 2.501552766 );
226 InitializeVectorPar1(parlow48Cd,0.025896549 , 0.008577291 , 2.485085985 );
227 InitializeVectorPar1(parlow49In,0.021622343 , 0.007503402 , 2.479363645 );
228 InitializeVectorPar1(parlow50Sn,0.019658368 , 0.005372974 , 2.448647025 );
229 InitializeVectorPar1(parlow51Sb,0.016726542 , 0.004540107 , 2.440841502 );
230 InitializeVectorPar1(parlow52Te,0.01210377 , 0.005793226 , 2.487112403 );
231 InitializeVectorPar1(parlow53I, 0.012068341 , 0.003340949 , 2.430469903 );
232 InitializeVectorPar1(parlow54Xe,0.011309929 , 0.002247705 , 2.39653105 );
233 InitializeVectorPar1(parlow55Cs,0.010322464 , 0.001584264 , 2.372751551 );
234 InitializeVectorPar1(parlow56Ba,0.00891645 , 0.001663984 , 2.365104652 );
235 InitializeVectorPar1(parlow57La,0.008600326 , 0.000973706 , 2.330302922 );
236 InitializeVectorPar1(parlow58Ce,0.00852455 , 0.000554424 , 2.286062759 );
237 InitializeVectorPar1(parlow59Pr,0.006464573 , 0.000751099 , 2.320217715 );
238 InitializeVectorPar1(parlow60Nd,0.005772215 , 0.000604636 , 2.308998547 );
239 InitializeVectorPar1(parlow61Pm,0.005439364 , 0.000395253 , 2.284991617 );
240 InitializeVectorPar1(parlow62Sm,0.004140155 , 0.000572248 , 2.319931149 );
241 InitializeVectorPar1(parlow63Eu,0.004391028 , 0.00026105 , 2.263027478 );
242 InitializeVectorPar1(parlow64Gd,0.004108655 , 0.00018564 , 2.242594759 );
243 InitializeVectorPar1(parlow65Tb,0.003335626 , 0.000214331 , 2.262351386 );
244 InitializeVectorPar1(parlow66Dy,0.003772419844 , 1.491354113E-4, 2.182513765 );
245 InitializeVectorPar1(parlow67Ho,0.003171261268 , 1.610171391E-4, 2.19445375 );
246 InitializeVectorPar1(parlow68Er,0.003511231198 , 6.431951513E-5, 2.131309482 );
247 InitializeVectorPar1(parlow69Tm,0.002951671258 , 7.185528546E-5, 2.14407165 );
248 InitializeVectorPar1(parlow70Yb,0.003090003974 , 3.232352598E-5, 2.10040042 );
249 InitializeVectorPar1(parlow71Lu,0.002504560425 , 4.55286561E-5 , 2.122593025 );
250 InitializeVectorPar1(parlow72Hf,0.002634273859 , 2.025018154E-5, 2.078762159 );
251 InitializeVectorPar1(parlow73Ta,0.002246920523 , 2.530957413E-5, 2.087036821 );
252 InitializeVectorPar1(parlow74W, 0.002138140746 , 1.685058881E-5, 2.071640226 );
253 InitializeVectorPar1(parlow75Re,0.002053720518 , 1.1699924E-5 , 2.052411395 );
254 InitializeVectorPar1(parlow76Os,0.001609160281 , 1.98304501E-5 , 2.07932865 );
255 InitializeVectorPar1(parlow77Ir,0.001451780286 , 1.89647459E-5 , 2.077028977 );
256 InitializeVectorPar1(parlow78Pt,0.1484928E-2 , 0.1015788E-4 , 2.043982 );
257 InitializeVectorPar1(parlow79Au,0.1401423E-2 , 0.7618721E-5 , 2.031677 );
258 InitializeVectorPar1(parlow80Hg,0.1279221E-2 , 0.7684558E-5 , 2.026922 );
259 InitializeVectorPar1(parlow81Tl,0.001116529691 , 8.665290032E-6, 2.035250841 );
260 InitializeVectorPar1(parlow82Pb,0.00123546161 , 3.264194012E-6, 1.984513455 );
261 InitializeVectorPar1(parlow83Bi,0.1168456E-2 , 0.2565728E-5 , 1.974040 );
262 InitializeVectorPar1(parlow84Po,0.001161495281 , 1.517549601E-6, 1.951280331 );
263 InitializeVectorPar1(parlow85At,0.9574954E-3 , 0.2619334E-5 , 1.974164 );
264 InitializeVectorPar1(parlow86Rn,0.9172934E-3 , 0.9124897E-6 , 1.966786 );
265 InitializeVectorPar1(parlow87Fr,0.7900983E-3 , 0.1254855E-5 , 1.978866 );
266 InitializeVectorPar1(parlow88Ra,0.7093558E-3 , 0.1436426E-5 , 1.981631 );
267 InitializeVectorPar1(parlow89Ac,0.0007906194724, 4.455647522E-7, 1.934485393 );
268 InitializeVectorPar1(parlow90Th,0.0007280666063, 4.668165706E-7, 1.931837161 );
269 InitializeVectorPar1(parlow91Pa,0.0006494787919, 5.572898409E-7, 1.937018718 );
270 InitializeVectorPar1(parlow92U, 0.0006712361299, 2.734932326E-7, 1.908976604 );
271}
272
273inline void G4hShellCrossSectionDoubleExpData::InitializeVectorPar2(std::vector<G4double> &vect, G4double value1, G4double value2, G4double value3, G4double value4, G4double value5) const
274{
275 vect.push_back(value1);
276 vect.push_back(value2);
277 vect.push_back(value3);
278 vect.push_back(value4);
279 vect.push_back(value5);
280}
281
282// Second parameters
283
284void G4hShellCrossSectionDoubleExpData::FillVectorValuesPar2()
285{
286 InitializeVectorPar2(parhigh6C, 0.2614505E+08, 0.1960089 , -0.9781201, 0.0, 0.0);
287 InitializeVectorPar2(parhigh7N, 0.1469968E+08, 0.9198590E-01, -0.9035943, 0.0, 0.0);
288 InitializeVectorPar2(parhigh8O, 9015526.0 , 0.4561001E-01, -0.8435532, 0.0, 0.0);
289 InitializeVectorPar2(parhigh9F, 5878802.0 , 0.2042609E-01, -0.8003597, 0.0, 0.0);
290 InitializeVectorPar2(parhigh10Ne,3951196.0 , 0.8858837E-02, -0.7637176, 0.0, 0.0);
291 InitializeVectorPar2(parhigh11Na,2838478.0 , 0.3098775E-02, -0.7451933, 0.0, 0.0);
292 InitializeVectorPar2(parhigh12Mg,2322929.0 , 0.6889491E-03, -0.7595924, 0.0, 0.0);
293 InitializeVectorPar2(parhigh13Al,1739866.0 , 0.1930171E-03, -0.7445110, 0.0, 0.0);
294 InitializeVectorPar2(parhigh14Si,1433611.0 , 0.3515504E-04, -0.7487031, 0.0, 0.0);
295 InitializeVectorPar2(parhigh15P, 1116443.0 , 0.7842201E-05, -0.7376342, 0.0, 0.0);
296 InitializeVectorPar2(parhigh16S, 1032940.0 , 0.7188930E-06, -0.7601197, 0.0, 0.0);
297 InitializeVectorPar2(parhigh17Cl,804101.7 , 0.1521286E-06, -0.7416714, 0.0, 0.0);
298 InitializeVectorPar2(parhigh18Ar,637702.8 , 0.2779553E-07, -0.7268317, 0.0, 0.0);
299 InitializeVectorPar2(parhigh19K, 519937.5 , 0.4463855E-08, -0.7158139, 0.0, 0.0);
300 InitializeVectorPar2(parhigh20Ca,436219.8 , 0.5878926E-09, -0.7102222, 0.0, 0.0);
301 InitializeVectorPar2(parhigh21Sc,368703.6 , 0.7282178E-10, -0.7046522, 0.0, 0.0);
302 InitializeVectorPar2(parhigh22Ti,397138.5 , 0.1330592E-11, -0.7196041, 0.0, 0.0);
303 InitializeVectorPar2(parhigh23V, 274427.9 , 0.6204048E-12, -0.6963667, 0.0, 0.0);
304 InitializeVectorPar2(parhigh24Cr,228919.3 , 0.7995359E-13, -0.6851181, 0.0, 0.0);
305 InitializeVectorPar2(parhigh25Mn,199919.3 , 0.4059536E-14, -0.6801181, 0.0, 0.0);
306 InitializeVectorPar2(parhigh26Fe,5888.956103 , -9427.905394 , 4519.130051 , -537.9587092, 0.0);
307 InitializeVectorPar2(parhigh27Co,5320.748659 , -8223.484184 , 3805.325189 , -441.4870379, 0.0);
308 InitializeVectorPar2(parhigh28Ni,4781.319646 , -7174.783041 , 3219.076628 , -364.8043497, 0.0);
309 InitializeVectorPar2(parhigh29Cu,5036.391497 , -7072.396111 , 2999.984178 , -330.4071107, 0.0);
310 InitializeVectorPar2(parhigh30Zn,4379.092796 , -6046.707901 , 2513.270363 , -271.5536118, 0.0);
311 InitializeVectorPar2(parhigh31Ga,4517.118137 , -5908.263246 , 2344.263972 , -247.0690309, 0.0);
312 InitializeVectorPar2(parhigh32Ge,3842.848695 , -4981.810605 , 1950.19123 , -202.2213944, 0.0);
313 InitializeVectorPar2(parhigh33As,3936.923647 , -4860.36463 , 1825.209324 , -184.9650549, 0.0);
314 InitializeVectorPar2(parhigh34Se,3405.652981 , -4159.784027 , 1540.917232 , -153.778945 , 0.0);
315 InitializeVectorPar2(parhigh35Br,3482.563108 , -4065.531877 , 1450.303751 , -141.6935145, 0.0);
316 InitializeVectorPar2(parhigh36Kr,27792.4379 , 2.49277E-26 , -0.500883086, 0.0, 0.0);
317 InitializeVectorPar2(parhigh37Rb,15982.48534 , 3.73258E-25 , -0.416741096, 0.0, 0.0);
318 InitializeVectorPar2(parhigh38Sr,17889.43232 , 1.13451E-27 , -0.454847146, 0.0, 0.0);
319 InitializeVectorPar2(parhigh39Y, 12451.29445 , 1.53532E-27 , -0.404893344, 0.0, 0.0);
320 InitializeVectorPar2(parhigh40Zr,10206.19067 , 1.77653E-28 , -0.383987982, 0.0, 0.0);
321 InitializeVectorPar2(parhigh41Nb,11017.03479 , 8.86594E-31 , -0.414463124, 0.0, 0.0);
322 InitializeVectorPar2(parhigh42Mo,8535.867912 , 3.16226E-31 , -0.383744046, 0.0, 0.0);
323 InitializeVectorPar2(parhigh43Tc,9820.134131 , 3.36129E-34 , -0.423630175, 0.0, 0.0);
324 InitializeVectorPar2(parhigh44Ru,7669.947859 , 1.20762E-34 , -0.39423283 , 0.0, 0.0);
325 InitializeVectorPar2(parhigh45Rh,8187.067923 , 3.36422E-37 , -0.420048481, 0.0, 0.0);
326 InitializeVectorPar2(parhigh46Pd,6449.690254 , 1.07461E-37 , -0.388615201, 0.0, 0.0);
327 InitializeVectorPar2(parhigh47Ag,7131.825038 , 1.13958E-40 , -0.419291952, 0.0, 0.0);
328 InitializeVectorPar2(parhigh48Cd,5622.898379 , 3.91721E-41 , -0.39138804 , 0.0, 0.0);
329 InitializeVectorPar2(parhigh49In,6154.018663 , 3.57242E-44 , -0.419285141, 0.0, 0.0);
330 InitializeVectorPar2(parhigh50Sn,5308.904577 , 1.83756E-45 , -0.405840662, 0.0, 0.0);
331 InitializeVectorPar2(parhigh51Sb,4205.562408 , 7.25156E-46 , -0.378710249, 0.0, 0.0);
332 InitializeVectorPar2(parhigh52Te,4590.865265 , 4.9712E-49 , -0.405349183, 0.0, 0.0);
333 InitializeVectorPar2(parhigh53I, 3930.781204 , 4.05142E-50 , -0.390685876, 0.0, 0.0);
334 InitializeVectorPar2(parhigh54Xe,3130.401154 , 1.74282E-50 , -0.364347309, 0.0, 0.0);
335 InitializeVectorPar2(parhigh55Cs,3744.705387 , 5.90267E-55 , -0.404336908, 0.0, 0.0);
336 InitializeVectorPar2(parhigh56Ba,677.6831684 , 1.8374E-42 , -0.129113534, 0.0, 0.0);
337 InitializeVectorPar2(parhigh57La,717.0531832 , 9.51565E-45 , -0.15227809 , 0.0, 0.0);
338 InitializeVectorPar2(parhigh58Ce,641.4062303 , 3.69483E-46 , -0.144632417, 0.0, 0.0);
339 InitializeVectorPar2(parhigh59Pr,507.1392883 , 3.50854E-46 , -0.117584038, 0.0, 0.0);
340 InitializeVectorPar2(parhigh60Nd,590.0643943 , 1.09726E-49 , -0.15428397 , 0.0, 0.0);
341 InitializeVectorPar2(parhigh61Pm,487.3334007 , 2.75498E-50 , -0.133698368, 0.0, 0.0);
342 InitializeVectorPar2(parhigh62Sm,358.3518729 , 2.70531E-49 , -0.095301068, 0.0, 0.0);
343 InitializeVectorPar2(parhigh63Eu,329.3455874 , 5.55138E-51 , -0.091729198, 0.0, 0.0);
344 InitializeVectorPar2(parhigh64Gd,413.9439156 , 1.23792E-55 , -0.139710531, 0.0, 0.0);
345 InitializeVectorPar2(parhigh65Tb,336.3177626 , 1.01038E-55 , -0.116827533, 0.0, 0.0);
346 InitializeVectorPar2(parhigh66Dy,-1424.412332 , 1495.300993 , -579.4404515, 97.27901953, -5.836831374 );
347 InitializeVectorPar2(parhigh67Ho,-1517.735661 , 1565.770524 , -596.9215188, 98.77020094, -5.861740123 );
348 InitializeVectorPar2(parhigh68Er,-1474.139215 , 1509.677298 , -571.276182 , 93.84192307, -5.533405798 );
349 InitializeVectorPar2(parhigh69Tm,-1468.335 , 1486.750 , -556.4769 , 90.48160 , -5.288727 );
350 InitializeVectorPar2(parhigh70Yb,-1224.051265 , 1244.798125 , -467.4078327, 76.09587249, -4.438193762 );
351 InitializeVectorPar2(parhigh71Lu,-1297.779 , 1295.939 , -478.6122 , 76.80235 , -4.430208 );
352 InitializeVectorPar2(parhigh72Hf,-1343.224 , 1320.527 , -480.7212 , 76.16011 , -4.348449 );
353 InitializeVectorPar2(parhigh73Ta,-1121.596 , 1113.986 , -408.6746 , 65.06258 , -3.719345 );
354 InitializeVectorPar2(parhigh74W, -1132.138 , 1109.449 , -402.0285 , 63.30364 , -3.586027 );
355 InitializeVectorPar2(parhigh75Re,-1099.159467 , 1069.010051 , -384.5257059, 60.12188988, -3.384193421 );
356 InitializeVectorPar2(parhigh76Os,-1190.387975 , 1140.34171 , -404.1551887, 62.32091335, -3.468141645 );
357 InitializeVectorPar2(parhigh77Ir,-1094.70162 , 1044.989734 , -368.8959923, 56.62618718, -3.134201442 );
358 InitializeVectorPar2(parhigh78Pt,-1163.104391 , 1092.439757 , -379.9958732, 57.58306885, -3.155663589 );
359 InitializeVectorPar2(parhigh79Au,-1144.500 , 1068.235 , -369.2752 , 55.62967 , -3.033460 );
360 InitializeVectorPar2(parhigh80Hg,-1003.592 , 937.0178 , -323.7748 , 48.69880 , -2.646195 );
361 InitializeVectorPar2(parhigh81Tl,-907.5804 , 850.2096 , -294.3571 , 44.29993 , -2.405155 );
362 InitializeVectorPar2(parhigh82Pb,-1015.548 , 929.8707 , -315.4740 , 46.66766 , -2.501258 );
363 InitializeVectorPar2(parhigh83Bi,-1076.950 , 971.7385 , -325.2844 , 47.55627 , -2.525684 );
364 InitializeVectorPar2(parhigh84Po,-1021.285 , 917.5267 , -305.7521 , 44.49191 , -2.351658 );
365 InitializeVectorPar2(parhigh85At,-967.7470 , 866.0588 , -287.4283 , 41.64792 , -2.191694 );
366 InitializeVectorPar2(parhigh86Rn,-338.6518 , 311.5210 , -105.8054 , 15.52788 , -0.8066356 );
367 InitializeVectorPar2(parhigh87Fr,-295.5938 , 271.6075 , -92.01641 , 13.43905 , -0.6914321 );
368 InitializeVectorPar2(parhigh88Ra,-348.9885768 , 317.6478799 , -106.8672925, 15.57672306, -0.8090625926);
369 InitializeVectorPar2(parhigh89Ac,-141.772905 , 134.4449818 , -46.38678324, 6.754305544, -0.3309070686);
370 InitializeVectorPar2(parhigh90Th,-265.8311235 , 240.8739384 , -80.49877194, 11.60745927, -0.5910321851);
371 InitializeVectorPar2(parhigh91Pa,-170.928272 , 157.4297802 , -53.12415098, 7.648148663, -0.379425533 );
372 InitializeVectorPar2(parhigh92U, -344.2548753 , 304.5108607 , -99.78215019, 14.21641181, -0.7269520491);
373}
374
375void G4hShellCrossSectionDoubleExpData::FillParameterMapEnergy()
376{
377 parameterMapEnergy [6] = &energy6C ;
378 parameterMapEnergy [7] = &energy7N ;
379 parameterMapEnergy [8] = &energy8O ;
380 parameterMapEnergy [9] = &energy9F ;
381 parameterMapEnergy [10] = &energy10Ne;
382 parameterMapEnergy [11] = &energy11Na;
383 parameterMapEnergy [12] = &energy12Mg;
384 parameterMapEnergy [13] = &energy13Al;
385 parameterMapEnergy [14] = &energy14Si;
386 parameterMapEnergy [15] = &energy15P ;
387 parameterMapEnergy [16] = &energy16S ;
388 parameterMapEnergy [17] = &energy17Cl;
389 parameterMapEnergy [18] = &energy18Ar;
390 parameterMapEnergy [19] = &energy19K ;
391 parameterMapEnergy [20] = &energy20Ca;
392 parameterMapEnergy [21] = &energy21Sc;
393 parameterMapEnergy [22] = &energy22Ti;
394 parameterMapEnergy [23] = &energy23V ;
395 parameterMapEnergy [24] = &energy24Cr;
396 parameterMapEnergy [25] = &energy25Mn;
397 parameterMapEnergy [26] = &energy26Fe;
398 parameterMapEnergy [27] = &energy27Co;
399 parameterMapEnergy [28] = &energy28Ni;
400 parameterMapEnergy [29] = &energy29Cu;
401 parameterMapEnergy [30] = &energy30Zn;
402 parameterMapEnergy [31] = &energy31Ga;
403 parameterMapEnergy [32] = &energy32Ge;
404 parameterMapEnergy [33] = &energy33As;
405 parameterMapEnergy [34] = &energy34Se;
406 parameterMapEnergy [35] = &energy35Br;
407 parameterMapEnergy [36] = &energy36Kr;
408 parameterMapEnergy [37] = &energy37Rb;
409 parameterMapEnergy [38] = &energy38Sr;
410 parameterMapEnergy [39] = &energy39Y ;
411 parameterMapEnergy [40] = &energy40Zr;
412 parameterMapEnergy [41] = &energy41Nb;
413 parameterMapEnergy [42] = &energy42Mo;
414 parameterMapEnergy [43] = &energy43Tc;
415 parameterMapEnergy [44] = &energy44Ru;
416 parameterMapEnergy [45] = &energy45Rh;
417 parameterMapEnergy [46] = &energy46Pd;
418 parameterMapEnergy [47] = &energy47Ag;
419 parameterMapEnergy [48] = &energy48Cd;
420 parameterMapEnergy [49] = &energy49In;
421 parameterMapEnergy [50] = &energy50Sn;
422 parameterMapEnergy [51] = &energy51Sb;
423 parameterMapEnergy [52] = &energy52Te;
424 parameterMapEnergy [53] = &energy53I ;
425 parameterMapEnergy [54] = &energy54Xe;
426 parameterMapEnergy [55] = &energy55Cs;
427 parameterMapEnergy [56] = &energy56Ba;
428 parameterMapEnergy [57] = &energy57La;
429 parameterMapEnergy [58] = &energy58Ce;
430 parameterMapEnergy [59] = &energy59Pr;
431 parameterMapEnergy [60] = &energy60Nd;
432 parameterMapEnergy [61] = &energy61Pm;
433 parameterMapEnergy [62] = &energy62Sm;
434 parameterMapEnergy [63] = &energy63Eu;
435 parameterMapEnergy [64] = &energy64Gd;
436 parameterMapEnergy [65] = &energy65Tb;
437 parameterMapEnergy [66] = &energy66Dy;
438 parameterMapEnergy [67] = &energy67Ho;
439 parameterMapEnergy [68] = &energy68Er;
440 parameterMapEnergy [69] = &energy69Tm;
441 parameterMapEnergy [70] = &energy70Yb;
442 parameterMapEnergy [71] = &energy71Lu;
443 parameterMapEnergy [72] = &energy72Hf;
444 parameterMapEnergy [73] = &energy73Ta;
445 parameterMapEnergy [74] = &energy74W ;
446 parameterMapEnergy [75] = &energy75Re;
447 parameterMapEnergy [76] = &energy76Os;
448 parameterMapEnergy [77] = &energy77Ir;
449 parameterMapEnergy [78] = &energy78Pt;
450 parameterMapEnergy [79] = &energy79Au;
451 parameterMapEnergy [80] = &energy80Hg;
452 parameterMapEnergy [81] = &energy81Tl;
453 parameterMapEnergy [82] = &energy82Pb;
454 parameterMapEnergy [83] = &energy83Bi;
455 parameterMapEnergy [84] = &energy84Po;
456 parameterMapEnergy [85] = &energy85At;
457 parameterMapEnergy [86] = &energy86Rn;
458 parameterMapEnergy [87] = &energy87Fr;
459 parameterMapEnergy [88] = &energy88Ra;
460 parameterMapEnergy [89] = &energy89Ac;
461 parameterMapEnergy [90] = &energy90Th;
462 parameterMapEnergy [91] = &energy91Pa;
463 parameterMapEnergy [92] = &energy92U ;
464}
465
466void G4hShellCrossSectionDoubleExpData::FillParameterMapPar1()
467{
468 parameterMapPar1 [6] = &parlow6C ;
469 parameterMapPar1 [7] = &parlow7N ;
470 parameterMapPar1 [8] = &parlow8O ;
471 parameterMapPar1 [9] = &parlow9F ;
472 parameterMapPar1 [10] = &parlow10Ne;
473 parameterMapPar1 [11] = &parlow11Na;
474 parameterMapPar1 [12] = &parlow12Mg;
475 parameterMapPar1 [13] = &parlow13Al;
476 parameterMapPar1 [14] = &parlow14Si;
477 parameterMapPar1 [15] = &parlow15P ;
478 parameterMapPar1 [16] = &parlow16S ;
479 parameterMapPar1 [17] = &parlow17Cl;
480 parameterMapPar1 [18] = &parlow18Ar;
481 parameterMapPar1 [19] = &parlow19K ;
482 parameterMapPar1 [20] = &parlow20Ca;
483 parameterMapPar1 [21] = &parlow21Sc;
484 parameterMapPar1 [22] = &parlow22Ti;
485 parameterMapPar1 [23] = &parlow23V ;
486 parameterMapPar1 [24] = &parlow24Cr;
487 parameterMapPar1 [25] = &parlow25Mn;
488 parameterMapPar1 [26] = &parlow26Fe;
489 parameterMapPar1 [27] = &parlow27Co;
490 parameterMapPar1 [28] = &parlow28Ni;
491 parameterMapPar1 [29] = &parlow29Cu;
492 parameterMapPar1 [30] = &parlow30Zn;
493 parameterMapPar1 [31] = &parlow31Ga;
494 parameterMapPar1 [32] = &parlow32Ge;
495 parameterMapPar1 [33] = &parlow33As;
496 parameterMapPar1 [34] = &parlow34Se;
497 parameterMapPar1 [35] = &parlow35Br;
498 parameterMapPar1 [36] = &parlow36Kr;
499 parameterMapPar1 [37] = &parlow37Rb;
500 parameterMapPar1 [38] = &parlow38Sr;
501 parameterMapPar1 [39] = &parlow39Y ;
502 parameterMapPar1 [40] = &parlow40Zr;
503 parameterMapPar1 [41] = &parlow41Nb;
504 parameterMapPar1 [42] = &parlow42Mo;
505 parameterMapPar1 [43] = &parlow43Tc;
506 parameterMapPar1 [44] = &parlow44Ru;
507 parameterMapPar1 [45] = &parlow45Rh;
508 parameterMapPar1 [46] = &parlow46Pd;
509 parameterMapPar1 [47] = &parlow47Ag;
510 parameterMapPar1 [48] = &parlow48Cd;
511 parameterMapPar1 [49] = &parlow49In;
512 parameterMapPar1 [50] = &parlow50Sn;
513 parameterMapPar1 [51] = &parlow51Sb;
514 parameterMapPar1 [52] = &parlow52Te;
515 parameterMapPar1 [53] = &parlow53I ;
516 parameterMapPar1 [54] = &parlow54Xe;
517 parameterMapPar1 [55] = &parlow55Cs;
518 parameterMapPar1 [56] = &parlow56Ba;
519 parameterMapPar1 [57] = &parlow57La;
520 parameterMapPar1 [58] = &parlow58Ce;
521 parameterMapPar1 [59] = &parlow59Pr;
522 parameterMapPar1 [60] = &parlow60Nd;
523 parameterMapPar1 [61] = &parlow61Pm;
524 parameterMapPar1 [62] = &parlow62Sm;
525 parameterMapPar1 [63] = &parlow63Eu;
526 parameterMapPar1 [64] = &parlow64Gd;
527 parameterMapPar1 [65] = &parlow65Tb;
528 parameterMapPar1 [66] = &parlow66Dy;
529 parameterMapPar1 [67] = &parlow67Ho;
530 parameterMapPar1 [68] = &parlow68Er;
531 parameterMapPar1 [69] = &parlow69Tm;
532 parameterMapPar1 [70] = &parlow70Yb;
533 parameterMapPar1 [71] = &parlow71Lu;
534 parameterMapPar1 [72] = &parlow72Hf;
535 parameterMapPar1 [73] = &parlow73Ta;
536 parameterMapPar1 [74] = &parlow74W ;
537 parameterMapPar1 [75] = &parlow75Re;
538 parameterMapPar1 [76] = &parlow76Os;
539 parameterMapPar1 [77] = &parlow77Ir;
540 parameterMapPar1 [78] = &parlow78Pt;
541 parameterMapPar1 [79] = &parlow79Au;
542 parameterMapPar1 [80] = &parlow80Hg;
543 parameterMapPar1 [81] = &parlow81Tl;
544 parameterMapPar1 [82] = &parlow82Pb;
545 parameterMapPar1 [83] = &parlow83Bi;
546 parameterMapPar1 [84] = &parlow84Po;
547 parameterMapPar1 [85] = &parlow85At;
548 parameterMapPar1 [86] = &parlow86Rn;
549 parameterMapPar1 [87] = &parlow87Fr;
550 parameterMapPar1 [88] = &parlow88Ra;
551 parameterMapPar1 [89] = &parlow89Ac;
552 parameterMapPar1 [90] = &parlow90Th;
553 parameterMapPar1 [91] = &parlow91Pa;
554 parameterMapPar1 [92] = &parlow92U ;
555}
556
557void G4hShellCrossSectionDoubleExpData::FillParameterMapPar2()
558{
559 parameterMapPar2 [6] = &parhigh6C ;
560 parameterMapPar2 [7] = &parhigh7N ;
561 parameterMapPar2 [8] = &parhigh8O ;
562 parameterMapPar2 [9] = &parhigh9F ;
563 parameterMapPar2 [10] = &parhigh10Ne;
564 parameterMapPar2 [11] = &parhigh11Na;
565 parameterMapPar2 [12] = &parhigh12Mg;
566 parameterMapPar2 [13] = &parhigh13Al;
567 parameterMapPar2 [14] = &parhigh14Si;
568 parameterMapPar2 [15] = &parhigh15P ;
569 parameterMapPar2 [16] = &parhigh16S ;
570 parameterMapPar2 [17] = &parhigh17Cl;
571 parameterMapPar2 [18] = &parhigh18Ar;
572 parameterMapPar2 [19] = &parhigh19K ;
573 parameterMapPar2 [20] = &parhigh20Ca;
574 parameterMapPar2 [21] = &parhigh21Sc;
575 parameterMapPar2 [22] = &parhigh22Ti;
576 parameterMapPar2 [23] = &parhigh23V ;
577 parameterMapPar2 [24] = &parhigh24Cr;
578 parameterMapPar2 [25] = &parhigh25Mn;
579 parameterMapPar2 [26] = &parhigh26Fe;
580 parameterMapPar2 [27] = &parhigh27Co;
581 parameterMapPar2 [28] = &parhigh28Ni;
582 parameterMapPar2 [29] = &parhigh29Cu;
583 parameterMapPar2 [30] = &parhigh30Zn;
584 parameterMapPar2 [31] = &parhigh31Ga;
585 parameterMapPar2 [32] = &parhigh32Ge;
586 parameterMapPar2 [33] = &parhigh33As;
587 parameterMapPar2 [34] = &parhigh34Se;
588 parameterMapPar2 [35] = &parhigh35Br;
589 parameterMapPar2 [36] = &parhigh36Kr;
590 parameterMapPar2 [37] = &parhigh37Rb;
591 parameterMapPar2 [38] = &parhigh38Sr;
592 parameterMapPar2 [39] = &parhigh39Y ;
593 parameterMapPar2 [40] = &parhigh40Zr;
594 parameterMapPar2 [41] = &parhigh41Nb;
595 parameterMapPar2 [42] = &parhigh42Mo;
596 parameterMapPar2 [43] = &parhigh43Tc;
597 parameterMapPar2 [44] = &parhigh44Ru;
598 parameterMapPar2 [45] = &parhigh45Rh;
599 parameterMapPar2 [46] = &parhigh46Pd;
600 parameterMapPar2 [47] = &parhigh47Ag;
601 parameterMapPar2 [48] = &parhigh48Cd;
602 parameterMapPar2 [49] = &parhigh49In;
603 parameterMapPar2 [50] = &parhigh50Sn;
604 parameterMapPar2 [51] = &parhigh51Sb;
605 parameterMapPar2 [52] = &parhigh52Te;
606 parameterMapPar2 [53] = &parhigh53I ;
607 parameterMapPar2 [54] = &parhigh54Xe;
608 parameterMapPar2 [55] = &parhigh55Cs;
609 parameterMapPar2 [56] = &parhigh56Ba;
610 parameterMapPar2 [57] = &parhigh57La;
611 parameterMapPar2 [58] = &parhigh58Ce;
612 parameterMapPar2 [59] = &parhigh59Pr;
613 parameterMapPar2 [60] = &parhigh60Nd;
614 parameterMapPar2 [61] = &parhigh61Pm;
615 parameterMapPar2 [62] = &parhigh62Sm;
616 parameterMapPar2 [63] = &parhigh63Eu;
617 parameterMapPar2 [64] = &parhigh64Gd;
618 parameterMapPar2 [65] = &parhigh65Tb;
619 parameterMapPar2 [66] = &parhigh66Dy;
620 parameterMapPar2 [67] = &parhigh67Ho;
621 parameterMapPar2 [68] = &parhigh68Er;
622 parameterMapPar2 [69] = &parhigh69Tm;
623 parameterMapPar2 [70] = &parhigh70Yb;
624 parameterMapPar2 [71] = &parhigh71Lu;
625 parameterMapPar2 [72] = &parhigh72Hf;
626 parameterMapPar2 [73] = &parhigh73Ta;
627 parameterMapPar2 [74] = &parhigh74W ;
628 parameterMapPar2 [75] = &parhigh75Re;
629 parameterMapPar2 [76] = &parhigh76Os;
630 parameterMapPar2 [77] = &parhigh77Ir;
631 parameterMapPar2 [78] = &parhigh78Pt;
632 parameterMapPar2 [79] = &parhigh79Au;
633 parameterMapPar2 [80] = &parhigh80Hg;
634 parameterMapPar2 [81] = &parhigh81Tl;
635 parameterMapPar2 [82] = &parhigh82Pb;
636 parameterMapPar2 [83] = &parhigh83Bi;
637 parameterMapPar2 [84] = &parhigh84Po;
638 parameterMapPar2 [85] = &parhigh85At;
639 parameterMapPar2 [86] = &parhigh86Rn;
640 parameterMapPar2 [87] = &parhigh87Fr;
641 parameterMapPar2 [88] = &parhigh88Ra;
642 parameterMapPar2 [89] = &parhigh89Ac;
643 parameterMapPar2 [90] = &parhigh90Th;
644 parameterMapPar2 [91] = &parhigh91Pa;
645 parameterMapPar2 [92] = &parhigh92U ;
646}
647
648std::vector<std::vector<G4double>*> G4hShellCrossSectionDoubleExpData::GetParam(G4int Z)
649{
650
651 std::vector<std::vector<G4double>*> energyPar1Par2;
652
653 // if (energyPar1Par2.size() != 0) {
654 // for (size_t i=1; i==energyPar1Par2.size(); i++)
655 // {
656 // energyPar1Par2.erase(energyPar1Par2.begin());
657 // }
658 // }
659
660 energyPar1Par2.push_back(parameterMapEnergy[Z]);
661 energyPar1Par2.push_back(parameterMapPar1[Z]);
662 energyPar1Par2.push_back(parameterMapPar2[Z]);
663 return energyPar1Par2;
664}
Note: See TracBrowser for help on using the repository browser.