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.4.p1)


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