source: trunk/source/materials/src/G4DensityEffectData.cc @ 1315

Last change on this file since 1315 was 1315, checked in by garnier, 14 years ago

update geant4-09-04-beta-cand-01 interfaces-V09-03-09 vis-V09-03-08

File size: 47.7 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// $Id: G4DensityEffectData.cc,v 1.12 2010/05/15 15:37:33 vnivanch Exp $
28// GEANT4 tag $Name: geant4-09-04-beta-cand-01 $
29
30//---------------------------------------------------------------------------
31//
32// GEANT4 Class file
33//
34// Description: Data on density effect
35//
36// Authors:    A.Bagulya, A.Ivanchenko 28.10.2009
37//
38// Modifications:
39//
40//
41//
42//----------------------------------------------------------------------------
43//
44
45//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
46
47#include "G4DensityEffectData.hh"
48#include <iomanip>
49
50//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
51
52G4DensityEffectData::G4DensityEffectData():index(0) 
53{
54  Initialize();
55}
56
57//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
58
59G4DensityEffectData::~G4DensityEffectData()
60{
61}
62
63void G4DensityEffectData::Initialize()
64{
65  // Data[9]:  Eplasma rho -C X_0 X_1 a m delta_0 DELTA_{max}       
66  // Eplasma - Plasma energy (in eV)
67  // rho     - Sternheimer adjustment factor for the atomic excitation energies
68  // -C
69  // X_0, X_1, m, a - parameters in fitting formulas
70  // delta_0 - Density-effect value delta(X_0)
71  // DELTA_{max} - Upper bound for the error inherent in the fitting procedure
72 
73  for (G4int i=0; i< NDENSELEM; ++i)
74    {
75      indexZ[i] = -1;
76      state[i]  = kStateSolid;
77    }
78
79//G4_lH2  index=0
80  G4double M1[NDENSARRAY]={7.031,1.546,3.2632,0.4759,1.9215,0.13483,5.6249,0,0.021};
81  AddMaterial(M1,"G4_lH2");
82
83//G4_H  index=1
84  G4double M0[NDENSARRAY]={0.263,1.412,9.5835,1.8639,3.2718,0.14092,5.7273,0.0,0.024};
85  AddMaterial(M0,"G4_H");
86  indexZ[1]=1;
87  state[1]=kStateGas;
88
89//G4_He  index=2
90  G4double M2[NDENSARRAY]={0.263,1.7,11.1393,2.2017,3.6122,0.13443,5.8347,0,0.024};
91  AddMaterial(M2,"G4_He");
92  indexZ[2]=2;
93  state[2]=kStateGas;
94
95//G4_Li  index=3
96  G4double M3[NDENSARRAY]={13.844,1.535,3.1221,0.1304,1.6397,0.95136,2.4993,0.14,0.062};
97  AddMaterial(M3,"G4_Li");
98  indexZ[3]=3;
99
100//G4_Be  index=4
101  G4double M4[NDENSARRAY]={26.096,1.908,2.7847,0.0392,1.6922,0.80392,2.4339,0.14,0.029};
102  AddMaterial(M4,"G4_Be");
103  indexZ[4]=4;
104
105//G4_B  index=5
106  G4double M5[NDENSARRAY]={30.17,2.32,2.8477,0.0305,1.9688,0.56224,2.4512,0.14,0.024};
107  AddMaterial(M5,"G4_B");
108  indexZ[5]=5;
109
110//G4_C  index=6
111  G4double M6[NDENSARRAY]={28.803,2.376,2.9925,-0.0351,2.486,0.2024,3.0036,0.1,0.038};
112  AddMaterial(M6,"G4_C");
113  indexZ[6]=6;
114
115//G4_N  index=7
116  G4double M7[NDENSARRAY]={0.695,1.984,10.5400,1.7378,4.1323,0.15349,3.2125,0.0,0.086};
117  AddMaterial(M7,"G4_N");
118  indexZ[7]=7;
119  state[7]=kStateGas;
120
121//G4_O  index=8
122  G4double M8[NDENSARRAY]={0.744,2.314,10.7004,1.7541,4.3213,0.11778,3.2913,0.0,0.101};
123  AddMaterial(M8,"G4_O");
124  indexZ[8]=8;
125  state[8]=kStateGas;
126
127//G4_F  index=9
128  G4double M9[NDENSARRAY]={0.788,2.450,10.9653,1.8433,4.4096,0.11083,3.2962,0.0,0.121};
129  AddMaterial(M9,"G4_F");
130  indexZ[9]=9;
131  state[9]=kStateGas;
132
133//G4_Ne  index=10
134  G4double M10[NDENSARRAY]={0.587,2.577,11.9041,2.0735,4.6421,0.08064,3.5771,0.0,0.110};
135  AddMaterial(M10,"G4_Ne");
136  indexZ[10]=10;
137  state[10]=kStateGas;
138
139//G4_Na  index=11
140  G4double M11[NDENSARRAY]={19.641,2.648,5.0526,0.2880,3.1962,0.07772,3.6452,0.08,0.098};
141  AddMaterial(M11,"G4_Na");
142  indexZ[11]=11;
143
144//G4_Mg  index=12
145  G4double M12[NDENSARRAY]={26.708,2.331,4.5297,0.1499,3.0668,0.081163,3.6166,0.08,0.073};
146  AddMaterial(M12,"G4_Mg");
147  indexZ[12]=12;
148
149//G4_Al  index=13
150  G4double M13[NDENSARRAY]={32.86,2.18,4.2395,0.1708,3.0127,0.08024,3.6345,0.12,0.061};
151  AddMaterial(M13,"G4_Al");
152  indexZ[13]=13;
153
154//G4_Si  index=14
155  G4double M14[NDENSARRAY]={31.055,2.103,4.4351,0.2014,2.8715,0.14921,3.2546,0.14,0.059};
156  AddMaterial(M14,"G4_Si");
157  indexZ[14]=14;
158
159//G4_P  index=15
160  G4double M15[NDENSARRAY]={29.743,2.056,4.5214,0.1696,2.7815,0.2361,2.9158,0.14,0.057};
161  AddMaterial(M15,"G4_P");
162  indexZ[15]=15;
163
164//G4_S  index=16
165  G4double M16[NDENSARRAY]={28.789,2.131,4.6659,0.158,2.7159,0.33992,2.6456,0.14,0.059};
166  AddMaterial(M16,"G4_S");
167  indexZ[16]=16;
168
169//G4_Cl  index=17
170  G4double M17[NDENSARRAY]={1.092,1.734,11.1421,1.5555,4.2994,0.19849,2.9702,0.0,0.041};
171  AddMaterial(M17,"G4_Cl");
172  indexZ[17]=17;
173  state[17]=kStateGas;
174
175//G4_Ar  index=18
176  G4double M18[NDENSARRAY]={0.789,1.753,11.9480,1.7635,4.4855,0.19714,2.9618,0.0,0.037};
177  AddMaterial(M18,"G4_Ar");
178  indexZ[18]=18;
179  state[18]=kStateGas;
180
181//G4_K  index=19
182  G4double M19[NDENSARRAY]={18.65,1.830,5.6423,0.3851,3.1724,0.19827,2.9233,0.10,0.035};
183  AddMaterial(M19,"G4_K");
184  indexZ[19]=19;
185
186//G4_Ca  index=20
187  G4double M20[NDENSARRAY]={25.342,1.666,5.0396,0.3228,3.1191,0.15643,3.0745,0.14,0.031};
188  AddMaterial(M20,"G4_Ca");
189  indexZ[20]=20;
190
191//G4_Sc  index=21
192  G4double M21[NDENSARRAY]={34.050,1.826,4.6949,0.1640,3.0593,0.15754,3.0517,0.10,0.027};
193  AddMaterial(M21,"G4_Sc");
194  indexZ[21]=21;
195
196//G4_Ti  index=22
197  G4double M22[NDENSARRAY]={41.619,1.969,4.4450,0.0957,3.0386,0.15662,3.0302,0.12,0.025};
198  AddMaterial(M22,"G4_Ti");
199  indexZ[22]=22;
200
201//G4_V  index=23
202  G4double M23[NDENSARRAY]={47.861,2.070,4.2659,0.0691,3.0322,0.15436,3.0163,0.14,0.024};
203  AddMaterial(M23,"G4_V");
204  indexZ[23]=23;
205
206//G4_Cr  index=24
207  G4double M24[NDENSARRAY]={52.458,2.181,4.1781,0.0340,3.0451,0.15419,2.9896,0.14,0.023};
208  AddMaterial(M24,"G4_Cr");
209  indexZ[24]=24;
210
211//G4_Mn  index=25
212  G4double M25[NDENSARRAY]={53.022,2.347,4.2702,0.0447,3.1074,0.14973,2.9796,0.14,0.021};
213  AddMaterial(M25,"G4_Mn");
214  indexZ[25]=25;
215
216//G4_Fe  index=26
217  G4double M26[NDENSARRAY]={55.172,2.504,4.2911,-0.0012,3.1531,0.146,2.9632,0.12,0.021};
218  AddMaterial(M26,"G4_Fe");
219  indexZ[26]=26;
220
221//G4_Co  index=27
222  G4double M27[NDENSARRAY]={58.188,2.626,4.2601,-0.0187,3.1790,0.14474,2.9502,0.12,0.019};
223  AddMaterial(M27,"G4_Co");
224  indexZ[27]=27;
225
226//G4_Ni  index=28
227  G4double M28[NDENSARRAY]={59.385,2.889,4.3115,-0.0566,3.1851,0.16496,2.843,0.10,0.020};
228  AddMaterial(M28,"G4_Ni");
229  indexZ[28]=28;
230
231//G4_Cu  index=29
232  G4double M29[NDENSARRAY]={58.270,2.956,4.4190,-0.0254,3.2792,0.14339,2.9044,0.08,0.019};
233  AddMaterial(M29,"G4_Cu");
234  indexZ[29]=29;
235
236//G4_Zn  index=30
237  G4double M30[NDENSARRAY]={52.132,3.142,4.6906,0.0049,3.3668,0.14714,2.8652,0.08,0.019};
238  AddMaterial(M30,"G4_Zn");
239  indexZ[30]=30;
240
241//G4_Ga  index=31
242  G4double M31[NDENSARRAY]={46.688,2.747,4.9353,0.2267,3.5434,0.09440,3.1314,0.14,0.019};
243  AddMaterial(M31,"G4_Ga");
244  indexZ[31]=31;
245
246//G4_Ge  index=32
247  G4double M32[NDENSARRAY]={44.141,2.461,5.1411,0.3376,3.6096,0.07188,3.3306,0.14,0.025};
248  AddMaterial(M32,"G4_Ge");
249  indexZ[32]=32;
250
251//G4_As  index=33
252  G4double M33[NDENSARRAY]={45.779,2.219,5.0510,0.1767,3.5702,0.06633,3.4176,0.00,0.030};
253  AddMaterial(M33,"G4_As");
254  indexZ[33]=33;
255
256//G4_Se  index=34
257  G4double M34[NDENSARRAY]={40.112,2.104,5.3210,0.2258,3.6264,0.06568,3.4317,0.10,0.024};
258  AddMaterial(M34,"G4_Se");
259  indexZ[34]=34;
260
261//G4_Br  index=35
262  G4double M35[NDENSARRAY]={1.604,1.845,11.7307,1.5262,4.9899,0.06335,3.467,0,0.022};
263  AddMaterial(M35,"G4_Br");
264  indexZ[35]=35;
265  state[35]=kStateGas;
266
267//G4_Kr  index=36
268  G4double M36[NDENSARRAY]={1.114,1.77,12.5115,1.7158,5.0748,0.07446,3.4051,0,0.025};
269  AddMaterial(M36,"G4_Kr");
270  indexZ[36]=36;
271  state[36]=kStateGas;
272
273//G4_Ru  index=37
274  G4double M37[NDENSARRAY]={23.467,1.823,6.4776,0.5737,3.7995,0.07261,3.4177,0.14,0.026};
275  AddMaterial(M37,"G4_Ru");
276  indexZ[37]=37;
277
278//G4_Sr  index=38
279  G4double M38[NDENSARRAY]={30.244,1.707,5.9867,0.4585,3.6778,0.07165,3.4435,0.14,0.026};
280  AddMaterial(M38,"G4_Sr");
281  indexZ[38]=38;
282
283//G4_Y  index=39
284  G4double M39[NDENSARRAY]={40.346,1.649,5.4801,0.3608,3.5542,0.07138,3.4565,0.14,0.027};
285  AddMaterial(M39,"G4_Y");
286  indexZ[39]=39;
287
288//G4_Zr  index=40
289  G4double M40[NDENSARRAY]={48.671,1.638,5.1774,0.2957,3.489,0.07177,3.4533,0.14,0.028};
290  AddMaterial(M40,"G4_Zr");
291  indexZ[40]=40;
292
293//G4_Nb  index=41
294  G4double M41[NDENSARRAY]={56.039,1.734,5.0141,0.1785,3.2201,0.13883,3.093,0.14,0.036};
295  AddMaterial(M41,"G4_Nb");
296  indexZ[41]=41;
297
298//G4_Mo  index=42
299  G4double M42[NDENSARRAY]={60.951,1.658,4.8793,0.2267,3.2784,0.10525,3.2549,0.14,0.03};
300  AddMaterial(M42,"G4_Mo");
301  indexZ[42]=42;
302
303//G4_Tc  index=43
304  G4double M43[NDENSARRAY]={64.760,1.727,4.7769,0.0949,3.1253,0.16572,2.9738,0.14,0.040};
305  AddMaterial(M43,"G4_Tc");
306  indexZ[43]=43;
307
308//G4_Ru  index=44
309  G4double M44[NDENSARRAY]={66.978,1.780,4.7694,0.0599,3.0834,0.19342,2.8707,0.14,0.046};
310  AddMaterial(M44,"G4_Ru");
311  indexZ[44]=44;
312
313//G4_Rh  index=45
314  G4double M45[NDENSARRAY]={67.128,1.804,4.8008,0.0576,3.1069,0.19205,2.8633,0.14,0.046};
315  AddMaterial(M45,"G4_Rh");
316  indexZ[45]=45;
317
318//G4_Pd  index=46
319  G4double M46[NDENSARRAY]={65.683,1.911,4.9358,0.0563,3.0555,0.24178,2.7239,0.14,0.047};
320  AddMaterial(M46,"G4_Pd");
321  indexZ[46]=46;
322
323//G4_Ag  index=47
324  G4double M47[NDENSARRAY]={61.635,1.933,5.0630,0.0657,3.1074,0.24585,2.6899,0.14,0.052};
325  AddMaterial(M47,"G4_Ag");
326  indexZ[47]=47;
327
328//G4_Cd  index=48
329  G4double M48[NDENSARRAY]={55.381,1.895,5.2727,0.1281,3.1667,0.24609,2.6772,0.14,0.051};
330  AddMaterial(M48,"G4_Cd");
331  indexZ[48]=48;
332
333//G4_In  index=49
334  G4double M49[NDENSARRAY]={50.896,1.851,5.5211,0.2406,3.2032,0.23879,2.7144,0.14,0.044};
335  AddMaterial(M49,"G4_In");
336  indexZ[49]=49;
337
338//G4_Sn  index=50
339  G4double M50[NDENSARRAY]={50.567,1.732,5.5340,0.2879,3.2959,0.18689,2.8576,0.14,0.037};
340  AddMaterial(M50,"G4_Sn");
341  indexZ[50]=50;
342
343//G4_Sb  index=51
344  G4double M51[NDENSARRAY]={48.242,1.645,5.6241,0.3189,3.3489,0.16652,2.9319,0.14,0.034};
345  AddMaterial(M51,"G4_Sb");
346  indexZ[51]=51;
347
348//G4_Te  index=52
349  G4double M52[NDENSARRAY]={45.952,1.577,5.7131,0.3296,3.4418,0.13815,3.0354,0.14,0.033};
350  AddMaterial(M52,"G4_Te");
351  indexZ[52]=52;
352
353//G4_I  index=53
354  G4double M53[NDENSARRAY]={41.348,1.498,5.9488,0.0549,3.2596,0.23766,2.7276,0.0,0.045};
355  AddMaterial(M53,"G4_I");
356  indexZ[53]=53;
357
358//G4_Xe  index=54
359  G4double M54[NDENSARRAY]={1.369,1.435,12.7281,1.563,4.7371,0.23314,2.7414,0,0.043};
360  AddMaterial(M54,"G4_Xe");
361  indexZ[54]=54;
362  state[54]=kStateGas;
363
364//G4_Cs  index=55
365  G4double M55[NDENSARRAY]={25.37,1.462,6.9135,0.5473,3.5914,0.18233,2.8866,0.14,0.035};
366  AddMaterial(M55,"G4_Cs");
367  indexZ[55]=55;
368
369//G4_Ba  index=56
370  G4double M56[NDENSARRAY]={34.425,1.410,6.3153,0.4190,3.4547,0.18268,2.8906,0.14,0.035};
371  AddMaterial(M56,"G4_Ba");
372  indexZ[56]=56;
373
374//G4_La  index=57
375  G4double M57[NDENSARRAY]={45.792,1.392,5.7850,0.3161,3.3293,0.18591,2.8828,0.14,0.036};
376  AddMaterial(M57,"G4_La");
377  indexZ[57]=57;
378
379//G4_Ce  index=58
380  G4double M58[NDENSARRAY]={47.834,1.461,5.7837,0.2713,3.3432,0.18885,2.8592,0.14,0.040};
381  AddMaterial(M58,"G4_Ce");
382  indexZ[58]=58;
383
384//G4_Pr  index=59
385  G4double M59[NDENSARRAY]={48.301,1.520,5.8096,0.2333,3.2773,0.23265,2.7331,0.14,0.041};
386  AddMaterial(M59,"G4_Pr");
387  indexZ[59]=59;
388
389//G4_Ne  index=60
390  G4double M60[NDENSARRAY]={48.819,1.588,5.8290,0.1984,3.3063,0.23530,2.7050,0.14,0.044};
391  AddMaterial(M60,"G4_Ne");
392  indexZ[60]=60;
393
394//G4_Pr  index=61
395  G4double M61[NDENSARRAY]={50.236,1.672,5.8224,0.1627,3.3199,0.24280,2.6674,0.14,0.048};
396  AddMaterial(M61,"G4_Pr");
397  indexZ[61]=61;
398
399//G4_Sa  index=62
400  G4double M62[NDENSARRAY]={50.540,1.749,5.8597,0.1520,3.3460,0.24698,2.6403,0.14,0.053};
401  AddMaterial(M62,"G4_Sa");
402  indexZ[62]=62;
403
404//G4_Eu  index=63
405  G4double M63[NDENSARRAY]={42.484,1.838,6.2278,0.1888,3.4633,0.24448,2.6245,0.14,0.06};
406  AddMaterial(M63,"G4_Eu");
407  indexZ[63]=63;
408
409//G4_Gd  index=64
410  G4double M64[NDENSARRAY]={51.672,1.882,5.8738,0.1058,3.3932,0.25109,2.5977,0.14,0.061};
411  AddMaterial(M64,"G4_Gd");
412  indexZ[64]=64;
413
414//G4_Tb  index=65
415  G4double M65[NDENSARRAY]={52.865,1.993,5.9045,0.0947,3.4224,0.24453,2.6056,0.14,0.063};
416  AddMaterial(M65,"G4_Tb");
417  indexZ[65]=65;
418
419//G4_Dy  index=66
420  G4double M66[NDENSARRAY]={53.698,2.081,5.9183,0.0822,3.4474,0.24665,2.5849,0.14,0.061};
421  AddMaterial(M66,"G4_Dy");
422  indexZ[66]=66;
423
424//G4_Ho  index=67
425  G4double M67[NDENSARRAY]={54.467,2.197,5.9587,0.0761,3.4782,0.24638,2.5726,0.14,0.062};
426  AddMaterial(M67,"G4_Ho");
427  indexZ[67]=67;
428
429//G4_Er  index=68
430  G4double M68[NDENSARRAY]={55.322,2.26,5.9521,0.0648,3.4922,0.24823,2.5573,0.14,0.061};
431  AddMaterial(M68,"G4_Er");
432  indexZ[68]=68;
433
434//G4_Tm  index=69
435  G4double M69[NDENSARRAY]={56.225,2.333,5.9677,0.0812,3.5085,0.24189,2.5469,0.14,0.062};
436  AddMaterial(M69,"G4_Tm");
437  indexZ[69]=69;
438
439//G4_Yb  index=70
440  G4double M70[NDENSARRAY]={47.546,2.505,6.3325,0.1199,3.6246,0.25295,2.5141,0.14,0.071};
441  AddMaterial(M70,"G4_Yb");
442  indexZ[70]=70;
443
444//G4_Lu  index=71
445  G4double M71[NDENSARRAY]={57.581,2.348,5.9785,0.1560,3.5218,0.24033,2.5643,0.14,0.054};
446  AddMaterial(M71,"G4_Lu");
447  indexZ[71]=71;
448
449//G4_Hf  index=72
450  G4double M72[NDENSARRAY]={66.770,2.174,5.7139,0.1965,3.4337,0.22918,2.6155,0.14,0.035};
451  AddMaterial(M72,"G4_Hf");
452  indexZ[72]=72;
453
454//G4_Ta  index=73
455  G4double M73[NDENSARRAY]={74.692,2.07,5.5262,0.2117,3.4805,0.17798,2.7623,0.14,0.03};
456  AddMaterial(M73,"G4_Ta");
457  indexZ[73]=73;
458
459//G4_W  index=74
460  G4double M74[NDENSARRAY]={80.315,1.997,5.4059,0.2167,3.496,0.15509,2.8447,0.14,0.027};
461  AddMaterial(M74,"G4_W");
462  indexZ[74]=74;
463
464//G4_Re  index=75
465  G4double M75[NDENSARRAY]={83.846,1.976,5.3445,0.0559,3.4845,0.15184,2.8627,0.08,0.026};
466  AddMaterial(M75,"G4_Re");
467  indexZ[75]=75;
468
469//G4_Os  index=76
470  G4double M76[NDENSARRAY]={86.537,1.947,5.3083,0.0891,3.5414,0.12751,2.9608,0.10,0.023};
471  AddMaterial(M76,"G4_Os");
472  indexZ[76]=76;
473
474//G4_Ir  index=77
475  G4double M77[NDENSARRAY]={86.357,1.927,5.3418,0.0819,3.5480,0.12690,2.9658,0.10,0.023};
476  AddMaterial(M77,"G4_Ir");
477  indexZ[77]=77;
478
479//G4_Pt  index=78
480  G4double M78[NDENSARRAY]={84.389,1.965,5.4732,0.1484,3.6212,0.11128,3.0417,0.12,0.021};
481  AddMaterial(M78,"G4_Pt");
482  indexZ[78]=78;
483
484//G4_Au  index=79
485  G4double M79[NDENSARRAY]={80.215,1.926,5.5747,0.2021,3.6979,0.09756,3.1101,0.14,0.020};
486  AddMaterial(M79,"G4_Au");
487  indexZ[79]=79;
488
489//G4_Hg  index=80
490  G4double M80[NDENSARRAY]={66.977,1.904,5.9605,0.2756,3.7275,0.11014,3.0519,0.14,0.021};
491  AddMaterial(M80,"G4_Hg");
492  indexZ[80]=80;
493
494//G4_Tl  index=81
495  G4double M81[NDENSARRAY]={62.104,1.814,6.1365,0.3491,3.8044,0.09455,3.1450,0.14,0.019};
496  AddMaterial(M81,"G4_Tl");
497  indexZ[81]=81;
498
499//G4_Pb  index=82
500  G4double M82[NDENSARRAY]={61.072,1.755,6.2018,0.3776,3.8073,0.09359,3.1608,0.14,0.019};
501  AddMaterial(M82,"G4_Pb");
502  indexZ[82]=82;
503
504//G4_Bi  index=83
505  G4double M83[NDENSARRAY]={56.696,1.684,6.3505,0.4152,3.8248,0.0941,3.1671,0.14,0.02};
506  AddMaterial(M83,"G4_Bi");
507  indexZ[83]=83;
508
509//G4_Po  index=84
510  G4double M84[NDENSARRAY]={55.773,1.637,6.4003,0.4267,3.8293,0.09282,3.183,0.14,0.02};
511  AddMaterial(M84,"G4_Po");
512  indexZ[84]=84;
513
514//G4_Rn  index=85
515  G4double M85[NDENSARRAY]={1.708,1.458,13.2839,1.5368,4.9889,0.20798,2.7409,0,0.057};
516  AddMaterial(M85,"G4_Rn");
517  indexZ[86]=85;
518  state[86]=kStateGas;
519
520//G4_Ra  index=86
521  G4double M86[NDENSARRAY]={40.205,1.403,7.0452,0.5991,3.9428,0.08804,3.2454,0.14,0.022};
522  AddMaterial(M86,"G4_Ra");
523  indexZ[88]=86;
524
525//G4_Ac  index=87
526  G4double M87[NDENSARRAY]={57.254,1.380,6.3742,0.4559,3.7966,0.08567,3.2683,0.14,0.023};
527  AddMaterial(M87,"G4_Ac");
528  indexZ[89]=87;
529
530//G4_Th  index=88
531  G4double M88[NDENSARRAY]={61.438,1.363,6.2473,0.4202,3.7681,0.08655,3.2610,0.14,0.025};
532  AddMaterial(M88,"G4_Th");
533  indexZ[90]=88;
534
535//G4_Pa  index=89
536  G4double M89[NDENSARRAY]={70.901,1.42,6.0327,0.3144,3.5079,.14770,2.9845,0.14,0.036};
537  AddMaterial(M89,"G4_Pa");
538  indexZ[91]=89;
539
540//G4_U  index=90
541  G4double M90[NDENSARRAY]={77.986,1.447,5.8694,0.2260,3.3721,.19677,2.8171,0.14,0.043};
542  AddMaterial(M90,"G4_U");
543  indexZ[92]=90;
544
545//G4_Np  index=91
546  G4double M91[NDENSARRAY]={81.221,1.468,5.8149,0.1869,3.369,0.19741,2.8082,0.14,0.043};
547  AddMaterial(M91,"G4_Np");
548  indexZ[93]=91;
549
550//G4_Pu  index=92
551  G4double M92[NDENSARRAY]={80.486,1.519,5.8748,0.1557,3.3981,0.20419,2.7679,0.14,0.057};
552  AddMaterial(M92,"G4_Pu");
553  indexZ[94]=92;
554
555//G4_Am  index=93
556  G4double M93[NDENSARRAY]={66.607,1.552,6.2813,0.2274,3.5021,0.20308,2.7615,0.14,0.056};
557  AddMaterial(M93,"G4_Am");
558  indexZ[95]=93;
559
560//G4_Cm  index=94
561  G4double M94[NDENSARRAY]={66.022,1.559,6.3097,0.2484,3.516,.20257,2.7579,0.14,0.056};
562  AddMaterial(M94,"G4_Cm");
563  indexZ[96]=94;
564
565//G4_Bk  index=95
566  G4double M95[NDENSARRAY]={67.557,1.574,6.2912,0.2378,3.5186,.20192,2.7560,0.14,0.062};
567  AddMaterial(M95,"G4_Bk");
568  indexZ[97]=95;
569
570//G4_A-150_TISSUE  index=96
571  G4double M96[NDENSARRAY]={22.667,1.950,3.1100,0.1329,2.6234,0.10783,3.4442,0,0.048};
572  AddMaterial(M96,"G4_A-150_TISSUE");
573
574//G4_ACETONE  index=97
575  G4double M97[NDENSARRAY]={19.010,1.976,3.4341,0.2197,2.6928,0.11100,3.4047,0,0.069};
576  AddMaterial(M97,"G4_ACETONE");
577
578//G4_ACETYLENE  index=98
579  G4double M98[NDENSARRAY]={0.700,1.784,9.8419,1.6017,4.0074,0.12167,3.4277,0,0.080};
580  AddMaterial(M98,"G4_ACETYLENE");
581
582//G4_ADENINE  index=99
583  G4double M99[NDENSARRAY]={24.098,1.892,3.1724,0.1295,2.4219,0.20908,3.0271,0,0.052};
584  AddMaterial(M99,"G4_ADENINE");
585
586//G4_ADIPOSE_TISSUE_ICRP  index=100
587  G4double M100[NDENSARRAY]={20.655,1.987,3.2367,0.1827,2.6530,0.10278,3.4817,0,0.060};
588  AddMaterial(M100,"G4_ADIPOSE_TISSUE_ICRP");
589
590//G4_AIR  index=101
591  G4double M101[NDENSARRAY]={0.707,2.054,10.5961,1.7418,4.2759,0.10914,3.3994,0,0.090};
592  AddMaterial(M101,"G4_AIR");
593
594//G4_ALANINE  index=102
595  G4double M102[NDENSARRAY]={25.204,2.074,3.0965,0.1354,2.6336,0.11484,3.3526,0,0.056};
596  AddMaterial(M102,"G4_ALANINE");
597
598//G4_ALUMINIM_OXIDE  index=103
599  G4double M103[NDENSARRAY]={40.206,2.394,3.5682,0.0402,2.8665,0.08500,3.5458,0,0.031};
600  AddMaterial(M103,"G4_ALUMINIM_OXIDE");
601
602//G4_AMBER  index=104
603  G4double M104[NDENSARRAY]={22.450,1.946,3.0701,0.1335,2.5610,0.11934,3.4098,0,0.053};
604  AddMaterial(M104,"G4_AMBER");
605
606//G4_AMMONIA  index=105
607  G4double M105[NDENSARRAY]={0.635,1.814,9.8763,1.6822,4.1158,0.08315,3.6464,0,0.102};
608  AddMaterial(M105,"G4_AMMONIA");
609
610//G4_ANILINE  index=106
611  G4double M106[NDENSARRAY]={21.361,1.938,3.2622,0.1618,2.5805,0.13134,3.3434,0,0.052};
612  AddMaterial(M106,"G4_ANILINE");
613
614//G4_ANTHRACENE  index=107
615  G4double M107[NDENSARRAY]={23.704,1.954,3.1514,0.1146,2.5213,0.14677,3.2831,0,0.042};
616  AddMaterial(M107,"G4_ANTHRACENE");
617
618//G4_B-100_BONE  index=108
619  G4double M108[NDENSARRAY]={25.199,2.013,3.4528,0.1252,3.042,0.05268,3.7365,0,0.043};
620  AddMaterial(M108,"G4_B-100_BONE");
621
622//G4_BAKELITE  index=109
623  G4double M109[NDENSARRAY]={23.408,2.046,3.2582,0.1471,2.6055,0.12713,3.347,0,0.052};
624  AddMaterial(M109,"G4_BAKELITE");
625
626//G4_BARIUM_FLUORIDE  index=110
627  G4double M110[NDENSARRAY]={41.398,1.727,5.4122,-0.0098,3.3871,0.15991,2.8867,0,0.034};
628  AddMaterial(M110,"G4_BARIUM_FLUORIDE");
629
630//G4_BARIUM_SULFATE  index=111
631  G4double M111[NDENSARRAY]={40.805,1.893,4.8923,-0.0128,3.4069,0.11747,3.0427,0,0.03};
632  AddMaterial(M111,"G4_BARIUM_SULFATE");
633
634//G4_BENZENE  index=112
635  G4double M112[NDENSARRAY]={19.806,1.873,3.3269,0.171,2.5091,0.16519,3.2174,0,0.052};
636  AddMaterial(M112,"G4_BENZENE");
637
638//G4_BERYLLIUM_OXIDE  index=113
639  G4double M113[NDENSARRAY]={34.629,2.296,2.9801,0.0241,2.5846,0.10755,3.4927,0,0.031};
640  AddMaterial(M113,"G4_BERYLLIUM_OXIDE");
641
642//G4_BGO  index=114
643  G4double M114[NDENSARRAY]={49.904,2.121,5.7409,0.0456,3.7816,0.09569,3.0781,0,0.023};
644  AddMaterial(M114,"G4_BGO");
645
646//G4_BLOOD_ICRP  index=115
647  G4double M115[NDENSARRAY]={22.001,2.184,3.4581,0.2239,2.8017,0.08492,3.5406,0,0.088};
648  AddMaterial(M115,"G4_BLOOD_ICRP");
649
650//G4_BONE_COMPACT_ICRU  index=116
651  G4double M116[NDENSARRAY]={28.536,2.091,3.3390,0.0944,3.0201,0.05822,3.6419,0,0.042};
652  AddMaterial(M116,"G4_BONE_COMPACT_ICRU");
653
654//G4_BONE_CORTICAL_ICRP  index=117
655  G4double M117[NDENSARRAY]={28.298,2.118,3.6488,0.1161,3.0919,0.06198,3.5919,0,0.04};
656  AddMaterial(M117,"G4_BONE_CORTICAL_ICRP");
657
658//G4_BORON_CARBIDE  index=118
659  G4double M118[NDENSARRAY]={31.38,2.14,2.9859,0.0093,2.1006,0.37087,2.8076,0,0.022};
660  AddMaterial(M118,"G4_BORON_CARBIDE");
661
662//G4_BORON_OXIDE  index=119
663  G4double M119[NDENSARRAY]={27.107,2.446,3.6027,0.1843,2.7379,0.11548,3.3832,0,0.053};
664  AddMaterial(M119,"G4_BORON_OXIDE");
665
666//G4_BRAIN_ICRP  index=120
667  G4double M120[NDENSARRAY]={21.772,2.162,3.4279,0.2206,2.8021,0.08255,3.5585,0,0.086};
668  AddMaterial(M120,"G4_BRAIN_ICRP");
669
670//G4_BUTANE  index=121
671  G4double M121[NDENSARRAY]={1.101,1.727,8.5633,1.3788,3.7524,0.10852,3.4184,0,0.1};
672  AddMaterial(M121,"G4_BUTANE");
673
674//G4_N-BUTYL_ALCOHOL  index=122
675  G4double M122[NDENSARRAY]={19.52,1.942,3.2425,0.1937,2.6439,0.10081,3.5139,0,0.065};
676  AddMaterial(M122,"G4_N-BUTYL_ALCOHOL");
677
678//G4_C-552  index=123
679  G4double M123[NDENSARRAY]={27.023,2.128,3.3338,0.151,2.7083,0.10492,3.4344,0,0.053};
680  AddMaterial(M123,"G4_C-552");
681
682//G4_CADMIUM_TELLURIDE  index=124
683  G4double M124[NDENSARRAY]={46.314,1.935,5.9096,0.0438,3.2836,0.2484,2.6665,0,0.057};
684  AddMaterial(M124,"G4_CADMIUM_TELLURIDE");
685
686//G4_CADMIUM_TUNGSTATE  index=125
687  G4double M125[NDENSARRAY]={52.954,2.289,5.3594,0.0123,3.5941,0.12861,2.915,0,0.027};
688  AddMaterial(M125,"G4_CADMIUM_TUNGSTATE");
689
690//G4_CALCIUM_CARBONATE  index=126
691  G4double M126[NDENSARRAY]={34.08,2.141,3.7738,0.0492,3.0549,0.08301,3.412,0,0.037};
692  AddMaterial(M126,"G4_CALCIUM_CARBONATE");
693
694//G4_CALCIUM_FLUORIDE  index=127
695  G4double M127[NDENSARRAY]={35.849,2.127,4.0653,0.0676,3.1683,0.06942,3.5263,0,0.044};
696  AddMaterial(M127,"G4_CALCIUM_FLUORIDE");
697
698//G4_CALCIUM_OXIDE  index=128
699  G4double M128[NDENSARRAY]={36.988,1.973,4.1209,-0.0172,3.0171,0.12128,3.1936,0,0.024};
700  AddMaterial(M128,"G4_CALCIUM_OXIDE");
701
702//G4_CALCIUM_SULFATE  index=129
703  G4double M129[NDENSARRAY]={35.038,2.179,3.9388,0.0587,3.1229,0.07708,3.4495,0,0.021};
704  AddMaterial(M129,"G4_CALCIUM_SULFATE");
705
706//G4_CALCIUM_TUNGSTATE  index=130
707  G4double M130[NDENSARRAY]={46.934,2.262,5.2603,0.0323,3.8932,0.0621,3.2649,0,0.021};
708  AddMaterial(M130,"G4_CALCIUM_TUNGSTATE");
709
710//G4_CARBON_DIOXIDE  index=131
711  G4double M131[NDENSARRAY]={0.874,2.118,10.1537,1.6294,4.1825,0.11768,3.3227,0,0.091};
712  AddMaterial(M131,"G4_CARBON_DIOXIDE");
713
714//G4_CARBON_TETRACHLORIDE  index=132
715  G4double M132[NDENSARRAY]={25.234,1.742,4.7712,0.1773,2.9165,0.19018,3.0116,0,0.041};
716  AddMaterial(M132,"G4_CARBON_TETRACHLORIDE");
717
718//G4_CELLULOSE_CELLOPHANE  index=133
719  G4double M133[NDENSARRAY]={25.008,2.17,3.2647,0.158,2.6778,0.11151,3.381,0,0.06};
720  AddMaterial(M133,"G4_CELLULOSE_CELLOPHANE");
721
722//G4_CELLULOSE_BUTYRATE  index=134
723  G4double M134[NDENSARRAY]={23.041,2.128,3.3497,0.1794,2.6809,0.11444,3.3738,0,0.056};
724  AddMaterial(M134,"G4_CELLULOSE_BUTYRATE");
725
726//G4_CELLULOSE_NITRATE  index=135
727  G4double M135[NDENSARRAY]={25.224,2.252,3.4762,0.1897,2.7253,0.11813,3.3237,0,0.063};
728  AddMaterial(M135,"G4_CELLULOSE_NITRATE");
729
730//G4_CERIC_SULFATE  index=136
731  G4double M136[NDENSARRAY]={21.743,2.205,3.5212,0.2363,2.8769,0.07666,3.5607,0,0.095};
732  AddMaterial(M136,"G4_CERIC_SULFATE");
733
734//G4_CESIUM_FLUORIDE  index=137
735  G4double M137[NDENSARRAY]={37.942,1.714,5.9046,0.0084,3.3374,0.22052,2.728,0,0.044};
736  AddMaterial(M137,"G4_CESIUM_FLUORIDE");
737
738//G4_CESIUM_IODIDE  index=138
739  G4double M138[NDENSARRAY]={39.455,1.672,6.2807,0.0395,3.3353,0.25381,2.6657,0,0.067};
740  AddMaterial(M138,"G4_CESIUM_IODIDE");
741
742//G4_CHLOROBENZENE  index=139
743  G4double M139[NDENSARRAY]={21.752,1.889,3.8201,0.1714,2.9272,0.09856,3.3797,0,0.031};
744  AddMaterial(M139,"G4_CHLOROBENZENE");
745
746//G4_CHLOROFORM  index=140
747  G4double M140[NDENSARRAY]={24.462,1.734,4.7055,0.1786,2.9581,0.16959,3.0627,0,0.038};
748  AddMaterial(M140,"G4_CHLOROFORM");
749
750//G4_CONCRETE  index=141
751  G4double M141[NDENSARRAY]={30.986,2.322,3.9464,0.1301,3.0466,0.07515,3.5467,0,0.024};
752  AddMaterial(M141,"G4_CONCRETE");
753
754//G4_CYCLOHEXANE  index=142
755  G4double M142[NDENSARRAY]={19.207,1.861,3.1544,0.1728,2.5549,0.12035,3.4278,0,0.057};
756  AddMaterial(M142,"G4_CYCLOHEXANE");
757
758//G4_1,2-DICHLOROBENZENE  index=143
759  G4double M143[NDENSARRAY]={23.354,1.862,4.0348,0.1587,2.8276,0.1601,3.0836,0,0.029};
760  AddMaterial(M143,"G4_1,2-DICHLOROBENZENE");
761
762//G4_DICHLORODIETHYL_ETHER  index=144
763  G4double M144[NDENSARRAY]={22.894,1.903,4.0135,0.1773,3.1586,0.06799,3.525,0,0.026};
764  AddMaterial(M144,"G4_DICHLORODIETHYL_ETHER");
765
766//G4_1,2-DICHLOROETHANE  index=145
767  G4double M145[NDENSARRAY]={22.764,1.618,4.1849,0.1375,2.9529,0.13383,3.1675,0,0.03};
768  AddMaterial(M145,"G4_1,2-DICHLOROETHANE");
769
770//G4_DIETHYL_ETHER  index=146
771  G4double M146[NDENSARRAY]={18.326,1.951,3.3721,0.2231,2.6745,0.1055,3.4586,0,0.07};
772  AddMaterial(M146,"G4_DIETHYL_ETHER");
773
774//G4_N,N-DIMETHYL_FORMAMIDE  index=147
775  G4double M147[NDENSARRAY]={20.763,2.005,3.3311,0.1977,2.6686,0.1147,3.371,0,0.065};
776  AddMaterial(M147,"G4_N,N-DIMETHYL_FORMAMIDE");
777
778//G4_DIMETHYL_SULFOXIDE  index=148
779  G4double M148[NDENSARRAY]={22.173,2.075,3.9844,0.2021,3.1263,0.06619,3.5708,0,0.03};
780  AddMaterial(M148,"G4_DIMETHYL_SULFOXIDE");
781
782//G4_ETHANE  index=149
783  G4double M149[NDENSARRAY]={0.789,1.69,9.1043,1.5107,3.8743,0.09627,3.6095,0,0.097};
784  AddMaterial(M149,"G4_ETHANE");
785
786//G4_ETHYL_ALCOHOL  index=150
787  G4double M150[NDENSARRAY]={19.232,2.013,3.3699,0.2218,2.7052,0.09878,3.4834,0,0.071};
788  AddMaterial(M150,"G4_ETHYL_ALCOHOL");
789
790//G4_ETHYL_CELLULOSE  index=151
791  G4double M151[NDENSARRAY]={22.594,2.065,3.2415,0.1683,2.6527,0.11077,3.4098,0,0.057};
792  AddMaterial(M151,"G4_ETHYL_CELLULOSE");
793
794//G4_ETHYLENE  index=152
795  G4double M152[NDENSARRAY]={0.746,1.733,9.438,1.5528,3.9327,0.10636,3.5387,0,0.085};
796  AddMaterial(M152,"G4_ETHYLENE");
797
798//G4_EYE_LENS_ICRP  index=153
799  G4double M153[NDENSARRAY]={22.388,2.154,3.372,0.207,2.7446,0.0969,3.455,0,0.077};
800  AddMaterial(M153,"G4_EYE_LENS_ICRP");
801
802//G4_FERRIC_OXIDE  index=154
803  G4double M154[NDENSARRAY]={45.331,2.747,4.2245,-0.0074,3.2573,0.10478,3.1313,0,0.026};
804  AddMaterial(M154,"G4_FERRIC_OXIDE");
805
806//G4_FERROBORIDE  index=155
807  G4double M155[NDENSARRAY]={52.546,2.726,4.2057,-0.0988,3.1749,0.12911,3.024,0,0.022};
808  AddMaterial(M155,"G4_FERROBORIDE");
809
810//G4_FERROUS_OXIDE  index=156
811  G4double M156[NDENSARRAY]={47.327,2.769,4.3175,-0.0279,3.2002,0.12959,3.0168,0,0.022};
812  AddMaterial(M156,"G4_FERROUS_OXIDE");
813
814//G4_FERROUS_SULFATE  index=157
815  G4double M157[NDENSARRAY]={21.69,2.208,3.5183,0.2378,2.8254,0.08759,3.4923,0,0.096};
816  AddMaterial(M157,"G4_FERROUS_SULFATE");
817
818//G4_FREON-12  index=158
819  G4double M158[NDENSARRAY]={21.121,1.974,4.8251,0.3035,3.2659,0.07978,3.4626,0,0.025};
820  AddMaterial(M158,"G4_FREON-12");
821
822//G4_FREON-12B2  index=159
823  G4double M159[NDENSARRAY]={25.877,2.195,5.7976,0.3406,3.7956,0.05144,3.5565,0,0.021};
824  AddMaterial(M159,"G4_FREON-12B2");
825
826//G4_FREON-13  index=160
827  G4double M160[NDENSARRAY]={19.432,2.116,4.7483,0.3659,3.2337,0.07238,3.5551,0,0.05};
828  AddMaterial(M160,"G4_FREON-13");
829
830//G4_FREON-13B1  index=161
831  G4double M161[NDENSARRAY]={23.849,2.233,5.3555,0.3522,3.7554,0.03925,3.7194,0,0.036};
832  AddMaterial(M161,"G4_FREON-13B1");
833
834//G4_FREON-13I1  index=162
835  G4double M162[NDENSARRAY]={25.615,1.924,5.8774,0.2847,3.728,0.09112,3.1658,0,0.025};
836  AddMaterial(M162,"G4_FREON-13I1");
837
838//G4_GADOLINIUM_OXYSULFIDE  index=163
839  G4double M163[NDENSARRAY]={51.099,2.179,5.5347,-0.1774,3.4045,0.22161,2.63,0,0.056};
840  AddMaterial(M163,"G4_GADOLINIUM_OXYSULFIDE");
841
842//G4_GALLIUM_ARSENIDE  index=164
843  G4double M164[NDENSARRAY]={44.17,2.652,5.3299,0.1764,3.642,0.07152,3.3356,0,0.027};
844  AddMaterial(M164,"G4_GALLIUM_ARSENIDE");
845
846//G4_GEL_PHOTO_EMULSION  index=165
847  G4double M165[NDENSARRAY]={24.058,2.156,3.2687,0.1709,2.7058,0.10102,3.4418,0,0.06};
848  AddMaterial(M165,"G4_GEL_PHOTO_EMULSION");
849
850//G4_Pyrex_Glass  index=166
851  G4double M166[NDENSARRAY]={30.339,2.369,3.9708,0.1479,2.9933,0.0827,3.5224,0,0.022};
852  AddMaterial(M166,"G4_Pyrex_Glass");
853
854//G4_GLASS_LEAD  index=167
855  G4double M167[NDENSARRAY]={46.631,2.085,5.8476,0.0614,3.8146,0.09544,3.074,0,0.025};
856  AddMaterial(M167,"G4_GLASS_LEAD");
857
858//G4_GLASS_PLATE  index=168
859  G4double M168[NDENSARRAY]={31.481,2.329,4.0602,0.1237,3.0649,0.07678,3.5381,0,0.025};
860  AddMaterial(M168,"G4_GLASS_PLATE");
861
862//G4_GLUCOSE  index=169
863  G4double M169[NDENSARRAY]={26.153,2.174,3.1649,0.1411,2.67,0.10783,3.3946,0,0.061};
864  AddMaterial(M169,"G4_GLUCOSE");
865
866//G4_GLUTAMINE  index=170
867  G4double M170[NDENSARRAY]={25.437,2.077,3.1167,0.1347,2.6301,0.11931,3.3254,0,0.055};
868  AddMaterial(M170,"G4_GLUTAMINE");
869
870//G4_GLYCEROL  index=171
871  G4double M171[NDENSARRAY]={23.846,2.12,3.2267,0.1653,2.6862,0.10168,3.4481,0,0.067};
872  AddMaterial(M171,"G4_GLYCEROL");
873
874//G4_GUANINE  index=172
875  G4double M172[NDENSARRAY]={26.022,1.97,3.1171,0.1163,2.4296,0.2053,3.0186,0,0.069};
876  AddMaterial(M172,"G4_GUANINE");
877
878//G4_GYPSUM  index=173
879  G4double M173[NDENSARRAY]={31.379,2.187,3.8382,0.0995,3.1206,0.06949,3.5134,0,0.038};
880  AddMaterial(M173,"G4_GYPSUM");
881
882//G4_N-HEPTANE  index=174
883  G4double M174[NDENSARRAY]={18.128,1.848,3.1978,0.1928,2.5706,0.11255,3.4885,0,0.059};
884  AddMaterial(M174,"G4_N-HEPTANE");
885
886//G4_N-HEXANE  index=175
887  G4double M175[NDENSARRAY]={17.836,1.843,3.2156,0.1984,2.5757,0.11085,3.5027,0,0.061};
888  AddMaterial(M175,"G4_N-HEXANE");
889
890//G4_KAPTON  index=176
891  G4double M176[NDENSARRAY]={24.586,2.109,3.3497,0.1509,2.5631,0.15972,3.1921,0,0.05};
892  AddMaterial(M176,"G4_KAPTON");
893
894//G4_LANTHANUM_OXYBROMIDE  index=177
895  G4double M177[NDENSARRAY]={47.125,1.831,5.4666,-0.035,3.3288,0.1783,2.8457,0,0.04};
896  AddMaterial(M177,"G4_LANTHANUM_OXYBROMIDE");
897
898//G4_LANTHANUM_OXYSULFIDE  index=178
899  G4double M178[NDENSARRAY]={45.394,1.681,5.6151,-0.0934,3.2741,0.22579,2.7075,0,0.065};
900  AddMaterial(M178,"G4_LANTHANUM_OXYSULFIDE");
901
902//G4_LEAD_OXIDE  index=179
903  G4double M179[NDENSARRAY]={56.488,2.012,6.2162,0.0356,3.5456,0.19645,2.7299,0,0.039};
904  AddMaterial(M179,"G4_LEAD_OXIDE");
905
906//G4_LITHIUM_AMIDE  index=180
907  G4double M180[NDENSARRAY]={22.609,1.74,2.7961,0.0198,2.5152,0.0874,3.7534,0,0.05};
908  AddMaterial(M180,"G4_LITHIUM_AMIDE");
909
910//G4_LITHIUM_CARBONATE index=181
911  G4double M181[NDENSARRAY]={29.217,2.246,3.2029,0.0551,2.6598,0.09936,3.5417,0,0.062};
912  AddMaterial(M181,"G4_LITHIUM_CARBONATE");
913
914//G4_LITHIUM_FLUORIDE  index=182
915  G4double M182[NDENSARRAY]={31.815,2.197,3.1667,0.0171,2.7049,0.07593,3.7478,0,0.084};
916  AddMaterial(M182,"G4_LITHIUM_FLUORIDE");
917
918//G4_LITHIUM_HYDRIDE  index=183
919  G4double M183[NDENSARRAY]={18.51,1.482,2.358,-0.0988,1.4515,0.90567,2.5849,0,0.035};
920  AddMaterial(M183,"G4_LITHIUM_HYDRIDE");
921
922//G4_LITHIUM_IODIDE  index=184
923  G4double M184[NDENSARRAY]={34.841,1.706,6.2671,0.0892,3.3702,0.23274,2.7146,0,0.043};
924  AddMaterial(M184,"G4_LITHIUM_IODIDE");
925
926//G4_LITHIUM_OXIDE  index=185
927  G4double M185[NDENSARRAY]={27.984,2.039,2.934,-0.0511,2.5874,0.08035,3.7878,0,0.043};
928  AddMaterial(M185,"G4_LITHIUM_OXIDE");
929
930//G4_LITHIUM_TETRABORATE  index=186
931  G4double M186[NDENSARRAY]={31.343,2.36,3.2093,0.0737,2.6502,0.11075,3.4389,0,0.048};
932  AddMaterial(M186,"G4_LITHIUM_TETRABORATE");
933
934//G4_LUNG_ICRP  index=187
935  G4double M187[NDENSARRAY]={21.891,2.184,3.4708,0.2261,2.8001,0.08588,3.5353,0,0.089};
936  AddMaterial(M187,"G4_LUNG_ICRP");
937
938//G4_M3_WAX  index=188
939  G4double M188[NDENSARRAY]={22,1.975,3.254,0.1523,2.7529,0.07864,3.6412,0,0.044};
940  AddMaterial(M188,"G4_M3_WAX");
941
942//G4_MAGNESIUM_CARBONATE  index=189
943  G4double M189[NDENSARRAY]={34.979,2.388,3.4319,0.086,2.7997,0.09219,3.5003,0,0.045};
944  AddMaterial(M189,"G4_MAGNESIUM_CARBONATE");
945
946//G4_MAGNESIUM_FLUORIDE  index=190
947  G4double M190[NDENSARRAY]={34.634,2.33,3.7105,0.1369,2.863,0.07934,3.6485,0,0.085};
948  AddMaterial(M190,"G4_MAGNESIUM_FLUORIDE");
949
950//G4_MAGNESIUM_OXIDE  index=191
951  G4double M191[NDENSARRAY]={38.407,2.412,3.6404,0.0575,2.858,0.08313,3.5968,0,0.055};
952  AddMaterial(M191,"G4_MAGNESIUM_OXIDE");
953
954//G4_MAGNESIUM_TETRABORATE  index=192
955  G4double M192[NDENSARRAY]={32.089,2.43,3.4328,0.1147,2.7635,0.09703,3.4893,0,0.044};
956  AddMaterial(M192,"G4_MAGNESIUM_TETRABORATE");
957
958//G4_MERCURIC_IODIDE  index=193
959  G4double M193[NDENSARRAY]={46.494,1.892,6.3787,0.104,3.4728,0.21513,2.7264,0,0.047};
960  AddMaterial(M193,"G4_MERCURIC_IODIDE");
961
962//G4_METHANE  index=194
963  G4double M194[NDENSARRAY]={0.588,1.662,9.5243,1.6263,3.9716,0.09253,3.6257,0,0.112};
964  AddMaterial(M194,"G4_METHANE");
965
966//G4_METHANOL  index=195
967  G4double M195[NDENSARRAY]={19.214,2.125,3.516,0.2529,2.7639,0.0897,3.5477,0,0.08};
968  AddMaterial(M195,"G4_METHANOL");
969
970//G4_MIX_D_WAX  index=196
971  G4double M196[NDENSARRAY]={21.547,1.905,3.078,0.1371,2.7145,0.0749,3.6823,0,0.047};
972  AddMaterial(M196,"G4_MIX_D_WAX");
973
974//G4_MS20_TISSUE  index=197
975  G4double M197[NDENSARRAY]={21.153,2.07,3.5241,0.1997,2.8033,0.08294,3.6061,0,0.053};
976  AddMaterial(M197,"G4_MS20_TISSUE");
977
978//G4_MUSCLE_SCELETAL_ICRP  index=198
979  G4double M198[NDENSARRAY]={21.781,2.185,3.4809,0.2282,2.7999,0.08636,3.533,0,0.089};
980  AddMaterial(M198,"G4_MUSCLE_SCELETAL_ICRP");
981
982//G4_MUSCLE_STRIATED_ICRU  index=199
983  G4double M199[NDENSARRAY]={21.795,2.174,3.4636,0.2249,2.8032,0.08507,3.5383,0,0.086};
984  AddMaterial(M199,"G4_MUSCLE_STRIATED_ICRU");
985
986//G4_MUSCLE_WITH_SUCROSE  index=200
987  G4double M200[NDENSARRAY]={22.48,2.169,3.391,0.2098,2.755,0.09481,3.4699,0,0.08};
988  AddMaterial(M200,"G4_MUSCLE_WITH_SUCROSE");
989
990//G4_MUSCLE_WITHOUT_SUCROSE  index=201
991  G4double M201[NDENSARRAY]={22.109,2.173,3.4216,0.2187,2.768,0.09143,3.4982,0,0.056};
992  AddMaterial(M201,"G4_MUSCLE_WITHOUT_SUCROSE");
993
994//G4_NAPHTALENE  index=202
995  G4double M202[NDENSARRAY]={22.459,1.956,3.2274,0.1374,2.5429,0.14766,3.2654,0,0.051};
996  AddMaterial(M202,"G4_NAPHTALENE");
997
998//G4_NITROBENZENE  index=203
999  G4double M203[NDENSARRAY]={22.747,2.065,3.4073,0.1777,2.663,0.12727,3.3091,0,0.051};
1000  AddMaterial(M203,"G4_NITROBENZENE");
1001
1002//G4_NITROUS_OXIDE  index=204
1003  G4double M204[NDENSARRAY]={0.872,2.059,10.1575,1.6477,4.1565,0.11992,3.3318,0,0.086};
1004  AddMaterial(M204,"G4_NITROUS_OXIDE");
1005
1006//G4_NYLON-8062  index=205
1007  G4double M205[NDENSARRAY]={22.221,1.967,3.125,0.1503,2.6004,0.11513,3.4044,0,0.054};
1008  AddMaterial(M205,"G4_NYLON-8062");
1009
1010//G4_NYLON-6/6  index=206
1011  G4double M206[NDENSARRAY]={22.774,1.931,3.0634,0.1336,2.5834,0.11818,3.3826,0,0.051};
1012  AddMaterial(M206,"G4_NYLON-6-6");
1013
1014//G4_NYLON-6/10  index=207
1015  G4double M207[NDENSARRAY]={22.866,1.942,3.0333,0.1304,2.5681,0.11852,3.3912,0,0.05};
1016  AddMaterial(M207,"G4_NYLON-6-10");
1017
1018//G4_NYLON-11_RILSAN  index=208
1019  G4double M208[NDENSARRAY]={25.661,1.902,2.7514,0.0678,2.4281,0.14868,3.2576,0,0.044};
1020  AddMaterial(M208,"G4_NYLON-11_RILSAN");
1021
1022//G4_OCTANE  index=209
1023  G4double M209[NDENSARRAY]={18.36,1.851,3.1834,0.1882,2.5664,0.11387,3.4776,0,0.057};
1024  AddMaterial(M209,"G4_OCTANE");
1025
1026//G4_PARAFFIN  index=210
1027  G4double M210[NDENSARRAY]={21.031,1.844,2.9551,0.1289,2.5084,0.12087,3.4288,0,0.052};
1028  AddMaterial(M210,"G4_PARAFFIN");
1029
1030//G4_N-PENTANE  index=211
1031  G4double M211[NDENSARRAY]={17.398,1.842,3.2504,0.2086,2.5855,0.10809,3.5265,0,0.064};
1032  AddMaterial(M211,"G4_N-PENTANE");
1033
1034//G4_PHOTO_EMULSION  index=212
1035  G4double M212[NDENSARRAY]={37.946,2.264,5.3319,0.1009,3.4866,0.12399,3.0094,0,0.028};
1036  AddMaterial(M212,"G4_PHOTO_EMULSION");
1037
1038//G4_PLASTIC_SC_VINYLTOLUENE  index=213
1039  G4double M213[NDENSARRAY]={21.54,1.929,3.1997,0.1464,2.4855,0.16101,3.2393,0,0.05};
1040  AddMaterial(M213,"G4_PLASTIC_SC_VINYLTOLUENE");
1041
1042//G4_PLUTONIUM_DIOXIDE  index=214
1043  G4double M214[NDENSARRAY]={62.143,1.846,5.9719,-0.2311,3.5554,0.20594,2.6522,0,0.111};
1044  AddMaterial(M214,"G4_PLUTONIUM_DIOXIDE");
1045
1046//G4_POLYACRYLONITRILE  index=215
1047  G4double M215[NDENSARRAY]={22.642,1.955,3.2459,0.1504,2.5159,0.16275,3.1975,0,0.05};
1048  AddMaterial(M215,"G4_POLYACRYLONITRILE");
1049
1050//G4_POLYCARBONATE  index=216
1051  G4double M216[NDENSARRAY]={22.915,2.06,3.3201,0.1606,2.6225,0.1286,3.3288,0,0.049};
1052  AddMaterial(M216,"G4_POLYCARBONATE");
1053
1054//G4_POLYCHLOROSTYRENE  index=217
1055  G4double M217[NDENSARRAY]={23.81,1.902,3.4659,0.1238,2.9241,0.0753,3.5441,0,0.029};
1056  AddMaterial(M217,"G4_POLYCHLOROSTYRENE");
1057
1058//G4_POLYETHYLENE  index=218
1059  G4double M218[NDENSARRAY]={21.099,1.882,3.0016,0.137,2.5177,0.12108,3.4292,0,0.051};
1060  AddMaterial(M218,"G4_POLYETHYLENE");
1061
1062//G4_MYLAR  index=219
1063  G4double M219[NDENSARRAY]={24.595,2.144,3.3262,0.1562,2.6507,0.12679,3.3076,0,0.052};
1064  AddMaterial(M219,"G4_MYLAR");
1065
1066//G4_LUCITE  index=220
1067  G4double M220[NDENSARRAY]={23.086,2.173,3.3297,0.1824,2.6681,0.11433,3.3836,0,0.056};
1068  AddMaterial(M220,"G4_LUCITE");
1069
1070//G4_POLYOXOMETHYLENE  index=221
1071  G4double M221[NDENSARRAY]={25.11,2.175,3.2514,0.1584,2.6838,0.10808,3.4002,0,0.063};
1072  AddMaterial(M221,"G4_POLYOXOMETHYLENE");
1073
1074//G4_POLYPROPYLENE  index=222
1075  G4double M222[NDENSARRAY]={20.457,1.884,3.1252,0.1534,2.4822,0.15045,3.2855,0,0.055};
1076  AddMaterial(M222,"G4_POLYPROPYLENE");
1077
1078//G4_POLYSTYRENE  index=223
1079  G4double M223[NDENSARRAY]={21.754,2.027,3.2999,0.1647,2.5031,0.16454,3.2224,0,0.051};
1080  AddMaterial(M223,"G4_POLYSTYRENE");
1081
1082//G4_TEFLON  index=224
1083  G4double M224[NDENSARRAY]={29.609,2.142,3.4161,0.1648,2.7404,0.10606,3.4046,0,0.073};
1084  AddMaterial(M224,"G4_TEFLON");
1085
1086//G4_POLYTRIFLUOROCHLOROETHYLENE  index=225
1087  G4double M225[NDENSARRAY]={28.955,2.094,3.8551,0.1714,3.0265,0.07727,3.5085,0,0.035};
1088  AddMaterial(M225,"G4_POLYTRIFLUOROCHLOROETHYLENE");
1089
1090//G4_POLYVINYL_ACETATE  index=226
1091  G4double M226[NDENSARRAY]={22.978,2.116,3.3309,0.1769,2.6747,0.11442,3.3762,0,0.055};
1092  AddMaterial(M226,"G4_POLYVINYL_ACETATE");
1093
1094//G4_PLOYVINYL_ALCOHOL  index=227
1095  G4double M227[NDENSARRAY]={24.251,2.071,3.1115,0.1401,2.6315,0.11178,3.3893,0,0.056};
1096  AddMaterial(M227,"G4_PLOYVINYL_ALCOHOL");
1097
1098//G4_POLYVINYL_BUTYRAL  index=228
1099  G4double M228[NDENSARRAY]={22.521,2.021,3.1865,0.1555,2.6186,0.11544,3.3983,0,0.054};
1100  AddMaterial(M228,"G4_POLYVINYL_BUTYRAL");
1101
1102//G4_POLYVINYL_CHLORIDE  index=229
1103  G4double M229[NDENSARRAY]={23.51,1.84,4.0532,0.1559,2.9415,0.12438,3.2104,0,0.027};
1104  AddMaterial(M229,"G4_POLYVINYL_CHLORIDE");
1105
1106//G4_POLYVINYLIDENE_CHLORIDE  index=230
1107  G4double M230[NDENSARRAY]={26.437,1.814,4.2506,0.1314,2.9009,0.15466,3.102,0,0.034};
1108  AddMaterial(M230,"G4_POLYVINYLIDENE_CHLORIDE");
1109
1110//G4_POLYVINYLIDENE_FLUORIDE  index=231
1111  G4double M231[NDENSARRAY]={27.024,2.16,3.3793,0.1717,2.7375,0.10316,3.42,0,0.067};
1112  AddMaterial(M231,"G4_POLYVINYLIDENE_FLUORIDE");
1113
1114//G4_POLYVINYL_PYRROLIDONE  index=232
1115  G4double M232[NDENSARRAY]={23.671,1.989,3.1017,0.1324,2.5867,0.12504,3.3326,0,0.031};
1116  AddMaterial(M232,"G4_POLYVINYL_PYRROLIDONE");
1117
1118//G4_POTASSIUM_IODIDE  index=233
1119  G4double M233[NDENSARRAY]={33.575,1.784,6.1088,0.1044,3.3442,0.22053,2.7558,0,0.042};
1120  AddMaterial(M233,"G4_POTASSIUM_IODIDE");
1121
1122//G4_POTASSIUM_OXIDE  index=234
1123  G4double M234[NDENSARRAY]={30.672,2.065,4.6463,0.048,3.011,0.16789,3.0121,0,0.027};
1124  AddMaterial(M234,"G4_POTASSIUM_OXIDE");
1125
1126//G4_PROPANE  index=235
1127  G4double M235[NDENSARRAY]={0.959,1.708,8.7878,1.4326,3.7998,0.09916,3.592,0,0.093};
1128  AddMaterial(M235,"G4_PROPANE");
1129
1130//G4_lPROPANE  index=236
1131  G4double M236[NDENSARRAY]={14.509,1.844,3.5529,0.2861,2.6568,0.10329,3.562,0,0.068};
1132  AddMaterial(M236,"G4_lPROPANE");
1133
1134//G4_N-PROPYL_ALCOHOL index=237
1135  G4double M237[NDENSARRAY]={19.429,1.972,3.2915,0.2046,2.6681,0.09644,3.5415,0,0.07};
1136  AddMaterial(M237,"N-PROPYL_ALCOHOL");
1137
1138//G4_PYRIDINE  index=238
1139  G4double M238[NDENSARRAY]={20.807,1.895,3.3148,0.167,2.5245,0.16399,3.1977,0,0.051};
1140  AddMaterial(M238,"G4_PYRIDINE");
1141
1142//G4_RUBBER_BUTYL  index=239
1143  G4double M239[NDENSARRAY]={20.873,1.852,2.9915,0.1347,2.5154,0.12108,3.4296,0,0.051};
1144  AddMaterial(M239,"G4_RUBBER_BUTYL");
1145
1146//G4_RUBBER_NATURAL  index=240
1147  G4double M240[NDENSARRAY]={20.644,1.889,3.1272,0.1512,2.4815,0.15058,3.2879,0,0.053};
1148  AddMaterial(M240,"G4_RUBBER_NATURAL");
1149
1150//G4_RUBBER_NEOPRENE  index=241
1151  G4double M241[NDENSARRAY]={23.036,1.874,3.7911,0.1501,2.9461,0.09763,3.3632,0,0.026};
1152  AddMaterial(M241,"G4_RUBBER_NEOPRENE");
1153
1154//G4_SILICON_DIOXIDE  index=242
1155  G4double M242[NDENSARRAY]={31.014,2.335,4.0029,0.1385,3.0025,0.08408,3.5064,0,0.018};
1156  AddMaterial(M242,"G4_SILICON_DIOXIDE");
1157
1158//G4_SILVER_BROMIDE  index=243
1159  G4double M243[NDENSARRAY]={48.448,2.271,5.6139,0.0352,3.2109,0.24582,2.682,0,0.043};
1160  AddMaterial(M243,"G4_SILVER_BROMIDE");
1161
1162//G4_SILVER_CHLORIDE  index=244
1163  G4double M244[NDENSARRAY]={45.405,2.096,5.3437,-0.0139,3.2022,0.22968,2.7041,0,0.062};
1164  AddMaterial(M244,"G4_SILVER_CHLORIDE");
1165
1166//G4_SILVER_HALIDES  index=245
1167  G4double M245[NDENSARRAY]={48.433,2.27,5.6166,0.0353,3.2117,0.24593,2.6814,0,0.043};
1168  AddMaterial(M245,"G4_SILVER_HALIDES");
1169
1170//G4_SILVER_IODIDE  index=246
1171  G4double M246[NDENSARRAY]={46.105,1.945,5.9342,0.0148,3.2908,0.25059,2.6572,0,0.071};
1172  AddMaterial(M246,"G4_SILVER_IODIDE");
1173
1174//G4_SKIN_ICRP  index=247
1175  G4double M247[NDENSARRAY]={22.4,2.14,3.3546,0.2019,2.7526,0.09459,3.4643,0,0.076};
1176  AddMaterial(M247,"G4_SKIN_ICRP");
1177
1178//G4_SODIUM_CARBONATE  index=248
1179  G4double M248[NDENSARRAY]={32.117,2.557,3.7178,0.1287,2.8591,0.08715,3.5638,0,0.074};
1180  AddMaterial(M248,"G4_SODIUM_CARBONATE");
1181
1182//G4_SODIUM_IODIDE  index=249
1183  G4double M249[NDENSARRAY]={36.057,1.857,6.0572,0.1203,3.592,0.12516,3.0398,0,0.031};
1184  AddMaterial(M249,"G4_SODIUM_IODIDE");
1185
1186//G4_SODIUM_MONOXIDE  index=250
1187  G4double M250[NDENSARRAY]={30.205,2.689,4.1892,0.1652,2.9793,0.07501,3.6943,0,0.097};
1188  AddMaterial(M250,"G4_SODIUM_MONOXIDE");
1189
1190//G4_SODIUM_NITRATE  index=251
1191  G4double M251[NDENSARRAY]={30.459,2.456,3.6502,0.1534,2.8221,0.09391,3.5097,0,0.081};
1192  AddMaterial(M251,"G4_SODIUM_NITRATE");
1193
1194//G4_STILBENE  index=252
1195  G4double M252[NDENSARRAY]={20.719,1.963,3.368,0.1734,2.5142,0.16659,3.2168,0,0.052};
1196  AddMaterial(M252,"G4_STILBENE");
1197
1198//G4_SUCROSE  index=253
1199  G4double M253[NDENSARRAY]={26.416,2.167,3.1526,0.1341,2.6558,0.11301,3.363,0,0.057};
1200  AddMaterial(M253,"G4_SUCROSE");
1201
1202//G4_TERPHENYL  index=254
1203  G4double M254[NDENSARRAY]={23.116,1.976,3.2639,0.1322,2.5429,0.14964,3.2685,0,0.043};
1204  AddMaterial(M254,"G4_TERPHENYL");
1205
1206//G4_TESTES_ICRP  index=255
1207  G4double M255[NDENSARRAY]={21.815,2.185,3.4698,0.2274,2.7988,0.08533,3.5428,0,0.091};
1208  AddMaterial(M255,"G4_TESTES_ICRP");
1209
1210//G4_TETRACHLOROETHYLENE  index=256
1211  G4double M256[NDENSARRAY]={25.513,1.79,4.6619,0.1713,2.9083,0.18595,3.0156,0,0.038};
1212  AddMaterial(M256,"G4_TETRACHLOROETHYLENE");
1213
1214//G4_THALIUM_CHLORIDE  index=257
1215  G4double M257[NDENSARRAY]={48.749,1.997,6.3009,0.0705,3.5716,0.18599,2.769,0,0.04};
1216  AddMaterial(M257,"G4_THALIUM_CHLORIDE");
1217
1218//G4_TISSUE_SOFT_ICRP  index=258
1219  G4double M258[NDENSARRAY]={21.394,2.144,3.4354,0.2211,2.7799,0.08926,3.511,0,0.077};
1220  AddMaterial(M258,"G4_TISSUE_SOFT_ICRP");
1221
1222//G4_TISSUE_SOFT_ICRU-4  index=259
1223  G4double M259[NDENSARRAY]={21.366,2.192,3.5087,0.2377,2.7908,0.09629,3.4371,0,0.092};
1224  AddMaterial(M259,"G4_TISSUE_SOFT_ICRU-4");
1225
1226//G4_TISSUE-METHANE  index=260
1227  G4double M260[NDENSARRAY]={0.697,1.89,9.95,1.6442,4.1399,0.09946,3.4708,0,0.098};
1228  AddMaterial(M260,"G4_TISSUE-METHANE");
1229
1230//G4_TISSUE-PROPANE  index=261
1231  G4double M261[NDENSARRAY]={0.913,1.856,9.3529,1.5139,3.9916,0.09802,3.5159,0,0.092};
1232  AddMaterial(M261,"G4_TISSUE-PROPANE");
1233
1234//G4_TITANIUM_DIOXIDE  index=262
1235  G4double M262[NDENSARRAY]={41.022,2.307,3.9522,-0.0119,3.1647,0.08569,3.3267,0,0.027};
1236  AddMaterial(M262,"G4_TITANIUM_DIOXIDE");
1237
1238//G4_TOLUENE  index=263
1239  G4double M263[NDENSARRAY]={19.764,1.88,3.3026,0.1722,2.5728,0.13284,3.3558,0,0.052};
1240  AddMaterial(M263,"G4_TOLUENE");
1241
1242//G4_TRICHLOROETHYLENE  index=264
1243  G4double M264[NDENSARRAY]={24.301,1.789,4.6148,0.1803,2.914,0.18272,3.0137,0,0.036};
1244  AddMaterial(M264,"G4_TRICHLOROETHYLENE");
1245
1246//G4_TRIETHYL_PHOSPHATE  index=265
1247  G4double M265[NDENSARRAY]={21.863,2.1,3.6242,0.2054,2.9428,0.06922,3.6302,0,0.049};
1248  AddMaterial(M265,"G4_TRIETHYL_PHOSPHATE");
1249
1250//G4_TUNGSTEN_HEXAFLUORIDE  index=266
1251  G4double M266[NDENSARRAY]={29.265,2.325,5.9881,0.302,4.2602,0.03658,3.5134,0,0.055};
1252  AddMaterial(M266,"G4_TUNGSTEN_HEXAFLUORIDE");
1253
1254//G4_URANIUM_DICARBIDE  index=267
1255  G4double M267[NDENSARRAY]={60.969,1.703,6.0247,-0.2191,3.5208,0.2112,2.6577,0,0.12};
1256  AddMaterial(M267,"G4_URANIUM_DICARBIDE");
1257
1258//G4_URANIUM_MONOCARBIDE  index=268
1259  G4double M268[NDENSARRAY]={66.602,1.68,6.121,-0.2524,3.4941,0.22972,2.6169,0,0.132};
1260  AddMaterial(M268,"G4_URANIUM_MONOCARBIDE");
1261
1262//G4_URANIUM_OXIDE  index=269
1263  G4double M269[NDENSARRAY]={60.332,1.76,5.9605,-0.1938,3.5292,0.20463,2.6711,0,0.098};
1264  AddMaterial(M269,"G4_URANIUM_OXIDE");
1265
1266//G4_UREA  index=270
1267  G4double M270[NDENSARRAY]={24.194,2.022,3.2032,0.1603,2.6225,0.11609,3.3461,0,0.06};
1268  AddMaterial(M270,"G4_UREA");
1269
1270//G4_VALINE  index=271
1271  G4double M271[NDENSARRAY]={23.622,2.024,3.1059,0.1441,2.6227,0.11386,3.3774,0,0.056};
1272  AddMaterial(M271,"G4_VALINE");
1273
1274//G4_VITON  index=272
1275  G4double M272[NDENSARRAY]={26.948,2.227,3.5943,0.2106,2.7874,0.09965,3.4556,0,0.07};
1276  AddMaterial(M272,"G4_VITON");
1277
1278//G4_WATER  index=273
1279  G4double M273[NDENSARRAY]={21.469,2.203,3.5017,0.24,2.8004,0.09116,3.4773,0,0.097};
1280  AddMaterial(M273,"G4_WATER");
1281
1282//G4_WATER_VAPOR  index=274
1283  G4double M274[NDENSARRAY]={0.59,2.175,10.5962,1.7952,4.3437,0.08101,3.5901,0,0.121};
1284  AddMaterial(M274,"G4_WATER_VAPOR");
1285
1286//G4_XYLENE  index=275
1287  G4double M275[NDENSARRAY]={19.866,1.882,3.2698,0.1695,2.5675,0.13216,3.3564,0,0.051};
1288  AddMaterial(M275,"G4_XYLENE");
1289
1290//G4_GRAPHITE  index=276
1291  G4double M276[NDENSARRAY]={30.652,2.29,2.868,-0.0178,2.3415,0.26142,2.8697,0.12,0.038};
1292  AddMaterial(M276,"G4_GRAPHITE");
1293
1294//G4_GRAPHITE_POROUS  index=277
1295  G4double  M277[NDENSARRAY]={26.555,2.49,3.155,0.048,2.5387,0.20762,2.9532,0.14,0.038};
1296  AddMaterial(M277,"G4_GRAPHITE_POROUS");
1297
1298}
1299
1300G4int G4DensityEffectData::GetElementIndex(G4int Z, G4State st)
1301{
1302  G4int idx = -1;
1303  if(Z > 0 && Z < NDENSELEM) { 
1304    if(st == state[Z])             { idx = indexZ[Z]; } 
1305    else if(st == kStateUndefined) { idx = indexZ[Z]; }
1306  } 
1307  return idx;
1308}
1309
1310G4int G4DensityEffectData::GetIndex(const G4String& matName)
1311{
1312  G4int idx = -1;
1313
1314  for (G4int i=0; i<NDENSDATA; ++i) {
1315    if ( names[i] == matName ) {
1316      idx = i;
1317      break;
1318    }
1319  }
1320
1321  return idx;
1322}
1323
1324void G4DensityEffectData::AddMaterial(G4double* val, const G4String& matName)
1325{
1326  data[index][0] = val[0]*eV; 
1327  for(G4int i=1; i<NDENSARRAY; ++i) { 
1328    data[index][i] = val[i]; 
1329  } 
1330  names.push_back( matName );
1331  ++index;
1332}
1333
1334void G4DensityEffectData::PrintData(const G4String& matName)
1335{
1336  if("" == matName || "all" == matName) {
1337    DumpData();
1338    return;
1339  }
1340  G4int idx = GetIndex(matName);
1341  if(idx >= 0) {
1342    G4cout << "G4DensityEffectData for <" << matName
1343           << "> index= " << idx << G4endl;
1344    G4cout << "Eplasma(eV)= " << data[idx][0]/eV
1345           << " rho= " << data[idx][1]
1346           << " -C= " << data[idx][2] 
1347           << " x0= " << data[idx][3]
1348           << " x1= " << data[idx][4]
1349           << " a= " << data[idx][5]
1350           << " m= " << data[idx][6] 
1351           << " d0= " << data[idx][7] 
1352           << " err= " << data[idx][8] 
1353           << G4endl;
1354  } else {
1355    G4cout << "G4DensityEffectData does not have <" << matName << ">" << G4endl; 
1356  }
1357}
1358
1359void G4DensityEffectData::DumpData()
1360{
1361  G4cout << "======================================================================"
1362         << G4endl; 
1363  G4cout << "     Material        Eplasma(eV)  rho  -C   x0   x1   a   m  d0  err"
1364         << G4endl;       
1365  G4cout << "======================================================================"
1366         << G4endl; 
1367  for(G4int i=0; i<NDENSDATA; ++i) {
1368    G4cout << std::setw(3) << i << ". " << std::setw(25) << names[i]
1369           << std::setw(8) << data[i][0]/eV;
1370    for(G4int j=1; j<NDENSARRAY; ++j) { G4cout << std::setw(8) << data[i][j]; }
1371    G4cout << G4endl;
1372  }
1373  G4cout << "======================================================================"
1374         << G4endl; 
1375}
Note: See TracBrowser for help on using the repository browser.