Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/materials/src/G4DensityEffectData.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /materials/src/G4DensityEffectData.cc (Version 11.3.0) and /materials/src/G4DensityEffectData.cc (Version 9.0.p2)


  1 //                                                  1 
  2 // *******************************************    
  3 // * License and Disclaimer                       
  4 // *                                              
  5 // * The  Geant4 software  is  copyright of th    
  6 // * the Geant4 Collaboration.  It is provided    
  7 // * conditions of the Geant4 Software License    
  8 // * LICENSE and available at  http://cern.ch/    
  9 // * include a list of copyright holders.         
 10 // *                                              
 11 // * Neither the authors of this software syst    
 12 // * institutes,nor the agencies providing fin    
 13 // * work  make  any representation or  warran    
 14 // * regarding  this  software system or assum    
 15 // * use.  Please see the license in the file     
 16 // * for the full disclaimer and the limitatio    
 17 // *                                              
 18 // * This  code  implementation is the result     
 19 // * technical work of the GEANT4 collaboratio    
 20 // * By using,  copying,  modifying or  distri    
 21 // * any work based  on the software)  you  ag    
 22 // * use  in  resulting  scientific  publicati    
 23 // * acceptance of all terms of the Geant4 Sof    
 24 // *******************************************    
 25 //                                                
 26 //--------------------------------------------    
 27 //                                                
 28 // GEANT4 Class file                              
 29 //                                                
 30 // Description: Data on density effect            
 31 //                                                
 32 // Authors:    A.Bagulya, A.Ivanchenko 28.10.2    
 33 //                                                
 34 //--------------------------------------------    
 35                                                   
 36 #include "G4DensityEffectData.hh"                 
 37                                                   
 38 #include "G4SystemOfUnits.hh"                     
 39                                                   
 40 #include <iomanip>                                
 41                                                   
 42 G4DensityEffectData::G4DensityEffectData() { I    
 43                                                   
 44 //....oooOO0OOooo........oooOO0OOooo........oo    
 45                                                   
 46 void G4DensityEffectData::Initialize()            
 47 {                                                 
 48   // R.M. Sternheimer et al. Density Effect fo    
 49   // Particles in Various Substances. Atom. Da    
 50   //                                              
 51   // Data[10]:  Eplasma rho -C X_0 X_1 a m del    
 52   // Eplasma - Plasma energy (in eV)              
 53   // rho     - Sternheimer adjustment factor f    
 54   // -C                                           
 55   // X_0, X_1, m, a - parameters in fitting fo    
 56   // delta_0 - Density-effect value delta(X_0)    
 57   // DELTA_{max} - Upper bound for the error i    
 58   // I - Mean ionisation potential in keV used    
 59                                                   
 60   for (G4int i = 0; i < NDENSELEM; ++i) {         
 61     indexZ[i] = -1;                               
 62     state[i] = kStateSolid;                       
 63   }                                               
 64                                                   
 65   // G4_lH2  index=0                              
 66   G4double M1[NDENSARRAY] = {                     
 67     7.031, 1.546, 3.2632, 0.4759, 1.9215, 0.13    
 68   AddMaterial(M1, "G4_lH2");                      
 69                                                   
 70   // G4_H  index=1                                
 71   G4double M0[NDENSARRAY] = {                     
 72     0.263, 1.412, 9.5835, 1.8639, 3.2718, 0.14    
 73   AddMaterial(M0, "G4_H");                        
 74   indexZ[1] = 1;                                  
 75   state[1] = kStateGas;                           
 76                                                   
 77   // G4_He  index=2                               
 78   G4double M2[NDENSARRAY] = {0.263, 1.7, 11.13    
 79   AddMaterial(M2, "G4_He");                       
 80   indexZ[2] = 2;                                  
 81   state[2] = kStateGas;                           
 82                                                   
 83   // G4_Li  index=3                               
 84   G4double M3[NDENSARRAY] = {                     
 85     13.844, 1.535, 3.1221, 0.1304, 1.6397, 0.9    
 86   AddMaterial(M3, "G4_Li");                       
 87   indexZ[3] = 3;                                  
 88                                                   
 89   // G4_Be  index=4                               
 90   G4double M4[NDENSARRAY] = {                     
 91     26.098, 1.908, 2.7847, 0.0392, 1.6922, 0.8    
 92   AddMaterial(M4, "G4_Be");                       
 93   indexZ[4] = 4;                                  
 94                                                   
 95   // G4_B  index=5                                
 96   G4double M5[NDENSARRAY] = {                     
 97     30.17, 2.32, 2.8477, 0.0305, 1.9688, 0.562    
 98   AddMaterial(M5, "G4_B");                        
 99   indexZ[5] = 5;                                  
100                                                   
101   // G4_C  index=6                                
102   G4double M6[NDENSARRAY] = {                     
103     28.803, 2.376, 2.9925, -0.0351, 2.486, 0.2    
104   AddMaterial(M6, "G4_C");                        
105   indexZ[6] = 6;                                  
106                                                   
107   // G4_N  index=7                                
108   G4double M7[NDENSARRAY] = {                     
109     0.695, 1.984, 10.5400, 1.7378, 4.1323, 0.1    
110   AddMaterial(M7, "G4_N");                        
111   indexZ[7] = 7;                                  
112   state[7] = kStateGas;                           
113                                                   
114   // G4_O  index=8                                
115   G4double M8[NDENSARRAY] = {                     
116     0.744, 2.314, 10.7004, 1.7541, 4.3213, 0.1    
117   AddMaterial(M8, "G4_O");                        
118   indexZ[8] = 8;                                  
119   state[8] = kStateGas;                           
120                                                   
121   // G4_F  index=9                                
122   G4double M9[NDENSARRAY] = {                     
123     0.788, 2.450, 10.9653, 1.8433, 4.4096, 0.1    
124   AddMaterial(M9, "G4_F");                        
125   indexZ[9] = 9;                                  
126   state[9] = kStateGas;                           
127                                                   
128   // G4_Ne  index=10                              
129   G4double M10[NDENSARRAY] = {                    
130     0.587, 2.577, 11.9041, 2.0735, 4.6421, 0.0    
131   AddMaterial(M10, "G4_Ne");                      
132   indexZ[10] = 10;                                
133   state[10] = kStateGas;                          
134                                                   
135   // G4_Na  index=11                              
136   G4double M11[NDENSARRAY] = {                    
137     19.641, 2.648, 5.0526, 0.2880, 3.1962, 0.0    
138   AddMaterial(M11, "G4_Na");                      
139   indexZ[11] = 11;                                
140                                                   
141   // G4_Mg  index=12                              
142   G4double M12[NDENSARRAY] = {                    
143     26.708, 2.331, 4.5297, 0.1499, 3.0668, 0.0    
144   AddMaterial(M12, "G4_Mg");                      
145   indexZ[12] = 12;                                
146                                                   
147   // G4_Al  index=13                              
148   G4double M13[NDENSARRAY] = {                    
149     32.86, 2.18, 4.2395, 0.1708, 3.0127, 0.080    
150   AddMaterial(M13, "G4_Al");                      
151   indexZ[13] = 13;                                
152                                                   
153   // G4_Si  index=14                              
154   G4double M14[NDENSARRAY] = {                    
155     31.055, 2.103, 4.4351, 0.2014, 2.8715, 0.1    
156   AddMaterial(M14, "G4_Si");                      
157   indexZ[14] = 14;                                
158                                                   
159   // G4_P  index=15                               
160   G4double M15[NDENSARRAY] = {                    
161     29.743, 2.056, 4.5214, 0.1696, 2.7815, 0.2    
162   AddMaterial(M15, "G4_P");                       
163   indexZ[15] = 15;                                
164                                                   
165   // G4_S  index=16                               
166   G4double M16[NDENSARRAY] = {                    
167     28.789, 2.131, 4.6659, 0.158, 2.7159, 0.33    
168   AddMaterial(M16, "G4_S");                       
169   indexZ[16] = 16;                                
170                                                   
171   // G4_Cl  index=17                              
172   G4double M17[NDENSARRAY] = {                    
173     1.092, 1.734, 11.1421, 1.5555, 4.2994, 0.1    
174   AddMaterial(M17, "G4_Cl");                      
175   indexZ[17] = 17;                                
176   state[17] = kStateGas;                          
177                                                   
178   // G4_Ar  index=18                              
179   G4double M18[NDENSARRAY] = {                    
180     0.789, 1.753, 11.9480, 1.7635, 4.4855, 0.1    
181   AddMaterial(M18, "G4_Ar");                      
182   indexZ[18] = 18;                                
183   state[18] = kStateGas;                          
184                                                   
185   // G4_K  index=19                               
186   G4double M19[NDENSARRAY] = {                    
187     18.65, 1.830, 5.6423, 0.3851, 3.1724, 0.19    
188   AddMaterial(M19, "G4_K");                       
189   indexZ[19] = 19;                                
190                                                   
191   // G4_Ca  index=20                              
192   G4double M20[NDENSARRAY] = {                    
193     25.342, 1.666, 5.0396, 0.3228, 3.1191, 0.1    
194   AddMaterial(M20, "G4_Ca");                      
195   indexZ[20] = 20;                                
196                                                   
197   // G4_Sc  index=21                              
198   G4double M21[NDENSARRAY] = {                    
199     34.050, 1.826, 4.6949, 0.1640, 3.0593, 0.1    
200   AddMaterial(M21, "G4_Sc");                      
201   indexZ[21] = 21;                                
202                                                   
203   // G4_Ti  index=22                              
204   G4double M22[NDENSARRAY] = {                    
205     41.619, 1.969, 4.4450, 0.0957, 3.0386, 0.1    
206   AddMaterial(M22, "G4_Ti");                      
207   indexZ[22] = 22;                                
208                                                   
209   // G4_V  index=23                               
210   G4double M23[NDENSARRAY] = {                    
211     47.861, 2.070, 4.2659, 0.0691, 3.0322, 0.1    
212   AddMaterial(M23, "G4_V");                       
213   indexZ[23] = 23;                                
214                                                   
215   // G4_Cr  index=24                              
216   G4double M24[NDENSARRAY] = {                    
217     52.458, 2.181, 4.1781, 0.0340, 3.0451, 0.1    
218   AddMaterial(M24, "G4_Cr");                      
219   indexZ[24] = 24;                                
220                                                   
221   // G4_Mn  index=25                              
222   G4double M25[NDENSARRAY] = {                    
223     53.022, 2.347, 4.2702, 0.0447, 3.1074, 0.1    
224   AddMaterial(M25, "G4_Mn");                      
225   indexZ[25] = 25;                                
226                                                   
227   // G4_Fe  index=26                              
228   G4double M26[NDENSARRAY] = {                    
229     55.172, 2.504, 4.2911, -0.0012, 3.1531, 0.    
230   AddMaterial(M26, "G4_Fe");                      
231   indexZ[26] = 26;                                
232                                                   
233   // G4_Co  index=27                              
234   G4double M27[NDENSARRAY] = {                    
235     58.188, 2.626, 4.2601, -0.0187, 3.1790, 0.    
236   AddMaterial(M27, "G4_Co");                      
237   indexZ[27] = 27;                                
238                                                   
239   // G4_Ni  index=28                              
240   G4double M28[NDENSARRAY] = {                    
241     59.385, 2.889, 4.3115, -0.0566, 3.1851, 0.    
242   AddMaterial(M28, "G4_Ni");                      
243   indexZ[28] = 28;                                
244                                                   
245   // G4_Cu  index=29                              
246   G4double M29[NDENSARRAY] = {                    
247     58.270, 2.956, 4.4190, -0.0254, 3.2792, 0.    
248   AddMaterial(M29, "G4_Cu");                      
249   indexZ[29] = 29;                                
250                                                   
251   // G4_Zn  index=30                              
252   G4double M30[NDENSARRAY] = {                    
253     52.132, 3.142, 4.6906, 0.0049, 3.3668, 0.1    
254   AddMaterial(M30, "G4_Zn");                      
255   indexZ[30] = 30;                                
256                                                   
257   // G4_Ga  index=31                              
258   G4double M31[NDENSARRAY] = {                    
259     46.688, 2.747, 4.9353, 0.2267, 3.5434, 0.0    
260   AddMaterial(M31, "G4_Ga");                      
261   indexZ[31] = 31;                                
262                                                   
263   // G4_Ge  index=32                              
264   G4double M32[NDENSARRAY] = {                    
265     44.141, 2.461, 5.1411, 0.3376, 3.6096, 0.0    
266   AddMaterial(M32, "G4_Ge");                      
267   indexZ[32] = 32;                                
268                                                   
269   // G4_As  index=33                              
270   G4double M33[NDENSARRAY] = {                    
271     45.779, 2.219, 5.0510, 0.1767, 3.5702, 0.0    
272   AddMaterial(M33, "G4_As");                      
273   indexZ[33] = 33;                                
274                                                   
275   // G4_Se  index=34                              
276   G4double M34[NDENSARRAY] = {                    
277     40.112, 2.104, 5.3210, 0.2258, 3.6264, 0.0    
278   AddMaterial(M34, "G4_Se");                      
279   indexZ[34] = 34;                                
280                                                   
281   // G4_Br  index=35                              
282   G4double M35[NDENSARRAY] = {                    
283     1.604, 1.845, 11.7307, 1.5262, 4.9899, 0.0    
284   AddMaterial(M35, "G4_Br");                      
285   indexZ[35] = 35;                                
286   state[35] = kStateGas;                          
287                                                   
288   // G4_Kr  index=36                              
289   G4double M36[NDENSARRAY] = {                    
290     1.114, 1.77, 12.5115, 1.7158, 5.0748, 0.07    
291   AddMaterial(M36, "G4_Kr");                      
292   indexZ[36] = 36;                                
293   state[36] = kStateGas;                          
294                                                   
295   // G4_Rb  index=37                              
296   G4double M37[NDENSARRAY] = {                    
297     23.467, 1.823, 6.4776, 0.5737, 3.7995, 0.0    
298   AddMaterial(M37, "G4_Rb");                      
299   indexZ[37] = 37;                                
300                                                   
301   // G4_Sr  index=38                              
302   G4double M38[NDENSARRAY] = {                    
303     30.244, 1.707, 5.9867, 0.4585, 3.6778, 0.0    
304   AddMaterial(M38, "G4_Sr");                      
305   indexZ[38] = 38;                                
306                                                   
307   // G4_Y  index=39                               
308   G4double M39[NDENSARRAY] = {                    
309     40.346, 1.649, 5.4801, 0.3608, 3.5542, 0.0    
310   AddMaterial(M39, "G4_Y");                       
311   indexZ[39] = 39;                                
312                                                   
313   // G4_Zr  index=40                              
314   G4double M40[NDENSARRAY] = {                    
315     48.671, 1.638, 5.1774, 0.2957, 3.489, 0.07    
316   AddMaterial(M40, "G4_Zr");                      
317   indexZ[40] = 40;                                
318                                                   
319   // G4_Nb  index=41                              
320   G4double M41[NDENSARRAY] = {                    
321     56.039, 1.734, 5.0141, 0.1785, 3.2201, 0.1    
322   AddMaterial(M41, "G4_Nb");                      
323   indexZ[41] = 41;                                
324                                                   
325   // G4_Mo  index=42                              
326   G4double M42[NDENSARRAY] = {                    
327     60.951, 1.658, 4.8793, 0.2267, 3.2784, 0.1    
328   AddMaterial(M42, "G4_Mo");                      
329   indexZ[42] = 42;                                
330                                                   
331   // G4_Tc  index=43                              
332   G4double M43[NDENSARRAY] = {                    
333     64.760, 1.727, 4.7769, 0.0949, 3.1253, 0.1    
334   AddMaterial(M43, "G4_Tc");                      
335   indexZ[43] = 43;                                
336                                                   
337   // G4_Ru  index=44                              
338   G4double M44[NDENSARRAY] = {                    
339     66.978, 1.780, 4.7694, 0.0599, 3.0834, 0.1    
340   AddMaterial(M44, "G4_Ru");                      
341   indexZ[44] = 44;                                
342                                                   
343   // G4_Rh  index=45                              
344   G4double M45[NDENSARRAY] = {                    
345     67.128, 1.804, 4.8008, 0.0576, 3.1069, 0.1    
346   AddMaterial(M45, "G4_Rh");                      
347   indexZ[45] = 45;                                
348                                                   
349   // G4_Pd  index=46                              
350   G4double M46[NDENSARRAY] = {                    
351     65.683, 1.911, 4.9358, 0.0563, 3.0555, 0.2    
352   AddMaterial(M46, "G4_Pd");                      
353   indexZ[46] = 46;                                
354                                                   
355   // G4_Ag  index=47                              
356   G4double M47[NDENSARRAY] = {                    
357     61.635, 1.933, 5.0630, 0.0657, 3.1074, 0.2    
358   AddMaterial(M47, "G4_Ag");                      
359   indexZ[47] = 47;                                
360                                                   
361   // G4_Cd  index=48                              
362   G4double M48[NDENSARRAY] = {                    
363     55.381, 1.895, 5.2727, 0.1281, 3.1667, 0.2    
364   AddMaterial(M48, "G4_Cd");                      
365   indexZ[48] = 48;                                
366                                                   
367   // G4_In  index=49                              
368   G4double M49[NDENSARRAY] = {                    
369     50.896, 1.851, 5.5211, 0.2406, 3.2032, 0.2    
370   AddMaterial(M49, "G4_In");                      
371   indexZ[49] = 49;                                
372                                                   
373   // G4_Sn  index=50                              
374   G4double M50[NDENSARRAY] = {                    
375     50.567, 1.732, 5.5340, 0.2879, 3.2959, 0.1    
376   AddMaterial(M50, "G4_Sn");                      
377   indexZ[50] = 50;                                
378                                                   
379   // G4_Sb  index=51                              
380   G4double M51[NDENSARRAY] = {                    
381     48.242, 1.645, 5.6241, 0.3189, 3.3489, 0.1    
382   AddMaterial(M51, "G4_Sb");                      
383   indexZ[51] = 51;                                
384                                                   
385   // G4_Te  index=52                              
386   G4double M52[NDENSARRAY] = {                    
387     45.952, 1.577, 5.7131, 0.3296, 3.4418, 0.1    
388   AddMaterial(M52, "G4_Te");                      
389   indexZ[52] = 52;                                
390                                                   
391   // G4_I  index=53                               
392   G4double M53[NDENSARRAY] = {                    
393     41.348, 1.498, 5.9488, 0.0549, 3.2596, 0.2    
394   AddMaterial(M53, "G4_I");                       
395   indexZ[53] = 53;                                
396                                                   
397   // G4_Xe  index=54                              
398   G4double M54[NDENSARRAY] = {                    
399     1.369, 1.435, 12.7281, 1.563, 4.7371, 0.23    
400   AddMaterial(M54, "G4_Xe");                      
401   indexZ[54] = 54;                                
402   state[54] = kStateGas;                          
403                                                   
404   // G4_Cs  index=55                              
405   G4double M55[NDENSARRAY] = {                    
406     25.37, 1.462, 6.9135, 0.5473, 3.5914, 0.18    
407   AddMaterial(M55, "G4_Cs");                      
408   indexZ[55] = 55;                                
409                                                   
410   // G4_Ba  index=56                              
411   G4double M56[NDENSARRAY] = {                    
412     34.425, 1.410, 6.3153, 0.4190, 3.4547, 0.1    
413   AddMaterial(M56, "G4_Ba");                      
414   indexZ[56] = 56;                                
415                                                   
416   // G4_La  index=57                              
417   G4double M57[NDENSARRAY] = {                    
418     45.792, 1.392, 5.7850, 0.3161, 3.3293, 0.1    
419   AddMaterial(M57, "G4_La");                      
420   indexZ[57] = 57;                                
421                                                   
422   // G4_Ce  index=58                              
423   G4double M58[NDENSARRAY] = {                    
424     47.834, 1.461, 5.7837, 0.2713, 3.3432, 0.1    
425   AddMaterial(M58, "G4_Ce");                      
426   indexZ[58] = 58;                                
427                                                   
428   // G4_Pr  index=59                              
429   G4double M59[NDENSARRAY] = {                    
430     48.301, 1.520, 5.8096, 0.2333, 3.2773, 0.2    
431   AddMaterial(M59, "G4_Pr");                      
432   indexZ[59] = 59;                                
433                                                   
434   // G4_Ne  index=60                              
435   G4double M60[NDENSARRAY] = {                    
436     48.819, 1.588, 5.8290, 0.1984, 3.3063, 0.2    
437   AddMaterial(M60, "G4_Ne");                      
438   indexZ[60] = 60;                                
439                                                   
440   // G4_Pr  index=61                              
441   G4double M61[NDENSARRAY] = {                    
442     50.236, 1.672, 5.8224, 0.1627, 3.3199, 0.2    
443   AddMaterial(M61, "G4_Pr");                      
444   indexZ[61] = 61;                                
445                                                   
446   // G4_Sa  index=62                              
447   G4double M62[NDENSARRAY] = {                    
448     50.540, 1.749, 5.8597, 0.1520, 3.3460, 0.2    
449   AddMaterial(M62, "G4_Sa");                      
450   indexZ[62] = 62;                                
451                                                   
452   // G4_Eu  index=63                              
453   G4double M63[NDENSARRAY] = {                    
454     42.484, 1.838, 6.2278, 0.1888, 3.4633, 0.2    
455   AddMaterial(M63, "G4_Eu");                      
456   indexZ[63] = 63;                                
457                                                   
458   // G4_Gd  index=64                              
459   G4double M64[NDENSARRAY] = {                    
460     51.672, 1.882, 5.8738, 0.1058, 3.3932, 0.2    
461   AddMaterial(M64, "G4_Gd");                      
462   indexZ[64] = 64;                                
463                                                   
464   // G4_Tb  index=65                              
465   G4double M65[NDENSARRAY] = {                    
466     52.865, 1.993, 5.9045, 0.0947, 3.4224, 0.2    
467   AddMaterial(M65, "G4_Tb");                      
468   indexZ[65] = 65;                                
469                                                   
470   // G4_Dy  index=66                              
471   G4double M66[NDENSARRAY] = {                    
472     53.698, 2.081, 5.9183, 0.0822, 3.4474, 0.2    
473   AddMaterial(M66, "G4_Dy");                      
474   indexZ[66] = 66;                                
475                                                   
476   // G4_Ho  index=67                              
477   G4double M67[NDENSARRAY] = {                    
478     54.467, 2.197, 5.9587, 0.0761, 3.4782, 0.2    
479   AddMaterial(M67, "G4_Ho");                      
480   indexZ[67] = 67;                                
481                                                   
482   // G4_Er  index=68                              
483   G4double M68[NDENSARRAY] = {                    
484     55.322, 2.26, 5.9521, 0.0648, 3.4922, 0.24    
485   AddMaterial(M68, "G4_Er");                      
486   indexZ[68] = 68;                                
487                                                   
488   // G4_Tm  index=69                              
489   G4double M69[NDENSARRAY] = {                    
490     56.225, 2.333, 5.9677, 0.0812, 3.5085, 0.2    
491   AddMaterial(M69, "G4_Tm");                      
492   indexZ[69] = 69;                                
493                                                   
494   // G4_Yb  index=70                              
495   G4double M70[NDENSARRAY] = {                    
496     47.546, 2.505, 6.3325, 0.1199, 3.6246, 0.2    
497   AddMaterial(M70, "G4_Yb");                      
498   indexZ[70] = 70;                                
499                                                   
500   // G4_Lu  index=71                              
501   G4double M71[NDENSARRAY] = {                    
502     57.581, 2.348, 5.9785, 0.1560, 3.5218, 0.2    
503   AddMaterial(M71, "G4_Lu");                      
504   indexZ[71] = 71;                                
505                                                   
506   // G4_Hf  index=72                              
507   G4double M72[NDENSARRAY] = {                    
508     66.770, 2.174, 5.7139, 0.1965, 3.4337, 0.2    
509   AddMaterial(M72, "G4_Hf");                      
510   indexZ[72] = 72;                                
511                                                   
512   // G4_Ta  index=73                              
513   G4double M73[NDENSARRAY] = {                    
514     74.692, 2.07, 5.5262, 0.2117, 3.4805, 0.17    
515   AddMaterial(M73, "G4_Ta");                      
516   indexZ[73] = 73;                                
517                                                   
518   // G4_W  index=74                               
519   G4double M74[NDENSARRAY] = {                    
520     80.315, 1.997, 5.4059, 0.2167, 3.496, 0.15    
521   AddMaterial(M74, "G4_W");                       
522   indexZ[74] = 74;                                
523                                                   
524   // G4_Re  index=75                              
525   G4double M75[NDENSARRAY] = {                    
526     83.846, 1.976, 5.3445, 0.0559, 3.4845, 0.1    
527   AddMaterial(M75, "G4_Re");                      
528   indexZ[75] = 75;                                
529                                                   
530   // G4_Os  index=76                              
531   G4double M76[NDENSARRAY] = {                    
532     86.537, 1.947, 5.3083, 0.0891, 3.5414, 0.1    
533   AddMaterial(M76, "G4_Os");                      
534   indexZ[76] = 76;                                
535                                                   
536   // G4_Ir  index=77                              
537   G4double M77[NDENSARRAY] = {                    
538     86.357, 1.927, 5.3418, 0.0819, 3.5480, 0.1    
539   AddMaterial(M77, "G4_Ir");                      
540   indexZ[77] = 77;                                
541                                                   
542   // G4_Pt  index=78                              
543   G4double M78[NDENSARRAY] = {                    
544     84.389, 1.965, 5.4732, 0.1484, 3.6212, 0.1    
545   AddMaterial(M78, "G4_Pt");                      
546   indexZ[78] = 78;                                
547                                                   
548   // G4_Au  index=79                              
549   G4double M79[NDENSARRAY] = {                    
550     80.215, 1.926, 5.5747, 0.2021, 3.6979, 0.0    
551   AddMaterial(M79, "G4_Au");                      
552   indexZ[79] = 79;                                
553                                                   
554   // G4_Hg  index=80                              
555   G4double M80[NDENSARRAY] = {                    
556     66.977, 1.904, 5.9605, 0.2756, 3.7275, 0.1    
557   AddMaterial(M80, "G4_Hg");                      
558   indexZ[80] = 80;                                
559                                                   
560   // G4_Tl  index=81                              
561   G4double M81[NDENSARRAY] = {                    
562     62.104, 1.814, 6.1365, 0.3491, 3.8044, 0.0    
563   AddMaterial(M81, "G4_Tl");                      
564   indexZ[81] = 81;                                
565                                                   
566   // G4_Pb  index=82                              
567   G4double M82[NDENSARRAY] = {                    
568     61.072, 1.755, 6.2018, 0.3776, 3.8073, 0.0    
569   AddMaterial(M82, "G4_Pb");                      
570   indexZ[82] = 82;                                
571                                                   
572   // G4_Bi  index=83                              
573   G4double M83[NDENSARRAY] = {                    
574     56.696, 1.684, 6.3505, 0.4152, 3.8248, 0.0    
575   AddMaterial(M83, "G4_Bi");                      
576   indexZ[83] = 83;                                
577                                                   
578   // G4_Po  index=84                              
579   G4double M84[NDENSARRAY] = {                    
580     55.773, 1.637, 6.4003, 0.4267, 3.8293, 0.0    
581   AddMaterial(M84, "G4_Po");                      
582   indexZ[84] = 84;                                
583                                                   
584   // G4_Rn  index=85                              
585   G4double M85[NDENSARRAY] = {                    
586     1.708, 1.458, 13.2839, 1.5368, 4.9889, 0.2    
587   AddMaterial(M85, "G4_Rn");                      
588   indexZ[86] = 85;                                
589   state[86] = kStateGas;                          
590                                                   
591   // G4_Ra  index=86                              
592   G4double M86[NDENSARRAY] = {                    
593     40.205, 1.403, 7.0452, 0.5991, 3.9428, 0.0    
594   AddMaterial(M86, "G4_Ra");                      
595   indexZ[88] = 86;                                
596                                                   
597   // G4_Ac  index=87                              
598   G4double M87[NDENSARRAY] = {                    
599     57.254, 1.380, 6.3742, 0.4559, 3.7966, 0.0    
600   AddMaterial(M87, "G4_Ac");                      
601   indexZ[89] = 87;                                
602                                                   
603   // G4_Th  index=88                              
604   G4double M88[NDENSARRAY] = {                    
605     61.438, 1.363, 6.2473, 0.4202, 3.7681, 0.0    
606   AddMaterial(M88, "G4_Th");                      
607   indexZ[90] = 88;                                
608                                                   
609   // G4_Pa  index=89                              
610   G4double M89[NDENSARRAY] = {                    
611     70.901, 1.42, 6.0327, 0.3144, 3.5079, .147    
612   AddMaterial(M89, "G4_Pa");                      
613   indexZ[91] = 89;                                
614                                                   
615   // G4_U  index=90                               
616   G4double M90[NDENSARRAY] = {                    
617     77.986, 1.447, 5.8694, 0.2260, 3.3721, .19    
618   AddMaterial(M90, "G4_U");                       
619   indexZ[92] = 90;                                
620                                                   
621   // G4_Np  index=91                              
622   G4double M91[NDENSARRAY] = {                    
623     81.221, 1.468, 5.8149, 0.1869, 3.369, 0.19    
624   AddMaterial(M91, "G4_Np");                      
625   indexZ[93] = 91;                                
626                                                   
627   // G4_Pu  index=92                              
628   G4double M92[NDENSARRAY] = {                    
629     80.486, 1.519, 5.8748, 0.1557, 3.3981, 0.2    
630   AddMaterial(M92, "G4_Pu");                      
631   indexZ[94] = 92;                                
632                                                   
633   // G4_Am  index=93                              
634   G4double M93[NDENSARRAY] = {                    
635     66.607, 1.552, 6.2813, 0.2274, 3.5021, 0.2    
636   AddMaterial(M93, "G4_Am");                      
637   indexZ[95] = 93;                                
638                                                   
639   // G4_Cm  index=94                              
640   G4double M94[NDENSARRAY] = {                    
641     66.022, 1.559, 6.3097, 0.2484, 3.516, .202    
642   AddMaterial(M94, "G4_Cm");                      
643   indexZ[96] = 94;                                
644                                                   
645   // G4_Bk  index=95                              
646   G4double M95[NDENSARRAY] = {                    
647     67.557, 1.574, 6.2912, 0.2378, 3.5186, .20    
648   AddMaterial(M95, "G4_Bk");                      
649   indexZ[97] = 95;                                
650                                                   
651   // G4_A-150_TISSUE  index=96                    
652   G4double M96[NDENSARRAY] = {                    
653     22.667, 1.950, 3.1100, 0.1329, 2.6234, 0.1    
654   AddMaterial(M96, "G4_A-150_TISSUE");            
655                                                   
656   // G4_ACETONE  index=97                         
657   G4double M97[NDENSARRAY] = {                    
658     19.010, 1.976, 3.4341, 0.2197, 2.6928, 0.1    
659   AddMaterial(M97, "G4_ACETONE");                 
660                                                   
661   // G4_ACETYLENE  index=98                       
662   G4double M98[NDENSARRAY] = {                    
663     0.700, 1.784, 9.8419, 1.6017, 4.0074, 0.12    
664   AddMaterial(M98, "G4_ACETYLENE");               
665                                                   
666   // G4_ADENINE  index=99                         
667   G4double M99[NDENSARRAY] = {                    
668     24.098, 1.892, 3.1724, 0.1295, 2.4219, 0.2    
669   AddMaterial(M99, "G4_ADENINE");                 
670                                                   
671   // G4_ADIPOSE_TISSUE_ICRP  index=100            
672   G4double M100[NDENSARRAY] = {                   
673     20.655, 1.987, 3.2367, 0.1827, 2.6530, 0.1    
674   AddMaterial(M100, "G4_ADIPOSE_TISSUE_ICRP");    
675                                                   
676   // G4_AIR  index=101                            
677   G4double M101[NDENSARRAY] = {                   
678     0.707, 2.054, 10.5961, 1.7418, 4.2759, 0.1    
679   AddMaterial(M101, "G4_AIR");                    
680                                                   
681   // G4_ALANINE  index=102                        
682   G4double M102[NDENSARRAY] = {                   
683     25.204, 2.074, 3.0965, 0.1354, 2.6336, 0.1    
684   AddMaterial(M102, "G4_ALANINE");                
685                                                   
686   // G4_ALUMINIM_OXIDE  index=103                 
687   G4double M103[NDENSARRAY] = {                   
688     40.206, 2.394, 3.5682, 0.0402, 2.8665, 0.0    
689   AddMaterial(M103, "G4_ALUMINIM_OXIDE");         
690                                                   
691   // G4_AMBER  index=104                          
692   G4double M104[NDENSARRAY] = {                   
693     22.450, 1.946, 3.0701, 0.1335, 2.5610, 0.1    
694   AddMaterial(M104, "G4_AMBER");                  
695                                                   
696   // G4_AMMONIA  index=105                        
697   G4double M105[NDENSARRAY] = {                   
698     0.635, 1.814, 9.8763, 1.6822, 4.1158, 0.08    
699   AddMaterial(M105, "G4_AMMONIA");                
700                                                   
701   // G4_ANILINE  index=106                        
702   G4double M106[NDENSARRAY] = {                   
703     21.361, 1.938, 3.2622, 0.1618, 2.5805, 0.1    
704   AddMaterial(M106, "G4_ANILINE");                
705                                                   
706   // G4_ANTHRACENE  index=107                     
707   G4double M107[NDENSARRAY] = {                   
708     23.704, 1.954, 3.1514, 0.1146, 2.5213, 0.1    
709   AddMaterial(M107, "G4_ANTHRACENE");             
710                                                   
711   // G4_B-100_BONE  index=108                     
712   G4double M108[NDENSARRAY] = {                   
713     25.199, 2.013, 3.4528, 0.1252, 3.042, 0.05    
714   AddMaterial(M108, "G4_B-100_BONE");             
715                                                   
716   // G4_BAKELITE  index=109                       
717   G4double M109[NDENSARRAY] = {                   
718     23.408, 2.046, 3.2582, 0.1471, 2.6055, 0.1    
719   AddMaterial(M109, "G4_BAKELITE");               
720                                                   
721   // G4_BARIUM_FLUORIDE  index=110                
722   G4double M110[NDENSARRAY] = {                   
723     41.398, 1.727, 5.4122, -0.0098, 3.3871, 0.    
724   AddMaterial(M110, "G4_BARIUM_FLUORIDE");        
725                                                   
726   // G4_BARIUM_SULFATE  index=111                 
727   G4double M111[NDENSARRAY] = {                   
728     40.805, 1.893, 4.8923, -0.0128, 3.4069, 0.    
729   AddMaterial(M111, "G4_BARIUM_SULFATE");         
730                                                   
731   // G4_BENZENE  index=112                        
732   G4double M112[NDENSARRAY] = {                   
733     19.806, 1.873, 3.3269, 0.171, 2.5091, 0.16    
734   AddMaterial(M112, "G4_BENZENE");                
735                                                   
736   // G4_BERYLLIUM_OXIDE  index=113                
737   G4double M113[NDENSARRAY] = {                   
738     34.629, 2.296, 2.9801, 0.0241, 2.5846, 0.1    
739   AddMaterial(M113, "G4_BERYLLIUM_OXIDE");        
740                                                   
741   // G4_BGO  index=114                            
742   G4double M114[NDENSARRAY] = {                   
743     49.904, 2.121, 5.7409, 0.0456, 3.7816, 0.0    
744   AddMaterial(M114, "G4_BGO");                    
745                                                   
746   // G4_BLOOD_ICRP  index=115                     
747   G4double M115[NDENSARRAY] = {                   
748     22.001, 2.184, 3.4581, 0.2239, 2.8017, 0.0    
749   AddMaterial(M115, "G4_BLOOD_ICRP");             
750                                                   
751   // G4_BONE_COMPACT_ICRU  index=116              
752   G4double M116[NDENSARRAY] = {                   
753     28.536, 2.091, 3.3390, 0.0944, 3.0201, 0.0    
754   AddMaterial(M116, "G4_BONE_COMPACT_ICRU");      
755                                                   
756   // G4_BONE_CORTICAL_ICRP  index=117             
757   G4double M117[NDENSARRAY] = {                   
758     28.298, 2.118, 3.6488, 0.1161, 3.0919, 0.0    
759   AddMaterial(M117, "G4_BONE_CORTICAL_ICRP");     
760                                                   
761   // G4_BORON_CARBIDE  index=118                  
762   G4double M118[NDENSARRAY] = {                   
763     31.38, 2.14, 2.9859, 0.0093, 2.1006, 0.370    
764   AddMaterial(M118, "G4_BORON_CARBIDE");          
765                                                   
766   // G4_BORON_OXIDE  index=119                    
767   G4double M119[NDENSARRAY] = {                   
768     27.107, 2.446, 3.6027, 0.1843, 2.7379, 0.1    
769   AddMaterial(M119, "G4_BORON_OXIDE");            
770                                                   
771   // G4_BRAIN_ICRP  index=120                     
772   G4double M120[NDENSARRAY] = {                   
773     21.772, 2.162, 3.4279, 0.2206, 2.8021, 0.0    
774   AddMaterial(M120, "G4_BRAIN_ICRP");             
775                                                   
776   // G4_BUTANE  index=121                         
777   G4double M121[NDENSARRAY] = {1.101, 1.727, 8    
778   AddMaterial(M121, "G4_BUTANE");                 
779                                                   
780   // G4_N-BUTYL_ALCOHOL  index=122                
781   G4double M122[NDENSARRAY] = {                   
782     19.52, 1.942, 3.2425, 0.1937, 2.6439, 0.10    
783   AddMaterial(M122, "G4_N-BUTYL_ALCOHOL");        
784                                                   
785   // G4_C-552  index=123                          
786   G4double M123[NDENSARRAY] = {                   
787     27.023, 2.128, 3.3338, 0.151, 2.7083, 0.10    
788   AddMaterial(M123, "G4_C-552");                  
789                                                   
790   // G4_CADMIUM_TELLURIDE  index=124              
791   G4double M124[NDENSARRAY] = {                   
792     46.314, 1.935, 5.9096, 0.0438, 3.2836, 0.2    
793   AddMaterial(M124, "G4_CADMIUM_TELLURIDE");      
794                                                   
795   // G4_CADMIUM_TUNGSTATE  index=125              
796   G4double M125[NDENSARRAY] = {                   
797     52.954, 2.289, 5.3594, 0.0123, 3.5941, 0.1    
798   AddMaterial(M125, "G4_CADMIUM_TUNGSTATE");      
799                                                   
800   // G4_CALCIUM_CARBONATE  index=126              
801   G4double M126[NDENSARRAY] = {                   
802     34.08, 2.141, 3.7738, 0.0492, 3.0549, 0.08    
803   AddMaterial(M126, "G4_CALCIUM_CARBONATE");      
804                                                   
805   // G4_CALCIUM_FLUORIDE  index=127               
806   G4double M127[NDENSARRAY] = {                   
807     35.849, 2.127, 4.0653, 0.0676, 3.1683, 0.0    
808   AddMaterial(M127, "G4_CALCIUM_FLUORIDE");       
809                                                   
810   // G4_CALCIUM_OXIDE  index=128                  
811   G4double M128[NDENSARRAY] = {                   
812     36.988, 1.973, 4.1209, -0.0172, 3.0171, 0.    
813   AddMaterial(M128, "G4_CALCIUM_OXIDE");          
814                                                   
815   // G4_CALCIUM_SULFATE  index=129                
816   G4double M129[NDENSARRAY] = {                   
817     35.038, 2.179, 3.9388, 0.0587, 3.1229, 0.0    
818   AddMaterial(M129, "G4_CALCIUM_SULFATE");        
819                                                   
820   // G4_CALCIUM_TUNGSTATE  index=130              
821   G4double M130[NDENSARRAY] = {                   
822     46.934, 2.262, 5.2603, 0.0323, 3.8932, 0.0    
823   AddMaterial(M130, "G4_CALCIUM_TUNGSTATE");      
824                                                   
825   // G4_CARBON_DIOXIDE  index=131                 
826   G4double M131[NDENSARRAY] = {                   
827     0.874, 2.118, 10.1537, 1.6294, 4.1825, 0.1    
828   AddMaterial(M131, "G4_CARBON_DIOXIDE");         
829                                                   
830   // G4_CARBON_TETRACHLORIDE  index=132           
831   G4double M132[NDENSARRAY] = {                   
832     25.234, 1.742, 4.7712, 0.1773, 2.9165, 0.1    
833   AddMaterial(M132, "G4_CARBON_TETRACHLORIDE")    
834                                                   
835   // G4_CELLULOSE_CELLOPHANE  index=133           
836   G4double M133[NDENSARRAY] = {25.008, 2.17, 3    
837   AddMaterial(M133, "G4_CELLULOSE_CELLOPHANE")    
838                                                   
839   // G4_CELLULOSE_BUTYRATE  index=134             
840   G4double M134[NDENSARRAY] = {                   
841     23.041, 2.128, 3.3497, 0.1794, 2.6809, 0.1    
842   AddMaterial(M134, "G4_CELLULOSE_BUTYRATE");     
843                                                   
844   // G4_CELLULOSE_NITRATE  index=135              
845   G4double M135[NDENSARRAY] = {                   
846     25.224, 2.252, 3.4762, 0.1897, 2.7253, 0.1    
847   AddMaterial(M135, "G4_CELLULOSE_NITRATE");      
848                                                   
849   // G4_CERIC_SULFATE  index=136                  
850   G4double M136[NDENSARRAY] = {                   
851     21.743, 2.205, 3.5212, 0.2363, 2.8769, 0.0    
852   AddMaterial(M136, "G4_CERIC_SULFATE");          
853                                                   
854   // G4_CESIUM_FLUORIDE  index=137                
855   G4double M137[NDENSARRAY] = {                   
856     37.942, 1.714, 5.9046, 0.0084, 3.3374, 0.2    
857   AddMaterial(M137, "G4_CESIUM_FLUORIDE");        
858                                                   
859   // G4_CESIUM_IODIDE  index=138                  
860   G4double M138[NDENSARRAY] = {                   
861     39.455, 1.672, 6.2807, 0.0395, 3.3353, 0.2    
862   AddMaterial(M138, "G4_CESIUM_IODIDE");          
863                                                   
864   // G4_CHLOROBENZENE  index=139                  
865   G4double M139[NDENSARRAY] = {                   
866     21.752, 1.889, 3.8201, 0.1714, 2.9272, 0.0    
867   AddMaterial(M139, "G4_CHLOROBENZENE");          
868                                                   
869   // G4_CHLOROFORM  index=140                     
870   G4double M140[NDENSARRAY] = {                   
871     24.462, 1.734, 4.7055, 0.1786, 2.9581, 0.1    
872   AddMaterial(M140, "G4_CHLOROFORM");             
873                                                   
874   // G4_CONCRETE  index=141                       
875   G4double M141[NDENSARRAY] = {                   
876     30.986, 2.322, 3.9464, 0.1301, 3.0466, 0.0    
877   AddMaterial(M141, "G4_CONCRETE");               
878                                                   
879   // G4_CYCLOHEXANE  index=142                    
880   G4double M142[NDENSARRAY] = {                   
881     19.207, 1.861, 3.1544, 0.1728, 2.5549, 0.1    
882   AddMaterial(M142, "G4_CYCLOHEXANE");            
883                                                   
884   // G4_1,2-DICHLOROBENZENE  index=143            
885   G4double M143[NDENSARRAY] = {                   
886     23.354, 1.862, 4.0348, 0.1587, 2.8276, 0.1    
887   AddMaterial(M143, "G4_1,2-DICHLOROBENZENE");    
888                                                   
889   // G4_DICHLORODIETHYL_ETHER  index=144          
890   G4double M144[NDENSARRAY] = {                   
891     22.894, 1.903, 4.0135, 0.1773, 3.1586, 0.0    
892   AddMaterial(M144, "G4_DICHLORODIETHYL_ETHER"    
893                                                   
894   // G4_1,2-DICHLOROETHANE  index=145             
895   G4double M145[NDENSARRAY] = {                   
896     22.764, 1.618, 4.1849, 0.1375, 2.9529, 0.1    
897   AddMaterial(M145, "G4_1,2-DICHLOROETHANE");     
898                                                   
899   // G4_DIETHYL_ETHER  index=146                  
900   G4double M146[NDENSARRAY] = {                   
901     18.326, 1.951, 3.3721, 0.2231, 2.6745, 0.1    
902   AddMaterial(M146, "G4_DIETHYL_ETHER");          
903                                                   
904   // G4_N,N-DIMETHYL_FORMAMIDE  index=147         
905   G4double M147[NDENSARRAY] = {                   
906     20.763, 2.005, 3.3311, 0.1977, 2.6686, 0.1    
907   AddMaterial(M147, "G4_N,N-DIMETHYL_FORMAMIDE    
908                                                   
909   // G4_DIMETHYL_SULFOXIDE  index=148             
910   G4double M148[NDENSARRAY] = {                   
911     22.173, 2.075, 3.9844, 0.2021, 3.1263, 0.0    
912   AddMaterial(M148, "G4_DIMETHYL_SULFOXIDE");     
913                                                   
914   // G4_ETHANE  index=149                         
915   G4double M149[NDENSARRAY] = {                   
916     0.789, 1.69, 9.1043, 1.5107, 3.8743, 0.096    
917   AddMaterial(M149, "G4_ETHANE");                 
918                                                   
919   // G4_ETHYL_ALCOHOL  index=150                  
920   G4double M150[NDENSARRAY] = {                   
921     19.232, 2.013, 3.3699, 0.2218, 2.7052, 0.0    
922   AddMaterial(M150, "G4_ETHYL_ALCOHOL");          
923                                                   
924   // G4_ETHYL_CELLULOSE  index=151                
925   G4double M151[NDENSARRAY] = {                   
926     22.594, 2.065, 3.2415, 0.1683, 2.6527, 0.1    
927   AddMaterial(M151, "G4_ETHYL_CELLULOSE");        
928                                                   
929   // G4_ETHYLENE  index=152                       
930   G4double M152[NDENSARRAY] = {                   
931     0.746, 1.733, 9.438, 1.5528, 3.9327, 0.106    
932   AddMaterial(M152, "G4_ETHYLENE");               
933                                                   
934   // G4_EYE_LENS_ICRP  index=153                  
935   G4double M153[NDENSARRAY] = {22.388, 2.154,     
936   AddMaterial(M153, "G4_EYE_LENS_ICRP");          
937                                                   
938   // G4_FERRIC_OXIDE  index=154                   
939   G4double M154[NDENSARRAY] = {                   
940     45.331, 2.747, 4.2245, -0.0074, 3.2573, 0.    
941   AddMaterial(M154, "G4_FERRIC_OXIDE");           
942                                                   
943   // G4_FERROBORIDE  index=155                    
944   G4double M155[NDENSARRAY] = {                   
945     52.546, 2.726, 4.2057, -0.0988, 3.1749, 0.    
946   AddMaterial(M155, "G4_FERROBORIDE");            
947                                                   
948   // G4_FERROUS_OXIDE  index=156                  
949   G4double M156[NDENSARRAY] = {                   
950     47.327, 2.769, 4.3175, -0.0279, 3.2002, 0.    
951   AddMaterial(M156, "G4_FERROUS_OXIDE");          
952                                                   
953   // G4_FERROUS_SULFATE  index=157                
954   G4double M157[NDENSARRAY] = {                   
955     21.69, 2.208, 3.5183, 0.2378, 2.8254, 0.08    
956   AddMaterial(M157, "G4_FERROUS_SULFATE");        
957                                                   
958   // G4_FREON-12  index=158                       
959   G4double M158[NDENSARRAY] = {                   
960     21.121, 1.974, 4.8251, 0.3035, 3.2659, 0.0    
961   AddMaterial(M158, "G4_FREON-12");               
962                                                   
963   // G4_FREON-12B2  index=159                     
964   G4double M159[NDENSARRAY] = {                   
965     25.877, 2.195, 5.7976, 0.3406, 3.7956, 0.0    
966   AddMaterial(M159, "G4_FREON-12B2");             
967                                                   
968   // G4_FREON-13  index=160                       
969   G4double M160[NDENSARRAY] = {                   
970     19.432, 2.116, 4.7483, 0.3659, 3.2337, 0.0    
971   AddMaterial(M160, "G4_FREON-13");               
972                                                   
973   // G4_FREON-13B1  index=161                     
974   G4double M161[NDENSARRAY] = {                   
975     23.849, 2.233, 5.3555, 0.3522, 3.7554, 0.0    
976   AddMaterial(M161, "G4_FREON-13B1");             
977                                                   
978   // G4_FREON-13I1  index=162                     
979   G4double M162[NDENSARRAY] = {                   
980     25.615, 1.924, 5.8774, 0.2847, 3.728, 0.09    
981   AddMaterial(M162, "G4_FREON-13I1");             
982                                                   
983   // G4_GADOLINIUM_OXYSULFIDE  index=163          
984   G4double M163[NDENSARRAY] = {                   
985     51.099, 2.179, 5.5347, -0.1774, 3.4045, 0.    
986   AddMaterial(M163, "G4_GADOLINIUM_OXYSULFIDE"    
987                                                   
988   // G4_GALLIUM_ARSENIDE  index=164               
989   G4double M164[NDENSARRAY] = {                   
990     44.17, 2.652, 5.3299, 0.1764, 3.642, 0.071    
991   AddMaterial(M164, "G4_GALLIUM_ARSENIDE");       
992                                                   
993   // G4_GEL_PHOTO_EMULSION  index=165             
994   G4double M165[NDENSARRAY] = {                   
995     24.058, 2.156, 3.2687, 0.1709, 2.7058, 0.1    
996   AddMaterial(M165, "G4_GEL_PHOTO_EMULSION");     
997                                                   
998   // G4_Pyrex_Glass  index=166                    
999   G4double M166[NDENSARRAY] = {                   
1000     30.339, 2.369, 3.9708, 0.1479, 2.9933, 0.    
1001   AddMaterial(M166, "G4_Pyrex_Glass");           
1002                                                  
1003   // G4_GLASS_LEAD  index=167                    
1004   G4double M167[NDENSARRAY] = {                  
1005     46.631, 2.085, 5.8476, 0.0614, 3.8146, 0.    
1006   AddMaterial(M167, "G4_GLASS_LEAD");            
1007                                                  
1008   // G4_GLASS_PLATE  index=168                   
1009   G4double M168[NDENSARRAY] = {                  
1010     31.481, 2.329, 4.0602, 0.1237, 3.0649, 0.    
1011   AddMaterial(M168, "G4_GLASS_PLATE");           
1012                                                  
1013   // G4_GLUCOSE  index=169                       
1014   G4double M169[NDENSARRAY] = {                  
1015     26.153, 2.174, 3.1649, 0.1411, 2.67, 0.10    
1016   AddMaterial(M169, "G4_GLUCOSE");               
1017                                                  
1018   // G4_GLUTAMINE  index=170                     
1019   G4double M170[NDENSARRAY] = {                  
1020     25.437, 2.077, 3.1167, 0.1347, 2.6301, 0.    
1021   AddMaterial(M170, "G4_GLUTAMINE");             
1022                                                  
1023   // G4_GLYCEROL  index=171                      
1024   G4double M171[NDENSARRAY] = {                  
1025     23.846, 2.12, 3.2267, 0.1653, 2.6862, 0.1    
1026   AddMaterial(M171, "G4_GLYCEROL");              
1027                                                  
1028   // G4_GUANINE  index=172                       
1029   G4double M172[NDENSARRAY] = {                  
1030     26.022, 1.97, 3.1171, 0.1163, 2.4296, 0.2    
1031   AddMaterial(M172, "G4_GUANINE");               
1032                                                  
1033   // G4_GYPSUM  index=173                        
1034   G4double M173[NDENSARRAY] = {                  
1035     31.379, 2.187, 3.8382, 0.0995, 3.1206, 0.    
1036   AddMaterial(M173, "G4_GYPSUM");                
1037                                                  
1038   // G4_N-HEPTANE  index=174                     
1039   G4double M174[NDENSARRAY] = {                  
1040     18.128, 1.848, 3.1978, 0.1928, 2.5706, 0.    
1041   AddMaterial(M174, "G4_N-HEPTANE");             
1042                                                  
1043   // G4_N-HEXANE  index=175                      
1044   G4double M175[NDENSARRAY] = {                  
1045     17.836, 1.843, 3.2156, 0.1984, 2.5757, 0.    
1046   AddMaterial(M175, "G4_N-HEXANE");              
1047                                                  
1048   // G4_KAPTON  index=176                        
1049   G4double M176[NDENSARRAY] = {                  
1050     24.586, 2.109, 3.3497, 0.1509, 2.5631, 0.    
1051   AddMaterial(M176, "G4_KAPTON");                
1052                                                  
1053   // G4_LANTHANUM_OXYBROMIDE  index=177          
1054   G4double M177[NDENSARRAY] = {                  
1055     47.125, 1.831, 5.4666, -0.035, 3.3288, 0.    
1056   AddMaterial(M177, "G4_LANTHANUM_OXYBROMIDE"    
1057                                                  
1058   // G4_LANTHANUM_OXYSULFIDE  index=178          
1059   G4double M178[NDENSARRAY] = {                  
1060     45.394, 1.681, 5.6151, -0.0934, 3.2741, 0    
1061   AddMaterial(M178, "G4_LANTHANUM_OXYSULFIDE"    
1062                                                  
1063   // G4_LEAD_OXIDE  index=179                    
1064   G4double M179[NDENSARRAY] = {                  
1065     56.488, 2.012, 6.2162, 0.0356, 3.5456, 0.    
1066   AddMaterial(M179, "G4_LEAD_OXIDE");            
1067                                                  
1068   // G4_LITHIUM_AMIDE  index=180                 
1069   G4double M180[NDENSARRAY] = {22.609, 1.74,     
1070   AddMaterial(M180, "G4_LITHIUM_AMIDE");         
1071                                                  
1072   // G4_LITHIUM_CARBONATE index=181              
1073   G4double M181[NDENSARRAY] = {                  
1074     29.217, 2.246, 3.2029, 0.0551, 2.6598, 0.    
1075   AddMaterial(M181, "G4_LITHIUM_CARBONATE");     
1076                                                  
1077   // G4_LITHIUM_FLUORIDE  index=182              
1078   G4double M182[NDENSARRAY] = {                  
1079     31.815, 2.197, 3.1667, 0.0171, 2.7049, 0.    
1080   AddMaterial(M182, "G4_LITHIUM_FLUORIDE");      
1081                                                  
1082   // G4_LITHIUM_HYDRIDE  index=183               
1083   G4double M183[NDENSARRAY] = {                  
1084     18.51, 1.482, 2.358, -0.0988, 1.4515, 0.9    
1085   AddMaterial(M183, "G4_LITHIUM_HYDRIDE");       
1086                                                  
1087   // G4_LITHIUM_IODIDE  index=184                
1088   G4double M184[NDENSARRAY] = {                  
1089     34.841, 1.706, 6.2671, 0.0892, 3.3702, 0.    
1090   AddMaterial(M184, "G4_LITHIUM_IODIDE");        
1091                                                  
1092   // G4_LITHIUM_OXIDE  index=185                 
1093   G4double M185[NDENSARRAY] = {                  
1094     27.984, 2.039, 2.934, -0.0511, 2.5874, 0.    
1095   AddMaterial(M185, "G4_LITHIUM_OXIDE");         
1096                                                  
1097   // G4_LITHIUM_TETRABORATE  index=186           
1098   G4double M186[NDENSARRAY] = {                  
1099     31.343, 2.36, 3.2093, 0.0737, 2.6502, 0.1    
1100   AddMaterial(M186, "G4_LITHIUM_TETRABORATE")    
1101                                                  
1102   // G4_LUNG_ICRP  index=187                     
1103   G4double M187[NDENSARRAY] = {                  
1104     21.891, 2.184, 3.4708, 0.2261, 2.8001, 0.    
1105   AddMaterial(M187, "G4_LUNG_ICRP");             
1106                                                  
1107   // G4_M3_WAX  index=188                        
1108   G4double M188[NDENSARRAY] = {22, 1.975, 3.2    
1109   AddMaterial(M188, "G4_M3_WAX");                
1110                                                  
1111   // G4_MAGNESIUM_CARBONATE  index=189           
1112   G4double M189[NDENSARRAY] = {                  
1113     34.979, 2.388, 3.4319, 0.086, 2.7997, 0.0    
1114   AddMaterial(M189, "G4_MAGNESIUM_CARBONATE")    
1115                                                  
1116   // G4_MAGNESIUM_FLUORIDE  index=190            
1117   G4double M190[NDENSARRAY] = {                  
1118     34.634, 2.33, 3.7105, 0.1369, 2.863, 0.07    
1119   AddMaterial(M190, "G4_MAGNESIUM_FLUORIDE");    
1120                                                  
1121   // G4_MAGNESIUM_OXIDE  index=191               
1122   G4double M191[NDENSARRAY] = {                  
1123     38.407, 2.412, 3.6404, 0.0575, 2.858, 0.0    
1124   AddMaterial(M191, "G4_MAGNESIUM_OXIDE");       
1125                                                  
1126   // G4_MAGNESIUM_TETRABORATE  index=192         
1127   G4double M192[NDENSARRAY] = {                  
1128     32.089, 2.43, 3.4328, 0.1147, 2.7635, 0.0    
1129   AddMaterial(M192, "G4_MAGNESIUM_TETRABORATE    
1130                                                  
1131   // G4_MERCURIC_IODIDE  index=193               
1132   G4double M193[NDENSARRAY] = {                  
1133     46.494, 1.892, 6.3787, 0.104, 3.4728, 0.2    
1134   AddMaterial(M193, "G4_MERCURIC_IODIDE");       
1135                                                  
1136   // G4_METHANE  index=194                       
1137   G4double M194[NDENSARRAY] = {                  
1138     0.588, 1.662, 9.5243, 1.6263, 3.9716, 0.0    
1139   AddMaterial(M194, "G4_METHANE");               
1140                                                  
1141   // G4_METHANOL  index=195                      
1142   G4double M195[NDENSARRAY] = {19.214, 2.125,    
1143   AddMaterial(M195, "G4_METHANOL");              
1144                                                  
1145   // G4_MIX_D_WAX  index=196                     
1146   G4double M196[NDENSARRAY] = {                  
1147     21.547, 1.905, 3.078, 0.1371, 2.7145, 0.0    
1148   AddMaterial(M196, "G4_MIX_D_WAX");             
1149                                                  
1150   // G4_MS20_TISSUE  index=197                   
1151   G4double M197[NDENSARRAY] = {                  
1152     21.153, 2.07, 3.5341, 0.1997, 2.8033, 0.0    
1153   AddMaterial(M197, "G4_MS20_TISSUE");           
1154                                                  
1155   // G4_MUSCLE_SCELETAL_ICRP  index=198          
1156   G4double M198[NDENSARRAY] = {                  
1157     21.781, 2.185, 3.4809, 0.2282, 2.7999, 0.    
1158   AddMaterial(M198, "G4_MUSCLE_SCELETAL_ICRP"    
1159                                                  
1160   // G4_MUSCLE_STRIATED_ICRU  index=199          
1161   G4double M199[NDENSARRAY] = {                  
1162     21.795, 2.174, 3.4636, 0.2249, 2.8032, 0.    
1163   AddMaterial(M199, "G4_MUSCLE_STRIATED_ICRU"    
1164                                                  
1165   // G4_MUSCLE_WITH_SUCROSE  index=200           
1166   G4double M200[NDENSARRAY] = {22.48, 2.169,     
1167   AddMaterial(M200, "G4_MUSCLE_WITH_SUCROSE")    
1168                                                  
1169   // G4_MUSCLE_WITHOUT_SUCROSE  index=201        
1170   G4double M201[NDENSARRAY] = {                  
1171     22.109, 2.173, 3.4216, 0.2187, 2.768, 0.0    
1172   AddMaterial(M201, "G4_MUSCLE_WITHOUT_SUCROS    
1173                                                  
1174   // G4_NAPHTALENE  index=202                    
1175   G4double M202[NDENSARRAY] = {                  
1176     22.459, 1.956, 3.2274, 0.1374, 2.5429, 0.    
1177   AddMaterial(M202, "G4_NAPHTALENE");            
1178                                                  
1179   // G4_NITROBENZENE  index=203                  
1180   G4double M203[NDENSARRAY] = {                  
1181     22.747, 2.065, 3.4073, 0.1777, 2.663, 0.1    
1182   AddMaterial(M203, "G4_NITROBENZENE");          
1183                                                  
1184   // G4_NITROUS_OXIDE  index=204                 
1185   G4double M204[NDENSARRAY] = {                  
1186     0.872, 2.059, 10.1575, 1.6477, 4.1565, 0.    
1187   AddMaterial(M204, "G4_NITROUS_OXIDE");         
1188                                                  
1189   // G4_NYLON-8062  index=205                    
1190   G4double M205[NDENSARRAY] = {                  
1191     22.221, 1.967, 3.125, 0.1503, 2.6004, 0.1    
1192   AddMaterial(M205, "G4_NYLON-8062");            
1193                                                  
1194   // G4_NYLON-6/6  index=206                     
1195   G4double M206[NDENSARRAY] = {                  
1196     22.774, 1.931, 3.0634, 0.1336, 2.5834, 0.    
1197   AddMaterial(M206, "G4_NYLON-6-6");             
1198                                                  
1199   // G4_NYLON-6/10  index=207                    
1200   G4double M207[NDENSARRAY] = {                  
1201     22.866, 1.942, 3.0333, 0.1304, 2.5681, 0.    
1202   AddMaterial(M207, "G4_NYLON-6-10");            
1203                                                  
1204   // G4_NYLON-11_RILSAN  index=208               
1205   G4double M208[NDENSARRAY] = {                  
1206     25.661, 1.902, 2.7514, 0.0678, 2.4281, 0.    
1207   AddMaterial(M208, "G4_NYLON-11_RILSAN");       
1208                                                  
1209   // G4_OCTANE  index=209                        
1210   G4double M209[NDENSARRAY] = {                  
1211     18.36, 1.851, 3.1834, 0.1882, 2.5664, 0.1    
1212   AddMaterial(M209, "G4_OCTANE");                
1213                                                  
1214   // G4_PARAFFIN  index=210                      
1215   G4double M210[NDENSARRAY] = {                  
1216     21.031, 1.844, 2.9551, 0.1289, 2.5084, 0.    
1217   AddMaterial(M210, "G4_PARAFFIN");              
1218                                                  
1219   // G4_N-PENTANE  index=211                     
1220   G4double M211[NDENSARRAY] = {                  
1221     17.398, 1.842, 3.2504, 0.2086, 2.5855, 0.    
1222   AddMaterial(M211, "G4_N-PENTANE");             
1223                                                  
1224   // G4_PHOTO_EMULSION  index=212                
1225   G4double M212[NDENSARRAY] = {                  
1226     37.946, 2.264, 5.3319, 0.1009, 3.4866, 0.    
1227   AddMaterial(M212, "G4_PHOTO_EMULSION");        
1228                                                  
1229   // G4_PLASTIC_SC_VINYLTOLUENE  index=213       
1230   G4double M213[NDENSARRAY] = {                  
1231     21.54, 1.929, 3.1997, 0.1464, 2.4855, 0.1    
1232   AddMaterial(M213, "G4_PLASTIC_SC_VINYLTOLUE    
1233                                                  
1234   // G4_PLUTONIUM_DIOXIDE  index=214             
1235   G4double M214[NDENSARRAY] = {                  
1236     62.143, 1.846, 5.9719, -0.2311, 3.5554, 0    
1237   AddMaterial(M214, "G4_PLUTONIUM_DIOXIDE");     
1238                                                  
1239   // G4_POLYACRYLONITRILE  index=215             
1240   G4double M215[NDENSARRAY] = {                  
1241     22.642, 1.955, 3.2459, 0.1504, 2.5159, 0.    
1242   AddMaterial(M215, "G4_POLYACRYLONITRILE");     
1243                                                  
1244   // G4_POLYCARBONATE  index=216                 
1245   G4double M216[NDENSARRAY] = {                  
1246     22.915, 2.06, 3.3201, 0.1606, 2.6225, 0.1    
1247   AddMaterial(M216, "G4_POLYCARBONATE");         
1248                                                  
1249   // G4_POLYCHLOROSTYRENE  index=217             
1250   G4double M217[NDENSARRAY] = {                  
1251     23.81, 1.902, 3.4659, 0.1238, 2.9241, 0.0    
1252   AddMaterial(M217, "G4_POLYCHLOROSTYRENE");     
1253                                                  
1254   // G4_POLYETHYLENE  index=218                  
1255   G4double M218[NDENSARRAY] = {                  
1256     21.099, 1.882, 3.0016, 0.137, 2.5177, 0.1    
1257   AddMaterial(M218, "G4_POLYETHYLENE");          
1258                                                  
1259   // G4_MYLAR  index=219                         
1260   G4double M219[NDENSARRAY] = {                  
1261     24.595, 2.144, 3.3262, 0.1562, 2.6507, 0.    
1262   AddMaterial(M219, "G4_MYLAR");                 
1263                                                  
1264   // G4_LUCITE  index=220                        
1265   G4double M220[NDENSARRAY] = {                  
1266     23.086, 2.173, 3.3297, 0.1824, 2.6681, 0.    
1267   AddMaterial(M220, "G4_LUCITE");                
1268                                                  
1269   // G4_POLYOXOMETHYLENE  index=221              
1270   G4double M221[NDENSARRAY] = {                  
1271     25.11, 2.175, 3.2514, 0.1584, 2.6838, 0.1    
1272   AddMaterial(M221, "G4_POLYOXOMETHYLENE");      
1273                                                  
1274   // G4_POLYPROPYLENE  index=222                 
1275   G4double M222[NDENSARRAY] = {                  
1276     20.457, 1.884, 3.1252, 0.1534, 2.4822, 0.    
1277   AddMaterial(M222, "G4_POLYPROPYLENE");         
1278                                                  
1279   // G4_POLYSTYRENE  index=223                   
1280   G4double M223[NDENSARRAY] = {                  
1281     21.754, 2.027, 3.2999, 0.1647, 2.5031, 0.    
1282   AddMaterial(M223, "G4_POLYSTYRENE");           
1283                                                  
1284   // G4_TEFLON  index=224                        
1285   G4double M224[NDENSARRAY] = {                  
1286     29.609, 2.142, 3.4161, 0.1648, 2.7404, 0.    
1287   AddMaterial(M224, "G4_TEFLON");                
1288                                                  
1289   // G4_POLYTRIFLUOROCHLOROETHYLENE  index=22    
1290   G4double M225[NDENSARRAY] = {                  
1291     28.955, 2.094, 3.8551, 0.1714, 3.0265, 0.    
1292   AddMaterial(M225, "G4_POLYTRIFLUOROCHLOROET    
1293                                                  
1294   // G4_POLYVINYL_ACETATE  index=226             
1295   G4double M226[NDENSARRAY] = {                  
1296     22.978, 2.116, 3.3309, 0.1769, 2.6747, 0.    
1297   AddMaterial(M226, "G4_POLYVINYL_ACETATE");     
1298                                                  
1299   // G4_PLOYVINYL_ALCOHOL  index=227             
1300   G4double M227[NDENSARRAY] = {                  
1301     24.251, 2.071, 3.1115, 0.1401, 2.6315, 0.    
1302   AddMaterial(M227, "G4_PLOYVINYL_ALCOHOL");     
1303                                                  
1304   // G4_POLYVINYL_BUTYRAL  index=228             
1305   G4double M228[NDENSARRAY] = {                  
1306     22.521, 2.021, 3.1865, 0.1555, 2.6186, 0.    
1307   AddMaterial(M228, "G4_POLYVINYL_BUTYRAL");     
1308                                                  
1309   // G4_POLYVINYL_CHLORIDE  index=229            
1310   G4double M229[NDENSARRAY] = {                  
1311     23.51, 1.84, 4.0532, 0.1559, 2.9415, 0.12    
1312   AddMaterial(M229, "G4_POLYVINYL_CHLORIDE");    
1313                                                  
1314   // G4_POLYVINYLIDENE_CHLORIDE  index=230       
1315   G4double M230[NDENSARRAY] = {                  
1316     26.437, 1.814, 4.2506, 0.1314, 2.9009, 0.    
1317   AddMaterial(M230, "G4_POLYVINYLIDENE_CHLORI    
1318                                                  
1319   // G4_POLYVINYLIDENE_FLUORIDE  index=231       
1320   G4double M231[NDENSARRAY] = {27.024, 2.16,     
1321   AddMaterial(M231, "G4_POLYVINYLIDENE_FLUORI    
1322                                                  
1323   // G4_POLYVINYL_PYRROLIDONE  index=232         
1324   G4double M232[NDENSARRAY] = {                  
1325     23.671, 1.989, 3.1017, 0.1324, 2.5867, 0.    
1326   AddMaterial(M232, "G4_POLYVINYL_PYRROLIDONE    
1327                                                  
1328   // G4_POTASSIUM_IODIDE  index=233              
1329   G4double M233[NDENSARRAY] = {                  
1330     33.575, 1.784, 6.1088, 0.1044, 3.3442, 0.    
1331   AddMaterial(M233, "G4_POTASSIUM_IODIDE");      
1332                                                  
1333   // G4_POTASSIUM_OXIDE  index=234               
1334   G4double M234[NDENSARRAY] = {                  
1335     30.672, 2.065, 4.6463, 0.048, 3.011, 0.16    
1336   AddMaterial(M234, "G4_POTASSIUM_OXIDE");       
1337                                                  
1338   // G4_PROPANE  index=235                       
1339   G4double M235[NDENSARRAY] = {                  
1340     0.959, 1.708, 8.7878, 1.4326, 3.7998, 0.0    
1341   AddMaterial(M235, "G4_PROPANE");               
1342                                                  
1343   // G4_lPROPANE  index=236                      
1344   G4double M236[NDENSARRAY] = {                  
1345     14.509, 1.844, 3.5529, 0.2861, 2.6568, 0.    
1346   AddMaterial(M236, "G4_lPROPANE");              
1347                                                  
1348   // G4_N-PROPYL_ALCOHOL index=237               
1349   G4double M237[NDENSARRAY] = {                  
1350     19.429, 1.972, 3.2915, 0.2046, 2.6681, 0.    
1351   AddMaterial(M237, "N-PROPYL_ALCOHOL");         
1352                                                  
1353   // G4_PYRIDINE  index=238                      
1354   G4double M238[NDENSARRAY] = {                  
1355     20.807, 1.895, 3.3148, 0.167, 2.5245, 0.1    
1356   AddMaterial(M238, "G4_PYRIDINE");              
1357                                                  
1358   // G4_RUBBER_BUTYL  index=239                  
1359   G4double M239[NDENSARRAY] = {                  
1360     20.873, 1.852, 2.9915, 0.1347, 2.5154, 0.    
1361   AddMaterial(M239, "G4_RUBBER_BUTYL");          
1362                                                  
1363   // G4_RUBBER_NATURAL  index=240                
1364   G4double M240[NDENSARRAY] = {                  
1365     20.644, 1.889, 3.1272, 0.1512, 2.4815, 0.    
1366   AddMaterial(M240, "G4_RUBBER_NATURAL");        
1367                                                  
1368   // G4_RUBBER_NEOPRENE  index=241               
1369   G4double M241[NDENSARRAY] = {                  
1370     23.036, 1.874, 3.7911, 0.1501, 2.9461, 0.    
1371   AddMaterial(M241, "G4_RUBBER_NEOPRENE");       
1372                                                  
1373   // G4_SILICON_DIOXIDE  index=242               
1374   G4double M242[NDENSARRAY] = {                  
1375     31.014, 2.335, 4.0029, 0.1385, 3.0025, 0.    
1376   AddMaterial(M242, "G4_SILICON_DIOXIDE");       
1377                                                  
1378   // G4_SILVER_BROMIDE  index=243                
1379   G4double M243[NDENSARRAY] = {                  
1380     48.448, 2.271, 5.6139, 0.0352, 3.2109, 0.    
1381   AddMaterial(M243, "G4_SILVER_BROMIDE");        
1382                                                  
1383   // G4_SILVER_CHLORIDE  index=244               
1384   G4double M244[NDENSARRAY] = {                  
1385     45.405, 2.096, 5.3437, -0.0139, 3.2022, 0    
1386   AddMaterial(M244, "G4_SILVER_CHLORIDE");       
1387                                                  
1388   // G4_SILVER_HALIDES  index=245                
1389   G4double M245[NDENSARRAY] = {                  
1390     48.433, 2.27, 5.6166, 0.0353, 3.2117, 0.2    
1391   AddMaterial(M245, "G4_SILVER_HALIDES");        
1392                                                  
1393   // G4_SILVER_IODIDE  index=246                 
1394   G4double M246[NDENSARRAY] = {                  
1395     46.105, 1.945, 5.9342, 0.0148, 3.2908, 0.    
1396   AddMaterial(M246, "G4_SILVER_IODIDE");         
1397                                                  
1398   // G4_SKIN_ICRP  index=247                     
1399   G4double M247[NDENSARRAY] = {22.4, 2.14, 3.    
1400   AddMaterial(M247, "G4_SKIN_ICRP");             
1401                                                  
1402   // G4_SODIUM_CARBONATE  index=248              
1403   G4double M248[NDENSARRAY] = {                  
1404     32.117, 2.557, 3.7178, 0.1287, 2.8591, 0.    
1405   AddMaterial(M248, "G4_SODIUM_CARBONATE");      
1406                                                  
1407   // G4_SODIUM_IODIDE  index=249                 
1408   G4double M249[NDENSARRAY] = {                  
1409     36.057, 1.857, 6.0572, 0.1203, 3.592, 0.1    
1410   AddMaterial(M249, "G4_SODIUM_IODIDE");         
1411                                                  
1412   // G4_SODIUM_MONOXIDE  index=250               
1413   G4double M250[NDENSARRAY] = {                  
1414     30.205, 2.689, 4.1892, 0.1652, 2.9793, 0.    
1415   AddMaterial(M250, "G4_SODIUM_MONOXIDE");       
1416                                                  
1417   // G4_SODIUM_NITRATE  index=251                
1418   G4double M251[NDENSARRAY] = {                  
1419     30.459, 2.456, 3.6502, 0.1534, 2.8221, 0.    
1420   AddMaterial(M251, "G4_SODIUM_NITRATE");        
1421                                                  
1422   // G4_STILBENE  index=252                      
1423   G4double M252[NDENSARRAY] = {                  
1424     20.719, 1.963, 3.368, 0.1734, 2.5142, 0.1    
1425   AddMaterial(M252, "G4_STILBENE");              
1426                                                  
1427   // G4_SUCROSE  index=253                       
1428   G4double M253[NDENSARRAY] = {                  
1429     26.416, 2.167, 3.1526, 0.1341, 2.6558, 0.    
1430   AddMaterial(M253, "G4_SUCROSE");               
1431                                                  
1432   // G4_TERPHENYL  index=254                     
1433   G4double M254[NDENSARRAY] = {                  
1434     23.116, 1.976, 3.2639, 0.1322, 2.5429, 0.    
1435   AddMaterial(M254, "G4_TERPHENYL");             
1436                                                  
1437   // G4_TESTES_ICRP  index=255                   
1438   G4double M255[NDENSARRAY] = {                  
1439     21.815, 2.185, 3.4698, 0.2274, 2.7988, 0.    
1440   AddMaterial(M255, "G4_TESTES_ICRP");           
1441                                                  
1442   // G4_TETRACHLOROETHYLENE  index=256           
1443   G4double M256[NDENSARRAY] = {                  
1444     25.513, 1.79, 4.6619, 0.1713, 2.9083, 0.1    
1445   AddMaterial(M256, "G4_TETRACHLOROETHYLENE")    
1446                                                  
1447   // G4_THALIUM_CHLORIDE  index=257              
1448   G4double M257[NDENSARRAY] = {                  
1449     48.749, 1.997, 6.3009, 0.0705, 3.5716, 0.    
1450   AddMaterial(M257, "G4_THALIUM_CHLORIDE");      
1451                                                  
1452   // G4_TISSUE_SOFT_ICRP  index=258              
1453   G4double M258[NDENSARRAY] = {                  
1454     21.394, 2.144, 3.4354, 0.2211, 2.7799, 0.    
1455   AddMaterial(M258, "G4_TISSUE_SOFT_ICRP");      
1456                                                  
1457   // G4_TISSUE_SOFT_ICRU-4  index=259            
1458   G4double M259[NDENSARRAY] = {                  
1459     21.366, 2.192, 3.5087, 0.2377, 2.7908, 0.    
1460   AddMaterial(M259, "G4_TISSUE_SOFT_ICRU-4");    
1461                                                  
1462   // G4_TISSUE-METHANE  index=260                
1463   G4double M260[NDENSARRAY] = {0.697, 1.89, 9    
1464   AddMaterial(M260, "G4_TISSUE-METHANE");        
1465                                                  
1466   // G4_TISSUE-PROPANE  index=261                
1467   G4double M261[NDENSARRAY] = {                  
1468     0.913, 1.856, 9.3529, 1.5139, 3.9916, 0.0    
1469   AddMaterial(M261, "G4_TISSUE-PROPANE");        
1470                                                  
1471   // G4_TITANIUM_DIOXIDE  index=262              
1472   G4double M262[NDENSARRAY] = {                  
1473     41.022, 2.307, 3.9522, -0.0119, 3.1647, 0    
1474   AddMaterial(M262, "G4_TITANIUM_DIOXIDE");      
1475                                                  
1476   // G4_TOLUENE  index=263                       
1477   G4double M263[NDENSARRAY] = {                  
1478     19.764, 1.88, 3.3026, 0.1722, 2.5728, 0.1    
1479   AddMaterial(M263, "G4_TOLUENE");               
1480                                                  
1481   // G4_TRICHLOROETHYLENE  index=264             
1482   G4double M264[NDENSARRAY] = {                  
1483     24.301, 1.789, 4.6148, 0.1803, 2.914, 0.1    
1484   AddMaterial(M264, "G4_TRICHLOROETHYLENE");     
1485                                                  
1486   // G4_TRIETHYL_PHOSPHATE  index=265            
1487   G4double M265[NDENSARRAY] = {                  
1488     21.863, 2.1, 3.6242, 0.2054, 2.9428, 0.06    
1489   AddMaterial(M265, "G4_TRIETHYL_PHOSPHATE");    
1490                                                  
1491   // G4_TUNGSTEN_HEXAFLUORIDE  index=266         
1492   G4double M266[NDENSARRAY] = {                  
1493     29.265, 2.325, 5.9881, 0.302, 4.2602, 0.0    
1494   AddMaterial(M266, "G4_TUNGSTEN_HEXAFLUORIDE    
1495                                                  
1496   // G4_URANIUM_DICARBIDE  index=267             
1497   G4double M267[NDENSARRAY] = {                  
1498     60.969, 1.703, 6.0247, -0.2191, 3.5208, 0    
1499   AddMaterial(M267, "G4_URANIUM_DICARBIDE");     
1500                                                  
1501   // G4_URANIUM_MONOCARBIDE  index=268           
1502   G4double M268[NDENSARRAY] = {                  
1503     66.602, 1.68, 6.121, -0.2524, 3.4941, 0.2    
1504   AddMaterial(M268, "G4_URANIUM_MONOCARBIDE")    
1505                                                  
1506   // G4_URANIUM_OXIDE  index=269                 
1507   G4double M269[NDENSARRAY] = {                  
1508     60.332, 1.76, 5.9605, -0.1938, 3.5292, 0.    
1509   AddMaterial(M269, "G4_URANIUM_OXIDE");         
1510                                                  
1511   // G4_UREA  index=270                          
1512   G4double M270[NDENSARRAY] = {                  
1513     24.194, 2.022, 3.2032, 0.1603, 2.6525, 0.    
1514   AddMaterial(M270, "G4_UREA");                  
1515                                                  
1516   // G4_VALINE  index=271                        
1517   G4double M271[NDENSARRAY] = {                  
1518     23.622, 2.024, 3.1059, 0.1441, 2.6227, 0.    
1519   AddMaterial(M271, "G4_VALINE");                
1520                                                  
1521   // G4_VITON  index=272                         
1522   G4double M272[NDENSARRAY] = {                  
1523     26.948, 2.227, 3.5943, 0.2106, 2.7874, 0.    
1524   AddMaterial(M272, "G4_VITON");                 
1525                                                  
1526   // G4_WATER  index=273                         
1527   G4double M273[NDENSARRAY] = {                  
1528     21.469, 2.203, 3.5017, 0.24, 2.8004, 0.09    
1529   AddMaterial(M273, "G4_WATER");                 
1530                                                  
1531   // G4_WATER_VAPOR  index=274                   
1532   G4double M274[NDENSARRAY] = {                  
1533     0.59, 2.175, 10.5962, 1.7952, 4.3437, 0.0    
1534   AddMaterial(M274, "G4_WATER_VAPOR");           
1535                                                  
1536   // G4_XYLENE  index=275                        
1537   G4double M275[NDENSARRAY] = {                  
1538     19.866, 1.882, 3.2698, 0.1695, 2.5675, 0.    
1539   AddMaterial(M275, "G4_XYLENE");                
1540                                                  
1541   // G4_GRAPHITE  index=276                      
1542   G4double M276[NDENSARRAY] = {                  
1543     30.652, 2.29, 2.868, -0.0178, 2.3415, 0.2    
1544   AddMaterial(M276, "G4_GRAPHITE");              
1545                                                  
1546   // G4_GRAPHITE_POROUS  index=277               
1547   G4double M277[NDENSARRAY] = {                  
1548     26.555, 2.49, 3.155, 0.048, 2.5387, 0.207    
1549   AddMaterial(M277, "G4_GRAPHITE_POROUS");       
1550 }                                                
1551                                                  
1552 G4int G4DensityEffectData::GetElementIndex(G4    
1553 {                                                
1554   return (Z >= 0 && Z < NDENSELEM) ? indexZ[Z    
1555 }                                                
1556                                                  
1557 G4int G4DensityEffectData::GetIndex(const G4S    
1558 {                                                
1559   G4int idx = -1;                                
1560                                                  
1561   for (G4int i = 0; i < NDENSDATA; ++i) {        
1562     if (names[i] == matName) {                   
1563       idx = i;                                   
1564       break;                                     
1565     }                                            
1566   }                                              
1567   return idx;                                    
1568 }                                                
1569                                                  
1570 void G4DensityEffectData::AddMaterial(G4doubl    
1571 {                                                
1572   for (G4int i = 0; i < NDENSARRAY; ++i) {       
1573     data[index][i] = val[i];                     
1574   }                                              
1575   data[index][0] *= CLHEP::eV;                   
1576   data[index][9] *= CLHEP::eV;                   
1577   names.push_back(matName);                      
1578   ++index;                                       
1579 }                                                
1580                                                  
1581 void G4DensityEffectData::PrintData(const G4S    
1582 {                                                
1583   if (matName.empty() || "all" == matName) {     
1584     DumpData();                                  
1585     return;                                      
1586   }                                              
1587   G4int idx = GetIndex(matName);                 
1588   if (idx >= 0) {                                
1589     G4cout << "G4DensityEffectData for <" <<     
1590     G4cout << "I(eV)= " << data[idx][9] / CLH    
1591            << " rho= " << data[idx][1] << " -    
1592            << " x1= " << data[idx][4] << " a=    
1593            << " d0= " << data[idx][7] << " er    
1594   }                                              
1595   else {                                         
1596     G4cout << "G4DensityEffectData does not h    
1597   }                                              
1598 }                                                
1599                                                  
1600 void G4DensityEffectData::DumpData() const       
1601 {                                                
1602   G4cout << "================================    
1603   G4cout << "     Material        Eplasma(eV)    
1604   G4cout << "================================    
1605   for (G4int i = 0; i < NDENSDATA; ++i) {        
1606     G4cout << std::setw(3) << i << ". " << st    
1607            << data[i][0] / eV;                   
1608     for (G4int j = 1; j < NDENSARRAY; ++j) {     
1609       G4cout << std::setw(8) << data[i][j];      
1610     }                                            
1611     G4cout << G4endl;                            
1612   }                                              
1613   G4cout << "================================    
1614 }                                                
1615