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 10.0.p4)


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