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