Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/materials/src/G4NistMaterialBuilder.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/G4NistMaterialBuilder.cc (Version 11.3.0) and /materials/src/G4NistMaterialBuilder.cc (Version 9.6.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 // $Id: G4NistMaterialBuilder.cc 67044 2013-01-30 08:50:06Z gcosmo $
                                                   >>  27 //
 26 // -------------------------------------------     28 // -------------------------------------------------------------------
 27 //                                                 29 //
 28 // GEANT4 Class file                               30 // GEANT4 Class file
 29 //                                                 31 //
 30 // File name:     G4NistMaterialBuilder            32 // File name:     G4NistMaterialBuilder
 31 //                                                 33 //
 32 // Author:        Vladimir Ivanchenko              34 // Author:        Vladimir Ivanchenko
 33 //                                                 35 //
 34 // Creation date: 23.12.2004                       36 // Creation date: 23.12.2004
 35 //                                                 37 //
 36 // Modifications:                                  38 // Modifications:
 37 // 31-10-05 Add chemical effect and gas proper     39 // 31-10-05 Add chemical effect and gas properties (V.Ivanchenko)
 38 // 27.02.06 V.Ivanchenko add ConstructNewGasMa <<  40 // 27.02.06 V.Ivanchneko add ConstructNewGasMaterial
 39 // 11.05.06 V.Ivanchenko add warning flag to F <<  41 // 11.05.06 V.Ivanchneko add warning flag to FindMaterial method
 40 // 27.06.06 V.Ivanchenko fix graphite descript <<  42 // 27.06.06 V.Ivanchneko fix graphite description
 41 // 27.07.07 V.Ivanchenko remove dependence on  <<  43 // 27.07.07 V.Ivanchneko remove dependence on NistManager
 42 // 30.10.09 V.Ivanchenko update density of G4_ <<  44 // 30.10.09 V.Ivanchneko update density of G4_GRAFITE from PDG'2008
 43 //                       added G4_GRAPHITE_POR     45 //                       added G4_GRAPHITE_POROUS
 44 // 03.11.09 A.Lechner changed following materi     46 // 03.11.09 A.Lechner changed following material names:
 45 //                    From G4_NYLON-6/6 to G4_     47 //                    From G4_NYLON-6/6 to G4_NYLON-6-6
 46 //                    From G4_NYLON-6/10 to G4     48 //                    From G4_NYLON-6/10 to G4_NYLON-6-10
 47 // 12.12.10 A.Ivantchenko added following mate     49 // 12.12.10 A.Ivantchenko added following materials methodes:
 48 //                    BioChemicalMaterials() a     50 //                    BioChemicalMaterials() and SpaceMaterials(),
 49 //                    where new materials are      51 //                    where new materials are introduced
 50 // 14.06.11 A.Ivantchenko updated body materia     52 // 14.06.11 A.Ivantchenko updated body materials (G4_....ICRP)
 51 //                    according ICRU Report 46 <<  53 //                    according ICRU Report 46 (1992) instead of 1975 
 52 //                    data from ICRU Report 37     54 //                    data from ICRU Report 37 used previously
 53 // 26.10.11 new scheme for G4Exception  (mma)      55 // 26.10.11 new scheme for G4Exception  (mma)
 54 // 09.02.12 P.Gumplinger add ConstructNewIdeal     56 // 09.02.12 P.Gumplinger add ConstructNewIdealGasMaterial
 55 // 30.04.13 M.Trocme & S.Seltzer:              <<  57 //
 56 //        - Replace AddElementByWeightFraction << 
 57 //          as much as possible                << 
 58 //        - Comment out ill-defined material G << 
 59 //        - Fixed density and atom composition << 
 60 //          POLYVINYL_BUTYRAL, TERPHENYL       << 
 61 // -------------------------------------------     58 // -------------------------------------------------------------------
 62 //                                                 59 //
 63 // Class Description:                              60 // Class Description:
 64 //                                                 61 //
 65 // Element data from the NIST DB on Atomic Wei     62 // Element data from the NIST DB on Atomic Weights and Isotope Compositions
 66 // http://physics.nist.gov/PhysRefData/Composi     63 // http://physics.nist.gov/PhysRefData/Compositions/index.html
                                                   >>  64 //
                                                   >>  65 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
                                                   >>  66 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 67                                                    67 
 68 #include "G4NistMaterialBuilder.hh"                68 #include "G4NistMaterialBuilder.hh"
 69                                                << 
 70 #include "G4ApplicationState.hh"               << 
 71 #include "G4AutoLock.hh"                       << 
 72 #include "G4Element.hh"                        << 
 73 #include "G4NistElementBuilder.hh"                 69 #include "G4NistElementBuilder.hh"
                                                   >>  70 #include "G4Element.hh"
 74 #include "G4PhysicalConstants.hh"                  71 #include "G4PhysicalConstants.hh"
 75 #include "G4StateManager.hh"                   << 
 76 #include "G4SystemOfUnits.hh"                      72 #include "G4SystemOfUnits.hh"
 77                                                    73 
 78 #include <iomanip>                             << 
 79                                                << 
 80 namespace                                      << 
 81 {                                              << 
 82 G4Mutex nistMaterialMutex = G4MUTEX_INITIALIZE << 
 83 }                                              << 
 84                                                << 
 85 //....oooOO0OOooo........oooOO0OOooo........oo     74 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 86                                                    75 
 87 G4NistMaterialBuilder::G4NistMaterialBuilder(G     76 G4NistMaterialBuilder::G4NistMaterialBuilder(G4NistElementBuilder* eb, G4int vb)
 88   : elmBuilder(eb), verbose(vb)                <<  77 : elmBuilder(eb),
                                                   >>  78   verbose(vb),
                                                   >>  79   nMaterials(0),
                                                   >>  80   nComponents(0),
                                                   >>  81   nCurrent(0),
                                                   >>  82   first(true)
 89 {                                                  83 {
 90   Initialise();                                    84   Initialise();
 91 }                                                  85 }
 92                                                    86 
 93 //....oooOO0OOooo........oooOO0OOooo........oo     87 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 94                                                    88 
 95 G4Material* G4NistMaterialBuilder::FindOrBuild <<  89 G4NistMaterialBuilder::~G4NistMaterialBuilder()
                                                   >>  90 {}
                                                   >>  91 
                                                   >>  92 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
                                                   >>  93 
                                                   >>  94 G4Material* G4NistMaterialBuilder::FindOrBuildMaterial(const G4String& matname,
                                                   >>  95                                                        G4bool isotopes,
                                                   >>  96                    G4bool warning)
 96 {                                                  97 {
 97   if (verbose > 1) {                           <<  98   if(first) {
 98     G4cout << "G4NistMaterialBuilder::FindOrBu <<  99     if(verbose > 0) {
                                                   >> 100       G4cout << "### NIST DataBase for Materials is used" << G4endl;
                                                   >> 101     }
                                                   >> 102     first = false;
 99   }                                               103   }
                                                   >> 104 
100   G4String name = matname;                        105   G4String name = matname;
101   if ("G4_NYLON-6/6" == matname) {             << 106   if("G4_NYLON-6/6" == matname)  { name = "G4_NYLON-6-6"; }
102     name = "G4_NYLON-6-6";                     << 107   if("G4_NYLON-6/10" == matname) { name = "G4_NYLON-6-10";}
103   }                                            << 
104   else if (name == "G4_NYLON-6/10") {          << 
105     name = "G4_NYLON-6-10";                    << 
106   }                                            << 
107                                                   108 
108   G4Material* mat = FindMaterial(name);        << 109   if (verbose > 1) {
109   if (mat != nullptr) {                        << 110     G4cout << "G4NistMaterialBuilder::FindOrBuildMaterial " << name << G4endl;
110     return mat;                                << 
111   }                                               111   }
                                                   >> 112   const G4MaterialTable* theMaterialTable = G4Material::GetMaterialTable();
                                                   >> 113   G4int nmat = theMaterialTable->size();
112                                                   114 
113   mat = BuildNistMaterial(name, warning);      << 115   // Check if name inside DB
114   return mat;                                  << 116   G4Material* mat = 0;
115 }                                              << 
116                                                   117 
117 //....oooOO0OOooo........oooOO0OOooo........oo << 118   for (G4int i=0; i<nMaterials; ++i) {
118                                                   119 
119 G4Material* G4NistMaterialBuilder::BuildNistMa << 
120 {                                              << 
121   G4Material* mat = nullptr;                   << 
122   // Check if name inside DB                   << 
123   for (G4int i = 0; i < nMaterials; ++i) {     << 
124     if (name == names[i]) {                       120     if (name == names[i]) {
125       if (matIndex[i] == -1) {                 << 121       // Build new Nist material 
126         // Build new Nist material             << 122       if(matIndex[i] == -1) { mat = BuildMaterial(i, isotopes); }
127         mat = BuildMaterial(i);                << 123       // Nist material was already built
128       }                                        << 124       else                  { mat = (*theMaterialTable)[matIndex[i]]; }
129       else {                                   << 
130         // Nist material was already built     << 
131         const G4MaterialTable* theMaterialTabl << 
132         mat = (*theMaterialTable)[matIndex[i]] << 
133       }                                        << 
134       return mat;                                 125       return mat;
135     }                                             126     }
136   }                                               127   }
137                                                   128 
138   if ((verbose == 1 && warning) || verbose > 1 << 129   // Check the list of all materials
139     G4cout << "G4NistMaterialBuilder::FindOrBu << 130   if (nmat > 0) {
140            << " material <" << name << "> is n << 131     for (G4int i=0; i<nmat; ++i) {
                                                   >> 132       if(name == ((*theMaterialTable)[i])->GetName()) {
                                                   >> 133         mat = (*theMaterialTable)[i];
                                                   >> 134   return mat;
                                                   >> 135       }
                                                   >> 136     }
141   }                                               137   }
142   return mat;                                  << 
143 }                                              << 
144                                                << 
145 //....oooOO0OOooo........oooOO0OOooo........oo << 
146                                                   138 
147 G4Material* G4NistMaterialBuilder::FindOrBuild << 139   if( (verbose == 1 && warning) || verbose > 1) {
148 {                                              << 140     G4cout << "G4NistMaterialBuilder::FindOrBuildMaterial WARNING:"
149   G4Material* mat = FindSimpleMaterial(Z);     << 141      << " material <" << name
150   if (mat == nullptr) {                        << 142      << "> is not found out" << G4endl;
151     mat = BuildNistMaterial(names[Z], warn);   << 
152   }                                               143   }
153   return mat;                                     144   return mat;
154 }                                                 145 }
155                                                   146 
156 //....oooOO0OOooo........oooOO0OOooo........oo    147 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
157                                                   148 
158 G4Material* G4NistMaterialBuilder::BuildMateri << 149 G4Material* G4NistMaterialBuilder::BuildMaterial(G4int i, G4bool isotopes)
159 {                                                 150 {
160   G4Material* mat = nullptr;                   << 151   if (verbose > 1) {
161   if (i >= nMaterials) {                       << 152     G4cout << "G4NistMaterialBuilder: BuildMaterial #" << i
162     return mat;                                << 153      << G4endl;
163   }                                               154   }
                                                   >> 155   G4Material* mat = 0;
                                                   >> 156   if (nMaterials == 0) { return mat; }
164                                                   157 
165   G4AutoLock l(&nistMaterialMutex);            << 158   G4int nc = components[i];
166   if (matIndex[i] >= 0) {                      << 
167     // Nist material was already built         << 
168     const G4MaterialTable* theMaterialTable =  << 
169     mat = (*theMaterialTable)[matIndex[i]];    << 
170   }                                            << 
171   else {                                       << 
172     if (verbose > 1) {                         << 
173       G4cout << "G4NistMaterialBuilder: BuildM << 
174     }                                          << 
175                                                << 
176     G4int nc = components[i];                  << 
177                                                   159 
178     // Check gas parameters:                   << 160   // Check gas parameters
179     // defaults may be changed via AddGas() me << 161   G4double t = STP_Temperature;
180     G4double t = NTP_Temperature;              << 162   G4double p = STP_Pressure;
181     G4double p = CLHEP::STP_Pressure;          << 163   if(kStateGas == states[i]) {
182     if (kStateGas == states[i]) {              << 164     size_t nn = idxGas.size();
183       size_t nn = idxGas.size();               << 165     if(nn > 0) {
184       if (nn > 0) {                            << 166       for(size_t j=0; j<nn; ++j) {
185         for (size_t j = 0; j < nn; ++j) {      << 167         if(i == idxGas[j]) {
186           if (i == idxGas[j]) {                << 168     t = gasTemperature[j];
187             t = gasTemperature[j];             << 169           p = gasPressure[j];
188             p = gasPressure[j];                << 170           break;
189             break;                             << 171   }
190           }                                    << 
191         }                                      << 
192       }                                           172       }
193     }                                             173     }
194     mat = new G4Material(names[i], densities[i << 174     // liquids
                                                   >> 175   } else if( !STP[i] ) { t = 0.0; }
195                                                   176 
196     if (verbose > 1) {                         << 177   mat = new G4Material(names[i],densities[i],nc,states[i],t,p);
197       G4cout << "New material nComponents= " < << 
198     }                                          << 
199     if (nc > 0) {                              << 
200       G4int idx = indexes[i];                  << 
201       for (G4int j = 0; j < nc; ++j) {         << 
202         G4int Z = elements[idx + j];           << 
203         G4Element* el = elmBuilder->FindOrBuil << 
204         if (el == nullptr) {                   << 
205           G4cout << "G4NistMaterialBuilder::Bu << 
206                  << "  ERROR: elements Z= " << << 
207                  << " for material " << names[ << 
208           G4Exception("G4NistMaterialBuilder:: << 
209             "Failed to construct material");   << 
210           return nullptr;                      << 
211         }                                      << 
212         if (atomCount[i]) {                    << 
213           mat->AddElement(el, G4lrint(fraction << 
214         }                                      << 
215         else {                                 << 
216           mat->AddElement(el, fractions[idx +  << 
217         }                                      << 
218       }                                        << 
219     }                                          << 
220                                                   178 
221     // Ionisation potential can be defined via << 179   if (verbose>1) { G4cout << "New material nComponents= " << nc << G4endl; }
222     // Chemical Formula (ICRU37 Report data)   << 180   if (nc > 0) {
223     G4IonisParamMat* ion = mat->GetIonisation( << 181     G4int idx = indexes[i];
224     G4double exc0 = ion->GetMeanExcitationEner << 182     for (G4int j=0; j<nc; ++j) {
225     G4double exc1 = exc0;                      << 183       G4int Z = elements[idx+j];
226     if (! chFormulas[i].empty()) {             << 184       G4Element* el = elmBuilder->FindOrBuildElement(Z, isotopes);
227       mat->SetChemicalFormula(chFormulas[i]);  << 185       if(!el) {
228       exc1 = ion->FindMeanExcitationEnergy(mat << 186   G4cout << "G4NistMaterialBuilder::BuildMaterial:"
229     }                                          << 187          << "  ERROR: elements Z= " << Z << " is not found "
230     // If exists, NIST DB data always overwrit << 188          << " for material " << names[i]
231     if (ionPotentials[i] > 0.0) {              << 189          << G4endl;
232       exc1 = ionPotentials[i];                 << 190   G4Exception("G4NistMaterialBuilder::BuildMaterial()", "mat103",
233     }                                          << 191                FatalException, "Fail to construct material");
234     if (exc0 != exc1) {                        << 192   return 0;
235       ion->SetMeanExcitationEnergy(exc1);      << 193       }
                                                   >> 194       if(atomCount[i]) {
                                                   >> 195   mat->AddElement(el,G4lrint(fractions[idx+j]));
                                                   >> 196       } else {
                                                   >> 197   mat->AddElement(el,fractions[idx+j]);
                                                   >> 198       }
236     }                                             199     }
                                                   >> 200   }
237                                                   201 
238     // Index in Material Table                 << 202   // Ionisation potential can be defined via NIST DB or 
239     matIndex[i] = (G4int)mat->GetIndex();      << 203   // Chemical Formula (ICRU37 Report data)
                                                   >> 204   G4IonisParamMat* ion = mat->GetIonisation();
                                                   >> 205   G4double exc0 = ion->GetMeanExcitationEnergy();
                                                   >> 206   G4double exc1 = exc0;
                                                   >> 207   if(chFormulas[i] != "") {
                                                   >> 208     mat->SetChemicalFormula(chFormulas[i]);
                                                   >> 209     exc1 = ion->FindMeanExcitationEnergy(chFormulas[i]);
240   }                                               210   }
241   l.unlock();                                  << 211   if(ionPotentials[i] > 0.0) { exc1 = ionPotentials[i]; }
                                                   >> 212   if(exc0 != exc1) { ion->SetMeanExcitationEnergy(exc1); }
                                                   >> 213 
                                                   >> 214   // Index in Material Table
                                                   >> 215   matIndex[i] = mat->GetIndex();
242   return mat;                                     216   return mat;
243 }                                                 217 }
244                                                   218 
245 //....oooOO0OOooo........oooOO0OOooo........oo    219 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
246                                                   220 
247 G4Material* G4NistMaterialBuilder::ConstructNe << 221 G4Material* G4NistMaterialBuilder::ConstructNewMaterial(
248   const std::vector<G4String>& elm, const std: << 222                                       const G4String& name,
249   G4double temp, G4double pres)                << 223                                       const std::vector<G4String>& elm,
                                                   >> 224                                       const std::vector<G4int>& nbAtoms,
                                                   >> 225               G4double dens, 
                                                   >> 226               G4bool isotopes,
                                                   >> 227               G4State state,     
                                                   >> 228               G4double temp,  
                                                   >> 229               G4double pres)
250 {                                                 230 {
251   // Material is in DB                            231   // Material is in DB
252   G4Material* mat = FindOrBuildMaterial(name);    232   G4Material* mat = FindOrBuildMaterial(name);
253   if (mat != nullptr) {                        << 233   if(mat) { 
254     G4cout << "G4NistMaterialBuilder::Construc    234     G4cout << "G4NistMaterialBuilder::ConstructNewMaterial:"
255            << "  WARNING: the material <" << n << 235            << "  WARNING: the material <" << name
256     G4cout << "      New material will NOT be  << 236      << "> is already exist" << G4endl;
257     return mat;                                << 237     G4cout << "      New material will NOT be built!"
                                                   >> 238      << G4endl;
                                                   >> 239     return mat; 
258   }                                               240   }
259                                                   241 
260   // Material not in DB                           242   // Material not in DB
261   auto els = (G4int)elm.size();                << 243   G4int els = elm.size();
262   if (els == 0) {                              << 244   if(els == 0) { 
263     G4cout << "G4NistMaterialBuilder::Construc    245     G4cout << "G4NistMaterialBuilder::ConstructNewMaterial:"
264            << "  WARNING: empty list of elemen << 246            << "  WARNING: empty list of elements for " << name
265     G4cout << "      New material will NOT be  << 247      << G4endl;
266     return nullptr;                            << 248     G4cout << "      New material will NOT be built!"
267   }                                            << 249      << G4endl;
                                                   >> 250     return 0;
                                                   >> 251   } 
268                                                   252 
269   // add parameters of material into internal     253   // add parameters of material into internal vectors
270   // density in g/cm3, mean ionisation potenti    254   // density in g/cm3, mean ionisation potential is not defined
271   G4bool stp = true;                              255   G4bool stp = true;
272   if (state == kStateGas && (temp != NTP_Tempe << 256   if(state == kStateGas && temp != STP_Temperature && pres != STP_Pressure)
273     stp = false;                               << 257     { stp = false; }
274   }                                            << 
275                                                   258 
276   AddMaterial(name, dens * CLHEP::cm3 / CLHEP: << 259   AddMaterial(name,dens*cm3/g,0,0.,els,state,stp);
277   if (! stp) {                                 << 260   if(!stp) { AddGas(name,temp,pres); }
278     AddGas(name, temp, pres);                  << 
279   }                                            << 
280                                                   261 
281   for (G4int i = 0; i < els; ++i) {            << 262   for (G4int i=0; i<els; ++i) {
282     AddElementByAtomCount(elmBuilder->GetZ(elm    263     AddElementByAtomCount(elmBuilder->GetZ(elm[i]), nbAtoms[i]);
283   }                                               264   }
284                                                   265 
285   return BuildMaterial(nMaterials - 1);        << 266   return BuildMaterial(nMaterials-1, isotopes);
286 }                                                 267 }
287                                                   268 
288 //....oooOO0OOooo........oooOO0OOooo........oo    269 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
289                                                   270 
290 G4Material* G4NistMaterialBuilder::ConstructNe << 271 G4Material* G4NistMaterialBuilder::ConstructNewMaterial(
291   const std::vector<G4String>& elm, const std: << 272                                       const G4String& name,
292   G4double temp, G4double pres)                << 273                                       const std::vector<G4String>& elm,
                                                   >> 274                                       const std::vector<G4double>& w,
                                                   >> 275               G4double dens, 
                                                   >> 276               G4bool isotopes,
                                                   >> 277               G4State state,     
                                                   >> 278               G4double temp,  
                                                   >> 279               G4double pres)
293 {                                                 280 {
294   // Material is in DB                            281   // Material is in DB
295   G4Material* mat = FindOrBuildMaterial(name);    282   G4Material* mat = FindOrBuildMaterial(name);
296   if (mat != nullptr) {                        << 283   if(mat) { 
297     G4cout << "G4NistMaterialBuilder::Construc    284     G4cout << "G4NistMaterialBuilder::ConstructNewMaterial:"
298            << "  WARNING: the material <" << n << 285            << "  WARNING: the material <" << name
299     G4cout << "      New material will NOT be  << 286      << "> is already exist" << G4endl;
300     return mat;                                << 287     G4cout << "      New material will NOT be built!"
                                                   >> 288      << G4endl;
                                                   >> 289     return mat; 
301   }                                               290   }
302                                                   291 
303   // Material not in DB                           292   // Material not in DB
304   auto els = (G4int)elm.size();                << 293   G4int els = elm.size();
305   if (els == 0) {                              << 294   if(els == 0) { 
306     G4cout << "G4NistMaterialBuilder::Construc    295     G4cout << "G4NistMaterialBuilder::ConstructNewMaterial:"
307            << "  WARNING: empty list of elemen << 296            << "  WARNING: empty list of elements for " << name
308     G4cout << "      New material will NOT be  << 297      << G4endl;
309     return nullptr;                            << 298     G4cout << "      New material will NOT be built!"
310   }                                            << 299      << G4endl;
                                                   >> 300     return 0;
                                                   >> 301   } 
311                                                   302 
312   // add parameters of material into internal     303   // add parameters of material into internal vectors
313   // density in g/cm3, mean ionisation potenti    304   // density in g/cm3, mean ionisation potential is not defined
314   G4bool stp = true;                              305   G4bool stp = true;
315   if (state == kStateGas && (temp != NTP_Tempe << 306   if(state == kStateGas && temp != STP_Temperature && pres != STP_Pressure)
316     stp = false;                               << 307     { stp = false; }
317   }                                            << 308   AddMaterial(name,dens*cm3/g,0,0.,els,state,stp);
318   AddMaterial(name, dens * CLHEP::cm3 / CLHEP: << 309   if(!stp) { AddGas(name,temp,pres); }
319   if (! stp) {                                 << 
320     AddGas(name, temp, pres);                  << 
321   }                                            << 
322                                                   310 
323   for (G4int i = 0; i < els; ++i) {            << 311   for (G4int i=0; i<els; ++i) {
324     AddElementByWeightFraction(elmBuilder->Get    312     AddElementByWeightFraction(elmBuilder->GetZ(elm[i]), w[i]);
325   }                                               313   }
326                                                << 314   
327   return BuildMaterial(nMaterials - 1);        << 315   return BuildMaterial(nMaterials-1, isotopes);    
328 }                                                 316 }
329                                                   317 
330 //....oooOO0OOooo........oooOO0OOooo........oo    318 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
331                                                   319 
332 G4Material* G4NistMaterialBuilder::ConstructNe    320 G4Material* G4NistMaterialBuilder::ConstructNewGasMaterial(
333   const G4String& name, const G4String& nameDB << 321               const G4String& name,
                                                   >> 322               const G4String& nameDB,
                                                   >> 323               G4double temp, 
                                                   >> 324               G4double pres, 
                                                   >> 325               G4bool)
334 {                                                 326 {
335   // Material name is in DB                       327   // Material name is in DB
336   G4Material* mat = FindOrBuildMaterial(name);    328   G4Material* mat = FindOrBuildMaterial(name);
337   if (mat != nullptr) {                        << 329   if(mat) { 
338     G4cout << "G4NistMaterialBuilder::Construc    330     G4cout << "G4NistMaterialBuilder::ConstructNewGasMaterial:"
339            << "  WARNING: the material <" << n << 331            << "  WARNING: the material <" << name
340     G4cout << "      New material will NOT be  << 332      << "> is already exist" << G4endl;
341     return mat;                                << 333     G4cout << "      New material will NOT be built!"
                                                   >> 334      << G4endl;
                                                   >> 335     return mat; 
342   }                                               336   }
343                                                   337 
344   G4Material* bmat = FindOrBuildMaterial(nameD    338   G4Material* bmat = FindOrBuildMaterial(nameDB);
345   if (bmat == nullptr) {                       << 339   if(!bmat) {
346     G4cout << "G4NistMaterialBuilder::Construc    340     G4cout << "G4NistMaterialBuilder::ConstructNewGasMaterial:"
347            << "  WARNING: the Name <" << nameD << 341      << "  WARNING: the Name <" << nameDB 
348            << "> is NOT in the database: no ne << 342      << "> is NOT in the DB: no new gas will be constructed"
349     return nullptr;                            << 343      << G4endl;
                                                   >> 344     return 0;
350   }                                               345   }
351   if (bmat->GetState() != kStateGas) {         << 346   if(bmat->GetState() != kStateGas) {
352     G4cout << "G4NistMaterialBuilder::Construc    347     G4cout << "G4NistMaterialBuilder::ConstructNewGasMaterial:"
353            << "  WARNING:  <" << nameDB << ">  << 348      << "  WARNING:  <" << nameDB 
354            << G4endl;                          << 349      << "> is NOT a gas -  no new gas will be constructed"
355     return nullptr;                            << 350      << G4endl;
                                                   >> 351     return 0;
356   }                                               352   }
357                                                   353 
358   G4double dens = bmat->GetDensity() * pres *  << 354   G4double dens = bmat->GetDensity()*pres*STP_Temperature/(temp*STP_Pressure);
359   mat = new G4Material(name, dens, bmat, kStat << 355   mat = new G4Material(name,dens,bmat,kStateGas,temp,pres);
360                                                   356 
361   if (verbose > 1) {                           << 357   if (verbose>1) {
362     G4cout << "G4NistMaterialBuilder::Construc    358     G4cout << "G4NistMaterialBuilder::ConstructNewGasMaterial: done" << G4endl;
363     G4cout << &mat << G4endl;                  << 359     G4cout << &mat << G4endl; 
364   }                                            << 360   } 
365   return mat;                                     361   return mat;
366 }                                                 362 }
367                                                   363 
368 //....oooOO0OOooo........oooOO0OOooo........oo    364 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
369                                                   365 
370 G4Material* G4NistMaterialBuilder::ConstructNe << 366 G4Material* G4NistMaterialBuilder::ConstructNewIdealGasMaterial(
371   const std::vector<G4String>& elm, const std: << 367                                       const G4String& name,
                                                   >> 368                                       const std::vector<G4String>& elm,
                                                   >> 369                                       const std::vector<G4int>& nbAtoms,
                                                   >> 370                                       G4bool isotopes,
                                                   >> 371                                       G4double temp,
                                                   >> 372                                       G4double pres)
372 {                                                 373 {
373   G4State state = kStateGas;                      374   G4State state = kStateGas;
374                                                   375 
375   // Material is in DB                            376   // Material is in DB
376   G4Material* mat = FindOrBuildMaterial(name);    377   G4Material* mat = FindOrBuildMaterial(name);
377   if (mat != nullptr) {                        << 378   if(mat) {
378     G4cout << "G4NistMaterialBuilder::Construc    379     G4cout << "G4NistMaterialBuilder::ConstructNewMaterial:"
379            << "  WARNING: the material <" << n << 380            << "  WARNING: the material <" << name
380     G4cout << "      New material will NOT be  << 381            << "> is already exist" << G4endl;
                                                   >> 382     G4cout << "      New material will NOT be built!"
                                                   >> 383            << G4endl;
381     return mat;                                   384     return mat;
382   }                                               385   }
383                                                   386 
384   // Material not in DB                           387   // Material not in DB
385   auto els = (G4int)elm.size();                << 388   G4int els = elm.size();
386   if (els == 0) {                              << 389   if(els == 0) {
387     G4cout << "G4NistMaterialBuilder::Construc    390     G4cout << "G4NistMaterialBuilder::ConstructNewMaterial:"
388            << "  WARNING: empty list of elemen << 391            << "  WARNING: empty list of elements for " << name
389     G4cout << "      New material will NOT be  << 392            << G4endl;
390     return nullptr;                            << 393     G4cout << "      New material will NOT be built!"
                                                   >> 394            << G4endl;
                                                   >> 395     return 0;
391   }                                               396   }
392                                                   397 
393   // add parameters of material into internal     398   // add parameters of material into internal vectors
394   // density in g/cm3, mean ionisation potenti    399   // density in g/cm3, mean ionisation potential is not defined
395   G4bool stp = true;                              400   G4bool stp = true;
396   if (temp != NTP_Temperature || pres != CLHEP << 401   if(temp != STP_Temperature && pres != STP_Pressure)
397     stp = false;                               << 402     { stp = false; }
398   }                                            << 
399                                                   403 
400   G4double massPerMole = 0.;                   << 404   G4double massPerMole = 0;
401                                                   405 
402   G4int Z = 0;                                    406   G4int Z = 0;
403   for (G4int i = 0; i < els; ++i) {            << 407   for (G4int i=0; i<els; ++i) {
404     Z = elmBuilder->GetZ(elm[i]);                 408     Z = elmBuilder->GetZ(elm[i]);
405     massPerMole += nbAtoms[i] * elmBuilder->Ge << 409     massPerMole += nbAtoms[i] * elmBuilder->GetAtomicMassAmu(Z) * amu_c2;
406   }                                               410   }
407                                                   411 
408   G4double dens = massPerMole / (CLHEP::Avogad << 412   G4double dens = massPerMole / (Avogadro*k_Boltzmann*temp/pres);
409                                                   413 
410   if (els == 1) {                              << 414   if (els == 1) { AddMaterial(name,dens,Z,0.,els,state,stp); }
411     AddMaterial(name, dens, Z, 0., els, state, << 
412   }                                            << 
413   else {                                          415   else {
414     AddMaterial(name, dens, 0, 0., els, state, << 416     AddMaterial(name,dens,0,0.,els,state,stp);
415     for (G4int i = 0; i < els; ++i) {          << 417     for (G4int i=0; i<els; ++i) {
416       AddElementByAtomCount(elmBuilder->GetZ(e    418       AddElementByAtomCount(elmBuilder->GetZ(elm[i]), nbAtoms[i]);
417     }                                             419     }
418   }                                               420   }
419                                                   421 
420   if (! stp) {                                 << 422   if(!stp) { AddGas(name,temp,pres); }
421     AddGas(name, temp, pres);                  << 
422   }                                            << 
423                                                   423 
424   return BuildMaterial(nMaterials - 1);        << 424   return BuildMaterial(nMaterials-1, isotopes);
425 }                                                 425 }
426                                                   426 
427 //....oooOO0OOooo........oooOO0OOooo........oo    427 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
428                                                   428 
429 void G4NistMaterialBuilder::AddMaterial(const  << 429 void G4NistMaterialBuilder::AddMaterial(const G4String& nameMat, G4double dens,
430   G4double pot, G4int ncomp, G4State state, G4 << 430           G4int Z, G4double pot, 
                                                   >> 431           G4int ncomp, G4State state, 
                                                   >> 432           G4bool stp)
431 {                                                 433 {
432   // add parameters of material into internal     434   // add parameters of material into internal vectors
433   // density in g/cm3, mean ionisation potenti    435   // density in g/cm3, mean ionisation potential in eV
434                                                   436 
435   // if ncomp == 1 then Z should be defined an << 437   // if ncomp == 1 then Z should be defined and 
436   // AddElement should not be applied             438   // AddElement should not be applied
437                                                   439 
438   if (nCurrent != 0) {                            440   if (nCurrent != 0) {
439     G4cout << "G4NistMaterialBuilder::AddMater    441     G4cout << "G4NistMaterialBuilder::AddMaterial WARNING: previous "
440            << "mixture " << nMaterials << " "  << 442      << "mixture " << nMaterials << " " << names[nMaterials] 
441            << G4endl;                          << 443      << " is not yet complete!"
442     G4cout << "         New material " << name << 444      << G4endl;
                                                   >> 445     G4cout << "         New material " << nameMat << " will not be added" 
                                                   >> 446      << G4endl;
443     return;                                       447     return;
444   }                                               448   }
445                                                   449 
446   // density in g/cm3, mean ionisation potenti    450   // density in g/cm3, mean ionisation potential in eV
447                                                   451 
448   names.push_back(nameMat);                       452   names.push_back(nameMat);
449   chFormulas.emplace_back("");                 << 453   chFormulas.push_back("");
450   densities.push_back(dens * CLHEP::g / CLHEP: << 454   densities.push_back(dens*g/cm3);
451   ionPotentials.push_back(pot * CLHEP::eV);    << 455   ionPotentials.push_back(pot*eV);
452   states.push_back(state);                        456   states.push_back(state);
453   components.push_back(ncomp);                    457   components.push_back(ncomp);
454   indexes.push_back(nComponents);                 458   indexes.push_back(nComponents);
455   STP.push_back(stp);                             459   STP.push_back(stp);
456   matIndex.push_back(-1);                         460   matIndex.push_back(-1);
457   atomCount.push_back(false);                     461   atomCount.push_back(false);
458                                                   462 
459   if (1 == ncomp && Z > 0) {                   << 463   if (ncomp == 1 && Z > 0) {
460     elements.push_back(Z);                        464     elements.push_back(Z);
461     fractions.push_back(1.0);                     465     fractions.push_back(1.0);
462     atomCount[nMaterials] = true;                 466     atomCount[nMaterials] = true;
463     ++nComponents;                                467     ++nComponents;
464     nCurrent = 0;                                 468     nCurrent = 0;
465   }                                            << 469   } else {
466   else {                                       << 
467     nCurrent = ncomp;                             470     nCurrent = ncomp;
468   }                                               471   }
469                                                   472 
470   ++nMaterials;                                   473   ++nMaterials;
471                                                   474 
472   if (verbose > 1) {                           << 475   if(verbose > 1) {
473     G4cout << "New material " << nameMat << "  << 476     G4cout << "New material " << nameMat << " is prepeared; "
474            << " nMaterials= " << nMaterials << << 477            << " nMaterials= " << nMaterials
475            << " nCurrent= " << nCurrent << G4e << 478            << " nComponents= " << nComponents
                                                   >> 479            << " nCurrent= " << nCurrent
                                                   >> 480            << G4endl;
476   }                                               481   }
477 }                                                 482 }
478                                                   483 
479 //....oooOO0OOooo........oooOO0OOooo........oo    484 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
480                                                   485 
481 void G4NistMaterialBuilder::SetVerbose(G4int v    486 void G4NistMaterialBuilder::SetVerbose(G4int val)
482 {                                                 487 {
483   verbose = val;                                  488   verbose = val;
484   elmBuilder->SetVerbose(verbose);                489   elmBuilder->SetVerbose(verbose);
485 }                                                 490 }
486                                                   491 
487 //....oooOO0OOooo........oooOO0OOooo........oo    492 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
488                                                   493 
489 void G4NistMaterialBuilder::ListMaterials(cons << 494 void G4NistMaterialBuilder::ListMaterials(const G4String& mnam)
490 {                                                 495 {
491   if (mnam == "simple") {                      << 496   if (mnam == "simple")           { ListNistSimpleMaterials(); }
492     ListNistSimpleMaterials();                 << 497   else if (mnam == "compound")    { ListNistCompoundMaterials(); }
493   }                                            << 498   else if (mnam == "hep")         { ListHepMaterials(); }
494   else if (mnam == "compound") {               << 499   else if (mnam == "space")       { ListSpaceMaterials(); }
495     ListNistCompoundMaterials();               << 500   else if (mnam == "biochemical") { ListBioChemicalMaterials(); }
496   }                                            << 
497   else if (mnam == "hep") {                    << 
498     ListHepMaterials();                        << 
499   }                                            << 
500   else if (mnam == "space") {                  << 
501     ListSpaceMaterials();                      << 
502   }                                            << 
503   else if (mnam == "bio") {                    << 
504     ListBioChemicalMaterials();                << 
505   }                                            << 
506                                                   501 
507   else if (mnam == "all") {                       502   else if (mnam == "all") {
508     ListNistSimpleMaterials();                    503     ListNistSimpleMaterials();
509     ListNistCompoundMaterials();                  504     ListNistCompoundMaterials();
510     ListHepMaterials();                           505     ListHepMaterials();
511     ListSpaceMaterials();                         506     ListSpaceMaterials();
512     ListBioChemicalMaterials();                   507     ListBioChemicalMaterials();
513   }                                            << 508 
514   else {                                       << 509   } else {
515     G4cout << "### G4NistMaterialBuilder::List << 510     G4cout << "### G4NistMaterialBuilder::ListMaterials: Warning " 
516            << G4endl;                          << 511      << mnam << " list is not known" << G4endl;
517   }                                               512   }
518 }                                                 513 }
519                                                   514 
520 //....oooOO0OOooo........oooOO0OOooo........oo    515 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
521                                                   516 
522 void G4NistMaterialBuilder::ListNistSimpleMate << 517 void G4NistMaterialBuilder::ListNistSimpleMaterials()
523 {                                                 518 {
524   G4cout << "=================================    519   G4cout << "=======================================================" << G4endl;
525   G4cout << "###   Simple Materials from the N << 520   G4cout << "###   Simple Materials from the NIST Data Base   ###" << G4endl;
526   G4cout << "=================================    521   G4cout << "=======================================================" << G4endl;
527   G4cout << " Z   Name   density(g/cm^3)  I(eV << 522   G4cout << " Z Name  ChFormula        density(g/cm^3)  I(eV)       " << G4endl;
                                                   >> 523   G4cout << "=======================================================" << G4endl;
                                                   >> 524   for (G4int i=0; i<nElementary; ++i) {DumpElm(i);}
528   G4cout << "=================================    525   G4cout << "=======================================================" << G4endl;
529   for (G4int i = 1; i < nElementary; ++i) {    << 
530     DumpElm(i);                                << 
531   }                                            << 
532 }                                                 526 }
533                                                   527 
534 //....oooOO0OOooo........oooOO0OOooo........oo    528 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
535                                                   529 
536 void G4NistMaterialBuilder::ListNistCompoundMa << 530 void G4NistMaterialBuilder::ListNistCompoundMaterials()
537 {                                                 531 {
538   G4cout << "================================= << 532   G4cout << "###    Compound Materials from the NIST Data Base    ##" << G4endl;
539   G4cout << "###    Compound Materials from th << 533   G4cout << "=======================================================" << G4endl;
540   G4cout << "================================= << 534   G4cout << " Ncomp Name  ChFormula        density(g/cm^3)  I(eV)   " << G4endl;
541   G4cout << " Ncomp             Name      dens << 535   G4cout << "=======================================================" << G4endl;
542   G4cout << "================================= << 536   for (G4int i=nElementary; i<nNIST; ++i) {DumpMix(i);}
543   for (G4int i = nElementary; i < nNIST; ++i)  << 537   G4cout << "=======================================================" << G4endl;
544     DumpMix(i);                                << 
545   }                                            << 
546   DumpMix(0);                                  << 
547 }                                                 538 }
548                                                   539 
549 //....oooOO0OOooo........oooOO0OOooo........oo    540 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
550                                                   541 
551 void G4NistMaterialBuilder::ListHepMaterials() << 542 void G4NistMaterialBuilder::ListHepMaterials()
552 {                                                 543 {
553   G4cout << "================================= << 544   G4cout << "=======================================================" << G4endl;
554   G4cout << "###           HEP & Nuclear Mater << 545   G4cout << "###           HEP & Nuclear Materials                ##" << G4endl;
555   G4cout << "================================= << 546   G4cout << "=======================================================" << G4endl;
556   G4cout << " Ncomp             Name      dens << 547   G4cout << " Ncomp Name  ChFormula        density(g/cm^3)  I(eV)   " << G4endl;
557   G4cout << "================================= << 548   G4cout << "=======================================================" << G4endl;
558   for (G4int i = nNIST; i < nHEP; ++i) {       << 549   for (G4int i=nNIST; i<nHEP; ++i) {DumpMix(i);}
559     DumpMix(i);                                << 550   G4cout << "=======================================================" << G4endl;
560   }                                            << 
561 }                                                 551 }
562                                                   552 
563 //....oooOO0OOooo........oooOO0OOooo........oo    553 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
564                                                   554 
565 void G4NistMaterialBuilder::ListSpaceMaterials << 555 void G4NistMaterialBuilder::ListSpaceMaterials()
566 {                                                 556 {
567   G4cout << "================================= << 557   G4cout << "=======================================================" << G4endl;
568   G4cout << "###           Space ISS Materials << 558   G4cout << "###           Space ISS Materials                    ##" << G4endl;
569   G4cout << "================================= << 559   G4cout << "=======================================================" << G4endl;
570   G4cout << " Ncomp             Name      dens << 560   G4cout << " Ncomp Name  ChFormula        density(g/cm^3)  I(eV)   " << G4endl;
571   G4cout << "================================= << 561   G4cout << "=======================================================" << G4endl;
572   for (G4int i = nHEP; i < nSpace; ++i) {      << 562   for (G4int i=nHEP; i<nSpace; ++i) {DumpMix(i);}
573     DumpMix(i);                                << 563   G4cout << "=======================================================" << G4endl;
574   }                                            << 
575 }                                                 564 }
576                                                   565 
577 //....oooOO0OOooo........oooOO0OOooo........oo    566 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
578                                                   567 
579 void G4NistMaterialBuilder::ListBioChemicalMat << 568 void G4NistMaterialBuilder::ListBioChemicalMaterials()
580 {                                                 569 {
581   G4cout << "================================= << 570   G4cout << "=======================================================" << G4endl;
582   G4cout << "###          Bio-Chemical Materia << 571   G4cout << "###          Bio-Chemical Materials                  ##" << G4endl;
583   G4cout << "================================= << 572   G4cout << "=======================================================" << G4endl;
584   G4cout << " Ncomp             Name      dens << 573   G4cout << " Ncomp Name  ChFormula        density(g/cm^3)  I(eV)   " << G4endl;
585   G4cout << "================================= << 574   G4cout << "=======================================================" << G4endl;
586   for (G4int i = nSpace; i < nMaterials; ++i)  << 575   for (G4int i=nSpace; i<nMaterials; ++i) {DumpMix(i);}
587     DumpMix(i);                                << 576   G4cout << "=======================================================" << G4endl;
588   }                                            << 
589   G4cout << "================================= << 
590 }                                                 577 }
591                                                   578 
592 //....oooOO0OOooo........oooOO0OOooo........oo    579 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
593                                                   580 
594 void G4NistMaterialBuilder::DumpElm(G4int i) c << 581 void G4NistMaterialBuilder::DumpElm(G4int i)
595 {                                                 582 {
596   G4cout << std::setw(2) << i << " " << std::s << 583   G4cout << i+1 << "  " << names[i] << "  " << chFormulas[i]
597          << densities[i] * cm3 / g << std::set << 584          << densities[i]*cm3/g << "     " << ionPotentials[i]/eV
                                                   >> 585    << G4endl;
598 }                                                 586 }
599                                                   587 
600 //....oooOO0OOooo........oooOO0OOooo........oo    588 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
601                                                   589 
602 void G4NistMaterialBuilder::DumpMix(G4int i) c << 590 void G4NistMaterialBuilder::DumpMix(G4int i)
603 {                                                 591 {
604   G4int nc = components[i];                       592   G4int nc = components[i];
605   G4cout << std::setw(2) << nc << " " << std:: << 593   G4cout << nc << "  " << names[i] << "  " << chFormulas[i]
606          << densities[i] * cm3 / g << std::set << 594          << densities[i]*cm3/g << "     " << ionPotentials[i]/eV
607          << chFormulas[i] << G4endl;           << 595    << G4endl;
608   if (nc > 1) {                                   596   if (nc > 1) {
609     G4int imin = indexes[i];                      597     G4int imin = indexes[i];
610     G4int imax = imin + nc;                       598     G4int imax = imin + nc;
611     for (G4int j = imin; j < imax; ++j) {      << 599     for (G4int j=imin; j<imax; ++j) {
612       G4cout << std::setw(10) << elements[j] < << 600       G4cout << "              " << elements[j] << "     " << fractions[j] 
                                                   >> 601              << G4endl;
613     }                                             602     }
614   }                                               603   }
615 }                                                 604 }
616                                                   605 
617 //....oooOO0OOooo........oooOO0OOooo........oo    606 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
618                                                   607 
619 void G4NistMaterialBuilder::AddGas(const G4Str << 608 void 
                                                   >> 609 G4NistMaterialBuilder::AddGas(const G4String& nameMat, G4double t, G4double p)
620 {                                                 610 {
621   for (G4int i = 0; i < nMaterials; ++i) {     << 611   G4int idx = nMaterials-1;
622     if (nameMat == names[i]) {                 << 612   if(nameMat != names[idx]) {
623       idxGas.push_back(i);                     << 613     idx = -1;
624       gasTemperature.push_back(t);             << 614     for(G4int i=0; i<nMaterials; ++i) {
625       gasPressure.push_back(p);                << 615       if(nameMat == names[i]) {
626       return;                                  << 616         idx = i; break;
                                                   >> 617       }
627     }                                             618     }
628   }                                               619   }
629   G4cout << "WARNING: G4NistMaterialBuilder::A << 620   if(idx >= 0) {
630          << " in the list of materials." << G4 << 621     idxGas.push_back(idx);
                                                   >> 622     gasTemperature.push_back(t);
                                                   >> 623     gasPressure.push_back(p);
                                                   >> 624   } else {
                                                   >> 625     G4cout << "WARNING: G4NistMaterialBuilder::AddGas problem: there is no "
                                                   >> 626      << nameMat << " in the list of materials;"
                                                   >> 627      << G4endl;
                                                   >> 628   }
631 }                                                 629 }
632                                                   630 
633 //....oooOO0OOooo........oooOO0OOooo........oo    631 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
634                                                   632 
635 void G4NistMaterialBuilder::AddElementByWeight    633 void G4NistMaterialBuilder::AddElementByWeightFraction(G4int Z, G4double w)
636 {                                                 634 {
637   elements.push_back(Z);                          635   elements.push_back(Z);
638   fractions.push_back(w);                         636   fractions.push_back(w);
639   --nCurrent;                                     637   --nCurrent;
640   ++nComponents;                                  638   ++nComponents;
641   if (nCurrent == 0) {                            639   if (nCurrent == 0) {
642     G4int n = nMaterials - 1;                     640     G4int n = nMaterials - 1;
643     G4double sum = 0.0;                           641     G4double sum = 0.0;
644     G4int imin = indexes[n];                      642     G4int imin = indexes[n];
645     G4int imax = imin + components[n];            643     G4int imax = imin + components[n];
646                                                   644 
647     if (! atomCount[n]) {                      << 645     if(!atomCount[n]) {
648       for (G4int i = imin; i < imax; ++i) {    << 646       for(G4int i=imin; i<imax; ++i) {sum += fractions[i];}
649         sum += fractions[i];                   << 647       if (sum > 0.0) for (G4int i=imin; i<imax; ++i) {fractions[i] /= sum;}
650       }                                        << 
651       if (sum > 0.0) {                         << 
652         for (G4int i = imin; i < imax; ++i) {  << 
653           fractions[i] /= sum;                 << 
654         }                                      << 
655       }                                        << 
656     }                                             648     }
657   }                                               649   }
658 }                                                 650 }
659                                                   651 
660 //....oooOO0OOooo........oooOO0OOooo........oo    652 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
661                                                   653 
662 void G4NistMaterialBuilder::AddElementByWeight << 654 void G4NistMaterialBuilder::AddElementByWeightFraction(const G4String& name,
                                                   >> 655                                                        G4double w)
663 {                                                 656 {
664   G4int Z = elmBuilder->GetZ(name);               657   G4int Z = elmBuilder->GetZ(name);
665   AddElementByWeightFraction(Z, w);               658   AddElementByWeightFraction(Z, w);
666 }                                                 659 }
667                                                   660 
668 //....oooOO0OOooo........oooOO0OOooo........oo    661 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
669                                                   662 
670 void G4NistMaterialBuilder::AddElementByAtomCo    663 void G4NistMaterialBuilder::AddElementByAtomCount(G4int Z, G4int nb)
671 {                                                 664 {
672   atomCount[nMaterials - 1] = true;            << 665   atomCount[nMaterials-1] = true;
673   auto w = (G4double)nb;                       << 666   G4double w = (G4double)nb;
674   AddElementByWeightFraction(Z, w);               667   AddElementByWeightFraction(Z, w);
675 }                                                 668 }
676                                                   669 
677 //....oooOO0OOooo........oooOO0OOooo........oo    670 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
678                                                   671 
679 void G4NistMaterialBuilder::AddElementByAtomCo << 672 void G4NistMaterialBuilder::AddElementByAtomCount(const G4String& name,
                                                   >> 673               G4int nb)
680 {                                                 674 {
681   atomCount[nMaterials - 1] = true;            << 675   atomCount[nMaterials-1] = true;
682   G4int Z = elmBuilder->GetZ(name);               676   G4int Z = elmBuilder->GetZ(name);
683   auto w = (G4double)nb;                       << 677   G4double w = (G4double)nb;
684   AddElementByWeightFraction(Z, w);               678   AddElementByWeightFraction(Z, w);
685 }                                                 679 }
686                                                   680 
                                                   >> 681 
687 //....oooOO0OOooo........oooOO0OOooo........oo    682 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
688                                                   683 
689 void G4NistMaterialBuilder::Initialise()          684 void G4NistMaterialBuilder::Initialise()
690 {                                                 685 {
691   if (verbose > 0) {                              686   if (verbose > 0) {
692     G4cout << "### G4NistMaterialBuilder::Init    687     G4cout << "### G4NistMaterialBuilder::Initialise()" << G4endl;
693   }                                               688   }
694   NistSimpleMaterials();                          689   NistSimpleMaterials();
695   NistCompoundMaterials();                        690   NistCompoundMaterials();
696   NistCompoundMaterials2();                    << 
697   HepAndNuclearMaterials();                       691   HepAndNuclearMaterials();
698   SpaceMaterials();                               692   SpaceMaterials();
699   BioChemicalMaterials();                         693   BioChemicalMaterials();
700                                                   694 
701   if (verbose > 1) {                           << 695   if (verbose > 1) { ListMaterials("all"); }
702     ListMaterials("all");                      << 
703   }                                            << 
704 }                                                 696 }
705                                                   697 
706 //....oooOO0OOooo........oooOO0OOooo........oo    698 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
707                                                   699 
708 void G4NistMaterialBuilder::NistSimpleMaterial    700 void G4NistMaterialBuilder::NistSimpleMaterials()
709 {                                                 701 {
710   // density in g/cm3, mean ionisation potenti    702   // density in g/cm3, mean ionisation potential in eV
711   // clang-format off                          << 
712   AddMaterial("G4_WATER", 1.0, 0, 78., 2, kSta << 
713   AddElementByAtomCount("H" ,  2);             << 
714   AddElementByAtomCount("O" ,  1);             << 
715   chFormulas[nMaterials-1] = "H_2O";           << 
716                                                   703 
717   AddMaterial("G4_H" ,  8.37480e-5,  1,  19.2,    704   AddMaterial("G4_H" ,  8.37480e-5,  1,  19.2, 1, kStateGas);
718   AddMaterial("G4_He",  1.66322e-4,  2,  41.8,    705   AddMaterial("G4_He",  1.66322e-4,  2,  41.8, 1, kStateGas);
719   AddMaterial("G4_Li",  0.534     ,  3,  40. )    706   AddMaterial("G4_Li",  0.534     ,  3,  40. );
720   AddMaterial("G4_Be",  1.848     ,  4,  63.7)    707   AddMaterial("G4_Be",  1.848     ,  4,  63.7);
721   AddMaterial("G4_B" ,  2.37      ,  5,  76. )    708   AddMaterial("G4_B" ,  2.37      ,  5,  76. );
722   AddMaterial("G4_C" ,  2.        ,  6,  81. )    709   AddMaterial("G4_C" ,  2.        ,  6,  81. );
723   AddMaterial("G4_N" ,  1.16520e-3,  7,  82. ,    710   AddMaterial("G4_N" ,  1.16520e-3,  7,  82. , 1, kStateGas);
724   AddMaterial("G4_O" ,  1.33151e-3,  8,  95. ,    711   AddMaterial("G4_O" ,  1.33151e-3,  8,  95. , 1, kStateGas);
725   AddMaterial("G4_F" ,  1.58029e-3,  9, 115. ,    712   AddMaterial("G4_F" ,  1.58029e-3,  9, 115. , 1, kStateGas);
726   AddMaterial("G4_Ne",  8.38505e-4, 10, 137. ,    713   AddMaterial("G4_Ne",  8.38505e-4, 10, 137. , 1, kStateGas);
727   AddMaterial("G4_Na",  0.971     , 11, 149. )    714   AddMaterial("G4_Na",  0.971     , 11, 149. );
728   AddMaterial("G4_Mg",  1.74      , 12, 156. )    715   AddMaterial("G4_Mg",  1.74      , 12, 156. );
729   AddMaterial("G4_Al",  2.699     , 13, 166. )    716   AddMaterial("G4_Al",  2.699     , 13, 166. );
730   AddMaterial("G4_Si",  2.33      , 14, 173. )    717   AddMaterial("G4_Si",  2.33      , 14, 173. );
731   AddMaterial("G4_P" ,  2.2       , 15, 173. )    718   AddMaterial("G4_P" ,  2.2       , 15, 173. );
732   AddMaterial("G4_S" ,  2.0       , 16, 180. )    719   AddMaterial("G4_S" ,  2.0       , 16, 180. );
733   AddMaterial("G4_Cl",  2.99473e-3, 17, 174. ,    720   AddMaterial("G4_Cl",  2.99473e-3, 17, 174. , 1, kStateGas);
734   AddMaterial("G4_Ar",  1.66201e-3, 18, 188.0,    721   AddMaterial("G4_Ar",  1.66201e-3, 18, 188.0, 1, kStateGas);
735   AddMaterial("G4_K" ,  0.862     , 19, 190. )    722   AddMaterial("G4_K" ,  0.862     , 19, 190. );
736   AddMaterial("G4_Ca",  1.55      , 20, 191. )    723   AddMaterial("G4_Ca",  1.55      , 20, 191. );
737   AddMaterial("G4_Sc",  2.989     , 21, 216. )    724   AddMaterial("G4_Sc",  2.989     , 21, 216. );
738   AddMaterial("G4_Ti",  4.54      , 22, 233. )    725   AddMaterial("G4_Ti",  4.54      , 22, 233. );
739   AddMaterial("G4_V" ,  6.11      , 23, 245. )    726   AddMaterial("G4_V" ,  6.11      , 23, 245. );
740   AddMaterial("G4_Cr",  7.18      , 24, 257. )    727   AddMaterial("G4_Cr",  7.18      , 24, 257. );
741   AddMaterial("G4_Mn",  7.44      , 25, 272. )    728   AddMaterial("G4_Mn",  7.44      , 25, 272. );
742   AddMaterial("G4_Fe",  7.874     , 26, 286. )    729   AddMaterial("G4_Fe",  7.874     , 26, 286. );
743   AddMaterial("G4_Co",  8.9       , 27, 297. )    730   AddMaterial("G4_Co",  8.9       , 27, 297. );
744   AddMaterial("G4_Ni",  8.902     , 28, 311. )    731   AddMaterial("G4_Ni",  8.902     , 28, 311. );
745   AddMaterial("G4_Cu",  8.96      , 29, 322. )    732   AddMaterial("G4_Cu",  8.96      , 29, 322. );
746   AddMaterial("G4_Zn",  7.133     , 30, 330. )    733   AddMaterial("G4_Zn",  7.133     , 30, 330. );
747   AddMaterial("G4_Ga",  5.904     , 31, 334. )    734   AddMaterial("G4_Ga",  5.904     , 31, 334. );
748   AddMaterial("G4_Ge",  5.323     , 32, 350. )    735   AddMaterial("G4_Ge",  5.323     , 32, 350. );
749   AddMaterial("G4_As",  5.73      , 33, 347. )    736   AddMaterial("G4_As",  5.73      , 33, 347. );
750   AddMaterial("G4_Se",  4.5       , 34, 348. )    737   AddMaterial("G4_Se",  4.5       , 34, 348. );
751   AddMaterial("G4_Br",  7.07210e-3, 35, 343. ,    738   AddMaterial("G4_Br",  7.07210e-3, 35, 343. , 1, kStateGas);
752   AddMaterial("G4_Kr",  3.47832e-3, 36, 352. ,    739   AddMaterial("G4_Kr",  3.47832e-3, 36, 352. , 1, kStateGas);
753   AddMaterial("G4_Rb",  1.532     , 37, 363. )    740   AddMaterial("G4_Rb",  1.532     , 37, 363. );
754   AddMaterial("G4_Sr",  2.54      , 38, 366. )    741   AddMaterial("G4_Sr",  2.54      , 38, 366. );
755   AddMaterial("G4_Y" ,  4.469     , 39, 379. )    742   AddMaterial("G4_Y" ,  4.469     , 39, 379. );
756   AddMaterial("G4_Zr",  6.506     , 40, 393. )    743   AddMaterial("G4_Zr",  6.506     , 40, 393. );
757   AddMaterial("G4_Nb",  8.57      , 41, 417. )    744   AddMaterial("G4_Nb",  8.57      , 41, 417. );
758   AddMaterial("G4_Mo", 10.22      , 42, 424. )    745   AddMaterial("G4_Mo", 10.22      , 42, 424. );
759   AddMaterial("G4_Tc", 11.50      , 43, 428. )    746   AddMaterial("G4_Tc", 11.50      , 43, 428. );
760   AddMaterial("G4_Ru", 12.41      , 44, 441. )    747   AddMaterial("G4_Ru", 12.41      , 44, 441. );
761   AddMaterial("G4_Rh", 12.41      , 45, 449. )    748   AddMaterial("G4_Rh", 12.41      , 45, 449. );
762   AddMaterial("G4_Pd", 12.02      , 46, 470. )    749   AddMaterial("G4_Pd", 12.02      , 46, 470. );
763   AddMaterial("G4_Ag", 10.5       , 47, 470. )    750   AddMaterial("G4_Ag", 10.5       , 47, 470. );
764   AddMaterial("G4_Cd",  8.65      , 48, 469. )    751   AddMaterial("G4_Cd",  8.65      , 48, 469. );
765   AddMaterial("G4_In",  7.31      , 49, 488. )    752   AddMaterial("G4_In",  7.31      , 49, 488. );
766   AddMaterial("G4_Sn",  7.31      , 50, 488. )    753   AddMaterial("G4_Sn",  7.31      , 50, 488. );
767   AddMaterial("G4_Sb",  6.691     , 51, 487. )    754   AddMaterial("G4_Sb",  6.691     , 51, 487. );
768   AddMaterial("G4_Te",  6.24      , 52, 485. )    755   AddMaterial("G4_Te",  6.24      , 52, 485. );
769   AddMaterial("G4_I" ,  4.93      , 53, 491. )    756   AddMaterial("G4_I" ,  4.93      , 53, 491. );
770   AddMaterial("G4_Xe",  5.48536e-3, 54, 482. ,    757   AddMaterial("G4_Xe",  5.48536e-3, 54, 482. , 1, kStateGas);
771   AddMaterial("G4_Cs",  1.873     , 55, 488. )    758   AddMaterial("G4_Cs",  1.873     , 55, 488. );
772   AddMaterial("G4_Ba",  3.5       , 56, 491. )    759   AddMaterial("G4_Ba",  3.5       , 56, 491. );
773   AddMaterial("G4_La",  6.154     , 57, 501. )    760   AddMaterial("G4_La",  6.154     , 57, 501. );
774   AddMaterial("G4_Ce",  6.657     , 58, 523. )    761   AddMaterial("G4_Ce",  6.657     , 58, 523. );
775   AddMaterial("G4_Pr",  6.71      , 59, 535. )    762   AddMaterial("G4_Pr",  6.71      , 59, 535. );
776   AddMaterial("G4_Nd",  6.9       , 60, 546. )    763   AddMaterial("G4_Nd",  6.9       , 60, 546. );
777   AddMaterial("G4_Pm",  7.22      , 61, 560. )    764   AddMaterial("G4_Pm",  7.22      , 61, 560. );
778   AddMaterial("G4_Sm",  7.46      , 62, 574. )    765   AddMaterial("G4_Sm",  7.46      , 62, 574. );
779   AddMaterial("G4_Eu",  5.243     , 63, 580. )    766   AddMaterial("G4_Eu",  5.243     , 63, 580. );
780   AddMaterial("G4_Gd",  7.9004    , 64, 591. )    767   AddMaterial("G4_Gd",  7.9004    , 64, 591. );
781   AddMaterial("G4_Tb",  8.229     , 65, 614. )    768   AddMaterial("G4_Tb",  8.229     , 65, 614. );
782   AddMaterial("G4_Dy",  8.55      , 66, 628. )    769   AddMaterial("G4_Dy",  8.55      , 66, 628. );
783   AddMaterial("G4_Ho",  8.795     , 67, 650. )    770   AddMaterial("G4_Ho",  8.795     , 67, 650. );
784   AddMaterial("G4_Er",  9.066     , 68, 658. )    771   AddMaterial("G4_Er",  9.066     , 68, 658. );
785   AddMaterial("G4_Tm",  9.321     , 69, 674. )    772   AddMaterial("G4_Tm",  9.321     , 69, 674. );
786   AddMaterial("G4_Yb",  6.73      , 70, 684. )    773   AddMaterial("G4_Yb",  6.73      , 70, 684. );
787   AddMaterial("G4_Lu",  9.84      , 71, 694. )    774   AddMaterial("G4_Lu",  9.84      , 71, 694. );
788   AddMaterial("G4_Hf", 13.31      , 72, 705. )    775   AddMaterial("G4_Hf", 13.31      , 72, 705. );
789   AddMaterial("G4_Ta", 16.654     , 73, 718. )    776   AddMaterial("G4_Ta", 16.654     , 73, 718. );
790   AddMaterial("G4_W" , 19.30      , 74, 727. )    777   AddMaterial("G4_W" , 19.30      , 74, 727. );
791   AddMaterial("G4_Re", 21.02      , 75, 736. )    778   AddMaterial("G4_Re", 21.02      , 75, 736. );
792   AddMaterial("G4_Os", 22.57      , 76, 746. )    779   AddMaterial("G4_Os", 22.57      , 76, 746. );
793   AddMaterial("G4_Ir", 22.42      , 77, 757. )    780   AddMaterial("G4_Ir", 22.42      , 77, 757. );
794   AddMaterial("G4_Pt", 21.45      , 78, 790. )    781   AddMaterial("G4_Pt", 21.45      , 78, 790. );
795   AddMaterial("G4_Au", 19.32      , 79, 790. )    782   AddMaterial("G4_Au", 19.32      , 79, 790. );
796   AddMaterial("G4_Hg", 13.546     , 80, 800. )    783   AddMaterial("G4_Hg", 13.546     , 80, 800. );
797   AddMaterial("G4_Tl", 11.72      , 81, 810. )    784   AddMaterial("G4_Tl", 11.72      , 81, 810. );
798   AddMaterial("G4_Pb", 11.35      , 82, 823. )    785   AddMaterial("G4_Pb", 11.35      , 82, 823. );
799   AddMaterial("G4_Bi",  9.747     , 83, 823. )    786   AddMaterial("G4_Bi",  9.747     , 83, 823. );
800   AddMaterial("G4_Po",  9.32      , 84, 830. )    787   AddMaterial("G4_Po",  9.32      , 84, 830. );
801   AddMaterial("G4_At",  9.32      , 85, 825. )    788   AddMaterial("G4_At",  9.32      , 85, 825. );
802   AddMaterial("G4_Rn",  9.00662e-3, 86, 794. ,    789   AddMaterial("G4_Rn",  9.00662e-3, 86, 794. , 1, kStateGas);
803   AddMaterial("G4_Fr",  1.00      , 87, 827. )    790   AddMaterial("G4_Fr",  1.00      , 87, 827. );
804   AddMaterial("G4_Ra",  5.00      , 88, 826. )    791   AddMaterial("G4_Ra",  5.00      , 88, 826. );
805   AddMaterial("G4_Ac", 10.07      , 89, 841. )    792   AddMaterial("G4_Ac", 10.07      , 89, 841. );
806   AddMaterial("G4_Th", 11.72      , 90, 847. )    793   AddMaterial("G4_Th", 11.72      , 90, 847. );
807   AddMaterial("G4_Pa", 15.37      , 91, 878. )    794   AddMaterial("G4_Pa", 15.37      , 91, 878. );
808   AddMaterial("G4_U" , 18.95      , 92, 890. )    795   AddMaterial("G4_U" , 18.95      , 92, 890. );
809   AddMaterial("G4_Np", 20.25      , 93, 902. )    796   AddMaterial("G4_Np", 20.25      , 93, 902. );
810   AddMaterial("G4_Pu", 19.84      , 94, 921. )    797   AddMaterial("G4_Pu", 19.84      , 94, 921. );
811   AddMaterial("G4_Am", 13.67      , 95, 934. )    798   AddMaterial("G4_Am", 13.67      , 95, 934. );
812   AddMaterial("G4_Cm", 13.51      , 96, 939. )    799   AddMaterial("G4_Cm", 13.51      , 96, 939. );
813   AddMaterial("G4_Bk", 14.00      , 97, 952. )    800   AddMaterial("G4_Bk", 14.00      , 97, 952. );
814   AddMaterial("G4_Cf", 10.00      , 98, 966. )    801   AddMaterial("G4_Cf", 10.00      , 98, 966. );
815                                                   802 
816   nElementary = nMaterials;                       803   nElementary = nMaterials;
817 }                                                 804 }
818                                                   805 
819 //....oooOO0OOooo........oooOO0OOooo........oo    806 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
820                                                   807 
821 void G4NistMaterialBuilder::NistCompoundMateri    808 void G4NistMaterialBuilder::NistCompoundMaterials()
822 {                                                 809 {
823   AddMaterial("G4_A-150_TISSUE", 1.127, 0, 65.    810   AddMaterial("G4_A-150_TISSUE", 1.127, 0, 65.1, 6);
824   AddElementByWeightFraction( 1, 0.101327);       811   AddElementByWeightFraction( 1, 0.101327);
825   AddElementByWeightFraction( 6, 0.775501);       812   AddElementByWeightFraction( 6, 0.775501);
826   AddElementByWeightFraction( 7, 0.035057);       813   AddElementByWeightFraction( 7, 0.035057);
827   AddElementByWeightFraction( 8, 0.052316);       814   AddElementByWeightFraction( 8, 0.052316);
828   AddElementByWeightFraction( 9, 0.017422);       815   AddElementByWeightFraction( 9, 0.017422);
829   AddElementByWeightFraction(20, 0.018378);       816   AddElementByWeightFraction(20, 0.018378);
830                                                   817 
831   AddMaterial("G4_ACETONE", 0.7899, 0, 64.2, 3 << 818   AddMaterial("G4_ACETONE", 0.7899, 0, 64.2, 3);
832   AddElementByAtomCount("C" ,  3);             << 819   AddElementByWeightFraction( 1, 0.104122);
833   AddElementByAtomCount("H" ,  6);             << 820   AddElementByWeightFraction( 6, 0.620405);
834   AddElementByAtomCount("O" ,  1);             << 821   AddElementByWeightFraction( 8, 0.275473);
835                                                   822 
836   AddMaterial("G4_ACETYLENE", 0.0010967, 0, 58    823   AddMaterial("G4_ACETYLENE", 0.0010967, 0, 58.2, 2, kStateGas);
837   AddElementByAtomCount("C" ,  2);             << 824   AddElementByWeightFraction( 1, 0.077418);
838   AddElementByAtomCount("H" ,  2);             << 825   AddElementByWeightFraction( 6, 0.922582);
839                                                   826 
840   //Tan, Z., et al. NIMB,2006(248)             << 827   AddMaterial("G4_ADENINE", 1.6/*1.35*/, 0, 71.4, 3);
841   AddMaterial("G4_ADENINE", 1.35, 0, 71.4, 3); << 828   AddElementByAtomCount("H",5 );
842   AddElementByAtomCount("C" ,  5);             << 829   AddElementByAtomCount("C",5 );
843   AddElementByAtomCount("H" ,  5);             << 830   AddElementByAtomCount("N",5 );
844   AddElementByAtomCount("N" ,  5);             << 
845                                                   831 
846   AddMaterial("G4_ADIPOSE_TISSUE_ICRP", 0.95,     832   AddMaterial("G4_ADIPOSE_TISSUE_ICRP", 0.95, 0, 63.2, 7);
847   AddElementByWeightFraction( 1, 0.114);          833   AddElementByWeightFraction( 1, 0.114);
848   AddElementByWeightFraction( 6, 0.598);          834   AddElementByWeightFraction( 6, 0.598);
849   AddElementByWeightFraction( 7, 0.007);          835   AddElementByWeightFraction( 7, 0.007);
850   AddElementByWeightFraction( 8, 0.278);          836   AddElementByWeightFraction( 8, 0.278);
851   AddElementByWeightFraction(11, 0.001);          837   AddElementByWeightFraction(11, 0.001);
852   AddElementByWeightFraction(16, 0.001);          838   AddElementByWeightFraction(16, 0.001);
853   AddElementByWeightFraction(17, 0.001);          839   AddElementByWeightFraction(17, 0.001);
854                                                   840 
855   AddMaterial("G4_AIR", 0.00120479, 0, 85.7, 4    841   AddMaterial("G4_AIR", 0.00120479, 0, 85.7, 4, kStateGas);
856   AddElementByWeightFraction( 6, 0.000124);       842   AddElementByWeightFraction( 6, 0.000124);
857   AddElementByWeightFraction( 7, 0.755267);       843   AddElementByWeightFraction( 7, 0.755267);
858   AddElementByWeightFraction( 8, 0.231781);       844   AddElementByWeightFraction( 8, 0.231781);
859   AddElementByWeightFraction(18, 0.012827);       845   AddElementByWeightFraction(18, 0.012827);
860                                                   846 
861   AddMaterial("G4_ALANINE", 1.42, 0, 71.9, 4);    847   AddMaterial("G4_ALANINE", 1.42, 0, 71.9, 4);
862   AddElementByAtomCount("C" ,  3);             << 848   AddElementByWeightFraction( 1, 0.07919 );
863   AddElementByAtomCount("H" ,  7);             << 849   AddElementByWeightFraction( 6, 0.404439);
864   AddElementByAtomCount("N" ,  1);             << 850   AddElementByWeightFraction( 7, 0.157213);
865   AddElementByAtomCount("O" ,  2);             << 851   AddElementByWeightFraction( 8, 0.359159);
866                                                   852 
867   AddMaterial("G4_ALUMINUM_OXIDE", 3.97, 0, 14    853   AddMaterial("G4_ALUMINUM_OXIDE", 3.97, 0, 145.2, 2);
868   AddElementByAtomCount("Al",  2);             << 854   AddElementByWeightFraction( 8, 0.470749);
869   AddElementByAtomCount("O" ,  3);             << 855   AddElementByWeightFraction(13, 0.529251);
870   chFormulas[nMaterials-1] = "Al_2O_3";           856   chFormulas[nMaterials-1] = "Al_2O_3";
871                                                   857 
872   AddMaterial("G4_AMBER", 1.1, 0, 63.2, 3);       858   AddMaterial("G4_AMBER", 1.1, 0, 63.2, 3);
873   AddElementByWeightFraction( 1, 0.10593 );       859   AddElementByWeightFraction( 1, 0.10593 );
874   AddElementByWeightFraction( 6, 0.788973);       860   AddElementByWeightFraction( 6, 0.788973);
875   AddElementByWeightFraction( 8, 0.105096);       861   AddElementByWeightFraction( 8, 0.105096);
876                                                   862 
877   AddMaterial("G4_AMMONIA", 0.000826019, 0, 53    863   AddMaterial("G4_AMMONIA", 0.000826019, 0, 53.7, 2, kStateGas);
878   AddElementByAtomCount("N" ,  1);             << 864   AddElementByWeightFraction( 1, 0.177547);
879   AddElementByAtomCount("H" ,  3);             << 865   AddElementByWeightFraction( 7, 0.822453);
880                                                   866 
881   AddMaterial("G4_ANILINE", 1.0235, 0, 66.2, 3    867   AddMaterial("G4_ANILINE", 1.0235, 0, 66.2, 3);
882   AddElementByAtomCount("C" ,  6);             << 868   AddElementByWeightFraction( 1, 0.075759);
883   AddElementByAtomCount("H" ,  7);             << 869   AddElementByWeightFraction( 6, 0.773838);
884   AddElementByAtomCount("N" ,  1);             << 870   AddElementByWeightFraction( 7, 0.150403);
885                                                   871 
886   AddMaterial("G4_ANTHRACENE", 1.283, 0, 69.5,    872   AddMaterial("G4_ANTHRACENE", 1.283, 0, 69.5, 2);
887   AddElementByAtomCount("C" , 14);             << 873   AddElementByWeightFraction( 1, 0.05655);
888   AddElementByAtomCount("H" , 10);             << 874   AddElementByWeightFraction( 6, 0.94345);
889                                                   875 
890   AddMaterial("G4_B-100_BONE", 1.45, 0, 85.9,     876   AddMaterial("G4_B-100_BONE", 1.45, 0, 85.9, 6);
891   AddElementByWeightFraction( 1, 0.065471);       877   AddElementByWeightFraction( 1, 0.065471);
892   AddElementByWeightFraction( 6, 0.536945);       878   AddElementByWeightFraction( 6, 0.536945);
893   AddElementByWeightFraction( 7, 0.0215  );       879   AddElementByWeightFraction( 7, 0.0215  );
894   AddElementByWeightFraction( 8, 0.032085);       880   AddElementByWeightFraction( 8, 0.032085);
895   AddElementByWeightFraction( 9, 0.167411);       881   AddElementByWeightFraction( 9, 0.167411);
896   AddElementByWeightFraction(20, 0.176589);       882   AddElementByWeightFraction(20, 0.176589);
897                                                   883 
898   AddMaterial("G4_BAKELITE", 1.25, 0, 72.4, 3)    884   AddMaterial("G4_BAKELITE", 1.25, 0, 72.4, 3);
899   AddElementByWeightFraction( 1, 0.057441);       885   AddElementByWeightFraction( 1, 0.057441);
900   AddElementByWeightFraction( 6, 0.774591);       886   AddElementByWeightFraction( 6, 0.774591);
901   AddElementByWeightFraction( 8, 0.167968);       887   AddElementByWeightFraction( 8, 0.167968);
902                                                   888 
903   AddMaterial("G4_BARIUM_FLUORIDE", 4.89 ,0, 3    889   AddMaterial("G4_BARIUM_FLUORIDE", 4.89 ,0, 375.9, 2);
904   AddElementByAtomCount("Ba",  1);             << 890   AddElementByWeightFraction( 9, 0.21672);
905   AddElementByAtomCount("F" ,  2);             << 891   AddElementByWeightFraction(56, 0.78328);
906                                                   892 
907   AddMaterial("G4_BARIUM_SULFATE", 4.5, 0, 285    893   AddMaterial("G4_BARIUM_SULFATE", 4.5, 0, 285.7, 3);
908   AddElementByAtomCount("Ba",  1);             << 894   AddElementByWeightFraction( 8,0.274212);
909   AddElementByAtomCount("S" ,  1);             << 895   AddElementByWeightFraction(16,0.137368);
910   AddElementByAtomCount("O" ,  4);             << 896   AddElementByWeightFraction(56,0.58842 );
911                                                << 897 
912   AddMaterial("G4_BENZENE", 0.87865, 0, 63.4,  << 898   AddMaterial("G4_BENZENE", 0.87865, 0, 63.4, 2);
913   AddElementByAtomCount("C" ,  6);             << 899   AddElementByWeightFraction( 1, 0.077418);
914   AddElementByAtomCount("H" ,  6);             << 900   AddElementByWeightFraction( 6, 0.922582);
915                                                   901 
916   AddMaterial("G4_BERYLLIUM_OXIDE", 3.01, 0, 9    902   AddMaterial("G4_BERYLLIUM_OXIDE", 3.01, 0, 93.2, 2);
917   AddElementByAtomCount("Be",  1);             << 903   AddElementByWeightFraction( 4, 0.36032);
918   AddElementByAtomCount("O" ,  1);             << 904   AddElementByWeightFraction( 8, 0.63968);
919                                                   905 
920   AddMaterial("G4_BGO", 7.13, 0, 534.1, 3);       906   AddMaterial("G4_BGO", 7.13, 0, 534.1, 3);
921   AddElementByAtomCount("Bi",  4);             << 907   AddElementByWeightFraction( 8, 0.154126);
922   AddElementByAtomCount("Ge",  3);             << 908   AddElementByWeightFraction(32, 0.17482 );
923   AddElementByAtomCount("O" , 12);             << 909   AddElementByWeightFraction(83, 0.671054);
924                                                   910 
925   AddMaterial("G4_BLOOD_ICRP", 1.06, 0, 75.2,  << 911   AddMaterial("G4_BLOOD_ICRP", 1.06, 0, 75.2, 10);
926   AddElementByWeightFraction( 1, 0.102);          912   AddElementByWeightFraction( 1, 0.102);
927   AddElementByWeightFraction( 6, 0.110);          913   AddElementByWeightFraction( 6, 0.110);
928   AddElementByWeightFraction( 7, 0.033);          914   AddElementByWeightFraction( 7, 0.033);
929   AddElementByWeightFraction( 8, 0.745);          915   AddElementByWeightFraction( 8, 0.745);
930   AddElementByWeightFraction(11, 0.001);          916   AddElementByWeightFraction(11, 0.001);
931   AddElementByWeightFraction(15, 0.001);          917   AddElementByWeightFraction(15, 0.001);
932   AddElementByWeightFraction(16, 0.002);          918   AddElementByWeightFraction(16, 0.002);
933   AddElementByWeightFraction(17, 0.003);          919   AddElementByWeightFraction(17, 0.003);
934   AddElementByWeightFraction(19, 0.002);          920   AddElementByWeightFraction(19, 0.002);
935   AddElementByWeightFraction(26, 0.001);          921   AddElementByWeightFraction(26, 0.001);
936                                                   922 
937   AddMaterial("G4_BONE_COMPACT_ICRU", 1.85, 0,    923   AddMaterial("G4_BONE_COMPACT_ICRU", 1.85, 0, 91.9, 8);
938   AddElementByWeightFraction( 1, 0.064);          924   AddElementByWeightFraction( 1, 0.064);
939   AddElementByWeightFraction( 6, 0.278);          925   AddElementByWeightFraction( 6, 0.278);
940   AddElementByWeightFraction( 7, 0.027);          926   AddElementByWeightFraction( 7, 0.027);
941   AddElementByWeightFraction( 8, 0.410);          927   AddElementByWeightFraction( 8, 0.410);
942   AddElementByWeightFraction(12, 0.002);          928   AddElementByWeightFraction(12, 0.002);
943   AddElementByWeightFraction(15, 0.07 );          929   AddElementByWeightFraction(15, 0.07 );
944   AddElementByWeightFraction(16, 0.002);          930   AddElementByWeightFraction(16, 0.002);
945   AddElementByWeightFraction(20, 0.147);          931   AddElementByWeightFraction(20, 0.147);
946                                                   932 
947   // Sceleton Cortical bone for Adult ICRU 46     933   // Sceleton Cortical bone for Adult ICRU 46
948   AddMaterial("G4_BONE_CORTICAL_ICRP", 1.92, 0    934   AddMaterial("G4_BONE_CORTICAL_ICRP", 1.92, 0, 110, 9);
949   AddElementByWeightFraction( 1, 0.034);          935   AddElementByWeightFraction( 1, 0.034);
950   AddElementByWeightFraction( 6, 0.155);          936   AddElementByWeightFraction( 6, 0.155);
951   AddElementByWeightFraction( 7, 0.042);          937   AddElementByWeightFraction( 7, 0.042);
952   AddElementByWeightFraction( 8, 0.435);          938   AddElementByWeightFraction( 8, 0.435);
953   AddElementByWeightFraction(11, 0.001);          939   AddElementByWeightFraction(11, 0.001);
954   AddElementByWeightFraction(12, 0.002);          940   AddElementByWeightFraction(12, 0.002);
955   AddElementByWeightFraction(15, 0.103);          941   AddElementByWeightFraction(15, 0.103);
956   AddElementByWeightFraction(16, 0.003);          942   AddElementByWeightFraction(16, 0.003);
957   AddElementByWeightFraction(20, 0.225);          943   AddElementByWeightFraction(20, 0.225);
958                                                   944 
959   AddMaterial("G4_BORON_CARBIDE", 2.52, 0, 84.    945   AddMaterial("G4_BORON_CARBIDE", 2.52, 0, 84.7, 2);
960   AddElementByAtomCount("B" ,  4);             << 946   AddElementByWeightFraction( 5, 0.78261);
961   AddElementByAtomCount("C" ,  1);             << 947   AddElementByWeightFraction( 6, 0.21739);
962                                                   948 
963   AddMaterial("G4_BORON_OXIDE", 1.812, 0, 99.6    949   AddMaterial("G4_BORON_OXIDE", 1.812, 0, 99.6, 2);
964   AddElementByAtomCount("B" ,  2);             << 950   AddElementByWeightFraction( 5, 0.310551);
965   AddElementByAtomCount("O" ,  3);             << 951   AddElementByWeightFraction( 8, 0.689449);
966                                                   952 
967   AddMaterial("G4_BRAIN_ICRP", 1.04, 0, 73.3,     953   AddMaterial("G4_BRAIN_ICRP", 1.04, 0, 73.3, 9);
968   AddElementByWeightFraction( 1, 0.107);          954   AddElementByWeightFraction( 1, 0.107);
969   AddElementByWeightFraction( 6, 0.145);          955   AddElementByWeightFraction( 6, 0.145);
970   AddElementByWeightFraction( 7, 0.022);          956   AddElementByWeightFraction( 7, 0.022);
971   AddElementByWeightFraction( 8, 0.712);          957   AddElementByWeightFraction( 8, 0.712);
972   AddElementByWeightFraction(11, 0.002);          958   AddElementByWeightFraction(11, 0.002);
973   AddElementByWeightFraction(15, 0.004);          959   AddElementByWeightFraction(15, 0.004);
974   AddElementByWeightFraction(16, 0.002);          960   AddElementByWeightFraction(16, 0.002);
975   AddElementByWeightFraction(17, 0.003);          961   AddElementByWeightFraction(17, 0.003);
976   AddElementByWeightFraction(19, 0.003);          962   AddElementByWeightFraction(19, 0.003);
977                                                   963 
978   AddMaterial("G4_BUTANE", 0.00249343, 0, 48.3    964   AddMaterial("G4_BUTANE", 0.00249343, 0, 48.3, 2, kStateGas);
979   AddElementByAtomCount("C" ,  4);             << 965   AddElementByWeightFraction( 1, 0.173408);
980   AddElementByAtomCount("H" , 10);             << 966   AddElementByWeightFraction( 6, 0.826592);
981                                                   967 
982   AddMaterial("G4_N-BUTYL_ALCOHOL", 0.8098, 0, << 968   AddMaterial("G4_N-BUTYL_ALCOHOL", 0.8098, 0, 59.9, 3);
983   AddElementByAtomCount("C" ,  4);             << 969   AddElementByWeightFraction( 1, 0.135978);
984   AddElementByAtomCount("H" , 10);             << 970   AddElementByWeightFraction( 6, 0.648171);
985   AddElementByAtomCount("O" ,  1);             << 971   AddElementByWeightFraction( 8, 0.215851);
986                                                   972 
987   AddMaterial("G4_C-552", 1.76, 0, 86.8, 5);      973   AddMaterial("G4_C-552", 1.76, 0, 86.8, 5);
988   AddElementByWeightFraction( 1, 0.02468 );       974   AddElementByWeightFraction( 1, 0.02468 );
989   AddElementByWeightFraction( 6, 0.50161 );       975   AddElementByWeightFraction( 6, 0.50161 );
990   AddElementByWeightFraction( 8, 0.004527);       976   AddElementByWeightFraction( 8, 0.004527);
991   AddElementByWeightFraction( 9, 0.465209);       977   AddElementByWeightFraction( 9, 0.465209);
992   AddElementByWeightFraction(14, 0.003973);       978   AddElementByWeightFraction(14, 0.003973);
993                                                   979 
994   AddMaterial("G4_CADMIUM_TELLURIDE", 6.2, 0,     980   AddMaterial("G4_CADMIUM_TELLURIDE", 6.2, 0, 539.3, 2);
995   AddElementByAtomCount("Cd",  1);             << 981   AddElementByWeightFraction(48, 0.468355);
996   AddElementByAtomCount("Te",  1);             << 982   AddElementByWeightFraction(52, 0.531645);
997                                                   983 
998   AddMaterial("G4_CADMIUM_TUNGSTATE", 7.9, 0,     984   AddMaterial("G4_CADMIUM_TUNGSTATE", 7.9, 0, 468.3, 3);
999   AddElementByAtomCount("Cd",  1);             << 985   AddElementByWeightFraction( 8, 0.177644);
1000   AddElementByAtomCount("W" ,  1);            << 986   AddElementByWeightFraction(48, 0.312027);
1001   AddElementByAtomCount("O" ,  4);            << 987   AddElementByWeightFraction(74, 0.510329);
1002                                                  988 
1003   AddMaterial("G4_CALCIUM_CARBONATE", 2.8, 0,    989   AddMaterial("G4_CALCIUM_CARBONATE", 2.8, 0, 136.4, 3);
1004   AddElementByAtomCount("Ca",  1);            << 990   AddElementByWeightFraction( 6, 0.120003);
1005   AddElementByAtomCount("C" ,  1);            << 991   AddElementByWeightFraction( 8, 0.479554);
1006   AddElementByAtomCount("O" ,  3);            << 992   AddElementByWeightFraction(20, 0.400443);
1007                                                  993 
1008   AddMaterial("G4_CALCIUM_FLUORIDE", 3.18, 0,    994   AddMaterial("G4_CALCIUM_FLUORIDE", 3.18, 0, 166., 2);
1009   AddElementByAtomCount("Ca",  1);            << 995   AddElementByWeightFraction( 9, 0.486659);
1010   AddElementByAtomCount("F" ,  2);            << 996   AddElementByWeightFraction(20, 0.513341);
1011                                                  997 
1012   AddMaterial("G4_CALCIUM_OXIDE", 3.3, 0, 176    998   AddMaterial("G4_CALCIUM_OXIDE", 3.3, 0, 176.1, 2);
1013   AddElementByAtomCount("Ca",  1);            << 999   AddElementByWeightFraction( 8, 0.285299);
1014   AddElementByAtomCount("O" ,  1);            << 1000   AddElementByWeightFraction(20, 0.714701);
1015                                                  1001 
1016   AddMaterial("G4_CALCIUM_SULFATE", 2.96, 0,     1002   AddMaterial("G4_CALCIUM_SULFATE", 2.96, 0, 152.3, 3);
1017   AddElementByAtomCount("Ca",  1);            << 1003   AddElementByWeightFraction( 8, 0.470095);
1018   AddElementByAtomCount("S" ,  1);            << 1004   AddElementByWeightFraction(16, 0.235497);
1019   AddElementByAtomCount("O" ,  4);            << 1005   AddElementByWeightFraction(20, 0.294408);
1020                                                  1006 
1021   AddMaterial("G4_CALCIUM_TUNGSTATE", 6.062,     1007   AddMaterial("G4_CALCIUM_TUNGSTATE", 6.062, 0, 395., 3);
1022   AddElementByAtomCount("Ca",  1);            << 1008   AddElementByWeightFraction( 8, 0.22227 );
1023   AddElementByAtomCount("W" ,  1);            << 1009   AddElementByWeightFraction(20, 0.139202);
1024   AddElementByAtomCount("O" ,  4);            << 1010   AddElementByWeightFraction(74, 0.638529);
1025                                                  1011 
1026   AddMaterial("G4_CARBON_DIOXIDE", 0.00184212    1012   AddMaterial("G4_CARBON_DIOXIDE", 0.00184212, 0, 85., 2, kStateGas);
1027   AddElementByAtomCount("C" ,  1);            << 1013   AddElementByWeightFraction( 6, 0.272916);
1028   AddElementByAtomCount("O" ,  2);            << 1014   AddElementByWeightFraction( 8, 0.727084);
1029   chFormulas[nMaterials-1] = "CO_2";             1015   chFormulas[nMaterials-1] = "CO_2";
1030                                                  1016 
1031   AddMaterial("G4_CARBON_TETRACHLORIDE", 1.59    1017   AddMaterial("G4_CARBON_TETRACHLORIDE", 1.594, 0, 166.3, 2);
1032   AddElementByAtomCount("C" ,  1);            << 1018   AddElementByWeightFraction( 6, 0.078083);
1033   AddElementByAtomCount("Cl",  4);            << 1019   AddElementByWeightFraction(17, 0.921917);
1034                                                  1020 
1035   AddMaterial("G4_CELLULOSE_CELLOPHANE", 1.42    1021   AddMaterial("G4_CELLULOSE_CELLOPHANE", 1.42, 0, 77.6, 3);
1036   AddElementByAtomCount("C" ,  6);            << 1022   AddElementByWeightFraction( 1, 0.062162);
1037   AddElementByAtomCount("H" , 10);            << 1023   AddElementByWeightFraction( 6, 0.444462);
1038   AddElementByAtomCount("O" ,  5);            << 1024   AddElementByWeightFraction( 8, 0.493376);
1039                                                  1025 
1040   AddMaterial("G4_CELLULOSE_BUTYRATE", 1.2, 0    1026   AddMaterial("G4_CELLULOSE_BUTYRATE", 1.2, 0, 74.6, 3);
1041   AddElementByWeightFraction( 1, 0.067125);      1027   AddElementByWeightFraction( 1, 0.067125);
1042   AddElementByWeightFraction( 6, 0.545403);      1028   AddElementByWeightFraction( 6, 0.545403);
1043   AddElementByWeightFraction( 8, 0.387472);      1029   AddElementByWeightFraction( 8, 0.387472);
1044                                                  1030 
1045   AddMaterial("G4_CELLULOSE_NITRATE", 1.49, 0    1031   AddMaterial("G4_CELLULOSE_NITRATE", 1.49, 0, 87., 4);
1046   AddElementByWeightFraction( 1, 0.029216);      1032   AddElementByWeightFraction( 1, 0.029216);
1047   AddElementByWeightFraction( 6, 0.271296);      1033   AddElementByWeightFraction( 6, 0.271296);
1048   AddElementByWeightFraction( 7, 0.121276);      1034   AddElementByWeightFraction( 7, 0.121276);
1049   AddElementByWeightFraction( 8, 0.578212);      1035   AddElementByWeightFraction( 8, 0.578212);
1050                                                  1036 
1051   AddMaterial("G4_CERIC_SULFATE", 1.03, 0, 76    1037   AddMaterial("G4_CERIC_SULFATE", 1.03, 0, 76.7, 5);
1052   AddElementByWeightFraction( 1, 0.107596);      1038   AddElementByWeightFraction( 1, 0.107596);
1053   AddElementByWeightFraction( 7, 0.0008  );      1039   AddElementByWeightFraction( 7, 0.0008  );
1054   AddElementByWeightFraction( 8, 0.874976);      1040   AddElementByWeightFraction( 8, 0.874976);
1055   AddElementByWeightFraction(16, 0.014627);      1041   AddElementByWeightFraction(16, 0.014627);
1056   AddElementByWeightFraction(58, 0.002001);      1042   AddElementByWeightFraction(58, 0.002001);
1057                                                  1043 
1058   AddMaterial("G4_CESIUM_FLUORIDE", 4.115, 0,    1044   AddMaterial("G4_CESIUM_FLUORIDE", 4.115, 0, 440.7, 2);
1059   AddElementByAtomCount("Cs",  1);            << 1045   AddElementByWeightFraction( 9, 0.125069);
1060   AddElementByAtomCount("F" ,  1);            << 1046   AddElementByWeightFraction(55, 0.874931);
1061                                                  1047 
1062   AddMaterial("G4_CESIUM_IODIDE", 4.51, 0, 55    1048   AddMaterial("G4_CESIUM_IODIDE", 4.51, 0, 553.1, 2);
1063   AddElementByAtomCount("Cs",  1);            << 1049   AddElementByWeightFraction(53, 0.488451);
1064   AddElementByAtomCount("I" ,  1);            << 1050   AddElementByWeightFraction(55, 0.511549);
1065                                                  1051 
1066   AddMaterial("G4_CHLOROBENZENE", 1.1058, 0,     1052   AddMaterial("G4_CHLOROBENZENE", 1.1058, 0, 89.1, 3);
1067   AddElementByAtomCount("C" ,  6);            << 1053   AddElementByWeightFraction( 1, 0.044772);
1068   AddElementByAtomCount("H" ,  5);            << 1054   AddElementByWeightFraction( 6, 0.640254);
1069   AddElementByAtomCount("Cl",  1);            << 1055   AddElementByWeightFraction(17, 0.314974);
1070                                                  1056 
1071   AddMaterial("G4_CHLOROFORM", 1.4832, 0, 156    1057   AddMaterial("G4_CHLOROFORM", 1.4832, 0, 156., 3);
1072   AddElementByAtomCount("C" ,  1);            << 1058   AddElementByWeightFraction( 1, 0.008443);
1073   AddElementByAtomCount("H" ,  1);            << 1059   AddElementByWeightFraction( 6, 0.100613);
1074   AddElementByAtomCount("Cl",  3);            << 1060   AddElementByWeightFraction(17, 0.890944);
1075                                                  1061 
1076   AddMaterial("G4_CONCRETE", 2.3, 0, 135.2, 1    1062   AddMaterial("G4_CONCRETE", 2.3, 0, 135.2, 10);
1077   AddElementByWeightFraction( 1, 0.01    );      1063   AddElementByWeightFraction( 1, 0.01    );
1078   AddElementByWeightFraction( 6, 0.001   );      1064   AddElementByWeightFraction( 6, 0.001   );
1079   AddElementByWeightFraction( 8, 0.529107);      1065   AddElementByWeightFraction( 8, 0.529107);
1080   AddElementByWeightFraction(11, 0.016   );      1066   AddElementByWeightFraction(11, 0.016   );
1081   AddElementByWeightFraction(12, 0.002   );      1067   AddElementByWeightFraction(12, 0.002   );
1082   AddElementByWeightFraction(13, 0.033872);      1068   AddElementByWeightFraction(13, 0.033872);
1083   AddElementByWeightFraction(14, 0.337021);      1069   AddElementByWeightFraction(14, 0.337021);
1084   AddElementByWeightFraction(19, 0.013   );      1070   AddElementByWeightFraction(19, 0.013   );
1085   AddElementByWeightFraction(20, 0.044   );      1071   AddElementByWeightFraction(20, 0.044   );
1086   AddElementByWeightFraction(26, 0.014   );      1072   AddElementByWeightFraction(26, 0.014   );
1087                                                  1073 
1088   AddMaterial("G4_CYCLOHEXANE", 0.779, 0, 56.    1074   AddMaterial("G4_CYCLOHEXANE", 0.779, 0, 56.4, 2);
1089   AddElementByAtomCount("C" ,  6);            << 1075   AddElementByWeightFraction( 1, 0.143711);
1090   AddElementByAtomCount("H" , 12);            << 1076   AddElementByWeightFraction( 6, 0.856289);
1091                                                  1077 
1092   AddMaterial("G4_1,2-DICHLOROBENZENE", 1.304    1078   AddMaterial("G4_1,2-DICHLOROBENZENE", 1.3048, 0, 106.5, 3);
1093   AddElementByAtomCount("C" ,  6);            << 1079   AddElementByWeightFraction( 1, 0.027425);
1094   AddElementByAtomCount("H" ,  4);            << 1080   AddElementByWeightFraction( 6, 0.490233);
1095   AddElementByAtomCount("Cl",  2);            << 1081   AddElementByWeightFraction(17, 0.482342);
1096                                                  1082 
1097   AddMaterial("G4_DICHLORODIETHYL_ETHER", 1.2    1083   AddMaterial("G4_DICHLORODIETHYL_ETHER", 1.2199, 0, 103.3, 4);
1098   AddElementByAtomCount("C" ,  4);            << 1084   AddElementByWeightFraction( 1, 0.056381);
1099   AddElementByAtomCount("H" ,  8);            << 1085   AddElementByWeightFraction( 6, 0.335942);
1100   AddElementByAtomCount("O" ,  1);            << 1086   AddElementByWeightFraction( 8, 0.111874);
1101   AddElementByAtomCount("Cl",  2);            << 1087   AddElementByWeightFraction(17, 0.495802);
1102                                               << 1088 
1103   AddMaterial("G4_1,2-DICHLOROETHANE", 1.2351 << 1089   AddMaterial("G4_1,2-DICHLOROETHANE", 1.2351, 0, 111.9, 3);
1104   AddElementByAtomCount("C" ,  2);            << 1090   AddElementByWeightFraction( 1, 0.04074 );
1105   AddElementByAtomCount("H" ,  4);            << 1091   AddElementByWeightFraction( 6, 0.242746);
1106   AddElementByAtomCount("Cl",  2);            << 1092   AddElementByWeightFraction(17, 0.716515);
1107                                               << 1093 
1108   AddMaterial("G4_DIETHYL_ETHER", 0.71378, 0, << 1094   AddMaterial("G4_DIETHYL_ETHER", 0.71378, 0, 60., 3);
1109   AddElementByAtomCount("C" ,  4);            << 1095   AddElementByWeightFraction( 1, 0.135978);
1110   AddElementByAtomCount("H" , 10);            << 1096   AddElementByWeightFraction( 6, 0.648171);
1111   AddElementByAtomCount("O" ,  1);            << 1097   AddElementByWeightFraction( 8, 0.215851);
1112                                               << 1098 
1113   AddMaterial("G4_N,N-DIMETHYL_FORMAMIDE", 0. << 1099   AddMaterial("G4_N,N-DIMETHYL_FORMAMIDE", 0.9487, 0, 66.6, 4);
1114   AddElementByAtomCount("C" ,  3);            << 1100   AddElementByWeightFraction( 1, 0.096523);
1115   AddElementByAtomCount("H" ,  7);            << 1101   AddElementByWeightFraction( 6, 0.492965);
1116   AddElementByAtomCount("N" ,  1);            << 1102   AddElementByWeightFraction( 7, 0.191625);
1117   AddElementByAtomCount("O" ,  1);            << 1103   AddElementByWeightFraction( 8, 0.218887);
1118                                               << 1104 
1119   AddMaterial("G4_DIMETHYL_SULFOXIDE", 1.1014 << 1105   AddMaterial("G4_DIMETHYL_SULFOXIDE", 1.1014, 0, 98.6, 4);
1120   AddElementByAtomCount("C" ,  2);            << 1106   AddElementByWeightFraction( 1, 0.077403);
1121   AddElementByAtomCount("H" ,  6);            << 1107   AddElementByWeightFraction( 6, 0.307467);
1122   AddElementByAtomCount("O" ,  1);            << 1108   AddElementByWeightFraction( 8, 0.204782);
1123   AddElementByAtomCount("S" ,  1);            << 1109   AddElementByWeightFraction(16, 0.410348);
1124                                                  1110 
1125   AddMaterial("G4_ETHANE", 0.00125324, 0, 45.    1111   AddMaterial("G4_ETHANE", 0.00125324, 0, 45.4, 2, kStateGas);
1126   AddElementByAtomCount("C" ,  2);            << 1112   AddElementByWeightFraction( 1, 0.201115);
1127   AddElementByAtomCount("H" ,  6);            << 1113   AddElementByWeightFraction( 6, 0.798885);
1128                                                  1114 
1129   AddMaterial("G4_ETHYL_ALCOHOL", 0.7893, 0,  << 1115   AddMaterial("G4_ETHYL_ALCOHOL", 0.7893, 0, 62.9, 3);
1130   AddElementByAtomCount("C" ,  2);            << 1116   AddElementByWeightFraction( 1, 0.131269);
1131   AddElementByAtomCount("H" ,  6);            << 1117   AddElementByWeightFraction( 6, 0.521438);
1132   AddElementByAtomCount("O" ,  1);            << 1118   AddElementByWeightFraction( 8, 0.347294);
1133                                                  1119 
1134   AddMaterial("G4_ETHYL_CELLULOSE", 1.13, 0,     1120   AddMaterial("G4_ETHYL_CELLULOSE", 1.13, 0, 69.3, 3);
1135   AddElementByWeightFraction( 1, 0.090027);      1121   AddElementByWeightFraction( 1, 0.090027);
1136   AddElementByWeightFraction( 6, 0.585182);      1122   AddElementByWeightFraction( 6, 0.585182);
1137   AddElementByWeightFraction( 8, 0.324791);      1123   AddElementByWeightFraction( 8, 0.324791);
1138                                                  1124 
1139   AddMaterial("G4_ETHYLENE", 0.00117497, 0, 5    1125   AddMaterial("G4_ETHYLENE", 0.00117497, 0, 50.7, 2, kStateGas);
1140   AddElementByAtomCount("C" ,  2);            << 1126   AddElementByWeightFraction( 1, 0.143711);
1141   AddElementByAtomCount("H" ,  4);            << 1127   AddElementByWeightFraction( 6, 0.856289);
1142                                                  1128 
1143   AddMaterial("G4_EYE_LENS_ICRP", 1.07, 0, 73    1129   AddMaterial("G4_EYE_LENS_ICRP", 1.07, 0, 73.3, 8);
1144   AddElementByWeightFraction( 1, 0.096);         1130   AddElementByWeightFraction( 1, 0.096);
1145   AddElementByWeightFraction( 6, 0.195);         1131   AddElementByWeightFraction( 6, 0.195);
1146   AddElementByWeightFraction( 7, 0.057);         1132   AddElementByWeightFraction( 7, 0.057);
1147   AddElementByWeightFraction( 8, 0.646);         1133   AddElementByWeightFraction( 8, 0.646);
1148   AddElementByWeightFraction(11, 0.001);         1134   AddElementByWeightFraction(11, 0.001);
1149   AddElementByWeightFraction(15, 0.001);         1135   AddElementByWeightFraction(15, 0.001);
1150   AddElementByWeightFraction(16, 0.003);         1136   AddElementByWeightFraction(16, 0.003);
1151   AddElementByWeightFraction(17, 0.001);         1137   AddElementByWeightFraction(17, 0.001);
1152                                                  1138 
1153   AddMaterial("G4_FERRIC_OXIDE", 5.2, 0, 227.    1139   AddMaterial("G4_FERRIC_OXIDE", 5.2, 0, 227.3, 2);
1154   AddElementByAtomCount("Fe",  2);            << 1140   AddElementByWeightFraction( 8, 0.300567);
1155   AddElementByAtomCount("O" ,  3);            << 1141   AddElementByWeightFraction(26, 0.699433);
1156                                                  1142 
1157   AddMaterial("G4_FERROBORIDE", 7.15, 0, 261.    1143   AddMaterial("G4_FERROBORIDE", 7.15, 0, 261., 2);
1158   AddElementByAtomCount("Fe",  1);            << 1144   AddElementByWeightFraction( 5, 0.162174);
1159   AddElementByAtomCount("B" ,  1);            << 1145   AddElementByWeightFraction(26, 0.837826);
1160                                                  1146 
1161   AddMaterial("G4_FERROUS_OXIDE", 5.7, 0, 248    1147   AddMaterial("G4_FERROUS_OXIDE", 5.7, 0, 248.6, 2);
1162   AddElementByAtomCount("Fe",  1);            << 1148   AddElementByWeightFraction( 8, 0.222689);
1163   AddElementByAtomCount("O" ,  1);            << 1149   AddElementByWeightFraction(26, 0.777311);
1164                                                  1150 
1165   AddMaterial("G4_FERROUS_SULFATE", 1.024, 0,    1151   AddMaterial("G4_FERROUS_SULFATE", 1.024, 0, 76.4, 7);
1166   AddElementByWeightFraction( 1, 0.108259);      1152   AddElementByWeightFraction( 1, 0.108259);
1167   AddElementByWeightFraction( 7, 2.7e-05 );      1153   AddElementByWeightFraction( 7, 2.7e-05 );
1168   AddElementByWeightFraction( 8, 0.878636);      1154   AddElementByWeightFraction( 8, 0.878636);
1169   AddElementByWeightFraction(11, 2.2e-05 );      1155   AddElementByWeightFraction(11, 2.2e-05 );
1170   AddElementByWeightFraction(16, 0.012968);      1156   AddElementByWeightFraction(16, 0.012968);
1171   AddElementByWeightFraction(17, 3.4e-05 );      1157   AddElementByWeightFraction(17, 3.4e-05 );
1172   AddElementByWeightFraction(26, 5.4e-05 );      1158   AddElementByWeightFraction(26, 5.4e-05 );
1173                                                  1159 
1174   AddMaterial("G4_FREON-12", 1.12, 0, 143., 3 << 1160   AddMaterial("G4_FREON-12", 1.12, 0, 143., 3);
1175   AddElementByWeightFraction( 6, 0.099335);      1161   AddElementByWeightFraction( 6, 0.099335);
1176   AddElementByWeightFraction( 9, 0.314247);      1162   AddElementByWeightFraction( 9, 0.314247);
1177   AddElementByWeightFraction(17, 0.586418);      1163   AddElementByWeightFraction(17, 0.586418);
1178                                                  1164 
1179   AddMaterial("G4_FREON-12B2", 1.8, 0, 284.9, << 1165   AddMaterial("G4_FREON-12B2", 1.8, 0, 284.9, 3);
1180   AddElementByWeightFraction( 6, 0.057245);      1166   AddElementByWeightFraction( 6, 0.057245);
1181   AddElementByWeightFraction( 9, 0.181096);      1167   AddElementByWeightFraction( 9, 0.181096);
1182   AddElementByWeightFraction(35, 0.761659);      1168   AddElementByWeightFraction(35, 0.761659);
1183                                                  1169 
1184   AddMaterial("G4_FREON-13", 0.95, 0, 126.6,  << 1170   AddMaterial("G4_FREON-13", 0.95, 0, 126.6, 3);
1185   AddElementByWeightFraction( 6, 0.114983);      1171   AddElementByWeightFraction( 6, 0.114983);
1186   AddElementByWeightFraction( 9, 0.545622);      1172   AddElementByWeightFraction( 9, 0.545622);
1187   AddElementByWeightFraction(17, 0.339396);      1173   AddElementByWeightFraction(17, 0.339396);
1188                                                  1174 
1189   AddMaterial("G4_FREON-13B1", 1.5, 0, 210.5, << 1175   AddMaterial("G4_FREON-13B1", 1.5, 0, 210.5, 3);
1190   AddElementByAtomCount("C" ,  1);            << 1176   AddElementByWeightFraction( 6, 0.080659);
1191   AddElementByAtomCount("F" ,  3);            << 1177   AddElementByWeightFraction( 9, 0.382749);
1192   AddElementByAtomCount("Br",  1);            << 1178   AddElementByWeightFraction(35, 0.536592);
1193                                                  1179 
1194   AddMaterial("G4_FREON-13I1", 1.8, 0, 293.5, << 1180   AddMaterial("G4_FREON-13I1", 1.8, 0, 293.5, 3);
1195   AddElementByWeightFraction( 6, 0.061309);      1181   AddElementByWeightFraction( 6, 0.061309);
1196   AddElementByWeightFraction( 9, 0.290924);      1182   AddElementByWeightFraction( 9, 0.290924);
1197   AddElementByWeightFraction(53, 0.647767);      1183   AddElementByWeightFraction(53, 0.647767);
1198                                                  1184 
1199   AddMaterial("G4_GADOLINIUM_OXYSULFIDE", 7.4    1185   AddMaterial("G4_GADOLINIUM_OXYSULFIDE", 7.44, 0, 493.3, 3);
1200   AddElementByAtomCount("Gd",  2);            << 1186   AddElementByWeightFraction( 8, 0.084528);
1201   AddElementByAtomCount("O" ,  2);            << 1187   AddElementByWeightFraction(16, 0.08469 );
1202   AddElementByAtomCount("S" ,  1);            << 1188   AddElementByWeightFraction(64, 0.830782);
1203                                                  1189 
1204   AddMaterial("G4_GALLIUM_ARSENIDE", 5.31, 0,    1190   AddMaterial("G4_GALLIUM_ARSENIDE", 5.31, 0, 384.9, 2);
1205   AddElementByAtomCount("Ga",  1);            << 1191   AddElementByWeightFraction(31, 0.482019);
1206   AddElementByAtomCount("As",  1);            << 1192   AddElementByWeightFraction(33, 0.517981);
1207                                                  1193 
1208   AddMaterial("G4_GEL_PHOTO_EMULSION", 1.2914    1194   AddMaterial("G4_GEL_PHOTO_EMULSION", 1.2914, 0, 74.8, 5);
1209   AddElementByWeightFraction( 1, 0.08118);       1195   AddElementByWeightFraction( 1, 0.08118);
1210   AddElementByWeightFraction( 6, 0.41606);       1196   AddElementByWeightFraction( 6, 0.41606);
1211   AddElementByWeightFraction( 7, 0.11124);       1197   AddElementByWeightFraction( 7, 0.11124);
1212   AddElementByWeightFraction( 8, 0.38064);       1198   AddElementByWeightFraction( 8, 0.38064);
1213   AddElementByWeightFraction(16, 0.01088);       1199   AddElementByWeightFraction(16, 0.01088);
1214                                                  1200 
1215   AddMaterial("G4_Pyrex_Glass", 2.23, 0, 134.    1201   AddMaterial("G4_Pyrex_Glass", 2.23, 0, 134., 6);
1216   AddElementByWeightFraction( 5, 0.040064);      1202   AddElementByWeightFraction( 5, 0.040064);
1217   AddElementByWeightFraction( 8, 0.539562);      1203   AddElementByWeightFraction( 8, 0.539562);
1218   AddElementByWeightFraction(11, 0.028191);      1204   AddElementByWeightFraction(11, 0.028191);
1219   AddElementByWeightFraction(13, 0.011644);      1205   AddElementByWeightFraction(13, 0.011644);
1220   AddElementByWeightFraction(14, 0.37722 );      1206   AddElementByWeightFraction(14, 0.37722 );
1221   AddElementByWeightFraction(19, 0.003321);      1207   AddElementByWeightFraction(19, 0.003321);
1222                                                  1208 
1223   AddMaterial("G4_GLASS_LEAD", 6.22, 0, 526.4    1209   AddMaterial("G4_GLASS_LEAD", 6.22, 0, 526.4, 5);
1224   AddElementByWeightFraction( 8, 0.156453);      1210   AddElementByWeightFraction( 8, 0.156453);
1225   AddElementByWeightFraction(14, 0.080866);      1211   AddElementByWeightFraction(14, 0.080866);
1226   AddElementByWeightFraction(22, 0.008092);      1212   AddElementByWeightFraction(22, 0.008092);
1227   AddElementByWeightFraction(33, 0.002651);      1213   AddElementByWeightFraction(33, 0.002651);
1228   AddElementByWeightFraction(82, 0.751938);      1214   AddElementByWeightFraction(82, 0.751938);
1229                                                  1215 
1230   AddMaterial("G4_GLASS_PLATE", 2.4, 0, 145.4    1216   AddMaterial("G4_GLASS_PLATE", 2.4, 0, 145.4, 4);
1231   AddElementByWeightFraction( 8, 0.4598  );      1217   AddElementByWeightFraction( 8, 0.4598  );
1232   AddElementByWeightFraction(11, 0.096441);      1218   AddElementByWeightFraction(11, 0.096441);
1233   AddElementByWeightFraction(14, 0.336553);      1219   AddElementByWeightFraction(14, 0.336553);
1234   AddElementByWeightFraction(20, 0.107205);      1220   AddElementByWeightFraction(20, 0.107205);
1235                                                  1221 
                                                   >> 1222   AddMaterial("G4_GLUCOSE", 1.54, 0, 77.2, 3);
                                                   >> 1223   AddElementByWeightFraction( 1, 0.071204);
                                                   >> 1224   AddElementByWeightFraction( 6, 0.363652);
                                                   >> 1225   AddElementByWeightFraction( 8, 0.565144);
                                                   >> 1226 
1236   AddMaterial("G4_GLUTAMINE", 1.46, 0, 73.3,     1227   AddMaterial("G4_GLUTAMINE", 1.46, 0, 73.3, 4);
1237   AddElementByAtomCount("C" ,  5);            << 1228   AddElementByWeightFraction( 1, 0.068965);
1238   AddElementByAtomCount("H" , 10);            << 1229   AddElementByWeightFraction( 6, 0.410926);
1239   AddElementByAtomCount("N" ,  2);            << 1230   AddElementByWeightFraction( 7, 0.191681);
1240   AddElementByAtomCount("O" ,  3);            << 1231   AddElementByWeightFraction( 8, 0.328427);
1241                                                  1232 
1242   AddMaterial("G4_GLYCEROL", 1.2613, 0, 72.6,    1233   AddMaterial("G4_GLYCEROL", 1.2613, 0, 72.6, 3);
1243   AddElementByAtomCount("C" ,  3);            << 1234   AddElementByWeightFraction( 1, 0.087554);
1244   AddElementByAtomCount("H" ,  8);            << 1235   AddElementByWeightFraction( 6, 0.391262);
1245   AddElementByAtomCount("O" ,  3);            << 1236   AddElementByWeightFraction( 8, 0.521185);
1246                                               << 1237 
1247   //Tan, Z., et al. NIMB,2006(248)            << 1238   AddMaterial("G4_GUANINE", 2.2/*1.58*/, 0, 75. ,4);
1248   AddMaterial("G4_GUANINE", 1.58, 0, 75. ,4); << 1239   AddElementByAtomCount("H",5 );
1249   AddElementByAtomCount("C" ,  5);            << 1240   AddElementByAtomCount("C",5 );
1250   AddElementByAtomCount("H" ,  5);            << 1241   AddElementByAtomCount("N",5 );
1251   AddElementByAtomCount("N" ,  5);            << 1242   AddElementByAtomCount("O",1 );
1252   AddElementByAtomCount("O" ,  1);            << 
1253                                                  1243 
1254   AddMaterial("G4_GYPSUM", 2.32, 0, 129.7, 4)    1244   AddMaterial("G4_GYPSUM", 2.32, 0, 129.7, 4);
1255   AddElementByAtomCount("Ca",  1);            << 1245   AddElementByWeightFraction( 1, 0.023416);
1256   AddElementByAtomCount("S" ,  1);            << 1246   AddElementByWeightFraction( 8, 0.557572);
1257   AddElementByAtomCount("O" ,  6);            << 1247   AddElementByWeightFraction(16, 0.186215);
1258   AddElementByAtomCount("H" ,  4);            << 1248   AddElementByWeightFraction(20, 0.232797);
1259                                                  1249 
1260   AddMaterial("G4_N-HEPTANE", 0.68376, 0, 54.    1250   AddMaterial("G4_N-HEPTANE", 0.68376, 0, 54.4, 2);
1261   AddElementByAtomCount("C" ,  7);            << 1251   AddElementByWeightFraction( 1, 0.160937);
1262   AddElementByAtomCount("H" , 16);            << 1252   AddElementByWeightFraction( 6, 0.839063);
1263                                                  1253 
1264   AddMaterial("G4_N-HEXANE", 0.6603, 0, 54.,     1254   AddMaterial("G4_N-HEXANE", 0.6603, 0, 54., 2);
1265   AddElementByAtomCount("C" ,  6);            << 1255   AddElementByWeightFraction( 1, 0.163741);
1266   AddElementByAtomCount("H" , 14);            << 1256   AddElementByWeightFraction( 6, 0.836259);
1267                                                  1257 
1268   AddMaterial("G4_KAPTON", 1.42, 0, 79.6, 4);    1258   AddMaterial("G4_KAPTON", 1.42, 0, 79.6, 4);
1269   AddElementByAtomCount("C" , 22);            << 1259   AddElementByWeightFraction( 1, 0.026362);
1270   AddElementByAtomCount("H" , 10);            << 1260   AddElementByWeightFraction( 6, 0.691133);
1271   AddElementByAtomCount("N" ,  2);            << 1261   AddElementByWeightFraction( 7, 0.07327 );
1272   AddElementByAtomCount("O" ,  5);            << 1262   AddElementByWeightFraction( 8, 0.209235);
1273                                                  1263 
1274   AddMaterial("G4_LANTHANUM_OXYBROMIDE", 6.28    1264   AddMaterial("G4_LANTHANUM_OXYBROMIDE", 6.28, 0, 439.7, 3);
1275   AddElementByAtomCount("La",  1);            << 1265   AddElementByWeightFraction( 8, 0.068138);
1276   AddElementByAtomCount("Br",  1);            << 1266   AddElementByWeightFraction(35, 0.340294);
1277   AddElementByAtomCount("O" ,  1);            << 1267   AddElementByWeightFraction(57, 0.591568);
1278                                                  1268 
1279   AddMaterial("G4_LANTHANUM_OXYSULFIDE", 5.86    1269   AddMaterial("G4_LANTHANUM_OXYSULFIDE", 5.86, 0, 421.2, 3);
1280   AddElementByAtomCount("La",  2);            << 1270   AddElementByWeightFraction( 8, 0.0936  );
1281   AddElementByAtomCount("O" ,  2);            << 1271   AddElementByWeightFraction(16, 0.093778);
1282   AddElementByAtomCount("S" ,  1);            << 1272   AddElementByWeightFraction(57, 0.812622);
1283                                                  1273 
1284   AddMaterial("G4_LEAD_OXIDE", 9.53, 0, 766.7    1274   AddMaterial("G4_LEAD_OXIDE", 9.53, 0, 766.7, 2);
1285   AddElementByWeightFraction( 8, 0.071682);      1275   AddElementByWeightFraction( 8, 0.071682);
1286   AddElementByWeightFraction(82, 0.928318);      1276   AddElementByWeightFraction(82, 0.928318);
1287                                                  1277 
1288   AddMaterial("G4_LITHIUM_AMIDE", 1.178, 0, 5    1278   AddMaterial("G4_LITHIUM_AMIDE", 1.178, 0, 55.5, 3);
1289   AddElementByAtomCount("Li",  1);            << 1279   AddElementByWeightFraction( 1, 0.087783);
1290   AddElementByAtomCount("N" ,  1);            << 1280   AddElementByWeightFraction( 3, 0.302262);
1291   AddElementByAtomCount("H" ,  2);            << 1281   AddElementByWeightFraction( 7, 0.609955);
1292                                                  1282 
1293   AddMaterial("G4_LITHIUM_CARBONATE", 2.11, 0    1283   AddMaterial("G4_LITHIUM_CARBONATE", 2.11, 0, 87.9, 3);
1294   AddElementByAtomCount("Li",  2);            << 1284   AddElementByWeightFraction( 3, 0.187871);
1295   AddElementByAtomCount("C" ,  1);            << 1285   AddElementByWeightFraction( 6, 0.16255 );
1296   AddElementByAtomCount("O" ,  3);            << 1286   AddElementByWeightFraction( 8, 0.649579);
1297                                                  1287 
1298   AddMaterial("G4_LITHIUM_FLUORIDE", 2.635, 0    1288   AddMaterial("G4_LITHIUM_FLUORIDE", 2.635, 0, 94., 2);
1299   AddElementByAtomCount("Li",  1);            << 1289   AddElementByWeightFraction( 3, 0.267585);
1300   AddElementByAtomCount("F" ,  1);            << 1290   AddElementByWeightFraction( 9, 0.732415);
1301                                                  1291 
1302   AddMaterial("G4_LITHIUM_HYDRIDE", 0.82, 0,     1292   AddMaterial("G4_LITHIUM_HYDRIDE", 0.82, 0, 36.5, 2);
1303   AddElementByAtomCount("Li",  1);            << 1293   AddElementByWeightFraction( 1, 0.126797);
1304   AddElementByAtomCount("H" ,  1);            << 1294   AddElementByWeightFraction( 3, 0.873203);
1305                                                  1295 
1306   AddMaterial("G4_LITHIUM_IODIDE", 3.494, 0,     1296   AddMaterial("G4_LITHIUM_IODIDE", 3.494, 0, 485.1, 2);
1307   AddElementByAtomCount("Li",  1);            << 1297   AddElementByWeightFraction( 3, 0.051858);
1308   AddElementByAtomCount("I" ,  1);            << 1298   AddElementByWeightFraction(53, 0.948142);
1309                                                  1299 
1310   AddMaterial("G4_LITHIUM_OXIDE", 2.013, 0, 7    1300   AddMaterial("G4_LITHIUM_OXIDE", 2.013, 0, 73.6, 2);
1311   AddElementByAtomCount("Li",  2);            << 1301   AddElementByWeightFraction( 3, 0.46457);
1312   AddElementByAtomCount("O" ,  1);            << 1302   AddElementByWeightFraction( 8, 0.53543);
1313                                                  1303 
1314   AddMaterial("G4_LITHIUM_TETRABORATE", 2.44,    1304   AddMaterial("G4_LITHIUM_TETRABORATE", 2.44, 0, 94.6, 3);
1315   AddElementByAtomCount("Li",  2);            << 1305   AddElementByWeightFraction( 3, 0.082085);
1316   AddElementByAtomCount("B" ,  4);            << 1306   AddElementByWeightFraction( 5, 0.25568 );
1317   AddElementByAtomCount("O" ,  7);            << 1307   AddElementByWeightFraction( 8, 0.662235);
1318 }                                             << 
1319                                                  1308 
1320 void G4NistMaterialBuilder::NistCompoundMater << 
1321 {                                             << 
1322   //Adult Lung congested                         1309   //Adult Lung congested
1323   AddMaterial("G4_LUNG_ICRP", 1.04, 0, 75.3,     1310   AddMaterial("G4_LUNG_ICRP", 1.04, 0, 75.3, 9);
1324   AddElementByWeightFraction( 1, 0.105);         1311   AddElementByWeightFraction( 1, 0.105);
1325   AddElementByWeightFraction( 6, 0.083);         1312   AddElementByWeightFraction( 6, 0.083);
1326   AddElementByWeightFraction( 7, 0.023);         1313   AddElementByWeightFraction( 7, 0.023);
1327   AddElementByWeightFraction( 8, 0.779);         1314   AddElementByWeightFraction( 8, 0.779);
1328   AddElementByWeightFraction(11, 0.002);         1315   AddElementByWeightFraction(11, 0.002);
1329   AddElementByWeightFraction(15, 0.001);         1316   AddElementByWeightFraction(15, 0.001);
1330   AddElementByWeightFraction(16, 0.002);         1317   AddElementByWeightFraction(16, 0.002);
1331   AddElementByWeightFraction(17, 0.003);         1318   AddElementByWeightFraction(17, 0.003);
1332   AddElementByWeightFraction(19, 0.002);         1319   AddElementByWeightFraction(19, 0.002);
1333                                                  1320 
1334   AddMaterial("G4_M3_WAX", 1.05, 0, 67.9, 5);    1321   AddMaterial("G4_M3_WAX", 1.05, 0, 67.9, 5);
1335   AddElementByWeightFraction( 1, 0.114318);      1322   AddElementByWeightFraction( 1, 0.114318);
1336   AddElementByWeightFraction( 6, 0.655823);      1323   AddElementByWeightFraction( 6, 0.655823);
1337   AddElementByWeightFraction( 8, 0.092183);      1324   AddElementByWeightFraction( 8, 0.092183);
1338   AddElementByWeightFraction(12, 0.134792);      1325   AddElementByWeightFraction(12, 0.134792);
1339   AddElementByWeightFraction(20, 0.002883);      1326   AddElementByWeightFraction(20, 0.002883);
1340                                                  1327 
1341   AddMaterial("G4_MAGNESIUM_CARBONATE", 2.958    1328   AddMaterial("G4_MAGNESIUM_CARBONATE", 2.958, 0, 118., 3);
1342   AddElementByAtomCount("Mg",  1);            << 1329   AddElementByWeightFraction( 6, 0.142455);
1343   AddElementByAtomCount("C" ,  1);            << 1330   AddElementByWeightFraction( 8, 0.569278);
1344   AddElementByAtomCount("O" ,  3);            << 1331   AddElementByWeightFraction(12, 0.288267);
1345                                                  1332 
1346   AddMaterial("G4_MAGNESIUM_FLUORIDE", 3.0, 0    1333   AddMaterial("G4_MAGNESIUM_FLUORIDE", 3.0, 0, 134.3, 2);
1347   AddElementByAtomCount("Mg",  1);            << 1334   AddElementByWeightFraction( 9, 0.609883);
1348   AddElementByAtomCount("F" ,  2);            << 1335   AddElementByWeightFraction(12, 0.390117);
1349                                                  1336 
1350   AddMaterial("G4_MAGNESIUM_OXIDE", 3.58, 0,     1337   AddMaterial("G4_MAGNESIUM_OXIDE", 3.58, 0, 143.8, 2);
1351   AddElementByAtomCount("Mg",  1);            << 1338   AddElementByWeightFraction( 8, 0.396964);
1352   AddElementByAtomCount("O" ,  1);            << 1339   AddElementByWeightFraction(12, 0.603036);
1353                                                  1340 
1354   AddMaterial("G4_MAGNESIUM_TETRABORATE", 2.5    1341   AddMaterial("G4_MAGNESIUM_TETRABORATE", 2.53, 0, 108.3, 3);
1355   AddElementByAtomCount("Mg",  1);            << 1342   AddElementByWeightFraction( 5, 0.240837);
1356   AddElementByAtomCount("B" ,  4);            << 1343   AddElementByWeightFraction( 8, 0.62379);
1357   AddElementByAtomCount("O" ,  7);            << 1344   AddElementByWeightFraction(12, 0.135373);
1358                                                  1345 
1359   AddMaterial("G4_MERCURIC_IODIDE", 6.36, 0,     1346   AddMaterial("G4_MERCURIC_IODIDE", 6.36, 0, 684.5, 2);
1360   AddElementByAtomCount("Hg",  1);            << 1347   AddElementByWeightFraction(53, 0.55856);
1361   AddElementByAtomCount("I" ,  2);            << 1348   AddElementByWeightFraction(80, 0.44144);
1362                                                  1349 
1363   AddMaterial("G4_METHANE", 0.000667151, 0, 4    1350   AddMaterial("G4_METHANE", 0.000667151, 0, 41.7, 2, kStateGas);
1364   AddElementByAtomCount("C" ,  1);            << 1351   AddElementByWeightFraction( 1, 0.251306);
1365   AddElementByAtomCount("H" ,  4);            << 1352   AddElementByWeightFraction( 6, 0.748694);
1366                                                  1353 
1367   AddMaterial("G4_METHANOL", 0.7914, 0, 67.6, << 1354   AddMaterial("G4_METHANOL", 0.7914, 0, 67.6, 3);
1368   AddElementByAtomCount("C" ,  1);            << 1355   AddElementByWeightFraction( 1, 0.125822);
1369   AddElementByAtomCount("H" ,  4);            << 1356   AddElementByWeightFraction( 6, 0.374852);
1370   AddElementByAtomCount("O" ,  1);            << 1357   AddElementByWeightFraction( 8, 0.499326);
1371                                                  1358 
1372   AddMaterial("G4_MIX_D_WAX", 0.99, 0, 60.9,     1359   AddMaterial("G4_MIX_D_WAX", 0.99, 0, 60.9, 5);
1373   AddElementByWeightFraction( 1, 0.13404 );      1360   AddElementByWeightFraction( 1, 0.13404 );
1374   AddElementByWeightFraction( 6, 0.77796 );      1361   AddElementByWeightFraction( 6, 0.77796 );
1375   AddElementByWeightFraction( 8, 0.03502 );      1362   AddElementByWeightFraction( 8, 0.03502 );
1376   AddElementByWeightFraction(12, 0.038594);      1363   AddElementByWeightFraction(12, 0.038594);
1377   AddElementByWeightFraction(22, 0.014386);      1364   AddElementByWeightFraction(22, 0.014386);
1378                                                  1365 
1379   AddMaterial("G4_MS20_TISSUE", 1.0, 0, 75.1,    1366   AddMaterial("G4_MS20_TISSUE", 1.0, 0, 75.1, 6);
1380   AddElementByWeightFraction( 1, 0.081192);      1367   AddElementByWeightFraction( 1, 0.081192);
1381   AddElementByWeightFraction( 6, 0.583442);      1368   AddElementByWeightFraction( 6, 0.583442);
1382   AddElementByWeightFraction( 7, 0.017798);      1369   AddElementByWeightFraction( 7, 0.017798);
1383   AddElementByWeightFraction( 8, 0.186381);      1370   AddElementByWeightFraction( 8, 0.186381);
1384   AddElementByWeightFraction(12, 0.130287);      1371   AddElementByWeightFraction(12, 0.130287);
1385   AddElementByWeightFraction(17, 0.0009  );      1372   AddElementByWeightFraction(17, 0.0009  );
1386                                                  1373 
1387   AddMaterial("G4_MUSCLE_SKELETAL_ICRP", 1.05    1374   AddMaterial("G4_MUSCLE_SKELETAL_ICRP", 1.05, 0, 75.3, 9);
1388   AddElementByWeightFraction( 1, 0.102);         1375   AddElementByWeightFraction( 1, 0.102);
1389   AddElementByWeightFraction( 6, 0.143);         1376   AddElementByWeightFraction( 6, 0.143);
1390   AddElementByWeightFraction( 7, 0.034);         1377   AddElementByWeightFraction( 7, 0.034);
1391   AddElementByWeightFraction( 8, 0.710);         1378   AddElementByWeightFraction( 8, 0.710);
1392   AddElementByWeightFraction(11, 0.001);         1379   AddElementByWeightFraction(11, 0.001);
1393   AddElementByWeightFraction(15, 0.002);         1380   AddElementByWeightFraction(15, 0.002);
1394   AddElementByWeightFraction(16, 0.003);         1381   AddElementByWeightFraction(16, 0.003);
1395   AddElementByWeightFraction(17, 0.001);         1382   AddElementByWeightFraction(17, 0.001);
1396   AddElementByWeightFraction(19, 0.004);         1383   AddElementByWeightFraction(19, 0.004);
1397                                                  1384 
1398   // from old ICRU report                        1385   // from old ICRU report
1399   AddMaterial("G4_MUSCLE_STRIATED_ICRU", 1.04    1386   AddMaterial("G4_MUSCLE_STRIATED_ICRU", 1.04, 0, 74.7, 8);
1400   AddElementByWeightFraction( 1, 0.102);         1387   AddElementByWeightFraction( 1, 0.102);
1401   AddElementByWeightFraction( 6, 0.123);         1388   AddElementByWeightFraction( 6, 0.123);
1402   AddElementByWeightFraction( 7, 0.035);         1389   AddElementByWeightFraction( 7, 0.035);
1403   AddElementByWeightFraction( 8, 0.729);         1390   AddElementByWeightFraction( 8, 0.729);
1404   AddElementByWeightFraction(11, 0.001);         1391   AddElementByWeightFraction(11, 0.001);
1405   AddElementByWeightFraction(15, 0.002);         1392   AddElementByWeightFraction(15, 0.002);
1406   AddElementByWeightFraction(16, 0.004);         1393   AddElementByWeightFraction(16, 0.004);
1407   AddElementByWeightFraction(19, 0.003);         1394   AddElementByWeightFraction(19, 0.003);
1408                                                  1395 
1409   AddMaterial("G4_MUSCLE_WITH_SUCROSE", 1.11,    1396   AddMaterial("G4_MUSCLE_WITH_SUCROSE", 1.11, 0, 74.3, 4);
1410   AddElementByWeightFraction( 1, 0.098234);      1397   AddElementByWeightFraction( 1, 0.098234);
1411   AddElementByWeightFraction( 6, 0.156214);      1398   AddElementByWeightFraction( 6, 0.156214);
1412   AddElementByWeightFraction( 7, 0.035451);      1399   AddElementByWeightFraction( 7, 0.035451);
1413   AddElementByWeightFraction( 8, 0.7101  );      1400   AddElementByWeightFraction( 8, 0.7101  );
1414                                               << 1401   
1415   AddMaterial("G4_MUSCLE_WITHOUT_SUCROSE", 1.    1402   AddMaterial("G4_MUSCLE_WITHOUT_SUCROSE", 1.07, 0, 74.2, 4);
1416   AddElementByWeightFraction( 1, 0.101969);      1403   AddElementByWeightFraction( 1, 0.101969);
1417   AddElementByWeightFraction( 6, 0.120058);      1404   AddElementByWeightFraction( 6, 0.120058);
1418   AddElementByWeightFraction( 7, 0.035451);      1405   AddElementByWeightFraction( 7, 0.035451);
1419   AddElementByWeightFraction( 8, 0.742522);      1406   AddElementByWeightFraction( 8, 0.742522);
1420                                                  1407 
1421   AddMaterial("G4_NAPHTHALENE", 1.145, 0, 68.    1408   AddMaterial("G4_NAPHTHALENE", 1.145, 0, 68.4, 2);
1422   AddElementByAtomCount("C" , 10);            << 1409   AddElementByWeightFraction( 1, 0.062909);
1423   AddElementByAtomCount("H" ,  8);            << 1410   AddElementByWeightFraction( 6, 0.937091);
1424                                                  1411 
1425   AddMaterial("G4_NITROBENZENE", 1.19867, 0,     1412   AddMaterial("G4_NITROBENZENE", 1.19867, 0, 75.8, 4);
1426   AddElementByAtomCount("C" ,  6);            << 1413   AddElementByWeightFraction( 1, 0.040935);
1427   AddElementByAtomCount("H" ,  5);            << 1414   AddElementByWeightFraction( 6, 0.585374);
1428   AddElementByAtomCount("N" ,  1);            << 1415   AddElementByWeightFraction( 7, 0.113773);
1429   AddElementByAtomCount("O" ,  2);            << 1416   AddElementByWeightFraction( 8, 0.259918);
1430                                                  1417 
1431   AddMaterial("G4_NITROUS_OXIDE", 0.00183094,    1418   AddMaterial("G4_NITROUS_OXIDE", 0.00183094, 0, 84.9, 2, kStateGas);
1432   AddElementByAtomCount("N" ,  2);            << 1419   AddElementByWeightFraction( 7, 0.636483);
1433   AddElementByAtomCount("O" ,  1);            << 1420   AddElementByWeightFraction( 8, 0.363517);
1434                                                  1421 
1435   AddMaterial("G4_NYLON-8062", 1.08, 0, 64.3,    1422   AddMaterial("G4_NYLON-8062", 1.08, 0, 64.3, 4);
1436   AddElementByWeightFraction( 1, 0.103509);      1423   AddElementByWeightFraction( 1, 0.103509);
1437   AddElementByWeightFraction( 6, 0.648415);      1424   AddElementByWeightFraction( 6, 0.648415);
1438   AddElementByWeightFraction( 7, 0.099536);      1425   AddElementByWeightFraction( 7, 0.099536);
1439   AddElementByWeightFraction( 8, 0.148539);      1426   AddElementByWeightFraction( 8, 0.148539);
1440                                                  1427 
1441   AddMaterial("G4_NYLON-6-6", 1.14, 0, 63.9,     1428   AddMaterial("G4_NYLON-6-6", 1.14, 0, 63.9, 4);
1442   AddElementByAtomCount("C" ,  6);            << 1429   AddElementByWeightFraction( 1, 0.097976);
1443   AddElementByAtomCount("H" , 11);            << 1430   AddElementByWeightFraction( 6, 0.636856);
1444   AddElementByAtomCount("N" ,  1);            << 1431   AddElementByWeightFraction( 7, 0.123779);
1445   AddElementByAtomCount("O" ,  1);            << 1432   AddElementByWeightFraction( 8, 0.141389);
1446                                                  1433 
1447   AddMaterial("G4_NYLON-6-10", 1.14, 0, 63.2,    1434   AddMaterial("G4_NYLON-6-10", 1.14, 0, 63.2, 4);
1448   AddElementByWeightFraction( 1, 0.107062);      1435   AddElementByWeightFraction( 1, 0.107062);
1449   AddElementByWeightFraction( 6, 0.680449);      1436   AddElementByWeightFraction( 6, 0.680449);
1450   AddElementByWeightFraction( 7, 0.099189);      1437   AddElementByWeightFraction( 7, 0.099189);
1451   AddElementByWeightFraction( 8, 0.1133  );      1438   AddElementByWeightFraction( 8, 0.1133  );
1452                                                  1439 
1453   AddMaterial("G4_NYLON-11_RILSAN", 1.425, 0,    1440   AddMaterial("G4_NYLON-11_RILSAN", 1.425, 0, 61.6, 4);
1454   AddElementByWeightFraction( 1, 0.115476);      1441   AddElementByWeightFraction( 1, 0.115476);
1455   AddElementByWeightFraction( 6, 0.720819);      1442   AddElementByWeightFraction( 6, 0.720819);
1456   AddElementByWeightFraction( 7, 0.076417);      1443   AddElementByWeightFraction( 7, 0.076417);
1457   AddElementByWeightFraction( 8, 0.087289);      1444   AddElementByWeightFraction( 8, 0.087289);
1458                                                  1445 
1459   AddMaterial("G4_OCTANE", 0.7026, 0, 54.7, 2 << 1446   AddMaterial("G4_OCTANE", 0.7026, 0, 54.7, 2);
1460   AddElementByAtomCount("C" ,  8);            << 1447   AddElementByWeightFraction( 1, 0.158821);
1461   AddElementByAtomCount("H" , 18);            << 1448   AddElementByWeightFraction( 6, 0.841179);
1462                                                  1449 
1463   AddMaterial("G4_PARAFFIN", 0.93, 0, 55.9, 2    1450   AddMaterial("G4_PARAFFIN", 0.93, 0, 55.9, 2);
1464   AddElementByAtomCount("C" , 25);            << 1451   AddElementByWeightFraction( 1, 0.148605);
1465   AddElementByAtomCount("H" , 52);            << 1452   AddElementByWeightFraction( 6, 0.851395);
1466                                                  1453 
1467   AddMaterial("G4_N-PENTANE", 0.6262, 0, 53.6 << 1454   AddMaterial("G4_N-PENTANE", 0.6262, 0, 53.6, 2);
1468   AddElementByAtomCount("C" ,  5);            << 1455   AddElementByWeightFraction( 1, 0.167635);
1469   AddElementByAtomCount("H" , 12);            << 1456   AddElementByWeightFraction (6, 0.832365);
1470                                                  1457 
1471   AddMaterial("G4_PHOTO_EMULSION", 3.815, 0,     1458   AddMaterial("G4_PHOTO_EMULSION", 3.815, 0, 331., 8);
1472   AddElementByWeightFraction( 1, 0.0141  );      1459   AddElementByWeightFraction( 1, 0.0141  );
1473   AddElementByWeightFraction( 6, 0.072261);      1460   AddElementByWeightFraction( 6, 0.072261);
1474   AddElementByWeightFraction( 7, 0.01932 );      1461   AddElementByWeightFraction( 7, 0.01932 );
1475   AddElementByWeightFraction( 8, 0.066101);      1462   AddElementByWeightFraction( 8, 0.066101);
1476   AddElementByWeightFraction(16, 0.00189 );      1463   AddElementByWeightFraction(16, 0.00189 );
1477   AddElementByWeightFraction(35, 0.349103);      1464   AddElementByWeightFraction(35, 0.349103);
1478   AddElementByWeightFraction(47, 0.474105);      1465   AddElementByWeightFraction(47, 0.474105);
1479   AddElementByWeightFraction(53, 0.00312 );      1466   AddElementByWeightFraction(53, 0.00312 );
1480                                                  1467 
1481   AddMaterial("G4_PLASTIC_SC_VINYLTOLUENE", 1    1468   AddMaterial("G4_PLASTIC_SC_VINYLTOLUENE", 1.032, 0, 64.7, 2);
1482   // AddElementByWeightFraction( 1, 0.085);   << 1469   AddElementByWeightFraction( 1, 0.085);
1483   // AddElementByWeightFraction( 6, 0.915);   << 1470   AddElementByWeightFraction( 6, 0.915);
1484   // Watch out! These weight fractions do not << 
1485   // (PolyVinylToluene, C_9H_10) but to an un << 
1486   // M.Trocme & S.Seltzer                     << 
1487   AddElementByAtomCount("C" ,  9);            << 
1488   AddElementByAtomCount("H" , 10);            << 
1489                                                  1471 
1490   AddMaterial("G4_PLUTONIUM_DIOXIDE", 11.46,     1472   AddMaterial("G4_PLUTONIUM_DIOXIDE", 11.46, 0, 746.5, 2);
1491   AddElementByAtomCount("Pu",  1);            << 1473   AddElementByWeightFraction( 8, 0.118055);
1492   AddElementByAtomCount("O" ,  2);            << 1474   AddElementByWeightFraction(94, 0.881945);
1493                                                  1475 
1494   AddMaterial("G4_POLYACRYLONITRILE", 1.17, 0    1476   AddMaterial("G4_POLYACRYLONITRILE", 1.17, 0, 69.6, 3);
1495   AddElementByAtomCount("C" ,  3);            << 1477   AddElementByWeightFraction( 1, 0.056983);
1496   AddElementByAtomCount("H" ,  3);            << 1478   AddElementByWeightFraction( 6, 0.679056);
1497   AddElementByAtomCount("N" ,  1);            << 1479   AddElementByWeightFraction( 7, 0.263962);
1498                                                  1480 
1499   AddMaterial("G4_POLYCARBONATE", 1.2, 0, 73.    1481   AddMaterial("G4_POLYCARBONATE", 1.2, 0, 73.1, 3);
1500   AddElementByAtomCount("C" , 16);            << 1482   AddElementByWeightFraction( 1, 0.055491);
1501   AddElementByAtomCount("H" , 14);            << 1483   AddElementByWeightFraction( 6, 0.755751);
1502   AddElementByAtomCount("O" ,  3);            << 1484   AddElementByWeightFraction( 8, 0.188758);
1503                                                  1485 
1504   AddMaterial("G4_POLYCHLOROSTYRENE", 1.3, 0,    1486   AddMaterial("G4_POLYCHLOROSTYRENE", 1.3, 0, 81.7, 3);
1505   //  AddElementByWeightFraction( 1, 0.061869 << 1487   AddElementByWeightFraction( 1, 0.061869);
1506   //  AddElementByWeightFraction( 6, 0.696325 << 1488   AddElementByWeightFraction( 6, 0.696325);
1507   //  AddElementByWeightFraction(17, 0.241806 << 1489   AddElementByWeightFraction(17, 0.241806);
1508   //  These weight fractions correspond to C_ << 
1509   //  POLYCHLOROSTYRENE. POLYCHLOROSTYRENE is << 
1510   //  M.Trocme & S.Seltzer                    << 
1511   AddElementByAtomCount("C" ,  8);            << 
1512   AddElementByAtomCount("H" ,  7);            << 
1513   AddElementByAtomCount("Cl",  1);            << 
1514                                                  1490 
1515   AddMaterial("G4_POLYETHYLENE", 0.94, 0, 57.    1491   AddMaterial("G4_POLYETHYLENE", 0.94, 0, 57.4, 2);
1516   AddElementByAtomCount("C" ,  1);            << 1492   AddElementByWeightFraction( 1, 0.143711);
1517   AddElementByAtomCount("H" ,  2);            << 1493   AddElementByWeightFraction( 6, 0.856289);
1518   chFormulas[nMaterials-1] = "(C_2H_4)_N-Poly    1494   chFormulas[nMaterials-1] = "(C_2H_4)_N-Polyethylene";
1519                                                  1495 
1520   AddMaterial("G4_MYLAR", 1.4, 0, 78.7, 3);      1496   AddMaterial("G4_MYLAR", 1.4, 0, 78.7, 3);
1521   AddElementByAtomCount("C" , 10);            << 1497   AddElementByWeightFraction( 1, 0.041959);
1522   AddElementByAtomCount("H" ,  8);            << 1498   AddElementByWeightFraction( 6, 0.625017);
1523   AddElementByAtomCount("O" ,  4);            << 1499   AddElementByWeightFraction( 8, 0.333025);
1524                                                  1500 
1525   AddMaterial("G4_PLEXIGLASS", 1.19, 0, 74.,     1501   AddMaterial("G4_PLEXIGLASS", 1.19, 0, 74., 3);
1526   AddElementByAtomCount("C" ,  5);            << 1502   AddElementByWeightFraction( 1, 0.080538);
1527   AddElementByAtomCount("H" ,  8);            << 1503   AddElementByWeightFraction( 6, 0.599848);
1528   AddElementByAtomCount("O" ,  2);            << 1504   AddElementByWeightFraction( 8, 0.319614);
1529                                                  1505 
1530   AddMaterial("G4_POLYOXYMETHYLENE", 1.425 ,0    1506   AddMaterial("G4_POLYOXYMETHYLENE", 1.425 ,0, 77.4, 3);
1531   AddElementByAtomCount("C" ,  1);            << 1507   AddElementByWeightFraction( 1, 0.067135);
1532   AddElementByAtomCount("H" ,  2);            << 1508   AddElementByWeightFraction( 6, 0.400017);
1533   AddElementByAtomCount("O" ,  1);            << 1509   AddElementByWeightFraction( 8, 0.532848);
1534                                                  1510 
1535   AddMaterial("G4_POLYPROPYLENE", 0.9, 0, 56.    1511   AddMaterial("G4_POLYPROPYLENE", 0.9, 0, 56.5, 2);
1536   AddElementByAtomCount("C" ,  2);            << 1512   AddElementByWeightFraction( 1, 0.143711);
1537   AddElementByAtomCount("H" ,  4);            << 1513   AddElementByWeightFraction( 6, 0.856289);
1538   chFormulas[nMaterials-1] = "(C_2H_4)_N-Poly    1514   chFormulas[nMaterials-1] = "(C_2H_4)_N-Polypropylene";
1539                                                  1515 
1540   AddMaterial("G4_POLYSTYRENE", 1.06, 0, 68.7    1516   AddMaterial("G4_POLYSTYRENE", 1.06, 0, 68.7, 2);
1541   AddElementByAtomCount("C" ,  8);            << 1517   AddElementByWeightFraction( 1, 0.077418);
1542   AddElementByAtomCount("H" ,  8);            << 1518   AddElementByWeightFraction( 6, 0.922582);
1543                                                  1519 
1544   AddMaterial("G4_TEFLON", 2.2, 0, 99.1, 2);     1520   AddMaterial("G4_TEFLON", 2.2, 0, 99.1, 2);
1545   AddElementByAtomCount("C" ,  2);            << 1521   AddElementByWeightFraction( 6, 0.240183);
1546   AddElementByAtomCount("F" ,  4);            << 1522   AddElementByWeightFraction( 9, 0.759817);
1547                                                  1523 
1548   AddMaterial("G4_POLYTRIFLUOROCHLOROETHYLENE    1524   AddMaterial("G4_POLYTRIFLUOROCHLOROETHYLENE", 2.1, 0, 120.7, 3);
1549   // correct chemical name Polychlorotrifluor    1525   // correct chemical name Polychlorotrifluoroethylene [CF2CClF]n, IvantchenkoA.
1550   AddElementByAtomCount("C" ,  2);            << 1526   AddElementByWeightFraction( 6, 0.20625 );
1551   AddElementByAtomCount("F" ,  3);            << 1527   AddElementByWeightFraction( 9, 0.489354);
1552   AddElementByAtomCount("Cl",  1);            << 1528   AddElementByWeightFraction(17, 0.304395);
1553                                                  1529 
1554   AddMaterial("G4_POLYVINYL_ACETATE", 1.19, 0    1530   AddMaterial("G4_POLYVINYL_ACETATE", 1.19, 0, 73.7, 3);
1555   AddElementByAtomCount("C" ,  4);            << 1531   AddElementByWeightFraction( 1, 0.070245);
1556   AddElementByAtomCount("H" ,  6);            << 1532   AddElementByWeightFraction( 6, 0.558066);
1557   AddElementByAtomCount("O" ,  2);            << 1533   AddElementByWeightFraction( 8, 0.371689);
1558                                               << 1534 
1559   AddMaterial("G4_POLYVINYL_ALCOHOL", 1.3, 0, << 1535   AddMaterial("G4_POLYVINYL_ALCOHOL", 1.3, 0, 69.7, 3);
1560   AddElementByAtomCount("C" ,  2);            << 1536   AddElementByWeightFraction( 1, 0.091517);
1561   AddElementByAtomCount("H" ,  4);            << 1537   AddElementByWeightFraction( 6, 0.545298);
1562   AddElementByAtomCount("O" ,  1);            << 1538   AddElementByWeightFraction( 8, 0.363185);
1563                                                  1539 
1564   AddMaterial("G4_POLYVINYL_BUTYRAL", 1.12, 0    1540   AddMaterial("G4_POLYVINYL_BUTYRAL", 1.12, 0, 67.2, 3);
1565   //  AddElementByWeightFraction( 1, 0.092802 << 1541   AddElementByWeightFraction( 1, 0.092802);
1566   //  AddElementByWeightFraction( 6, 0.680561 << 1542   AddElementByWeightFraction( 6, 0.680561);
1567   //  AddElementByWeightFraction( 8, 0.226637 << 1543   AddElementByWeightFraction( 8, 0.226637);
1568   //  These weight fractions correspond to C_ << 
1569   //  POLYVINYL_BUTYRAL. POLYVINYL_BUTYRAL is << 
1570   //  M.Trocme & S.Seltzer                    << 
1571   AddElementByAtomCount("C" ,  8);            << 
1572   AddElementByAtomCount("H" , 14);            << 
1573   AddElementByAtomCount("O" ,  2);            << 
1574                                                  1544 
1575   AddMaterial("G4_POLYVINYL_CHLORIDE", 1.3, 0    1545   AddMaterial("G4_POLYVINYL_CHLORIDE", 1.3, 0, 108.2, 3);
1576   AddElementByAtomCount("C" ,  2);            << 1546   AddElementByWeightFraction( 1, 0.04838);
1577   AddElementByAtomCount("H" ,  3);            << 1547   AddElementByWeightFraction( 6, 0.38436);
1578   AddElementByAtomCount("Cl",  1);            << 1548   AddElementByWeightFraction(17, 0.56726);
1579                                                  1549 
1580   AddMaterial("G4_POLYVINYLIDENE_CHLORIDE", 1    1550   AddMaterial("G4_POLYVINYLIDENE_CHLORIDE", 1.7, 0, 134.3, 3);
1581   AddElementByAtomCount("C" ,  2);            << 1551   AddElementByWeightFraction( 1, 0.020793);
1582   AddElementByAtomCount("H" ,  2);            << 1552   AddElementByWeightFraction( 6, 0.247793);
1583   AddElementByAtomCount("Cl",  2);            << 1553   AddElementByWeightFraction(17, 0.731413);
1584                                                  1554 
1585   AddMaterial("G4_POLYVINYLIDENE_FLUORIDE", 1    1555   AddMaterial("G4_POLYVINYLIDENE_FLUORIDE", 1.76, 0, 88.8, 3);
1586   AddElementByAtomCount("C" ,  2);            << 1556   AddElementByWeightFraction( 1, 0.03148 );
1587   AddElementByAtomCount("H" ,  2);            << 1557   AddElementByWeightFraction( 6, 0.375141);
1588   AddElementByAtomCount("F" ,  2);            << 1558   AddElementByWeightFraction( 9, 0.593379);
1589                                                  1559 
1590   AddMaterial("G4_POLYVINYL_PYRROLIDONE", 1.2    1560   AddMaterial("G4_POLYVINYL_PYRROLIDONE", 1.25, 0, 67.7, 4);
1591   AddElementByAtomCount("C" ,  6);            << 1561   AddElementByWeightFraction( 1, 0.081616);
1592   AddElementByAtomCount("H" ,  9);            << 1562   AddElementByWeightFraction( 6, 0.648407);
1593   AddElementByAtomCount("N" ,  1);            << 1563   AddElementByWeightFraction( 7, 0.126024);
1594   AddElementByAtomCount("O" ,  1);            << 1564   AddElementByWeightFraction( 8, 0.143953);
1595                                                  1565 
1596   AddMaterial("G4_POTASSIUM_IODIDE", 3.13, 0,    1566   AddMaterial("G4_POTASSIUM_IODIDE", 3.13, 0, 431.9, 2);
1597   AddElementByAtomCount("K" ,  1);            << 1567   AddElementByWeightFraction(19, 0.235528);
1598   AddElementByAtomCount("I" ,  1);            << 1568   AddElementByWeightFraction(53, 0.764472);
1599                                                  1569 
1600   AddMaterial("G4_POTASSIUM_OXIDE", 2.32, 0,     1570   AddMaterial("G4_POTASSIUM_OXIDE", 2.32, 0, 189.9, 2);
1601   AddElementByAtomCount("K" ,  2);            << 1571   AddElementByWeightFraction( 8, 0.169852);
1602   AddElementByAtomCount("O" ,  1);            << 1572   AddElementByWeightFraction(19, 0.830148);
1603                                                  1573 
1604   AddMaterial("G4_PROPANE", 0.00187939, 0, 47    1574   AddMaterial("G4_PROPANE", 0.00187939, 0, 47.1, 2, kStateGas);
1605   AddElementByAtomCount("C" ,  3);            << 1575   AddElementByWeightFraction( 1, 0.182855);
1606   AddElementByAtomCount("H" ,  8);            << 1576   AddElementByWeightFraction( 6, 0.817145);
1607                                                  1577 
1608   AddMaterial("G4_lPROPANE", 0.43, 0, 52., 2, << 1578   AddMaterial("G4_lPROPANE", 0.43, 0, 52., 2);
1609   AddElementByAtomCount("C" ,  3);            << 1579   AddElementByWeightFraction( 1, 0.182855);
1610   AddElementByAtomCount("H" ,  8);            << 1580   AddElementByWeightFraction( 6, 0.817145);
1611                                               << 1581 
1612   AddMaterial("G4_N-PROPYL_ALCOHOL", 0.8035,  << 1582   AddMaterial("G4_N-PROPYL_ALCOHOL", 0.8035, 0, 61.1, 3);
1613   AddElementByAtomCount("C" ,  3);            << 1583   AddElementByWeightFraction( 1, 0.134173);
1614   AddElementByAtomCount("H" ,  8);            << 1584   AddElementByWeightFraction( 6, 0.599595);
1615   AddElementByAtomCount("O" ,  1);            << 1585   AddElementByWeightFraction( 8, 0.266232);
1616                                                  1586 
1617   AddMaterial("G4_PYRIDINE", 0.9819, 0, 66.2,    1587   AddMaterial("G4_PYRIDINE", 0.9819, 0, 66.2, 3);
1618   AddElementByAtomCount("C" ,  5);            << 1588   AddElementByWeightFraction( 1, 0.06371 );
1619   AddElementByAtomCount("H" ,  5);            << 1589   AddElementByWeightFraction( 6, 0.759217);
1620   AddElementByAtomCount("N" ,  1);            << 1590   AddElementByWeightFraction( 7, 0.177073);
1621                                                  1591 
1622   AddMaterial("G4_RUBBER_BUTYL", 0.92, 0, 56.    1592   AddMaterial("G4_RUBBER_BUTYL", 0.92, 0, 56.5, 2);
1623   AddElementByWeightFraction( 1, 0.143711);      1593   AddElementByWeightFraction( 1, 0.143711);
1624   AddElementByWeightFraction( 6, 0.856289);      1594   AddElementByWeightFraction( 6, 0.856289);
1625                                                  1595 
1626   AddMaterial("G4_RUBBER_NATURAL", 0.92, 0, 5    1596   AddMaterial("G4_RUBBER_NATURAL", 0.92, 0, 59.8, 2);
1627   AddElementByWeightFraction( 1, 0.118371);      1597   AddElementByWeightFraction( 1, 0.118371);
1628   AddElementByWeightFraction( 6, 0.881629);      1598   AddElementByWeightFraction( 6, 0.881629);
1629                                                  1599 
1630   AddMaterial("G4_RUBBER_NEOPRENE", 1.23, 0,     1600   AddMaterial("G4_RUBBER_NEOPRENE", 1.23, 0, 93., 3);
1631   AddElementByWeightFraction( 1, 0.05692 );      1601   AddElementByWeightFraction( 1, 0.05692 );
1632   AddElementByWeightFraction( 6, 0.542646);      1602   AddElementByWeightFraction( 6, 0.542646);
1633   AddElementByWeightFraction(17, 0.400434);      1603   AddElementByWeightFraction(17, 0.400434);
1634                                                  1604 
1635   AddMaterial("G4_SILICON_DIOXIDE", 2.32, 0,     1605   AddMaterial("G4_SILICON_DIOXIDE", 2.32, 0, 139.2, 2);
1636   AddElementByAtomCount("Si",  1);            << 1606   AddElementByWeightFraction( 8, 0.532565);
1637   AddElementByAtomCount("O" ,  2);            << 1607   AddElementByWeightFraction(14, 0.467435);
1638   chFormulas[nMaterials-1] = "SiO_2";            1608   chFormulas[nMaterials-1] = "SiO_2";
1639                                                  1609 
1640   AddMaterial("G4_SILVER_BROMIDE", 6.473, 0,     1610   AddMaterial("G4_SILVER_BROMIDE", 6.473, 0, 486.6, 2);
1641   AddElementByAtomCount("Ag",  1);            << 1611   AddElementByWeightFraction(35, 0.425537);
1642   AddElementByAtomCount("Br",  1);            << 1612   AddElementByWeightFraction(47, 0.574463);
1643                                                  1613 
1644   AddMaterial("G4_SILVER_CHLORIDE", 5.56, 0,     1614   AddMaterial("G4_SILVER_CHLORIDE", 5.56, 0, 398.4, 2);
1645   AddElementByAtomCount("Ag",  1);            << 1615   AddElementByWeightFraction(17, 0.247368);
1646   AddElementByAtomCount("Cl",  1);            << 1616   AddElementByWeightFraction(47, 0.752632);
1647                                                  1617 
1648   AddMaterial("G4_SILVER_HALIDES", 6.47, 0, 4    1618   AddMaterial("G4_SILVER_HALIDES", 6.47, 0, 487.1, 3);
1649   AddElementByWeightFraction(35, 0.422895);      1619   AddElementByWeightFraction(35, 0.422895);
1650   AddElementByWeightFraction(47, 0.573748);      1620   AddElementByWeightFraction(47, 0.573748);
1651   AddElementByWeightFraction(53, 0.003357);      1621   AddElementByWeightFraction(53, 0.003357);
1652                                                  1622 
1653   AddMaterial("G4_SILVER_IODIDE", 6.01, 0, 54    1623   AddMaterial("G4_SILVER_IODIDE", 6.01, 0, 543.5, 2);
1654   AddElementByAtomCount("Ag",  1);            << 1624   AddElementByWeightFraction(47, 0.459458);
1655   AddElementByAtomCount("I" ,  1);            << 1625   AddElementByWeightFraction(53, 0.540542);
1656                                                  1626 
1657   AddMaterial("G4_SKIN_ICRP", 1.09, 0, 72.7,     1627   AddMaterial("G4_SKIN_ICRP", 1.09, 0, 72.7, 9);
1658   AddElementByWeightFraction( 1, 0.100);         1628   AddElementByWeightFraction( 1, 0.100);
1659   AddElementByWeightFraction( 6, 0.204);         1629   AddElementByWeightFraction( 6, 0.204);
1660   AddElementByWeightFraction( 7, 0.042);         1630   AddElementByWeightFraction( 7, 0.042);
1661   AddElementByWeightFraction( 8, 0.645);         1631   AddElementByWeightFraction( 8, 0.645);
1662   AddElementByWeightFraction(11, 0.002);         1632   AddElementByWeightFraction(11, 0.002);
1663   AddElementByWeightFraction(15, 0.001);         1633   AddElementByWeightFraction(15, 0.001);
1664   AddElementByWeightFraction(16, 0.002);         1634   AddElementByWeightFraction(16, 0.002);
1665   AddElementByWeightFraction(17, 0.003);         1635   AddElementByWeightFraction(17, 0.003);
1666   AddElementByWeightFraction(19, 0.001);         1636   AddElementByWeightFraction(19, 0.001);
1667                                                  1637 
1668   AddMaterial("G4_SODIUM_CARBONATE", 2.532, 0    1638   AddMaterial("G4_SODIUM_CARBONATE", 2.532, 0, 125., 3);
1669   AddElementByAtomCount("Na",  2);            << 1639   AddElementByWeightFraction( 6, 0.113323);
1670   AddElementByAtomCount("C" ,  1);            << 1640   AddElementByWeightFraction( 8, 0.452861);
1671   AddElementByAtomCount("O" ,  3);            << 1641   AddElementByWeightFraction(11, 0.433815);
1672                                                  1642 
1673   AddMaterial("G4_SODIUM_IODIDE", 3.667, 0, 4    1643   AddMaterial("G4_SODIUM_IODIDE", 3.667, 0, 452., 2);
1674   AddElementByAtomCount("Na",  1);            << 1644   AddElementByWeightFraction(11, 0.153373);
1675   AddElementByAtomCount("I" ,  1);            << 1645   AddElementByWeightFraction(53, 0.846627);
1676                                                  1646 
1677   AddMaterial("G4_SODIUM_MONOXIDE", 2.27, 0,     1647   AddMaterial("G4_SODIUM_MONOXIDE", 2.27, 0, 148.8, 2);
1678   AddElementByAtomCount("Na",  2);            << 1648   AddElementByWeightFraction( 8, 0.258143);
1679   AddElementByAtomCount("O" ,  1);            << 1649   AddElementByWeightFraction(11, 0.741857);
1680                                                  1650 
1681   AddMaterial("G4_SODIUM_NITRATE", 2.261, 0,     1651   AddMaterial("G4_SODIUM_NITRATE", 2.261, 0, 114.6, 3);
1682   AddElementByAtomCount("Na",  1);            << 1652   AddElementByWeightFraction( 7, 0.164795);
1683   AddElementByAtomCount("N" ,  1);            << 1653   AddElementByWeightFraction( 8, 0.56472 );
1684   AddElementByAtomCount("O" ,  3);            << 1654   AddElementByWeightFraction(11, 0.270485);
1685                                                  1655 
1686   AddMaterial("G4_STILBENE", 0.9707, 0, 67.7,    1656   AddMaterial("G4_STILBENE", 0.9707, 0, 67.7, 2);
1687   AddElementByAtomCount("C" , 14);            << 1657   AddElementByWeightFraction( 1, 0.067101);
1688   AddElementByAtomCount("H" , 12);            << 1658   AddElementByWeightFraction( 6, 0.932899);
1689                                                  1659 
1690   AddMaterial("G4_SUCROSE", 1.5805, 0, 77.5,     1660   AddMaterial("G4_SUCROSE", 1.5805, 0, 77.5, 3);
1691   AddElementByAtomCount("C" , 12);            << 1661   AddElementByWeightFraction( 1, 0.064779);
1692   AddElementByAtomCount("H" , 22);            << 1662   AddElementByWeightFraction( 6, 0.42107);
1693   AddElementByAtomCount("O" , 11);            << 1663   AddElementByWeightFraction( 8, 0.514151);
1694                                               << 1664 
1695   AddMaterial("G4_TERPHENYL", 1.24 /*1.234*/, << 1665   AddMaterial("G4_TERPHENYL", 1.234, 0, 71.7, 2);
1696   //  AddElementByWeightFraction( 1, 0.044543 << 1666   AddElementByWeightFraction( 1, 0.044543);
1697   //  AddElementByWeightFraction( 6, 0.955457 << 1667   AddElementByWeightFraction( 6, 0.955457);
1698   //  These weight fractions correspond to C_ << 
1699   //  TERPHENYL is C_18H_14! The current dens << 
1700   //  M.Trocme & S.Seltzer                    << 
1701   AddElementByAtomCount("C" , 18);            << 
1702   AddElementByAtomCount("H" , 14);            << 
1703                                                  1668 
1704   AddMaterial("G4_TESTIS_ICRP", 1.04, 0, 75.,    1669   AddMaterial("G4_TESTIS_ICRP", 1.04, 0, 75., 9);
1705   AddElementByWeightFraction( 1, 0.106);         1670   AddElementByWeightFraction( 1, 0.106);
1706   AddElementByWeightFraction( 6, 0.099);         1671   AddElementByWeightFraction( 6, 0.099);
1707   AddElementByWeightFraction( 7, 0.020);         1672   AddElementByWeightFraction( 7, 0.020);
1708   AddElementByWeightFraction( 8, 0.766);         1673   AddElementByWeightFraction( 8, 0.766);
1709   AddElementByWeightFraction(11, 0.002);         1674   AddElementByWeightFraction(11, 0.002);
1710   AddElementByWeightFraction(15, 0.001);         1675   AddElementByWeightFraction(15, 0.001);
1711   AddElementByWeightFraction(16, 0.002);         1676   AddElementByWeightFraction(16, 0.002);
1712   AddElementByWeightFraction(17, 0.002);         1677   AddElementByWeightFraction(17, 0.002);
1713   AddElementByWeightFraction(19, 0.002);         1678   AddElementByWeightFraction(19, 0.002);
1714                                                  1679 
1715   AddMaterial("G4_TETRACHLOROETHYLENE", 1.625    1680   AddMaterial("G4_TETRACHLOROETHYLENE", 1.625, 0, 159.2, 2);
1716   AddElementByAtomCount("C" ,  2);            << 1681   AddElementByWeightFraction( 6, 0.144856);
1717   AddElementByAtomCount("Cl",  4);            << 1682   AddElementByWeightFraction(17, 0.855144);
1718                                                  1683 
1719   AddMaterial("G4_THALLIUM_CHLORIDE", 7.004,     1684   AddMaterial("G4_THALLIUM_CHLORIDE", 7.004, 0, 690.3, 2);
1720   AddElementByAtomCount("Tl",  1);            << 1685   AddElementByWeightFraction(17, 0.147822);
1721   AddElementByAtomCount("Cl",  1);            << 1686   AddElementByWeightFraction(81, 0.852178);
1722                                                  1687 
1723   // TISSUE_SOFT_MALE ICRU-44/46 (1989)          1688   // TISSUE_SOFT_MALE ICRU-44/46 (1989)
1724   AddMaterial("G4_TISSUE_SOFT_ICRP", 1.03, 0,    1689   AddMaterial("G4_TISSUE_SOFT_ICRP", 1.03, 0, 72.3, 9);
1725   AddElementByWeightFraction( 1, 0.105);         1690   AddElementByWeightFraction( 1, 0.105);
1726   AddElementByWeightFraction( 6, 0.256);         1691   AddElementByWeightFraction( 6, 0.256);
1727   AddElementByWeightFraction( 7, 0.027);         1692   AddElementByWeightFraction( 7, 0.027);
1728   AddElementByWeightFraction( 8, 0.602);         1693   AddElementByWeightFraction( 8, 0.602);
1729   AddElementByWeightFraction(11, 0.001);         1694   AddElementByWeightFraction(11, 0.001);
1730   AddElementByWeightFraction(15, 0.002);         1695   AddElementByWeightFraction(15, 0.002);
1731   AddElementByWeightFraction(16, 0.003);         1696   AddElementByWeightFraction(16, 0.003);
1732   AddElementByWeightFraction(17, 0.002);         1697   AddElementByWeightFraction(17, 0.002);
1733   AddElementByWeightFraction(19, 0.002);         1698   AddElementByWeightFraction(19, 0.002);
1734                                                  1699 
1735   // Tissue soft adult ICRU-33 (1980)            1700   // Tissue soft adult ICRU-33 (1980)
1736   AddMaterial("G4_TISSUE_SOFT_ICRU-4", 1.0, 0    1701   AddMaterial("G4_TISSUE_SOFT_ICRU-4", 1.0, 0, 74.9, 4);
1737   AddElementByWeightFraction( 1, 0.101);         1702   AddElementByWeightFraction( 1, 0.101);
1738   AddElementByWeightFraction( 6, 0.111);         1703   AddElementByWeightFraction( 6, 0.111);
1739   AddElementByWeightFraction( 7, 0.026);         1704   AddElementByWeightFraction( 7, 0.026);
1740   AddElementByWeightFraction( 8, 0.762);         1705   AddElementByWeightFraction( 8, 0.762);
1741                                                  1706 
1742   AddMaterial("G4_TISSUE-METHANE", 0.00106409    1707   AddMaterial("G4_TISSUE-METHANE", 0.00106409, 0, 61.2, 4, kStateGas);
1743   AddElementByWeightFraction( 1, 0.101869);      1708   AddElementByWeightFraction( 1, 0.101869);
1744   AddElementByWeightFraction( 6, 0.456179);      1709   AddElementByWeightFraction( 6, 0.456179);
1745   AddElementByWeightFraction( 7, 0.035172);      1710   AddElementByWeightFraction( 7, 0.035172);
1746   AddElementByWeightFraction( 8, 0.40678 );      1711   AddElementByWeightFraction( 8, 0.40678 );
1747                                                  1712 
1748   AddMaterial("G4_TISSUE-PROPANE", 0.00182628    1713   AddMaterial("G4_TISSUE-PROPANE", 0.00182628, 0, 59.5, 4, kStateGas);
1749   AddElementByWeightFraction( 1, 0.102672);      1714   AddElementByWeightFraction( 1, 0.102672);
1750   AddElementByWeightFraction( 6, 0.56894 );      1715   AddElementByWeightFraction( 6, 0.56894 );
1751   AddElementByWeightFraction( 7, 0.035022);      1716   AddElementByWeightFraction( 7, 0.035022);
1752   AddElementByWeightFraction( 8, 0.293366);      1717   AddElementByWeightFraction( 8, 0.293366);
1753                                                  1718 
1754   AddMaterial("G4_TITANIUM_DIOXIDE", 4.26, 0,    1719   AddMaterial("G4_TITANIUM_DIOXIDE", 4.26, 0, 179.5, 2);
1755   AddElementByAtomCount("Ti",  1);            << 1720   AddElementByWeightFraction( 8, 0.400592);
1756   AddElementByAtomCount("O" ,  2);            << 1721   AddElementByWeightFraction(22, 0.599408);
1757                                                  1722 
1758   AddMaterial("G4_TOLUENE", 0.8669, 0, 62.5,     1723   AddMaterial("G4_TOLUENE", 0.8669, 0, 62.5, 2);
1759   AddElementByAtomCount("C" ,  7);            << 1724   AddElementByWeightFraction( 1, 0.08751);
1760   AddElementByAtomCount("H" ,  8);            << 1725   AddElementByWeightFraction( 6, 0.91249);
1761                                                  1726 
1762   AddMaterial("G4_TRICHLOROETHYLENE", 1.46, 0    1727   AddMaterial("G4_TRICHLOROETHYLENE", 1.46, 0, 148.1, 3);
1763   AddElementByAtomCount("C" ,  2);            << 1728   AddElementByWeightFraction( 1, 0.007671);
1764   AddElementByAtomCount("H" ,  1);            << 1729   AddElementByWeightFraction( 6, 0.182831);
1765   AddElementByAtomCount("Cl",  3);            << 1730   AddElementByWeightFraction(17, 0.809498);
1766                                                  1731 
1767   AddMaterial("G4_TRIETHYL_PHOSPHATE", 1.07,     1732   AddMaterial("G4_TRIETHYL_PHOSPHATE", 1.07, 0, 81.2, 4);
1768   AddElementByAtomCount("C" ,  6);            << 1733   AddElementByWeightFraction( 1, 0.082998);
1769   AddElementByAtomCount("H" , 15);            << 1734   AddElementByWeightFraction( 6, 0.395628);
1770   AddElementByAtomCount("O" ,  4);            << 1735   AddElementByWeightFraction( 8, 0.351334);
1771   AddElementByAtomCount("P" ,  1);            << 1736   AddElementByWeightFraction(15, 0.17004 );
1772                                                  1737 
1773   AddMaterial("G4_TUNGSTEN_HEXAFLUORIDE", 2.4    1738   AddMaterial("G4_TUNGSTEN_HEXAFLUORIDE", 2.4, 0, 354.4, 2);
1774   AddElementByAtomCount("W" ,  1);            << 1739   AddElementByWeightFraction( 9, 0.382723);
1775   AddElementByAtomCount("F" ,  6);            << 1740   AddElementByWeightFraction(74, 0.617277);
1776                                                  1741 
1777   AddMaterial("G4_URANIUM_DICARBIDE", 11.28,     1742   AddMaterial("G4_URANIUM_DICARBIDE", 11.28, 0, 752., 2);
1778   AddElementByAtomCount("U" ,  1);            << 1743   AddElementByWeightFraction( 6, 0.091669);
1779   AddElementByAtomCount("C" ,  2);            << 1744   AddElementByWeightFraction(92, 0.908331);
1780                                                  1745 
1781   AddMaterial("G4_URANIUM_MONOCARBIDE", 13.63    1746   AddMaterial("G4_URANIUM_MONOCARBIDE", 13.63, 0, 862., 2);
1782   AddElementByAtomCount("U" ,  1);            << 1747   AddElementByWeightFraction( 6, 0.048036);
1783   AddElementByAtomCount("C" ,  1);            << 1748   AddElementByWeightFraction(92, 0.951964);
1784                                                  1749 
1785   AddMaterial("G4_URANIUM_OXIDE", 10.96, 0, 7    1750   AddMaterial("G4_URANIUM_OXIDE", 10.96, 0, 720.6, 2);
1786   AddElementByAtomCount("U" ,  1);            << 1751   AddElementByWeightFraction( 8, 0.118502);
1787   AddElementByAtomCount("O" ,  2);            << 1752   AddElementByWeightFraction(92, 0.881498);
1788                                                  1753 
1789   AddMaterial("G4_UREA", 1.323, 0, 72.8, 4);     1754   AddMaterial("G4_UREA", 1.323, 0, 72.8, 4);
1790   AddElementByAtomCount("C" ,  1);            << 1755   AddElementByWeightFraction( 1, 0.067131);
1791   AddElementByAtomCount("H" ,  4);            << 1756   AddElementByWeightFraction( 6, 0.199999);
1792   AddElementByAtomCount("N" ,  2);            << 1757   AddElementByWeightFraction( 7, 0.466459);
1793   AddElementByAtomCount("O" ,  1);            << 1758   AddElementByWeightFraction( 8, 0.266411);
1794                                                  1759 
1795   AddMaterial("G4_VALINE", 1.23, 0, 67.7, 4);    1760   AddMaterial("G4_VALINE", 1.23, 0, 67.7, 4);
1796   AddElementByAtomCount("C" ,  5);            << 1761   AddElementByWeightFraction( 1, 0.094641);
1797   AddElementByAtomCount("H" , 11);            << 1762   AddElementByWeightFraction( 6, 0.512645);
1798   AddElementByAtomCount("N" ,  1);            << 1763   AddElementByWeightFraction( 7, 0.119565);
1799   AddElementByAtomCount("O" ,  2);            << 1764   AddElementByWeightFraction( 8, 0.27315 );
1800                                                  1765 
1801   AddMaterial("G4_VITON", 1.8, 0, 98.6, 3);      1766   AddMaterial("G4_VITON", 1.8, 0, 98.6, 3);
1802   AddElementByWeightFraction( 1, 0.009417);      1767   AddElementByWeightFraction( 1, 0.009417);
1803   AddElementByWeightFraction( 6, 0.280555);      1768   AddElementByWeightFraction( 6, 0.280555);
1804   AddElementByWeightFraction( 9, 0.710028);      1769   AddElementByWeightFraction( 9, 0.710028);
1805                                                  1770 
                                                   >> 1771   AddMaterial("G4_WATER", 1.0,0, 78., 2);
                                                   >> 1772   AddElementByAtomCount("H", 2);
                                                   >> 1773   AddElementByAtomCount("O", 1);
                                                   >> 1774   chFormulas[nMaterials-1] = "H_2O";
                                                   >> 1775 
1806   AddMaterial("G4_WATER_VAPOR", 0.000756182,     1776   AddMaterial("G4_WATER_VAPOR", 0.000756182, 0, 71.6, 2, kStateGas);
1807   AddElementByAtomCount("H" ,  2);            << 1777   AddElementByAtomCount("H", 2);
1808   AddElementByAtomCount("O" ,  1);            << 1778   AddElementByAtomCount("O", 1);
1809   chFormulas[nMaterials-1] = "H_2O-Gas";         1779   chFormulas[nMaterials-1] = "H_2O-Gas";
1810                                                  1780 
1811   AddMaterial("G4_XYLENE", 0.87, 0, 61.8, 2);    1781   AddMaterial("G4_XYLENE", 0.87, 0, 61.8, 2);
1812   AddElementByAtomCount("C" ,  8);            << 1782   AddElementByWeightFraction( 1, 0.094935);
1813   AddElementByAtomCount("H" , 10);            << 1783   AddElementByWeightFraction( 6, 0.905065);
1814                                                  1784 
1815   AddMaterial("G4_GRAPHITE", 2.21, 6, 81.);   << 1785   AddMaterial("G4_GRAPHITE", 2.21, 6, 78.);
1816   chFormulas[nMaterials-1] = "Graphite";         1786   chFormulas[nMaterials-1] = "Graphite";
1817                                                  1787 
1818   nNIST = nMaterials;                            1788   nNIST = nMaterials;
1819 }                                                1789 }
1820                                                  1790 
1821 //....oooOO0OOooo........oooOO0OOooo........o    1791 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
1822                                                  1792 
1823 void G4NistMaterialBuilder::HepAndNuclearMate    1793 void G4NistMaterialBuilder::HepAndNuclearMaterials()
1824 {                                                1794 {
1825   AddMaterial("G4_lH2", 0.0708,  1,  21.8, 1,    1795   AddMaterial("G4_lH2", 0.0708,  1,  21.8, 1, kStateLiquid, false);
1826   AddMaterial("G4_lN2", 0.807,   7,  82.,  1,    1796   AddMaterial("G4_lN2", 0.807,   7,  82.,  1, kStateLiquid, false);
1827   AddMaterial("G4_lO2", 1.141,   8,  95.,  1,    1797   AddMaterial("G4_lO2", 1.141,   8,  95.,  1, kStateLiquid, false);
1828   AddMaterial("G4_lAr", 1.396 , 18, 188. , 1,    1798   AddMaterial("G4_lAr", 1.396 , 18, 188. , 1, kStateLiquid, false);
1829   AddMaterial("G4_lBr", 3.1028, 35, 343. , 1, << 
1830   AddMaterial("G4_lKr", 2.418 , 36, 352. , 1,    1799   AddMaterial("G4_lKr", 2.418 , 36, 352. , 1, kStateLiquid, false);
1831   AddMaterial("G4_lXe", 2.953 , 54, 482. , 1,    1800   AddMaterial("G4_lXe", 2.953 , 54, 482. , 1, kStateLiquid, false);
1832                                                  1801 
1833   AddMaterial("G4_PbWO4", 8.28, 0, 0.0, 3);      1802   AddMaterial("G4_PbWO4", 8.28, 0, 0.0, 3);
1834   AddElementByAtomCount("O" , 4);                1803   AddElementByAtomCount("O" , 4);
1835   AddElementByAtomCount("Pb", 1);                1804   AddElementByAtomCount("Pb", 1);
1836   AddElementByAtomCount("W" , 1);                1805   AddElementByAtomCount("W" , 1);
1837                                                  1806 
1838   G4double density = universe_mean_density*cm    1807   G4double density = universe_mean_density*cm3/g;
1839   AddMaterial("G4_Galactic", density, 1, 21.8    1808   AddMaterial("G4_Galactic", density, 1, 21.8, 1, kStateGas);
1840   AddGas("G4_Galactic",2.73*kelvin, 3.e-18*he << 1809   AddGas("G4_Galactic",2.73*kelvin, 3.e-18*pascal);
1841                                                  1810 
1842   AddMaterial("G4_GRAPHITE_POROUS", 1.7, 6, 8 << 1811   AddMaterial("G4_GRAPHITE_POROUS", 1.7, 6, 78.);
1843   chFormulas[nMaterials-1] = "Graphite";         1812   chFormulas[nMaterials-1] = "Graphite";
1844                                                  1813 
1845   // LUCITE is equal to plustiglass              1814   // LUCITE is equal to plustiglass
1846   AddMaterial("G4_LUCITE", 1.19, 0, 74., 3);     1815   AddMaterial("G4_LUCITE", 1.19, 0, 74., 3);
1847   AddElementByWeightFraction( 1, 0.080538);      1816   AddElementByWeightFraction( 1, 0.080538);
1848   AddElementByWeightFraction( 6, 0.599848);      1817   AddElementByWeightFraction( 6, 0.599848);
1849   AddElementByWeightFraction( 8, 0.319614);      1818   AddElementByWeightFraction( 8, 0.319614);
1850                                                  1819 
1851   // SRIM-2008 materials                         1820   // SRIM-2008 materials
1852   AddMaterial("G4_BRASS", 8.52, 0, 0.0, 3);      1821   AddMaterial("G4_BRASS", 8.52, 0, 0.0, 3);
1853   AddElementByWeightFraction( 29, 0.57515);   << 1822   AddElementByAtomCount("Cu", 62);
1854   AddElementByWeightFraction( 30, 0.33415);   << 1823   AddElementByAtomCount("Zn", 35);
1855   AddElementByWeightFraction( 82, 0.0907);    << 1824   AddElementByAtomCount("Pb" , 3);
1856                                                  1825 
1857   AddMaterial("G4_BRONZE", 8.82, 0, 0.0, 3);     1826   AddMaterial("G4_BRONZE", 8.82, 0, 0.0, 3);
1858   AddElementByWeightFraction( 29, 0.8494);    << 1827   AddElementByAtomCount("Cu", 89);
1859   AddElementByWeightFraction( 30, 0.0884);    << 1828   AddElementByAtomCount("Zn",  9);
1860   AddElementByWeightFraction( 82, 0.0622);    << 1829   AddElementByAtomCount("Pb" , 2);
1861                                                  1830 
1862   // parameters are taken from                   1831   // parameters are taken from
1863   //  http://www.azom.com/article.aspx?Articl    1832   //  http://www.azom.com/article.aspx?ArticleID=965
1864   AddMaterial("G4_STAINLESS-STEEL", 8.00, 0,     1833   AddMaterial("G4_STAINLESS-STEEL", 8.00, 0, 0.0, 3);
1865   AddElementByWeightFraction( 26, 0.7462);    << 1834   AddElementByAtomCount("Fe", 74);
1866   AddElementByWeightFraction( 24, 0.1690);    << 1835   AddElementByAtomCount("Cr", 18);
1867   AddElementByWeightFraction( 28, 0.0848);    << 1836   AddElementByAtomCount("Ni" , 8);
1868                                                  1837 
1869   AddMaterial("G4_CR39", 1.32, 0, 0.0, 3);       1838   AddMaterial("G4_CR39", 1.32, 0, 0.0, 3);
1870   AddElementByAtomCount("H", 18);                1839   AddElementByAtomCount("H", 18);
1871   AddElementByAtomCount("C", 12);                1840   AddElementByAtomCount("C", 12);
1872   AddElementByAtomCount("O", 7);                 1841   AddElementByAtomCount("O", 7);
1873                                                  1842 
1874   AddMaterial("G4_OCTADECANOL", 0.812, 0, 0.0    1843   AddMaterial("G4_OCTADECANOL", 0.812, 0, 0.0, 3);
1875   AddElementByAtomCount("H", 38);                1844   AddElementByAtomCount("H", 38);
1876   AddElementByAtomCount("C", 18);                1845   AddElementByAtomCount("C", 18);
1877   AddElementByAtomCount("O", 1);                 1846   AddElementByAtomCount("O", 1);
1878                                                  1847 
1879   nHEP = nMaterials;                             1848   nHEP = nMaterials;
1880 }                                                1849 }
1881                                                  1850 
1882 //....oooOO0OOooo........oooOO0OOooo........o    1851 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
1883                                                  1852 
1884 void G4NistMaterialBuilder::SpaceMaterials()     1853 void G4NistMaterialBuilder::SpaceMaterials()
1885 {                                                1854 {
1886   // density in g/cm3                            1855   // density in g/cm3
1887   AddMaterial("G4_KEVLAR" , 1.44, 0, 0.0, 4);    1856   AddMaterial("G4_KEVLAR" , 1.44, 0, 0.0, 4);
1888   AddElementByAtomCount("C", 14);                1857   AddElementByAtomCount("C", 14);
1889   AddElementByAtomCount("H", 10);                1858   AddElementByAtomCount("H", 10);
1890   AddElementByAtomCount("O", 2);                 1859   AddElementByAtomCount("O", 2);
1891   AddElementByAtomCount("N", 2);                 1860   AddElementByAtomCount("N", 2);
1892                                                  1861 
1893   AddMaterial("G4_DACRON" , 1.40, 0, 0.0, 3);    1862   AddMaterial("G4_DACRON" , 1.40, 0, 0.0, 3);   // G4_POLYETHYLENE_TEREPHTALATE
1894   AddElementByAtomCount("C", 10);                1863   AddElementByAtomCount("C", 10);
1895   AddElementByAtomCount("H", 8);                 1864   AddElementByAtomCount("H", 8);
1896   AddElementByAtomCount("O", 4);                 1865   AddElementByAtomCount("O", 4);
1897                                                  1866 
1898   AddMaterial("G4_NEOPRENE" , 1.23, 0, 0.0, 3    1867   AddMaterial("G4_NEOPRENE" , 1.23, 0, 0.0, 3);   // POLYCLOROPRENE
1899   AddElementByAtomCount("C", 4);                 1868   AddElementByAtomCount("C", 4);
1900   AddElementByAtomCount("H", 5);                 1869   AddElementByAtomCount("H", 5);
1901   AddElementByAtomCount("Cl",1);              << 1870   AddElementByAtomCount("Cl", 1);
1902                                                  1871 
1903   nSpace = nMaterials;                           1872   nSpace = nMaterials;
1904 }                                                1873 }
1905                                                  1874 
1906 //....oooOO0OOooo........oooOO0OOooo........o    1875 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
1907                                                  1876 
1908 /**                                           << 
1909   Build biochemical materials used in G4DNA A << 
1910   Materials are defined in bonded and unbonde << 
1911   following schema:                           << 
1912   G4_MATERIAL: Molecule in its free state     << 
1913   G4_DNA_MATERIAL: Molecule, considering atom << 
1914 */                                            << 
1915                                               << 
1916 void G4NistMaterialBuilder::BioChemicalMateri    1877 void G4NistMaterialBuilder::BioChemicalMaterials()
1917 {                                                1878 {
1918   // BEGIN MATERIALS IN THEIR UNBONDED FORM   << 1879   AddMaterial("G4_CYTOSINE", 1.55, 0, 72., 4);
1919                                               << 
1920   // G4_ADENINE, G4_GUANINE are defined in    << 
1921   // G4NistMaterialBuilder::NistCompoundMater << 
1922                                               << 
1923   //Tan, Z., et al. NIMB,2006(248)            << 
1924   AddMaterial("G4_CYTOSINE", 1.3, 0, 72., 4); << 
1925   AddElementByAtomCount("H", 5);                 1880   AddElementByAtomCount("H", 5);
1926   AddElementByAtomCount("C", 4);                 1881   AddElementByAtomCount("C", 4);
1927   AddElementByAtomCount("N", 3);                 1882   AddElementByAtomCount("N", 3);
1928   AddElementByAtomCount("O", 1);                 1883   AddElementByAtomCount("O", 1);
1929                                                  1884 
1930   AddMaterial("G4_THYMINE", 1.48, 0, 72., 4); << 1885   AddMaterial("G4_THYMINE", 1.23, 0, 72., 4);
1931   AddElementByAtomCount("H", 6);                 1886   AddElementByAtomCount("H", 6);
1932   AddElementByAtomCount("C", 5);                 1887   AddElementByAtomCount("C", 5);
1933   AddElementByAtomCount("N", 2);                 1888   AddElementByAtomCount("N", 2);
1934   AddElementByAtomCount("O", 2);                 1889   AddElementByAtomCount("O", 2);
1935                                                  1890 
1936   AddMaterial("G4_URACIL", 1.32, 0, 72., 4);     1891   AddMaterial("G4_URACIL", 1.32, 0, 72., 4);
1937   AddElementByAtomCount("H", 4);                 1892   AddElementByAtomCount("H", 4);
1938   AddElementByAtomCount("C", 4);                 1893   AddElementByAtomCount("C", 4);
1939   AddElementByAtomCount("N", 2);                 1894   AddElementByAtomCount("N", 2);
1940   AddElementByAtomCount("O", 2);                 1895   AddElementByAtomCount("O", 2);
1941                                                  1896 
1942   //ACD Labs Percepta Plateform - PhysChem Mo << 1897   // DNA_Nucleobase (Nucleobase-1H)
1943   //(https://www.acdlabs.com/products/percept << 
1944   AddMaterial("G4_DEOXYRIBOSE", 1.5, 0, 72, 3 << 
1945   AddElementByAtomCount("H", 10);             << 
1946   AddElementByAtomCount("C", 5);              << 
1947   AddElementByAtomCount("O", 4);              << 
1948                                               << 
1949   //Egan, E.P. and B.B. Luff,                 << 
1950   //Industrial & Engineering Chemistry, 1955. << 
1951   AddMaterial("G4_PHOSPHORIC_ACID", 1.87, 0,  << 
1952   AddElementByAtomCount("H", 3);              << 
1953   AddElementByAtomCount("P", 1);              << 
1954   AddElementByAtomCount("O", 4);              << 
1955                                               << 
1956   // END UNBONDED MATERIALS / BEGIN BONDED MA << 
1957                                               << 
1958   // Deoxyribose loses 3 OH groups in bonding << 
1959   AddMaterial("G4_DNA_DEOXYRIBOSE", 1, 0, 72. << 
1960   AddElementByAtomCount("H", 7);              << 
1961   AddElementByAtomCount("C", 5);              << 
1962   AddElementByAtomCount("O", 1);              << 
1963                                               << 
1964   // Typically there are no H atoms considere << 
1965   AddMaterial("G4_DNA_PHOSPHATE", 1, 0, 72.,  << 
1966   AddElementByAtomCount("P", 1);              << 
1967   AddElementByAtomCount("O", 4);              << 
1968                                               << 
1969   // GATCU bases bonded to a deoxyribose (the << 
1970   AddMaterial("G4_DNA_ADENINE", 1, 0, 72., 3)    1898   AddMaterial("G4_DNA_ADENINE", 1, 0, 72., 3);
1971   AddElementByAtomCount("H",4 );                 1899   AddElementByAtomCount("H",4 );
1972   AddElementByAtomCount("C",5 );                 1900   AddElementByAtomCount("C",5 );
1973   AddElementByAtomCount("N",5 );                 1901   AddElementByAtomCount("N",5 );
1974                                                  1902 
1975   AddMaterial("G4_DNA_GUANINE", 1, 0, 72., 4) << 1903   AddMaterial("G4_DNA_GUANINE", 1, 0, 72. ,4);
1976   AddElementByAtomCount("H",4 );                 1904   AddElementByAtomCount("H",4 );
1977   AddElementByAtomCount("C",5 );                 1905   AddElementByAtomCount("C",5 );
1978   AddElementByAtomCount("N",5 );                 1906   AddElementByAtomCount("N",5 );
1979   AddElementByAtomCount("O",1 );                 1907   AddElementByAtomCount("O",1 );
1980                                                  1908 
1981   AddMaterial("G4_DNA_CYTOSINE", 1, 0, 72., 4    1909   AddMaterial("G4_DNA_CYTOSINE", 1, 0, 72., 4);
1982   AddElementByAtomCount("H", 4);                 1910   AddElementByAtomCount("H", 4);
1983   AddElementByAtomCount("C", 4);                 1911   AddElementByAtomCount("C", 4);
1984   AddElementByAtomCount("N", 3);                 1912   AddElementByAtomCount("N", 3);
1985   AddElementByAtomCount("O", 1);                 1913   AddElementByAtomCount("O", 1);
1986                                                  1914 
1987   AddMaterial("G4_DNA_THYMINE", 1, 0, 72., 4)    1915   AddMaterial("G4_DNA_THYMINE", 1, 0, 72., 4);
1988   AddElementByAtomCount("H", 5);                 1916   AddElementByAtomCount("H", 5);
1989   AddElementByAtomCount("C", 5);                 1917   AddElementByAtomCount("C", 5);
1990   AddElementByAtomCount("N", 2);                 1918   AddElementByAtomCount("N", 2);
1991   AddElementByAtomCount("O", 2);                 1919   AddElementByAtomCount("O", 2);
1992                                                  1920 
1993   AddMaterial("G4_DNA_URACIL", 1, 0, 72., 4);    1921   AddMaterial("G4_DNA_URACIL", 1, 0, 72., 4);
1994   AddElementByAtomCount("H", 3);                 1922   AddElementByAtomCount("H", 3);
1995   AddElementByAtomCount("C", 4);                 1923   AddElementByAtomCount("C", 4);
1996   AddElementByAtomCount("N", 2);                 1924   AddElementByAtomCount("N", 2);
1997   AddElementByAtomCount("O", 2);                 1925   AddElementByAtomCount("O", 2);
1998                                                  1926 
1999   // END BONDED MATERIALS                     << 1927   // DNA_Nucleoside (Nucleoside-3H)
                                                   >> 1928   AddMaterial("G4_DNA_ADENOSINE", 1, 0, 72., 4);
                                                   >> 1929   AddElementByAtomCount("H", 10);
                                                   >> 1930   AddElementByAtomCount("C", 10);
                                                   >> 1931   AddElementByAtomCount("N", 5);
                                                   >> 1932   AddElementByAtomCount("O", 4);
                                                   >> 1933 
                                                   >> 1934   AddMaterial("G4_DNA_GUANOSINE", 1, 0, 72. ,4);
                                                   >> 1935   AddElementByAtomCount("H", 10);
                                                   >> 1936   AddElementByAtomCount("C", 10);
                                                   >> 1937   AddElementByAtomCount("N", 5);
                                                   >> 1938   AddElementByAtomCount("O", 5);
                                                   >> 1939 
                                                   >> 1940   AddMaterial("G4_DNA_CYTIDINE", 1, 0, 72., 4);
                                                   >> 1941   AddElementByAtomCount("H", 10);
                                                   >> 1942   AddElementByAtomCount("C", 9);
                                                   >> 1943   AddElementByAtomCount("N", 3);
                                                   >> 1944   AddElementByAtomCount("O", 5);
2000                                                  1945 
2001   // clang-format on                          << 1946   AddMaterial("G4_DNA_URIDINE", 1, 0, 72., 4);
                                                   >> 1947   AddElementByAtomCount("H", 9);
                                                   >> 1948   AddElementByAtomCount("C", 9);
                                                   >> 1949   AddElementByAtomCount("N", 2);
                                                   >> 1950   AddElementByAtomCount("O", 6);
                                                   >> 1951 
                                                   >> 1952   AddMaterial("G4_DNA_METHYLURIDINE", 1, 0, 72., 4);
                                                   >> 1953   AddElementByAtomCount("H", 11);
                                                   >> 1954   AddElementByAtomCount("C", 10);
                                                   >> 1955   AddElementByAtomCount("N", 2);
                                                   >> 1956   AddElementByAtomCount("O", 6);
                                                   >> 1957 
                                                   >> 1958   AddMaterial("G4_DNA_MONOPHOSPHATE", 1, 0, 72., 2);
                                                   >> 1959   AddElementByAtomCount("P", 1);
                                                   >> 1960   AddElementByAtomCount("O", 3);
                                                   >> 1961 
                                                   >> 1962   AddMaterial("G4_DNA_A", 1, 0, 72., 5);  //Adenine base
                                                   >> 1963   AddElementByAtomCount("H", 10);
                                                   >> 1964   AddElementByAtomCount("C", 10);
                                                   >> 1965   AddElementByAtomCount("N", 5);
                                                   >> 1966   AddElementByAtomCount("O", 7);
                                                   >> 1967   AddElementByAtomCount("P", 1);
                                                   >> 1968 
                                                   >> 1969   AddMaterial("G4_DNA_G", 1, 0, 72. ,5); //Guanine base 
                                                   >> 1970   AddElementByAtomCount("H", 10);
                                                   >> 1971   AddElementByAtomCount("C", 10);
                                                   >> 1972   AddElementByAtomCount("N", 5);
                                                   >> 1973   AddElementByAtomCount("O", 8);
                                                   >> 1974   AddElementByAtomCount("P", 1);
                                                   >> 1975 
                                                   >> 1976   AddMaterial("G4_DNA_C", 1, 0, 72., 5); // Cytosine base
                                                   >> 1977   AddElementByAtomCount("H", 10);
                                                   >> 1978   AddElementByAtomCount("C", 9);
                                                   >> 1979   AddElementByAtomCount("N", 3);
                                                   >> 1980   AddElementByAtomCount("O", 8);
                                                   >> 1981   AddElementByAtomCount("P", 1);
                                                   >> 1982 
                                                   >> 1983   AddMaterial("G4_DNA_U", 1, 0, 72., 5); // Uracil base
                                                   >> 1984   AddElementByAtomCount("H", 9);
                                                   >> 1985   AddElementByAtomCount("C", 9);
                                                   >> 1986   AddElementByAtomCount("N", 2);
                                                   >> 1987   AddElementByAtomCount("O", 9);
                                                   >> 1988   AddElementByAtomCount("P", 1);
                                                   >> 1989 
                                                   >> 1990   AddMaterial("G4_DNA_MU", 1, 0, 72., 5);  // MethaUracil base
                                                   >> 1991   AddElementByAtomCount("H", 11);
                                                   >> 1992   AddElementByAtomCount("C", 10);
                                                   >> 1993   AddElementByAtomCount("N", 2);
                                                   >> 1994   AddElementByAtomCount("O", 9);
                                                   >> 1995   AddElementByAtomCount("P", 1);
                                                   >> 1996   /*
                                                   >> 1997   // Complete 70 kg body of adult men from en.wikipedia.org/ see References there
                                                   >> 1998   AddMaterial("G4_BODY", 1.8, 0, 78, 12);
                                                   >> 1999   AddElementByWeightFraction( 8, 0.650);
                                                   >> 2000   AddElementByWeightFraction( 6, 0.180);
                                                   >> 2001   AddElementByWeightFraction( 1, 0.100);
                                                   >> 2002   AddElementByWeightFraction( 7, 0.030);
                                                   >> 2003   AddElementByWeightFraction(20, 0.015);
                                                   >> 2004   AddElementByWeightFraction(15, 0.010);
                                                   >> 2005   AddElementByWeightFraction(19, 0.0025);
                                                   >> 2006   AddElementByWeightFraction(16, 0.0025);
                                                   >> 2007   AddElementByWeightFraction(11, 0.0015);
                                                   >> 2008   AddElementByWeightFraction(17, 0.0015);
                                                   >> 2009   AddElementByWeightFraction(12, 0.0005);
                                                   >> 2010   AddElementByWeightFraction(26, 0.00006);
                                                   >> 2011   */
2002 }                                                2012 }
                                                   >> 2013 
                                                   >> 2014 
                                                   >> 2015 
                                                   >> 2016 
2003                                                  2017