| Geant4 Cross Reference |
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 // Code developed by:
27 // S.Guatelli, M. Large and A. Malaroda, University of Wollongong
28 //
29 //
30 #include "ICRP110PhantomMaterial_Female.hh"
31 #include "globals.hh"
32 #include "G4SystemOfUnits.hh"
33 #include "G4MaterialPropertiesTable.hh"
34 #include "G4MaterialPropertyVector.hh"
35 #include "G4MaterialTable.hh"
36 #include "Randomize.hh"
37 #include "G4RunManager.hh"
38 #include "G4Element.hh"
39 #include "G4ElementTable.hh"
40
41 ICRP110PhantomMaterial_Female::ICRP110PhantomMaterial_Female():
42 fLung(nullptr), fTeeth(nullptr), fBone(nullptr), fHumeri_upper(nullptr), fHumeri_lower(nullptr),
43 fArm_lower(nullptr), fHand(nullptr), fClavicle(nullptr), fCranium(nullptr), fFemora_upper(nullptr),
44 fFemora_lower(nullptr), fLeg_lower(nullptr), fFoot(nullptr), fMandible(nullptr),
45 fPelvis(nullptr), fRibs(nullptr), fScapulae(nullptr), fSpine_cervical(nullptr),
46 fSpine_lumbar(nullptr), fSpine_thoratic(nullptr), fSacrum(nullptr),
47 fSternum(nullptr), fHf_upper(nullptr), fHf_lower(nullptr), fMed_lowerleg(nullptr),
48 fMed_lowerarm(nullptr), fCartilage(nullptr), fSkin(nullptr), fBlood(nullptr),
49 fMuscle(nullptr), fLiver(nullptr), fPancreas(nullptr), fBrain(nullptr), fHeart(nullptr), fEye(nullptr),
50 fKidney(nullptr), fStomach(nullptr), fIntestine_sml(nullptr), fIntestine_lrg(nullptr),
51 fSpleen(nullptr), fThyroid(nullptr), fBladder(nullptr), fOvaries_testes(nullptr), fAdrenals(nullptr),
52 fOesophagus(nullptr), fMisc(nullptr), fUterus_prostate(nullptr), fLymph(nullptr),
53 fBreast_glandular(nullptr), fBreast_adipose(nullptr), fGastro_content(nullptr),
54 fUrine(nullptr)
55 {;}
56
57 ICRP110PhantomMaterial_Female::~ICRP110PhantomMaterial_Female()
58 {;}
59
60 void ICRP110PhantomMaterial_Female::DefineMaterials()
61 {
62 // Define required materials
63
64 G4double A; // atomic mass
65 G4double Z; // atomic number
66 G4double d; // density
67
68 // General elements
69
70 A = 1.01*g/mole;
71 auto elH = new G4Element ("Hydrogen","H",Z = 1.,A);
72
73 A = 12.011*g/mole;
74 auto elC = new G4Element("Carbon","C",Z = 6.,A);
75
76 A = 14.01*g/mole;
77 auto elN = new G4Element("Nitrogen","N",Z = 7.,A);
78
79 A = 16.00*g/mole;
80 auto elO = new G4Element("Oxygen","O",Z = 8.,A);
81
82 A = 22.99*g/mole;
83 auto elNa = new G4Element("Sodium","Na",Z = 11.,A);
84
85 A = 24.305*g/mole;
86 auto elMg = new G4Element("Magnesium","Mg",Z = 12.,A);
87
88 A = 30.974*g/mole;
89 auto elP = new G4Element("Phosphorus","P",Z = 15.,A);
90
91 A = 32.064*g/mole;
92 auto elS = new G4Element("Sulfur","S",Z = 16.,A);
93
94 A = 35.453*g/mole;
95 auto elCl = new G4Element("Chlorine","Cl",Z = 17.,A);
96
97 A = 39.098*g/mole;
98 auto elK = new G4Element("Potassium","K",Z = 19.,A);
99
100 A = 40.08*g/mole;
101 auto elCa = new G4Element("Calcium","Ca",Z = 20.,A);
102
103 A = 55.85*g/mole;
104 auto elFe = new G4Element("Iron","Fe",Z = 26.,A);
105
106 A = 126.90447 *g/mole;
107 auto elI = new G4Element("Iodine","I", Z = 53.,A);
108
109 //Added tissues for phantom including their tissue ID
110 //Teeth -> Tissue ID 1
111 d = 2.750 *g/cm3;
112 fTeeth = new G4Material("teeth",d,7);
113 fTeeth -> AddElement(elH,0.022);
114 fTeeth -> AddElement(elC,0.095);
115 fTeeth -> AddElement(elN,0.029);
116 fTeeth -> AddElement(elO,0.421);
117 fTeeth -> AddElement(elMg,0.007);
118 fTeeth -> AddElement(elP,0.137);
119 fTeeth -> AddElement(elCa,0.289);
120
121 //Mineral Bone -> Tissue ID 2
122 d = 1.920 *g/cm3;
123 fBone = new G4Material("bone",d,9);
124 fBone -> AddElement(elH,0.036);
125 fBone -> AddElement(elC,0.159);
126 fBone -> AddElement(elN,0.042);
127 fBone -> AddElement(elO,0.448);
128 fBone -> AddElement(elNa,0.003);
129 fBone -> AddElement(elMg,0.002);
130 fBone -> AddElement(elP,0.094);
131 fBone -> AddElement(elS,0.003);
132 fBone -> AddElement(elCa,0.213);
133
134 //Humeri, upper half, spongiosa -> ID 3
135 d = 1.185 *g/cm3;
136 fHumeri_upper = new G4Material("humeri_upper",d,11);
137 fHumeri_upper -> AddElement(elH,0.087) ;
138 fHumeri_upper -> AddElement(elC,0.366);
139 fHumeri_upper -> AddElement(elN,0.025);
140 fHumeri_upper -> AddElement(elO,0.422);
141 fHumeri_upper -> AddElement(elNa,0.002);
142 fHumeri_upper -> AddElement(elMg,0.001);
143 fHumeri_upper -> AddElement(elP,0.030);
144 fHumeri_upper -> AddElement(elS,0.003);
145 fHumeri_upper -> AddElement(elCl,0.001);
146 fHumeri_upper -> AddElement(elK, 0.001);
147 fHumeri_upper -> AddElement(elCa,0.062);
148
149 //Humeri, lower half, spongiosa -> ID 4
150 d = 1.117 *g/cm3;
151 fHumeri_lower = new G4Material("humeri_lower",d,9);
152 fHumeri_lower -> AddElement(elH,0.096);
153 fHumeri_lower -> AddElement(elC,0.473);
154 fHumeri_lower -> AddElement(elN,0.017);
155 fHumeri_lower -> AddElement(elO,0.341);
156 fHumeri_lower -> AddElement(elNa,0.002);
157 fHumeri_lower -> AddElement(elP,0.022);
158 fHumeri_lower -> AddElement(elS,0.002);
159 fHumeri_lower -> AddElement(elCl,0.001);
160 fHumeri_lower -> AddElement(elCa,0.046);
161
162 //Lower arm bones, spongiosa -> ID 5
163 d = 1.117 *g/cm3;
164 fArm_lower = new G4Material("arm_lower",d,9);
165 fArm_lower -> AddElement(elH,0.096);
166 fArm_lower -> AddElement(elC,0.473);
167 fArm_lower -> AddElement(elN,0.017);
168 fArm_lower -> AddElement(elO,0.341);
169 fArm_lower -> AddElement(elNa,0.002);
170 fArm_lower -> AddElement(elP,0.022);
171 fArm_lower -> AddElement(elS,0.002);
172 fArm_lower -> AddElement(elCl,0.001);
173 fArm_lower -> AddElement(elCa,0.046);
174
175 //Hand Bones, Spongiosa ->ID 6
176 d = 1.117 *g/cm3;
177 fHand = new G4Material("hand",d,9);
178 fHand -> AddElement(elH,0.096);
179 fHand -> AddElement(elC,0.473);
180 fHand -> AddElement(elN,0.017);
181 fHand -> AddElement(elO,0.341);
182 fHand -> AddElement(elNa,0.002);
183 fHand -> AddElement(elP,0.022);
184 fHand -> AddElement(elS,0.002);
185 fHand -> AddElement(elCl,0.001);
186 fHand -> AddElement(elCa,0.046);
187
188 //Clavicles, spongiosa -> ID 7
189 d = 1.191 *g/cm3;
190 fClavicle = new G4Material("clavicle",d,11);
191 fClavicle -> AddElement(elH,0.087);
192 fClavicle -> AddElement(elC,0.361);
193 fClavicle -> AddElement(elN,0.025);
194 fClavicle -> AddElement(elO,0.424);
195 fClavicle -> AddElement(elNa,0.002);
196 fClavicle -> AddElement(elMg,0.001);
197 fClavicle -> AddElement(elP,0.031);
198 fClavicle -> AddElement(elS,0.003);
199 fClavicle -> AddElement(elCl,0.001);
200 fClavicle -> AddElement(elK,0.001);
201 fClavicle -> AddElement(elCa,0.064);
202
203 // Cranium, spongiosa -> ID 8
204 d = 1.245 *g/cm3;
205 fCranium = new G4Material("cranium",d,11);
206 fCranium -> AddElement(elH,0.081);
207 fCranium -> AddElement(elC,0.317);
208 fCranium -> AddElement(elN,0.028);
209 fCranium -> AddElement(elO,0.451);
210 fCranium -> AddElement(elNa,0.002);
211 fCranium -> AddElement(elMg, 0.001);
212 fCranium -> AddElement(elP,0.037);
213 fCranium -> AddElement(elS,0.003);
214 fCranium -> AddElement(elCl,0.001);
215 fCranium -> AddElement(elK,0.001);
216 fCranium -> AddElement(elCa,0.078);
217
218 //Femora, upper half, spongiosa -> ID 9
219 d = 1.046 *g/cm3;
220 fFemora_upper = new G4Material("femora_upper",d,10);
221 fFemora_upper -> AddElement(elH,0.104);
222 fFemora_upper -> AddElement(elC,0.496);
223 fFemora_upper -> AddElement(elN,0.018);
224 fFemora_upper -> AddElement(elO,0.349);
225 fFemora_upper -> AddElement(elNa,0.001);
226 fFemora_upper -> AddElement(elP,0.009);
227 fFemora_upper -> AddElement(elS,0.002);
228 fFemora_upper -> AddElement(elCl,0.001);
229 fFemora_upper -> AddElement(elK,0.001);
230 fFemora_upper -> AddElement(elCa,0.019);
231
232 //Femora, lower half, spongiosa -> ID 10
233 d = 1.117 *g/cm3;
234 fFemora_lower = new G4Material("femora_lower",d,9);
235 fFemora_lower -> AddElement(elH,0.096);
236 fFemora_lower -> AddElement(elC,0.473);
237 fFemora_lower -> AddElement(elN,0.017);
238 fFemora_lower -> AddElement(elO,0.341);
239 fFemora_lower -> AddElement(elNa,0.002);
240 fFemora_lower -> AddElement(elP,0.022);
241 fFemora_lower -> AddElement(elS,0.002);
242 fFemora_lower -> AddElement(elCl,0.001);
243 fFemora_lower -> AddElement(elCa,0.046);
244
245 //Lower leg bones, spongiosa -> ID 11
246 d = 1.117 *g/cm3;
247 fLeg_lower = new G4Material("leg_lower",d,9);
248 fLeg_lower -> AddElement(elH,0.096);
249 fLeg_lower -> AddElement(elC,0.473);
250 fLeg_lower -> AddElement(elN,0.017);
251 fLeg_lower -> AddElement(elO,0.341);
252 fLeg_lower -> AddElement(elNa,0.002);
253 fLeg_lower -> AddElement(elP,0.022);
254 fLeg_lower -> AddElement(elS,0.002);
255 fLeg_lower -> AddElement(elCl,0.001);
256 fLeg_lower -> AddElement(elCa,0.046);
257
258 //Foot bones, spongiosa ->ID 12
259 d = 1.117 *g/cm3;
260 fFoot = new G4Material("foot",d,9);
261 fFoot -> AddElement(elH,0.096);
262 fFoot -> AddElement(elC,0.473);
263 fFoot -> AddElement(elN,0.017);
264 fFoot -> AddElement(elO,0.341);
265 fFoot -> AddElement(elNa,0.002);
266 fFoot -> AddElement(elP,0.022);
267 fFoot -> AddElement(elS,0.002);
268 fFoot -> AddElement(elCl,0.001);
269 fFoot -> AddElement(elCa,0.046);
270
271 //Mandible, spongiosa -> ID 13
272 d = 1.189 *g/cm3;
273 fMandible = new G4Material("mandible",d,11);
274 fMandible -> AddElement(elH,0.087);
275 fMandible -> AddElement(elC,0.357);
276 fMandible -> AddElement(elN,0.026);
277 fMandible -> AddElement(elO,0.429);
278 fMandible -> AddElement(elNa,0.002);
279 fMandible -> AddElement(elMg,0.001);
280 fMandible -> AddElement(elP,0.030);
281 fMandible -> AddElement(elS,0.003);
282 fMandible -> AddElement(elCl,0.001);
283 fMandible -> AddElement(elK,0.001);
284 fMandible -> AddElement(elCa,0.063);
285
286 //Pelvis, Spongiosa -> ID 14
287 d = 1.109 *g/cm3;
288 fPelvis = new G4Material("pelvis",d,10);
289 fPelvis -> AddElement(elH,0.096);
290 fPelvis -> AddElement(elC,0.406);
291 fPelvis -> AddElement(elN,0.025);
292 fPelvis -> AddElement(elO,0.412);
293 fPelvis -> AddElement(elNa,0.001);
294 fPelvis -> AddElement(elP,0.018);
295 fPelvis -> AddElement(elS,0.002);
296 fPelvis -> AddElement(elCl,0.001);
297 fPelvis -> AddElement(elK,0.001);
298 fPelvis -> AddElement(elCa,0.038);
299
300 //Ribs, spongiosa -> ID 15
301 d = 1.092 *g/cm3;
302 fRibs = new G4Material("ribs",d,11);
303 fRibs -> AddElement(elH,0.097);
304 fRibs -> AddElement(elC,0.381);
305 fRibs -> AddElement(elN,0.028);
306 fRibs -> AddElement(elO,0.445);
307 fRibs -> AddElement(elNa,0.001);
308 fRibs -> AddElement(elP,0.014);
309 fRibs -> AddElement(elS,0.002);
310 fRibs -> AddElement(elCl,0.002);
311 fRibs -> AddElement(elK,0.001);
312 fRibs -> AddElement(elCa,0.028);
313 fRibs -> AddElement(elFe,0.001);
314
315 //Scapulae, spongiosa -> ID 16
316 d = 1.128 *g/cm3;
317 fScapulae = new G4Material("scapulae",d,10);
318 fScapulae -> AddElement(elH,0.094);
319 fScapulae -> AddElement(elC,0.406);
320 fScapulae -> AddElement(elN,0.024);
321 fScapulae -> AddElement(elO,0.404);
322 fScapulae -> AddElement(elNa,0.001);
323 fScapulae -> AddElement(elP,0.022);
324 fScapulae -> AddElement(elS,0.002);
325 fScapulae -> AddElement(elCl,0.001);
326 fScapulae -> AddElement(elK,0.001);
327 fScapulae -> AddElement(elCa,0.045);
328
329 //Cervical Spine, spongiosa -> ID 17
330 d = 1.135 *g/cm3;
331 fSpine_cervical = new G4Material("spine_cervical",d,10);
332 fSpine_cervical -> AddElement(elH,0.092);
333 fSpine_cervical -> AddElement(elC,0.351);
334 fSpine_cervical -> AddElement(elN,0.029);
335 fSpine_cervical -> AddElement(elO,0.458);
336 fSpine_cervical -> AddElement(elNa,0.001);
337 fSpine_cervical -> AddElement(elP,0.021);
338 fSpine_cervical -> AddElement(elS,0.002);
339 fSpine_cervical -> AddElement(elCl,0.002);
340 fSpine_cervical -> AddElement(elK,0.002);
341 fSpine_cervical -> AddElement(elCa,0.043);
342
343 //Thoratic Spine, spongiosa -> ID 18
344 d = 1.084 *g/cm3;
345 fSpine_thoratic = new G4Material("spine_thoratic",d,11);
346 fSpine_thoratic -> AddElement(elH,0.098);
347 fSpine_thoratic -> AddElement(elC,0.386);
348 fSpine_thoratic -> AddElement(elN,0.028);
349 fSpine_thoratic -> AddElement(elO,0.442);
350 fSpine_thoratic -> AddElement(elNa,0.001);
351 fSpine_thoratic -> AddElement(elP,0.013);
352 fSpine_thoratic -> AddElement(elS,0.002);
353 fSpine_thoratic -> AddElement(elCl,0.002);
354 fSpine_thoratic -> AddElement(elK,0.001);
355 fSpine_thoratic -> AddElement(elCa,0.026);
356 fSpine_thoratic -> AddElement(elFe,0.001);
357
358 //Lumbar Spine, spongiosa -> ID 19
359 d = 1.171 *g/cm3;
360 fSpine_lumbar = new G4Material("spine_lumbar",d,11);
361 fSpine_lumbar -> AddElement(elH,0.088);
362 fSpine_lumbar -> AddElement(elC,0.329);
363 fSpine_lumbar -> AddElement(elN,0.030);
364 fSpine_lumbar -> AddElement(elO,0.466);
365 fSpine_lumbar -> AddElement(elNa,0.001);
366 fSpine_lumbar -> AddElement(elMg,0.001);
367 fSpine_lumbar -> AddElement(elP,0.026);
368 fSpine_lumbar -> AddElement(elS,0.003);
369 fSpine_lumbar -> AddElement(elCl,0.001);
370 fSpine_lumbar -> AddElement(elK,0.001);
371 fSpine_lumbar -> AddElement(elCa,0.054);
372
373 //Sacrum, spongiosa -> ID 20
374 d = 1.052 *g/cm3;
375 fSacrum = new G4Material("sacrum",d,11);
376 fSacrum -> AddElement(elH,0.102);
377 fSacrum -> AddElement(elC,0.410);
378 fSacrum -> AddElement(elN,0.027);
379 fSacrum -> AddElement(elO,0.433);
380 fSacrum -> AddElement(elNa,0.001);
381 fSacrum -> AddElement(elP,0.007);
382 fSacrum -> AddElement(elS,0.002);
383 fSacrum -> AddElement(elCl,0.002);
384 fSacrum -> AddElement(elK,0.001);
385 fSacrum -> AddElement(elCa,0.014);
386 fSacrum -> AddElement(elFe,0.001);
387
388 //Sternum, spongiosa -> ID 21
389 d = 1.076 *g/cm3;
390 fSternum = new G4Material("sternum",d,11);
391 fSternum -> AddElement(elH,0.099);
392 fSternum -> AddElement(elC,0.392);
393 fSternum -> AddElement(elN,0.028);
394 fSternum -> AddElement(elO,0.439);
395 fSternum -> AddElement(elNa,0.001);
396 fSternum -> AddElement(elP,0.012);
397 fSternum -> AddElement(elS,0.002);
398 fSternum -> AddElement(elCl,0.002);
399 fSternum -> AddElement(elK,0.001);
400 fSternum -> AddElement(elCa,0.023);
401 fSternum -> AddElement(elFe,0.001);
402
403 //Humeri and femora, upper halves, medullary cavity -> ID 22
404 d = 0.980 *g/cm3;
405 fHf_upper = new G4Material("hf_upper",d,7);
406 fHf_upper -> AddElement(elH,0.115);
407 fHf_upper -> AddElement(elC,0.637);
408 fHf_upper -> AddElement(elN,0.007);
409 fHf_upper -> AddElement(elO,0.238);
410 fHf_upper -> AddElement(elNa,0.001);
411 fHf_upper -> AddElement(elS,0.001);
412 fHf_upper -> AddElement(elCl,0.001);
413
414 //Humeri and femora, lower halves, medullary cavity -> ID 23
415 d = 0.980 *g/cm3;
416 fHf_lower = new G4Material("hf_lower",d,7);
417 fHf_lower -> AddElement(elH,0.115);
418 fHf_lower -> AddElement(elC,0.637);
419 fHf_lower -> AddElement(elN,0.007);
420 fHf_lower -> AddElement(elO,0.238);
421 fHf_lower -> AddElement(elNa,0.001);
422 fHf_lower -> AddElement(elS,0.001);
423 fHf_lower -> AddElement(elCl,0.001);
424
425 //Lower arm bones, medullary cavity -> ID 24
426 d = 0.980 *g/cm3;
427 fMed_lowerarm = new G4Material("med_lowerarm",d,7);
428 fMed_lowerarm -> AddElement(elH,0.115);
429 fMed_lowerarm -> AddElement(elC,0.637);
430 fMed_lowerarm -> AddElement(elN,0.007);
431 fMed_lowerarm -> AddElement(elO,0.238);
432 fMed_lowerarm -> AddElement(elNa,0.001);
433 fMed_lowerarm -> AddElement(elS,0.001);
434 fMed_lowerarm -> AddElement(elCl,0.001);
435
436 //Lower leg bones, medullary cavity -> ID 25
437 d = 0.980 *g/cm3;
438 fMed_lowerleg = new G4Material("med_lowerleg",d,7);
439 fMed_lowerleg -> AddElement(elH,0.115);
440 fMed_lowerleg -> AddElement(elC,0.637);
441 fMed_lowerleg -> AddElement(elN,0.007);
442 fMed_lowerleg -> AddElement(elO,0.238);
443 fMed_lowerleg -> AddElement(elNa,0.001);
444 fMed_lowerleg -> AddElement(elS,0.001);
445 fMed_lowerleg -> AddElement(elCl,0.001);
446
447 //Cartilage -> ID 26
448 d = 1.100 *g/cm3;
449 fCartilage = new G4Material("cartilage",d,8);
450 fCartilage -> AddElement(elH,0.096);
451 fCartilage -> AddElement(elC,0.099);
452 fCartilage -> AddElement(elN,0.022);
453 fCartilage -> AddElement(elO,0.744);
454 fCartilage -> AddElement(elNa,0.005);
455 fCartilage -> AddElement(elP,0.022);
456 fCartilage -> AddElement(elS,0.009);
457 fCartilage -> AddElement(elCl,0.003);
458
459 //Skin -> Id 27
460 d = 1.090 *g/cm3;
461 fSkin = new G4Material("skin",d,9);
462 fSkin -> AddElement(elH,0.100);
463 fSkin -> AddElement(elC,0.199);
464 fSkin -> AddElement(elN,0.042);
465 fSkin -> AddElement(elO,0.650);
466 fSkin -> AddElement(elNa,0.002);
467 fSkin -> AddElement(elP,0.001);
468 fSkin -> AddElement(elS,0.002);
469 fSkin -> AddElement(elCl,0.003);
470 fSkin -> AddElement(elK,0.001);
471
472 //Blood -> ID 28
473 d = 1.060 *g/cm3;
474 fBlood = new G4Material("blood",d,10);
475 fBlood -> AddElement(elH,0.102);
476 fBlood -> AddElement(elC,0.110);
477 fBlood -> AddElement(elN,0.033);
478 fBlood -> AddElement(elO,0.745);
479 fBlood -> AddElement(elNa,0.001);
480 fBlood -> AddElement(elP,0.001);
481 fBlood -> AddElement(elS,0.002);
482 fBlood -> AddElement(elCl,0.003);
483 fBlood -> AddElement(elK,0.002);
484 fBlood -> AddElement(elFe,0.001);
485
486 //Muscular Tissue -> ID 29
487 d = 1.050 *g/cm3;
488 fMuscle = new G4Material("muscle",d,9);
489 fMuscle -> AddElement(elH,0.102);
490 fMuscle -> AddElement(elC,0.142);
491 fMuscle -> AddElement(elN,0.034);
492 fMuscle -> AddElement(elO,0.711);
493 fMuscle -> AddElement(elNa,0.001);
494 fMuscle -> AddElement(elP,0.002);
495 fMuscle -> AddElement(elS,0.003);
496 fMuscle -> AddElement(elCl,0.001);
497 fMuscle -> AddElement(elK,0.004);
498
499 //Liver -> ID 30
500 d = 1.050 *g/cm3;
501 fLiver = new G4Material("liver",d,9);
502 fLiver -> AddElement(elH,0.102);
503 fLiver -> AddElement(elC,0.131);
504 fLiver -> AddElement(elN,0.031);
505 fLiver -> AddElement(elO,0.724);
506 fLiver -> AddElement(elNa,0.002);
507 fLiver -> AddElement(elP,0.002);
508 fLiver -> AddElement(elS,0.003);
509 fLiver -> AddElement(elCl,0.002);
510 fLiver -> AddElement(elK,0.003);
511
512 //Pancreas ->ID 31
513 d = 1.050 *g/cm3;
514 fPancreas = new G4Material("pancreas",d,9);
515 fPancreas -> AddElement(elH,0.105);
516 fPancreas -> AddElement(elC,0.157);
517 fPancreas -> AddElement(elN,0.024);
518 fPancreas -> AddElement(elO,0.705);
519 fPancreas -> AddElement(elNa,0.002);
520 fPancreas -> AddElement(elP,0.002);
521 fPancreas -> AddElement(elS,0.001);
522 fPancreas -> AddElement(elCl,0.002);
523 fPancreas -> AddElement(elK,0.002);
524
525 //Brain -> ID 32
526 d = 1.050 *g/cm3;
527 fBrain = new G4Material("brain",d,9);
528 fBrain -> AddElement(elH,0.107);
529 fBrain -> AddElement(elC,0.144);
530 fBrain -> AddElement(elN,0.022);
531 fBrain -> AddElement(elO,0.713);
532 fBrain -> AddElement(elNa,0.002);
533 fBrain -> AddElement(elP,0.004);
534 fBrain -> AddElement(elS,0.002);
535 fBrain -> AddElement(elCl,0.003);
536 fBrain -> AddElement(elK,0.003);
537
538 //Heart -> ID 33
539 d = 1.050 *g/cm3;
540 fHeart = new G4Material("heart",d,9);
541 fHeart -> AddElement(elH,0.104);
542 fHeart -> AddElement(elC,0.138);
543 fHeart -> AddElement(elN,0.029);
544 fHeart -> AddElement(elO,0.719);
545 fHeart -> AddElement(elNa,0.001);
546 fHeart -> AddElement(elP,0.002);
547 fHeart -> AddElement(elS,0.002);
548 fHeart -> AddElement(elCl,0.002);
549 fHeart -> AddElement(elK,0.003);
550
551 //Eye ->ID 34
552 d = 1.050 *g/cm3;
553 fEye = new G4Material("eye",d,8);
554 fEye -> AddElement(elH,0.097);
555 fEye -> AddElement(elC,0.183);
556 fEye -> AddElement(elN,0.054);
557 fEye -> AddElement(elO,0.660);
558 fEye -> AddElement(elCa,0.001);
559 fEye -> AddElement(elP,0.001);
560 fEye -> AddElement(elS,0.003);
561 fEye -> AddElement(elCl,0.001);
562
563 //Kidneys -> ID 35
564 d = 1.050 *g/cm3;
565 fKidney = new G4Material("kidney",d,10);
566 fKidney -> AddElement(elH,0.103);
567 fKidney -> AddElement(elC,0.125);
568 fKidney -> AddElement(elN,0.031);
569 fKidney -> AddElement(elO,0.730);
570 fKidney -> AddElement(elNa,0.002);
571 fKidney -> AddElement(elP,0.002);
572 fKidney -> AddElement(elS,0.002);
573 fKidney -> AddElement(elCl,0.002);
574 fKidney -> AddElement(elK,0.002);
575 fKidney -> AddElement(elCa,0.001);
576
577 //Stomach ->ID 36
578 d = 1.040 *g/cm3;
579 fStomach = new G4Material("stomach",d,9);
580 fStomach -> AddElement(elH,0.105);
581 fStomach -> AddElement(elC,0.114);
582 fStomach -> AddElement(elN,0.025);
583 fStomach -> AddElement(elO,0.750);
584 fStomach -> AddElement(elNa,0.001);
585 fStomach -> AddElement(elP,0.001);
586 fStomach -> AddElement(elS,0.001);
587 fStomach -> AddElement(elCl,0.002);
588 fStomach -> AddElement(elK,0.001);
589
590 //Small intestine ->ID 37
591 d = 1.040 *g/cm3;
592 fIntestine_sml = new G4Material("intestine_sml",d,9);
593 fIntestine_sml -> AddElement(elH,0.105);
594 fIntestine_sml -> AddElement(elC,0.114);
595 fIntestine_sml -> AddElement(elN,0.025);
596 fIntestine_sml -> AddElement(elO,0.750);
597 fIntestine_sml -> AddElement(elNa,0.001);
598 fIntestine_sml -> AddElement(elP,0.001);
599 fIntestine_sml -> AddElement(elS,0.001);
600 fIntestine_sml -> AddElement(elCl,0.002);
601 fIntestine_sml -> AddElement(elK,0.001);
602
603 //Large intestine ->ID 38
604 d = 1.040 *g/cm3;
605 fIntestine_lrg = new G4Material("intestine_lrg",d,9);
606 fIntestine_lrg -> AddElement(elH,0.105);
607 fIntestine_lrg -> AddElement(elC,0.114);
608 fIntestine_lrg -> AddElement(elN,0.025);
609 fIntestine_lrg -> AddElement(elO,0.750);
610 fIntestine_lrg -> AddElement(elNa,0.001);
611 fIntestine_lrg -> AddElement(elP,0.001);
612 fIntestine_lrg -> AddElement(elS,0.001);
613 fIntestine_lrg -> AddElement(elCl,0.002);
614 fIntestine_lrg -> AddElement(elK,0.001);
615
616 //Spleen -> ID 39
617 d = 1.040 *g/cm3;
618 fSpleen = new G4Material("spleen",d,9);
619 fSpleen -> AddElement(elH,0.103);
620 fSpleen -> AddElement(elC,0.112);
621 fSpleen -> AddElement(elN,0.032);
622 fSpleen -> AddElement(elO,0.743);
623 fSpleen -> AddElement(elNa,0.001);
624 fSpleen -> AddElement(elP,0.002);
625 fSpleen -> AddElement(elS,0.002);
626 fSpleen -> AddElement(elCl,0.002);
627 fSpleen -> AddElement(elK,0.003);
628
629 //Thyroid -> ID 40
630 d = 1.040 *g/cm3;
631 fThyroid = new G4Material("thyroid",d,10);
632 fThyroid -> AddElement(elH,0.104);
633 fThyroid -> AddElement(elC,0.118);
634 fThyroid -> AddElement(elN,0.025);
635 fThyroid -> AddElement(elO,0.745);
636 fThyroid -> AddElement(elNa,0.002);
637 fThyroid -> AddElement(elP,0.001);
638 fThyroid -> AddElement(elS,0.001);
639 fThyroid -> AddElement(elCl,0.002);
640 fThyroid -> AddElement(elK,0.001);
641 fThyroid -> AddElement(elI,0.001);
642
643 //Urinary Bladder -> ID 41
644 d = 1.040 *g/cm3;
645 fBladder = new G4Material("bladder",d,9);
646 fBladder -> AddElement(elH,0.105);
647 fBladder -> AddElement(elC,0.096);
648 fBladder -> AddElement(elN,0.026);
649 fBladder -> AddElement(elO,0.761);
650 fBladder -> AddElement(elNa,0.002);
651 fBladder -> AddElement(elP,0.002);
652 fBladder -> AddElement(elS,0.002);
653 fBladder -> AddElement(elCl,0.003);
654 fBladder -> AddElement(elK,0.003);
655
656 // Ovaries (Defined as ovaries_testes to be called upon by ColourMap.dat
657 //for visualisation purposes) -> ID 42
658 d = 1.040 *g/cm3;
659 fOvaries_testes = new G4Material("ovaries_testes",d,9);
660 fOvaries_testes -> AddElement(elH,0.105);
661 fOvaries_testes -> AddElement(elC,0.094);
662 fOvaries_testes -> AddElement(elN,0.025);
663 fOvaries_testes -> AddElement(elO,0.766);
664 fOvaries_testes -> AddElement(elNa,0.002);
665 fOvaries_testes -> AddElement(elP,0.002);
666 fOvaries_testes -> AddElement(elS,0.002);
667 fOvaries_testes -> AddElement(elCl,0.002);
668 fOvaries_testes -> AddElement(elK,0.002);
669
670 //Adrenals -> ID 43
671 d = 1.030 *g/cm3;
672 fAdrenals = new G4Material("adrenals",d,9);
673 fAdrenals -> AddElement(elH,0.104);
674 fAdrenals -> AddElement(elC,0.228);
675 fAdrenals -> AddElement(elN,0.028);
676 fAdrenals -> AddElement(elO,0.630);
677 fAdrenals -> AddElement(elNa,0.001);
678 fAdrenals -> AddElement(elP,0.002);
679 fAdrenals -> AddElement(elS,0.003);
680 fAdrenals -> AddElement(elCl,0.002);
681 fAdrenals -> AddElement(elK,0.002);
682
683 //Oesophagus -> ID 44
684 d = 1.030 *g/cm3;
685 fOesophagus = new G4Material("oesophagus",d,9);
686 fOesophagus -> AddElement(elH,0.104);
687 fOesophagus -> AddElement(elC,0.222);
688 fOesophagus -> AddElement(elN,0.028);
689 fOesophagus -> AddElement(elO,0.636);
690 fOesophagus -> AddElement(elNa,0.001);
691 fOesophagus -> AddElement(elP,0.002);
692 fOesophagus -> AddElement(elS,0.003);
693 fOesophagus -> AddElement(elCl,0.002);
694 fOesophagus -> AddElement(elK,0.002);
695
696 //Miscillaneous (Gallbladder, Trachea, Thymus, Tonsils, Ureters, ...) -> ID 45
697 d = 1.030 *g/cm3;
698 fMisc = new G4Material("misc",d,9);
699 fMisc -> AddElement(elH,0.105);
700 fMisc -> AddElement(elC,0.235);
701 fMisc -> AddElement(elN,0.028);
702 fMisc -> AddElement(elO,0.622);
703 fMisc -> AddElement(elNa,0.001);
704 fMisc -> AddElement(elP,0.002);
705 fMisc -> AddElement(elS,0.003);
706 fMisc -> AddElement(elCl,0.002);
707 fMisc -> AddElement(elK,0.002);
708
709 //Uterus (Defined as fUterus_prostate to be called upon by ColourMap.dat
710 // for visualisation purposes) -> ID 46
711 d = 1.030 *g/cm3;
712 fUterus_prostate = new G4Material("uterus_prostate",d,9);
713 fUterus_prostate -> AddElement(elH,0.105);
714 fUterus_prostate -> AddElement(elC,0.286);
715 fUterus_prostate -> AddElement(elN,0.025);
716 fUterus_prostate -> AddElement(elO,0.576);
717 fUterus_prostate -> AddElement(elNa,0.001);
718 fUterus_prostate -> AddElement(elP,0.002);
719 fUterus_prostate -> AddElement(elS,0.002);
720 fUterus_prostate -> AddElement(elCl,0.001);
721 fUterus_prostate -> AddElement(elK,0.002);
722
723 //Lymph -> ID 47
724 d = 1.030 *g/cm3;
725 fLymph = new G4Material("lymph",d,7);
726 fLymph -> AddElement(elH,0.108);
727 fLymph -> AddElement(elC,0.042);
728 fLymph -> AddElement(elN,0.011);
729 fLymph -> AddElement(elO,0.831);
730 fLymph -> AddElement(elNa,0.003);
731 fLymph -> AddElement(elS,0.001);
732 fLymph -> AddElement(elCl,0.004);
733
734 //Breast (Mammary Gland)-> ID 48
735 d = 1.020 *g/cm3;
736 fBreast_glandular = new G4Material("breast_glandular",d,4);
737 fBreast_glandular -> AddElement(elH,0.114);
738 fBreast_glandular -> AddElement(elC,0.461);
739 fBreast_glandular -> AddElement(elN,0.005);
740 fBreast_glandular -> AddElement(elO,0.420);
741
742 //Adipose tissue (breast) -> ID 49
743 d = 0.950 *g/cm3;
744 fBreast_adipose = new G4Material("breast_adipose",d,7);
745 fBreast_adipose -> AddElement(elH,0.114);
746 fBreast_adipose -> AddElement(elC,0.589);
747 fBreast_adipose -> AddElement(elN,0.007);
748 fBreast_adipose -> AddElement(elO,0.287);
749 fBreast_adipose -> AddElement(elNa,0.001);
750 fBreast_adipose -> AddElement(elS,0.001);
751 fBreast_adipose -> AddElement(elCl,0.001);
752
753 //Lung Tissue (Compressed Lung) -> ID 50
754 d = 0.385 *g/cm3;
755 fLung = new G4Material("lung",d,9);
756 fLung -> AddElement(elH,0.103);
757 fLung -> AddElement(elC,0.107);
758 fLung -> AddElement(elN,0.032);
759 fLung -> AddElement(elO,0.746);
760 fLung -> AddElement(elNa,0.002);
761 fLung -> AddElement(elP,0.002);
762 fLung -> AddElement(elS,0.003);
763 fLung -> AddElement(elCl,0.003);
764 fLung -> AddElement(elK,0.002);
765
766 //Contents of gastro-intestinal tract -> ID 51
767 d = 1.040 *g/cm3;
768 fGastro_content = new G4Material("gastro_content",d,10);
769 fGastro_content -> AddElement(elH,0.100);
770 fGastro_content -> AddElement(elC,0.222);
771 fGastro_content -> AddElement(elN,0.022);
772 fGastro_content -> AddElement(elO,0.644);
773 fGastro_content -> AddElement(elNa,0.001);
774 fGastro_content -> AddElement(elP,0.002);
775 fGastro_content -> AddElement(elS,0.003);
776 fGastro_content -> AddElement(elCl,0.001);
777 fGastro_content -> AddElement(elK,0.004);
778 fGastro_content -> AddElement(elCa,0.001);
779
780 //Urine -> ID 52
781 d = 1.040 *g/cm3;
782 fUrine = new G4Material("urine",d,7);
783 fUrine -> AddElement(elH,0.107);
784 fUrine -> AddElement(elC,0.003);
785 fUrine -> AddElement(elN,0.010);
786 fUrine -> AddElement(elO,0.873);
787 fUrine -> AddElement(elNa,0.004);
788 fUrine -> AddElement(elP,0.001);
789 fUrine -> AddElement(elK,0.002);
790
791 }
792
793 G4Material* ICRP110PhantomMaterial_Female::GetMaterial(G4String material)
794 {
795 // Returns a material
796 G4Material* pttoMaterial = G4Material::GetMaterial(material);
797 if (!pttoMaterial) G4cout << "WARNING: material '" << material << "' is not defined!" << G4endl;
798 return pttoMaterial;
799 }
800