141{
142
143 G4double z,
a, density;
145
147 symbol = "C",
148 z = 6.0,
a = 12.011 * g/mole );
150 symbol = "H",
151 z = 1.0,
a = 1.008 * g/mole );
153 symbol = "N",
154 z = 7.0,
a = 14.007 * g/mole );
156 symbol = "O",
157 z = 8.0,
a = 16.00 * g/mole );
159 symbol = "Na",
160 z= 11.0,
a = 22.98977* g/mole );
162 symbol = "Mg",
163 z = 12.0,
a = 24.3050* g/mole );
165 symbol = "P",
166 z = 15.0,
a = 30.973976* g/mole );
168 symbol = "S",
169 z = 16.0,
a = 32.065* g/mole );
171 symbol = "P",
172 z = 17.0,
a = 35.453* g/mole );
174 symbol = "P",
175 z = 19.0,
a = 30.0983* g/mole );
176
178 symbol = "Fe",
179 z = 26,
a = 56.845* g/mole );
180
182 symbol = "Ca",
183 z = 20.0,
a = 40.078* g/mole );
184
186 symbol = "Zn",
187 z = 30.0,
a = 65.382* g/mole );
188
189
190 G4int numberofElements;
191
192
194 1.290*mg/cm3,
195 numberofElements = 2 );
196 fAir->AddElement(elN, 0.7);
197 fAir->AddElement(elO, 0.3);
198
199
200
202 numberofElements = 13);
203 softTissue->AddElement(elH, 10.4472*perCent);
204 softTissue->AddElement(elC, 23.219*perCent);
205 softTissue->AddElement(elN, 2.488*perCent);
206 softTissue->AddElement(elO, 63.0238*perCent);
207 softTissue->AddElement(elNa, 0.113*perCent);
208 softTissue->AddElement(elMg, 0.0113*perCent);
209 softTissue->AddElement(elP, 0.113*perCent);
210 softTissue->AddElement(elS, 0.199*perCent);
211 softTissue->AddElement(elCl, 0.134*perCent);
212 softTissue->AddElement(elK, 0.199*perCent);
213 softTissue->AddElement(elCa, 0.023*perCent);
214 softTissue->AddElement(elFe, 0.005*perCent);
215 softTissue->AddElement(elZn, 0.003*perCent);
216
217
219 density = 0.217*g/cm3,
220 numberofElements = 9);
221 lunginhale->AddElement(elH,0.103);
222 lunginhale->AddElement(elC,0.105);
223 lunginhale->AddElement(elN,0.031);
224 lunginhale->AddElement(elO,0.749);
225 lunginhale->AddElement(elNa,0.002);
226 lunginhale->AddElement(elP,0.002);
227 lunginhale->AddElement(elS,0.003);
228 lunginhale->AddElement(elCl,0.002);
229 lunginhale->AddElement(elK,0.003);
230
231
233 density = 0.508*g/cm3,
234 numberofElements = 9 );
235 lungexhale->AddElement(elH,0.103);
236 lungexhale->AddElement(elC,0.105);
237 lungexhale->AddElement(elN,0.031);
238 lungexhale->AddElement(elO,0.749);
239 lungexhale->AddElement(elNa,0.002);
240 lungexhale->AddElement(elP,0.002);
241 lungexhale->AddElement(elS,0.003);
242 lungexhale->AddElement(elCl,0.002);
243 lungexhale->AddElement(elK,0.003);
244
245
247 density = 0.967*g/cm3,
248 numberofElements = 7);
249 adiposeTissue->AddElement(elH,0.114);
250 adiposeTissue->AddElement(elC,0.598);
251 adiposeTissue->AddElement(elN,0.007);
252 adiposeTissue->AddElement(elO,0.278);
253 adiposeTissue->AddElement(elNa,0.001);
254 adiposeTissue->AddElement(elS,0.001);
255 adiposeTissue->AddElement(elCl,0.001);
256
257
259 numberofElements = 13);
260 brainTissue->AddElement(elH, 11.0667*perCent);
261 brainTissue->AddElement(elC, 12.542*perCent);
262 brainTissue->AddElement(elN, 1.328*perCent);
263 brainTissue->AddElement(elO, 73.7723*perCent);
264 brainTissue->AddElement(elNa, 0.1840*perCent);
265 brainTissue->AddElement(elMg, 0.015*perCent);
266 brainTissue->AddElement(elP, 0.356*perCent);
267 brainTissue->AddElement(elS, 0.177*perCent);
268 brainTissue->AddElement(elCl, 0.236*perCent);
269 brainTissue->AddElement(elK, 0.31*perCent);
270 brainTissue->AddElement(elCa, 0.009*perCent);
271 brainTissue->AddElement(elFe, 0.005*perCent);
272 brainTissue->AddElement(elZn, 0.001*perCent);
273
274
275
277 density = 0.990*g/cm3,
278 numberofElements = 8 );
279 breast->AddElement(elH,0.109);
280 breast->AddElement(elC,0.506);
281 breast->AddElement(elN,0.023);
282 breast->AddElement(elO,0.358);
283 breast->AddElement(elNa,0.001);
284 breast->AddElement(elP,0.001);
285 breast->AddElement(elS,0.001);
286 breast->AddElement(elCl,0.001);
287
288
290 numberofElements = 8);
291 spinalDisc->AddElement(elH, 9.60*perCent);
292 spinalDisc->AddElement(elC, 9.90*perCent);
293 spinalDisc->AddElement(elN, 2.20*perCent);
294 spinalDisc->AddElement(elO, 74.40*perCent);
295 spinalDisc->AddElement(elNa, 0.50*perCent);
296 spinalDisc->AddElement(elP, 2.20*perCent);
297 spinalDisc->AddElement(elS, 0.90*perCent);
298 spinalDisc->AddElement(elCl, 0.30*perCent);
299
300
301
303 density = 1.0*g/cm3,
304 numberofElements = 2 );
305 water->AddElement(elH,0.112);
306 water->AddElement(elO,0.888);
307
308
310 density = 1.061*g/cm3,
311 numberofElements = 9 );
312 muscle->AddElement(elH,0.102);
313 muscle->AddElement(elC,0.143);
314 muscle->AddElement(elN,0.034);
315 muscle->AddElement(elO,0.710);
316 muscle->AddElement(elNa,0.001);
317 muscle->AddElement(elP,0.002);
318 muscle->AddElement(elS,0.003);
319 muscle->AddElement(elCl,0.001);
320 muscle->AddElement(elK,0.004);
321
322
324 density = 1.071*g/cm3,
325 numberofElements = 9);
326 liver->AddElement(elH,0.102);
327 liver->AddElement(elC,0.139);
328 liver->AddElement(elN,0.030);
329 liver->AddElement(elO,0.716);
330 liver->AddElement(elNa,0.002);
331 liver->AddElement(elP,0.003);
332 liver->AddElement(elS,0.003);
333 liver->AddElement(elCl,0.002);
334 liver->AddElement(elK,0.003);
335
336
338 numberofElements = 10);
339 toothDentin->AddElement(elH, 2.67*perCent);
340 toothDentin->AddElement(elC, 12.77*perCent);
341 toothDentin->AddElement(elN, 4.27*perCent);
342 toothDentin->AddElement(elO, 40.40*perCent);
343 toothDentin->AddElement(elNa, 0.65*perCent);
344 toothDentin->AddElement(elMg, 0.59*perCent);
345 toothDentin->AddElement(elP, 11.86*perCent);
346 toothDentin->AddElement(elCl, 0.04*perCent);
347 toothDentin->AddElement(elCa, 26.74*perCent);
348 toothDentin->AddElement(elZn, 0.01*perCent);
349
350
351
353 density = 1.159*g/cm3,
354 numberofElements = 12 );
355 trabecularBone->AddElement(elH,0.085);
356 trabecularBone->AddElement(elC,0.404);
357 trabecularBone->AddElement(elN,0.058);
358 trabecularBone->AddElement(elO,0.367);
359 trabecularBone->AddElement(elNa,0.001);
360 trabecularBone->AddElement(elMg,0.001);
361 trabecularBone->AddElement(elP,0.034);
362 trabecularBone->AddElement(elS,0.002);
363 trabecularBone->AddElement(elCl,0.002);
364 trabecularBone->AddElement(elK,0.001);
365 trabecularBone->AddElement(elCa,0.044);
366 trabecularBone->AddElement(elFe,0.001);
367
368
369
371 1.18 * g/cm3,
372 numberofElements = 12);
373 trabecularBone_head->AddElement(elH, 8.50*perCent);
374 trabecularBone_head->AddElement(elC, 40.40*perCent);
375 trabecularBone_head->AddElement(elN, 2.80*perCent);
376 trabecularBone_head->AddElement(elO, 36.70*perCent);
377 trabecularBone_head->AddElement(elNa, 0.10*perCent);
378 trabecularBone_head->AddElement(elMg, 0.10*perCent);
379 trabecularBone_head->AddElement(elP, 3.40*perCent);
380 trabecularBone_head->AddElement(elS, 0.20*perCent);
381 trabecularBone_head->AddElement(elCl, 0.20*perCent);
382 trabecularBone_head->AddElement(elK, 0.10*perCent);
383 trabecularBone_head->AddElement(elCa, 7.40*perCent);
384 trabecularBone_head->AddElement(elFe, 0.10*perCent);
385
386
388 density = 1.575*g/cm3,
389 numberofElements = 11 );
390 denseBone->AddElement(elH,0.056);
391 denseBone->AddElement(elC,0.235);
392 denseBone->AddElement(elN,0.050);
393 denseBone->AddElement(elO,0.434);
394 denseBone->AddElement(elNa,0.001);
395 denseBone->AddElement(elMg,0.001);
396 denseBone->AddElement(elP,0.072);
397 denseBone->AddElement(elS,0.003);
398 denseBone->AddElement(elCl,0.001);
399 denseBone->AddElement(elK,0.001);
400 denseBone->AddElement(elCa,0.146);
401
402
404 numberofElements = 9);
405 corticalBone->AddElement(elH, 4.7234*perCent);
406 corticalBone->AddElement(elC, 14.4330*perCent);
407 corticalBone->AddElement(elN, 4.199*perCent);
408 corticalBone->AddElement(elO, 44.6096*perCent);
409 corticalBone->AddElement(elMg, 0.22*perCent);
410 corticalBone->AddElement(elP, 10.497*perCent);
411 corticalBone->AddElement(elS, 0.315*perCent);
412 corticalBone->AddElement(elCa, 20.993*perCent);
413 corticalBone->AddElement(elZn, 0.01*perCent);
414
415
416
418 numberofElements = 10);
419 toothEnamel->AddElement(elH, 0.95*perCent);
420 toothEnamel->AddElement(elC, 1.11*perCent);
421 toothEnamel->AddElement(elN, 0.23*perCent);
422 toothEnamel->AddElement(elO,41.66*perCent);
423 toothEnamel->AddElement(elNa, 0.79*perCent);
424 toothEnamel->AddElement(elMg, 0.23*perCent);
425 toothEnamel->AddElement(elP, 18.71*perCent);
426 toothEnamel->AddElement(elCl, 0.34*perCent);
427 toothEnamel->AddElement(elCa, 35.97*perCent);
428 toothEnamel->AddElement(elZn, 0.02*perCent);
429
430#ifdef DICOM_USE_HEAD
431
440 G4cout << "The materials of the DICOM Head have been used" << G4endl;
441#else
452 G4cout << "Default materials of the DICOM Extended examples have been used"
453 << G4endl;
454#endif
455}
std::vector< ExP01TrackerHit * > a
std::vector< G4Material * > fOriginalMaterials