ICRP110_HumanPhantoms/src/ICRP110PhantomMaterial_Male.cc

0001 //
0002 // ********************************************************************
0003 // * License and Disclaimer                                           *
0004 // *                                                                  *
0005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
0006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
0007 // * conditions of the Geant4 Software License,  included in the file *
0008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
0009 // * include a list of copyright holders.                             *
0010 // *                                                                  *
0011 // * Neither the authors of this software system, nor their employing *
0012 // * institutes,nor the agencies providing financial support for this *
0013 // * work  make  any representation or  warranty, express or implied, *
0014 // * regarding  this  software system or assume any liability for its *
0015 // * use.  Please see the license in the file  LICENSE  and URL above *
0016 // * for the full disclaimer and the limitation of liability.         *
0017 // *                                                                  *
0018 // * This  code  implementation is the result of  the  scientific and *
0019 // * technical work of the GEANT4 collaboration.                      *
0020 // * By using,  copying,  modifying or  distributing the software (or *
0021 // * any work based  on the software)  you  agree  to acknowledge its *
0022 // * use  in  resulting  scientific  publications,  and indicate your *
0023 // * acceptance of all terms of the Geant4 Software license.          *
0024 // ********************************************************************
0025 //
0026 // Code developed by:
0027 // S.Guatelli, M. Large and A. Malaroda, University of Wollongong
0028 //
0029 #include "ICRP110PhantomMaterial_Male.hh"
0030 #include "globals.hh"
0031 #include "G4SystemOfUnits.hh"
0032 #include "G4MaterialPropertiesTable.hh"
0033 #include "G4MaterialPropertyVector.hh"
0034 #include "G4MaterialTable.hh"
0035 #include "Randomize.hh"
0036 #include "G4RunManager.hh"
0037 #include "G4Element.hh"
0038 #include "G4ElementTable.hh"
0039
0040 ICRP110PhantomMaterial_Male::ICRP110PhantomMaterial_Male():
0041   fLung(nullptr), fTeeth(nullptr), fBone(nullptr), fHumeri_upper(nullptr), fHumeri_lower(nullptr),
0042   fArm_lower(nullptr), fHand(nullptr), fClavicle(nullptr), fCranium(nullptr), fFemora_upper(nullptr),
0043   fFemora_lower(nullptr), fLeg_lower(nullptr), fFoot(nullptr), fMandible(nullptr),
0044   fPelvis(nullptr), fRibs(nullptr), fScapulae(nullptr), fSpine_cervical(nullptr),
0045   fSpine_lumbar(nullptr), fSpine_thoratic(nullptr), fSacrum(nullptr),
0046   fSternum(nullptr), fHf_upper(nullptr), fHf_lower(nullptr), fMed_lowerleg(nullptr),
0047   fMed_lowerarm(nullptr), fCartilage(nullptr), fSkin(nullptr), fBlood(nullptr),
0048   fMuscle(nullptr), fLiver(nullptr), fPancreas(nullptr), fBrain(nullptr), fHeart(nullptr), fEye(nullptr),
0049   fKidney(nullptr), fStomach(nullptr), fIntestine_sml(nullptr), fIntestine_lrg(nullptr),
0050   fSpleen(nullptr), fThyroid(nullptr), fBladder(nullptr), fOvaries_testes(nullptr), fAdrenals(nullptr),
0051   fOesophagus(nullptr), fMisc(nullptr), fUterus_prostate(nullptr), fLymph(nullptr),
0052   fBreast_glandular(nullptr), fBreast_adipose(nullptr), fGastro_content(nullptr),
0053   fUrine(nullptr)
0054 {;}
0055
0056 ICRP110PhantomMaterial_Male::~ICRP110PhantomMaterial_Male()
0057 {;}
0058
0059 void ICRP110PhantomMaterial_Male::DefineMaterials()
0060 {
0061   // Define required materials
0062
0063   G4double A;  // atomic mass
0064   G4double Z;  // atomic number
0065   G4double d;  // density
0066
0067   // General elements
0068
0069   A = 1.01*g/mole;
0070   auto elH = new G4Element ("Hydrogen","H",Z = 1.,A);
0071
0072   A = 12.011*g/mole;
0073   auto elC = new G4Element("Carbon","C",Z = 6.,A);
0074
0075   A = 14.01*g/mole;
0076   auto elN = new G4Element("Nitrogen","N",Z = 7.,A);
0077
0078   A = 16.00*g/mole;
0079   auto elO = new G4Element("Oxygen","O",Z = 8.,A);
0080
0081   A = 22.99*g/mole;
0082   auto elNa = new G4Element("Sodium","Na",Z = 11.,A);
0083
0084   A = 24.305*g/mole;
0085   auto elMg = new G4Element("Magnesium","Mg",Z = 12.,A);
0086
0087   A = 30.974*g/mole;
0088   auto elP = new G4Element("Phosphorus","P",Z = 15.,A);
0089
0090   A = 32.064*g/mole;
0091   auto elS = new G4Element("Sulfur","S",Z = 16.,A);
0092
0093   A = 35.453*g/mole;
0094   auto elCl = new G4Element("Chlorine","Cl",Z = 17.,A);
0095
0096   A = 39.098*g/mole;
0097   auto elK = new G4Element("Potassium","K",Z = 19.,A);
0098
0099   A = 40.08*g/mole;
0100   auto elCa = new G4Element("Calcium","Ca",Z = 20.,A);
0101
0102   A = 55.85*g/mole;
0103   auto elFe  = new G4Element("Iron","Fe",Z = 26.,A);
0104
0105   A = 126.90447 *g/mole;
0106   auto elI = new G4Element("Iodine","I", Z = 53.,A);
0107
0108
0109   //Added tissues for phantom including their tissue ID
0110   //Teeth -> Tissue ID 1
0111   d = 2.750 *g/cm3;
0112   fTeeth = new G4Material("teeth",d,7);
0113   fTeeth -> AddElement(elH,0.022);
0114   fTeeth -> AddElement(elC,0.095);
0115   fTeeth -> AddElement(elN,0.029);
0116   fTeeth -> AddElement(elO,0.421);
0117   fTeeth -> AddElement(elMg,0.007);
0118   fTeeth -> AddElement(elP,0.137);
0119   fTeeth -> AddElement(elCa,0.289);
0120
0121   //Mineral Bone -> Tissue ID 2
0122   d = 1.920 *g/cm3;
0123   fBone = new G4Material("bone",d,9);
0124   fBone -> AddElement(elH,0.036);
0125   fBone -> AddElement(elC,0.159);
0126   fBone -> AddElement(elN,0.042);
0127   fBone -> AddElement(elO,0.448);
0128   fBone -> AddElement(elNa,0.003);
0129   fBone -> AddElement(elMg,0.002);
0130   fBone -> AddElement(elP,0.094);
0131   fBone -> AddElement(elS,0.003);
0132   fBone -> AddElement(elCa,0.213);
0133
0134   //Humeri, upper half, spongiosa -> ID 3
0135   d = 1.205 *g/cm3;
0136   fHumeri_upper = new G4Material("humeri_upper",d,10);
0137   fHumeri_upper -> AddElement(elH,0.085) ;
0138   fHumeri_upper -> AddElement(elC,0.288);
0139   fHumeri_upper -> AddElement(elN,0.026);
0140   fHumeri_upper -> AddElement(elO,0.498);
0141   fHumeri_upper -> AddElement(elNa,0.002);
0142   fHumeri_upper -> AddElement(elMg,0.001);
0143   fHumeri_upper -> AddElement(elP,0.033);
0144   fHumeri_upper -> AddElement(elS,0.004);
0145   fHumeri_upper -> AddElement(elCl,0.002);
0146   fHumeri_upper -> AddElement(elCa,0.061);
0147
0148   //Humeri, lower half, spongiosa -> ID 4
0149   d = 1.108 *g/cm3;
0150   fHumeri_lower = new G4Material("humeri_lower",d,9);
0151   fHumeri_lower -> AddElement(elH,0.097);
0152   fHumeri_lower -> AddElement(elC,0.439);
0153   fHumeri_lower -> AddElement(elN,0.017);
0154   fHumeri_lower -> AddElement(elO,0.381);
0155   fHumeri_lower -> AddElement(elNa,0.002);
0156   fHumeri_lower -> AddElement(elP,0.021);
0157   fHumeri_lower -> AddElement(elS,0.003);
0158   fHumeri_lower -> AddElement(elCl,0.001);
0159   fHumeri_lower -> AddElement(elCa,0.039);
0160
0161   //Lower arm bones, spongiosa -> ID 5
0162   d = 1.108 *g/cm3;
0163   fArm_lower = new G4Material("arm_lower",d,9);
0164   fArm_lower -> AddElement(elH,0.097);
0165   fArm_lower -> AddElement(elC,0.439);
0166   fArm_lower -> AddElement(elN,0.017);
0167   fArm_lower -> AddElement(elO,0.381);
0168   fArm_lower -> AddElement(elNa,0.002);
0169   fArm_lower -> AddElement(elP,0.021);
0170   fArm_lower -> AddElement(elS,0.003);
0171   fArm_lower -> AddElement(elCl,0.001);
0172   fArm_lower -> AddElement(elCa,0.039);
0173
0174   //Hand Bones, Spongiosa ->ID 6
0175   d = 1.108 *g/cm3;
0176   fHand = new G4Material("hand",d,9);
0177   fHand -> AddElement(elH,0.097);
0178   fHand -> AddElement(elC,0.439);
0179   fHand -> AddElement(elN,0.017);
0180   fHand -> AddElement(elO,0.381);
0181   fHand -> AddElement(elNa,0.002);
0182   fHand -> AddElement(elP,0.021);
0183   fHand -> AddElement(elS,0.003);
0184   fHand -> AddElement(elCl,0.001);
0185   fHand -> AddElement(elCa,0.039);
0186
0187   //Clavicles, spongiosa -> ID 7
0188   d = 1.151 *g/cm3;
0189   fClavicle = new G4Material("clavicle",d,9);
0190   fClavicle -> AddElement(elH,0.091);
0191   fClavicle -> AddElement(elC,0.348);
0192   fClavicle -> AddElement(elN,0.024);
0193   fClavicle -> AddElement(elO,0.457);
0194   fClavicle -> AddElement(elNa,0.002);
0195   fClavicle -> AddElement(elP,0.026);
0196   fClavicle -> AddElement(elS,0.003);
0197   fClavicle -> AddElement(elCl,0.001);
0198   fClavicle -> AddElement(elCa,0.048);
0199
0200   //cranium, spongiosa -> ID 8
0201   d = 1.157 *g/cm3;
0202   fCranium = new G4Material("cranium",d,10);
0203   fCranium -> AddElement(elH,0.090);
0204   fCranium -> AddElement(elC,0.335);
0205   fCranium -> AddElement(elN,0.025);
0206   fCranium -> AddElement(elO,0.467);
0207   fCranium -> AddElement(elNa,0.002);
0208   fCranium -> AddElement(elP,0.026);
0209   fCranium -> AddElement(elS,0.003);
0210   fCranium -> AddElement(elCl,0.002);
0211   fCranium -> AddElement(elK,0.001);
0212   fCranium -> AddElement(elCa,0.049);
0213
0214   //femora, upper half, spongiosa -> ID 9
0215   d = 1.124 *g/cm3;
0216   fFemora_upper = new G4Material("femora_upper",d,9);
0217   fFemora_upper->AddElement(elH,0.094);
0218   fFemora_upper->AddElement(elC,0.385);
0219   fFemora_upper->AddElement(elN,0.022);
0220   fFemora_upper->AddElement(elO,0.430);
0221   fFemora_upper->AddElement(elNa,0.002);
0222   fFemora_upper->AddElement(elP,0.022);
0223   fFemora_upper->AddElement(elS,0.003);
0224   fFemora_upper->AddElement(elCl,0.001);
0225   fFemora_upper->AddElement(elCa,0.041);
0226
0227   //femora, lower half, spongiosa -> ID 10
0228   d = 1.108 *g/cm3;
0229   fFemora_lower = new G4Material("femora_lower",d,9);
0230   fFemora_lower->AddElement(elH,0.097);
0231   fFemora_lower->AddElement(elC,0.439);
0232   fFemora_lower->AddElement(elN,0.017);
0233   fFemora_lower->AddElement(elO,0.381);
0234   fFemora_lower->AddElement(elNa,0.002);
0235   fFemora_lower->AddElement(elP,0.021);
0236   fFemora_lower->AddElement(elS,0.003);
0237   fFemora_lower->AddElement(elCl,0.001);
0238   fFemora_lower->AddElement(elCa,0.039);
0239
0240   //Lower leg bones, spongiosa -> ID 11
0241   d = 1.108 *g/cm3;
0242   fLeg_lower = new G4Material("leg_lower",d,9);
0243   fLeg_lower -> AddElement(elH,0.097);
0244   fLeg_lower -> AddElement(elC,0.439);
0245   fLeg_lower -> AddElement(elN,0.017);
0246   fLeg_lower -> AddElement(elO,0.381);
0247   fLeg_lower -> AddElement(elNa,0.002);
0248   fLeg_lower -> AddElement(elP,0.021);
0249   fLeg_lower -> AddElement(elS,0.003);
0250   fLeg_lower -> AddElement(elCl,0.001);
0251   fLeg_lower -> AddElement(elCa,0.039);
0252
0253   //Foot bones, spongiosa ->ID 12
0254   d = 1.108 *g/cm3;
0255   fFoot = new G4Material("foot",d,9);
0256   fFoot -> AddElement(elH,0.097);
0257   fFoot -> AddElement(elC,0.439);
0258   fFoot -> AddElement(elN,0.017);
0259   fFoot -> AddElement(elO,0.381);
0260   fFoot -> AddElement(elNa,0.002);
0261   fFoot -> AddElement(elP,0.021);
0262   fFoot -> AddElement(elS,0.003);
0263   fFoot -> AddElement(elCl,0.001);
0264   fFoot -> AddElement(elCa,0.039);
0265
0266   //Mandible, spongiosa -> ID 13
0267   d = 1.228 *g/cm3;
0268   fMandible = new G4Material("mandible",d,10);
0269   fMandible -> AddElement(elH,0.083);
0270   fMandible -> AddElement(elC,0.266);
0271   fMandible -> AddElement(elN,0.027);
0272   fMandible -> AddElement(elO,0.511);
0273   fMandible -> AddElement(elNa,0.003);
0274   fMandible -> AddElement(elMg,0.001);
0275   fMandible -> AddElement(elP,0.036);
0276   fMandible -> AddElement(elS,0.004);
0277   fMandible -> AddElement(elCl,0.002);
0278   fMandible -> AddElement(elCa,0.067);
0279
0280   //Pelvis, Spongiosa -> ID 14
0281   d = 1.123 *g/cm3;
0282   fPelvis = new G4Material("pelvis",d,10);
0283   fPelvis -> AddElement(elH,0.094);
0284   fPelvis -> AddElement(elC,0.360);
0285   fPelvis -> AddElement(elN,0.025);
0286   fPelvis -> AddElement(elO,0.454);
0287   fPelvis -> AddElement(elNa,0.002);
0288   fPelvis -> AddElement(elP,0.021);
0289   fPelvis -> AddElement(elS,0.003);
0290   fPelvis -> AddElement(elCl,0.002);
0291   fPelvis -> AddElement(elK,0.001);
0292   fPelvis -> AddElement(elCa,0.038);
0293
0294   //Ribs, spongiosa -> ID 15
0295   d = 1.165 *g/cm3;
0296   fRibs = new G4Material("ribs",d,10);
0297   fRibs -> AddElement(elH,0.089);
0298   fRibs -> AddElement(elC,0.292);
0299   fRibs -> AddElement(elN,0.029);
0300   fRibs -> AddElement(elO,0.507);
0301   fRibs -> AddElement(elNa,0.002);
0302   fRibs -> AddElement(elP,0.026);
0303   fRibs -> AddElement(elS,0.004);
0304   fRibs -> AddElement(elCl,0.002);
0305   fRibs -> AddElement(elK,0.001);
0306   fRibs -> AddElement(elCa,0.048);
0307
0308   //scapulae, spongiosa -> ID 16
0309   d = 1.183 *g/cm3;
0310   fScapulae = new G4Material("scapulae",d,10);
0311   fScapulae -> AddElement(elH,0.087);
0312   fScapulae -> AddElement(elC,0.309);
0313   fScapulae -> AddElement(elN,0.026);
0314   fScapulae -> AddElement(elO,0.483);
0315   fScapulae -> AddElement(elNa,0.002);
0316   fScapulae -> AddElement(elMg,0.001);
0317   fScapulae -> AddElement(elP,0.030);
0318   fScapulae -> AddElement(elS,0.004);
0319   fScapulae -> AddElement(elCl,0.002);
0320   fScapulae -> AddElement(elCa,0.056);
0321
0322   //Cervical fSpine, spongiosa -> ID 17
0323   d = 1.050 *g/cm3;
0324   fSpine_cervical = new G4Material("spine_cervical",d,11);
0325   fSpine_cervical -> AddElement(elH,0.103);
0326   fSpine_cervical -> AddElement(elC,0.400);
0327   fSpine_cervical -> AddElement(elN,0.027);
0328   fSpine_cervical -> AddElement(elO,0.444);
0329   fSpine_cervical -> AddElement(elNa,0.001);
0330   fSpine_cervical -> AddElement(elP,0.007);
0331   fSpine_cervical -> AddElement(elS,0.002);
0332   fSpine_cervical -> AddElement(elCl,0.002);
0333   fSpine_cervical -> AddElement(elK,0.001);
0334   fSpine_cervical -> AddElement(elCa,0.012);
0335   fSpine_cervical -> AddElement(elFe,0.001);
0336
0337   //Thoratic spine, spongiosa -> ID 18
0338   d = 1.074 *g/cm3;
0339   fSpine_thoratic = new G4Material("spine_thoratic",d,11);
0340   fSpine_thoratic -> AddElement(elH,0.099);
0341   fSpine_thoratic -> AddElement(elC,0.376);
0342   fSpine_thoratic -> AddElement(elN,0.027);
0343   fSpine_thoratic -> AddElement(elO,0.459);
0344   fSpine_thoratic -> AddElement(elNa,0.001);
0345   fSpine_thoratic -> AddElement(elP,0.012);
0346   fSpine_thoratic -> AddElement(elS,0.002);
0347   fSpine_thoratic -> AddElement(elCl,0.002);
0348   fSpine_thoratic -> AddElement(elK,0.001);
0349   fSpine_thoratic -> AddElement(elCa,0.020);
0350   fSpine_thoratic -> AddElement(elFe,0.001);
0351
0352   //Lumbar spine, spongiosa -> ID 19
0353   d = 1.112 *g/cm3;
0354   fSpine_lumbar = new G4Material("spine_lumbar",d,10);
0355   fSpine_lumbar -> AddElement(elH,0.095);
0356   fSpine_lumbar -> AddElement(elC,0.340);
0357   fSpine_lumbar -> AddElement(elN,0.028);
0358   fSpine_lumbar -> AddElement(elO,0.480);
0359   fSpine_lumbar -> AddElement(elNa,0.001);
0360   fSpine_lumbar -> AddElement(elP,0.018);
0361   fSpine_lumbar -> AddElement(elS,0.003);
0362   fSpine_lumbar -> AddElement(elCl,0.002);
0363   fSpine_lumbar -> AddElement(elK,0.001);
0364   fSpine_lumbar -> AddElement(elCa,0.032);
0365
0366   //sacrum, spongiosa -> ID 20
0367   d = 1.031 *g/cm3;
0368   fSacrum = new G4Material("sacrum",d,11);
0369   fSacrum -> AddElement(elH,0.105);
0370   fSacrum -> AddElement(elC,0.419);
0371   fSacrum -> AddElement(elN,0.027);
0372   fSacrum -> AddElement(elO,0.432);
0373   fSacrum -> AddElement(elNa,0.001);
0374   fSacrum -> AddElement(elP,0.004);
0375   fSacrum -> AddElement(elS,0.002);
0376   fSacrum -> AddElement(elCl,0.002);
0377   fSacrum -> AddElement(elK,0.001);
0378   fSacrum -> AddElement(elCa,0.006);
0379   fSacrum -> AddElement(elFe,0.001);
0380
0381   //sternum, spongiosa -> ID 21
0382   d = 1.041 *g/cm3;
0383   fSternum = new G4Material("sternum",d,11);
0384   fSternum->AddElement(elH,0.104);
0385   fSternum->AddElement(elC,0.409);
0386   fSternum->AddElement(elN,0.027);
0387   fSternum->AddElement(elO,0.438);
0388   fSternum->AddElement(elNa,0.001);
0389   fSternum->AddElement(elP,0.006);
0390   fSternum->AddElement(elS,0.002);
0391   fSternum->AddElement(elCl,0.002);
0392   fSternum->AddElement(elK,0.001);
0393   fSternum->AddElement(elCa,0.009);
0394   fSternum->AddElement(elFe,0.001);
0395
0396   //Humeri and femora, upper halves, medullary cavity -> ID 22
0397   d = 0.980 *g/cm3;
0398   fHf_upper = new G4Material("hf_upper",d,7);
0399   fHf_upper -> AddElement(elH,0.115);
0400   fHf_upper -> AddElement(elC,0.636);
0401   fHf_upper -> AddElement(elN,0.007);
0402   fHf_upper -> AddElement(elO,0.239);
0403   fHf_upper -> AddElement(elNa,0.001);
0404   fHf_upper -> AddElement(elS,0.001);
0405   fHf_upper -> AddElement(elCl,0.001);
0406
0407   //Humeri and femora, lower halves, medullary cavity -> ID 23
0408   d = 0.980 *g/cm3;
0409   fHf_lower = new G4Material("hf_lower",d,7);
0410   fHf_lower -> AddElement(elH,0.115);
0411   fHf_lower -> AddElement(elC,0.636);
0412   fHf_lower -> AddElement(elN,0.007);
0413   fHf_lower -> AddElement(elO,0.239);
0414   fHf_lower -> AddElement(elNa,0.001);
0415   fHf_lower -> AddElement(elS,0.001);
0416   fHf_lower -> AddElement(elCl,0.001);
0417
0418   //Lower arm bones, medullary cavity -> ID 24
0419   d = 0.980 *g/cm3;
0420   fMed_lowerarm = new G4Material("med_lowerarm",d,7);
0421   fMed_lowerarm -> AddElement(elH,0.115);
0422   fMed_lowerarm -> AddElement(elC,0.636);
0423   fMed_lowerarm -> AddElement(elN,0.007);
0424   fMed_lowerarm -> AddElement(elO,0.239);
0425   fMed_lowerarm -> AddElement(elNa,0.001);
0426   fMed_lowerarm -> AddElement(elS,0.001);
0427   fMed_lowerarm -> AddElement(elCl,0.001);
0428
0429   //Lower leg bones, medullary cavity -> ID 25
0430   d = 0.980 *g/cm3;
0431   fMed_lowerleg = new G4Material("med_lowerleg",d,7);
0432   fMed_lowerleg -> AddElement(elH,0.115);
0433   fMed_lowerleg -> AddElement(elC,0.636);
0434   fMed_lowerleg -> AddElement(elN,0.007);
0435   fMed_lowerleg -> AddElement(elO,0.239);
0436   fMed_lowerleg -> AddElement(elNa,0.001);
0437   fMed_lowerleg -> AddElement(elS,0.001);
0438   fMed_lowerleg -> AddElement(elCl,0.001);
0439
0440   //Cartilage -> ID 26
0441   d = 1.100 *g/cm3;
0442   fCartilage = new G4Material("cartilage",d,8);
0443   fCartilage -> AddElement(elH,0.096);
0444   fCartilage -> AddElement(elC,0.099);
0445   fCartilage -> AddElement(elN,0.022);
0446   fCartilage -> AddElement(elO,0.744);
0447   fCartilage -> AddElement(elNa,0.005);
0448   fCartilage -> AddElement(elP,0.022);
0449   fCartilage -> AddElement(elS,0.009);
0450   fCartilage -> AddElement(elCl,0.003);
0451
0452   //Skin -> Id 27
0453   d = 1.090 *g/cm3;
0454   fSkin = new G4Material("skin",d,9);
0455   fSkin -> AddElement(elH,0.100);
0456   fSkin -> AddElement(elC,0.199);
0457   fSkin -> AddElement(elN,0.042);
0458   fSkin -> AddElement(elO,0.650);
0459   fSkin -> AddElement(elNa,0.002);
0460   fSkin -> AddElement(elP,0.001);
0461   fSkin -> AddElement(elS,0.002);
0462   fSkin -> AddElement(elCl,0.003);
0463   fSkin -> AddElement(elK,0.001);
0464
0465   //Blood -> ID 28
0466   d = 1.060 *g/cm3;
0467   fBlood = new G4Material("blood",d,10);
0468   fBlood -> AddElement(elH,0.102);
0469   fBlood -> AddElement(elC,0.110);
0470   fBlood -> AddElement(elN,0.033);
0471   fBlood -> AddElement(elO,0.745);
0472   fBlood -> AddElement(elNa,0.001);
0473   fBlood -> AddElement(elP,0.001);
0474   fBlood -> AddElement(elS,0.002);
0475   fBlood -> AddElement(elCl,0.003);
0476   fBlood -> AddElement(elK,0.002);
0477   fBlood -> AddElement(elFe,0.001);
0478
0479   //Muscular Tissue -> ID 29
0480   d = 1.050 *g/cm3;
0481   fMuscle = new G4Material("muscle",d,9);
0482   fMuscle -> AddElement(elH,0.102);
0483   fMuscle -> AddElement(elC,0.142);
0484   fMuscle -> AddElement(elN,0.034);
0485   fMuscle -> AddElement(elO,0.711);
0486   fMuscle -> AddElement(elNa,0.001);
0487   fMuscle -> AddElement(elP,0.002);
0488   fMuscle -> AddElement(elS,0.003);
0489   fMuscle -> AddElement(elCl,0.001);
0490   fMuscle -> AddElement(elK,0.004);
0491
0492   //Liver -> ID 30
0493   d = 1.050 *g/cm3;
0494   fLiver = new G4Material("liver",d,9);
0495   fLiver -> AddElement(elH,0.102);
0496   fLiver -> AddElement(elC,0.130);
0497   fLiver -> AddElement(elN,0.031);
0498   fLiver -> AddElement(elO,0.725);
0499   fLiver -> AddElement(elNa,0.002);
0500   fLiver -> AddElement(elP,0.002);
0501   fLiver -> AddElement(elS,0.003);
0502   fLiver -> AddElement(elCl,0.002);
0503   fLiver -> AddElement(elK,0.003);
0504
0505   //Pancreas ->ID 31
0506   d = 1.050 *g/cm3;
0507   fPancreas = new G4Material("pancreas",d,9);
0508   fPancreas -> AddElement(elH,0.105);
0509   fPancreas -> AddElement(elC,0.155);
0510   fPancreas -> AddElement(elN,0.025);
0511   fPancreas -> AddElement(elO,0.706);
0512   fPancreas -> AddElement(elNa,0.002);
0513   fPancreas -> AddElement(elP,0.002);
0514   fPancreas -> AddElement(elS,0.001);
0515   fPancreas -> AddElement(elCl,0.002);
0516   fPancreas -> AddElement(elK,0.002);
0517
0518   //Brain -> ID 32
0519   d = 1.050 *g/cm3;
0520   fBrain = new G4Material("brain",d,9);
0521   fBrain -> AddElement(elH,0.107);
0522   fBrain -> AddElement(elC,0.143);
0523   fBrain -> AddElement(elN,0.023);
0524   fBrain -> AddElement(elO,0.713);
0525   fBrain -> AddElement(elNa,0.002);
0526   fBrain -> AddElement(elP,0.004);
0527   fBrain -> AddElement(elS,0.002);
0528   fBrain -> AddElement(elCl,0.003);
0529   fBrain -> AddElement(elK,0.003);
0530
0531   //Heart -> ID 33
0532   d = 1.050 *g/cm3;
0533   fHeart = new G4Material("heart",d,9);
0534   fHeart -> AddElement(elH,0.104);
0535   fHeart -> AddElement(elC,0.138);
0536   fHeart -> AddElement(elN,0.029);
0537   fHeart -> AddElement(elO,0.719);
0538   fHeart -> AddElement(elNa,0.001);
0539   fHeart -> AddElement(elP,0.002);
0540   fHeart -> AddElement(elS,0.002);
0541   fHeart -> AddElement(elCl,0.002);
0542   fHeart -> AddElement(elK,0.003);
0543
0544   //Eye ->ID 34
0545   d = 1.050 *g/cm3;
0546   fEye = new G4Material("eye",d,8);
0547   fEye -> AddElement(elH,0.097);
0548   fEye -> AddElement(elC,0.181);
0549   fEye -> AddElement(elN,0.053);
0550   fEye -> AddElement(elO,0.663);
0551   fEye -> AddElement(elCa,0.001);
0552   fEye -> AddElement(elP,0.001);
0553   fEye -> AddElement(elS,0.003);
0554   fEye -> AddElement(elCl,0.001);
0555
0556   //Kidneys -> ID 35
0557   d = 1.050 *g/cm3;
0558   fKidney = new G4Material("kidney",d,10);
0559   fKidney -> AddElement(elH,0.103);
0560   fKidney -> AddElement(elC,0.124);
0561   fKidney -> AddElement(elN,0.031);
0562   fKidney -> AddElement(elO,0.731);
0563   fKidney -> AddElement(elNa,0.002);
0564   fKidney -> AddElement(elP,0.002);
0565   fKidney -> AddElement(elS,0.002);
0566   fKidney -> AddElement(elCl,0.002);
0567   fKidney -> AddElement(elK,0.002);
0568   fKidney -> AddElement(elCa,0.001);
0569
0570   //Stomach ->ID 36
0571   d = 1.040 *g/cm3;
0572   fStomach = new G4Material("stomach",d,9);
0573   fStomach -> AddElement(elH,0.105);
0574   fStomach -> AddElement(elC,0.114);
0575   fStomach -> AddElement(elN,0.025);
0576   fStomach -> AddElement(elO,0.750);
0577   fStomach -> AddElement(elNa,0.001);
0578   fStomach -> AddElement(elP,0.001);
0579   fStomach -> AddElement(elS,0.001);
0580   fStomach -> AddElement(elCl,0.002);
0581   fStomach -> AddElement(elK,0.001);
0582
0583   //Small intestine ->ID 37
0584   d = 1.040 *g/cm3;
0585   fIntestine_sml = new G4Material("intestine_sml",d,9);
0586   fIntestine_sml -> AddElement(elH,0.105);
0587   fIntestine_sml -> AddElement(elC,0.113);
0588   fIntestine_sml -> AddElement(elN,0.026);
0589   fIntestine_sml -> AddElement(elO,0.750);
0590   fIntestine_sml -> AddElement(elNa,0.001);
0591   fIntestine_sml -> AddElement(elP,0.001);
0592   fIntestine_sml -> AddElement(elS,0.001);
0593   fIntestine_sml -> AddElement(elCl,0.002);
0594   fIntestine_sml -> AddElement(elK,0.001);
0595
0596   //Large intestine ->ID 38
0597   d = 1.040 *g/cm3;
0598   fIntestine_lrg = new G4Material("intestine_lrg",d,9);
0599   fIntestine_lrg -> AddElement(elH,0.105);
0600   fIntestine_lrg -> AddElement(elC,0.113);
0601   fIntestine_lrg -> AddElement(elN,0.026);
0602   fIntestine_lrg -> AddElement(elO,0.750);
0603   fIntestine_lrg -> AddElement(elNa,0.001);
0604   fIntestine_lrg -> AddElement(elP,0.001);
0605   fIntestine_lrg -> AddElement(elS,0.001);
0606   fIntestine_lrg -> AddElement(elCl,0.002);
0607   fIntestine_lrg -> AddElement(elK,0.001);
0608
0609   //Spleen -> ID 39
0610   d = 1.040 *g/cm3;
0611   fSpleen = new G4Material("spleen",d,10);
0612   fSpleen -> AddElement(elH,0.102);
0613   fSpleen -> AddElement(elC,0.111);
0614   fSpleen -> AddElement(elN,0.033);
0615   fSpleen -> AddElement(elO,0.743);
0616   fSpleen -> AddElement(elNa,0.001);
0617   fSpleen -> AddElement(elP,0.002);
0618   fSpleen -> AddElement(elS,0.002);
0619   fSpleen -> AddElement(elCl,0.003);
0620   fSpleen -> AddElement(elK,0.002);
0621   fSpleen -> AddElement(elFe,0.001);
0622
0623   //Thyroid -> ID 40
0624   d = 1.040 *g/cm3;
0625   fThyroid = new G4Material("thyroid",d,10);
0626   fThyroid -> AddElement(elH,0.104);
0627   fThyroid -> AddElement(elC,0.117);
0628   fThyroid -> AddElement(elN,0.026);
0629   fThyroid -> AddElement(elO,0.745);
0630   fThyroid -> AddElement(elNa,0.002);
0631   fThyroid -> AddElement(elP,0.001);
0632   fThyroid -> AddElement(elS,0.001);
0633   fThyroid -> AddElement(elCl,0.002);
0634   fThyroid -> AddElement(elK,0.001);
0635   fThyroid -> AddElement(elI,0.001);
0636
0637   //Urinary Bladder -> ID 41
0638   d = 1.040 *g/cm3;
0639   fBladder = new G4Material("bladder",d,9);
0640   fBladder -> AddElement(elH,0.105);
0641   fBladder -> AddElement(elC,0.096);
0642   fBladder -> AddElement(elN,0.026);
0643   fBladder -> AddElement(elO,0.761);
0644   fBladder -> AddElement(elNa,0.002);
0645   fBladder -> AddElement(elP,0.002);
0646   fBladder -> AddElement(elS,0.002);
0647   fBladder -> AddElement(elCl,0.003);
0648   fBladder -> AddElement(elK,0.003);
0649
0650   //  Testes (Defined as ovaries_testes for visualisation purposes) -> ID 42
0651   d = 1.040 *g/cm3;
0652   fOvaries_testes = new G4Material("ovaries_testes",d,9);
0653   fOvaries_testes -> AddElement(elH,0.106);
0654   fOvaries_testes -> AddElement(elC,0.100);
0655   fOvaries_testes -> AddElement(elN,0.021);
0656   fOvaries_testes -> AddElement(elO,0.764);
0657   fOvaries_testes -> AddElement(elNa,0.002);
0658   fOvaries_testes -> AddElement(elP,0.001);
0659   fOvaries_testes -> AddElement(elS,0.002);
0660   fOvaries_testes -> AddElement(elCl,0.002);
0661   fOvaries_testes -> AddElement(elK,0.002);
0662
0663   //Adrenals -> ID 43
0664   d = 1.030 *g/cm3;
0665   fAdrenals = new G4Material("adrenals",d,9);
0666   fAdrenals -> AddElement(elH,0.104);
0667   fAdrenals -> AddElement(elC,0.221);
0668   fAdrenals -> AddElement(elN,0.028);
0669   fAdrenals -> AddElement(elO,0.637);
0670   fAdrenals -> AddElement(elNa,0.001);
0671   fAdrenals -> AddElement(elP,0.002);
0672   fAdrenals -> AddElement(elS,0.003);
0673   fAdrenals -> AddElement(elCl,0.002);
0674   fAdrenals -> AddElement(elK,0.002);
0675
0676   //Oesophagus -> ID 44
0677   d = 1.030 *g/cm3;
0678   fOesophagus = new G4Material("oesophagus",d,9);
0679   fOesophagus -> AddElement(elH,0.104);
0680   fOesophagus -> AddElement(elC,0.213);
0681   fOesophagus -> AddElement(elN,0.029);
0682   fOesophagus -> AddElement(elO,0.644);
0683   fOesophagus -> AddElement(elNa,0.001);
0684   fOesophagus -> AddElement(elP,0.002);
0685   fOesophagus -> AddElement(elS,0.003);
0686   fOesophagus -> AddElement(elCl,0.002);
0687   fOesophagus -> AddElement(elK,0.002);
0688
0689   //Miscillaneous (Gallbladder, Trachea, Thymus, Tonsils, Ureters, ...) -> ID 45
0690   d = 1.030 *g/cm3;
0691   fMisc = new G4Material("misc",d,9);
0692   fMisc -> AddElement(elH,0.104);
0693   fMisc -> AddElement(elC,0.231);
0694   fMisc -> AddElement(elN,0.028);
0695   fMisc -> AddElement(elO,0.627);
0696   fMisc -> AddElement(elNa,0.001);
0697   fMisc -> AddElement(elP,0.002);
0698   fMisc -> AddElement(elS,0.003);
0699   fMisc -> AddElement(elCl,0.002);
0700   fMisc -> AddElement(elK,0.002);
0701
0702   //Prostate (Defined as Uterus_Prostate for visualisation purposes) -> ID 46
0703   d = 1.030 *g/cm3;
0704   fUterus_prostate = new G4Material("uterus_prostate",d,9);
0705   fUterus_prostate -> AddElement(elH,0.104);
0706   fUterus_prostate -> AddElement(elC,0.231);
0707   fUterus_prostate -> AddElement(elN,0.028);
0708   fUterus_prostate -> AddElement(elO,0.627);
0709   fUterus_prostate -> AddElement(elNa,0.001);
0710   fUterus_prostate -> AddElement(elP,0.002);
0711   fUterus_prostate -> AddElement(elS,0.003);
0712   fUterus_prostate -> AddElement(elCl,0.002);
0713   fUterus_prostate -> AddElement(elK,0.002);
0714
0715   //Lymph -> ID 47
0716   d = 1.030 *g/cm3;
0717   fLymph = new G4Material("lymph",d,7);
0718   fLymph -> AddElement(elH,0.108);
0719   fLymph -> AddElement(elC,0.042);
0720   fLymph -> AddElement(elN,0.011);
0721   fLymph -> AddElement(elO,0.831);
0722   fLymph -> AddElement(elNa,0.003);
0723   fLymph -> AddElement(elS,0.001);
0724   fLymph -> AddElement(elCl,0.004);
0725
0726   //Breast (Mammary Gland)-> ID 48
0727   d = 1.020 *g/cm3;
0728   fBreast_glandular = new G4Material("breast_glandular",d,7);
0729   fBreast_glandular -> AddElement(elH,0.112);
0730   fBreast_glandular -> AddElement(elC,0.516);
0731   fBreast_glandular -> AddElement(elN,0.011);
0732   fBreast_glandular -> AddElement(elO,0.358);
0733   fBreast_glandular -> AddElement(elNa,0.001);
0734   fBreast_glandular -> AddElement(elS,0.001);
0735   fBreast_glandular -> AddElement(elCl,0.001);
0736
0737   //Adipose tissue (fBreast) -> ID 49
0738   d = 0.950 *g/cm3;
0739   fBreast_adipose = new G4Material("breast_adipose",d,7);
0740   fBreast_adipose -> AddElement(elH,0.114);
0741   fBreast_adipose -> AddElement(elC,0.588);
0742   fBreast_adipose -> AddElement(elN,0.008);
0743   fBreast_adipose -> AddElement(elO,0.287);
0744   fBreast_adipose -> AddElement(elNa,0.001);
0745   fBreast_adipose -> AddElement(elS,0.001);
0746   fBreast_adipose -> AddElement(elCl,0.001);
0747
0748   //Lung Tissue (Compressed Lung) -> ID 50
0749   d = 0.382 *g/cm3;
0750   fLung = new G4Material("lung",d,9);
0751   fLung -> AddElement(elH,0.103);
0752   fLung -> AddElement(elC,0.107);
0753   fLung -> AddElement(elN,0.032);
0754   fLung -> AddElement(elO,0.746);
0755   fLung -> AddElement(elNa,0.002);
0756   fLung -> AddElement(elP,0.002);
0757   fLung -> AddElement(elS,0.003);
0758   fLung -> AddElement(elCl,0.003);
0759   fLung -> AddElement(elK,0.002);
0760
0761   //Contents of Gastro-intestinal tract -> ID 51
0762   d = 1.040 *g/cm3;
0763   fGastro_content = new G4Material("gastro_content",d,10);
0764   fGastro_content -> AddElement(elH,0.100);
0765   fGastro_content -> AddElement(elC,0.222);
0766   fGastro_content -> AddElement(elN,0.022);
0767   fGastro_content -> AddElement(elO,0.644);
0768   fGastro_content -> AddElement(elNa,0.001);
0769   fGastro_content -> AddElement(elP,0.002);
0770   fGastro_content -> AddElement(elS,0.003);
0771   fGastro_content -> AddElement(elCl,0.001);
0772   fGastro_content -> AddElement(elK,0.004);
0773   fGastro_content -> AddElement(elCa,0.001);
0774
0775   //Urine -> ID 52
0776   d = 1.040 *g/cm3;
0777   fUrine = new G4Material("urine",d,7);
0778   fUrine -> AddElement(elH,0.107);
0779   fUrine -> AddElement(elC,0.003);
0780   fUrine -> AddElement(elN,0.010);
0781   fUrine -> AddElement(elO,0.873);
0782   fUrine -> AddElement(elNa,0.004);
0783   fUrine -> AddElement(elP,0.001);
0784   fUrine -> AddElement(elK,0.002);
0785
0786 }
0787
0788 G4Material* ICRP110PhantomMaterial_Male::GetMaterial(G4String material)
0789 {
0790   // Returns a material
0791   G4Material* pttoMaterial = G4Material::GetMaterial(material);
0792
0793   if (!pttoMaterial) G4cout << "WARNING: material '" << material << "' is not defined!" << G4endl;
0794
0795   return pttoMaterial;
0796 }