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