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

Last change on this file since 1347 was 1347, checked in by garnier, 13 years ago

geant4 tag 9.4
File size: 34.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: G4hShellCrossSectionDoubleExpData.cc,v 1.5 2009/06/10 13:32:36 mantero Exp $
36// GEANT4 tag $Name: geant4-09-04-ref-00 $
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.