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