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