ICRP110_HumanPhantoms/src/ICRP110PhantomMaterial_Female.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 //
0030 #include "ICRP110PhantomMaterial_Female.hh"
0031 #include "globals.hh"
0032 #include "G4SystemOfUnits.hh"
0033 #include "G4MaterialPropertiesTable.hh"
0034 #include "G4MaterialPropertyVector.hh"
0035 #include "G4MaterialTable.hh"
0036 #include "Randomize.hh"
0037 #include "G4RunManager.hh"
0038 #include "G4Element.hh"
0039 #include "G4ElementTable.hh"
0040
0041 ICRP110PhantomMaterial_Female::ICRP110PhantomMaterial_Female():
0042   fLung(nullptr), fTeeth(nullptr), fBone(nullptr), fHumeri_upper(nullptr), fHumeri_lower(nullptr),
0043   fArm_lower(nullptr), fHand(nullptr), fClavicle(nullptr), fCranium(nullptr), fFemora_upper(nullptr),
0044   fFemora_lower(nullptr), fLeg_lower(nullptr), fFoot(nullptr), fMandible(nullptr),
0045   fPelvis(nullptr), fRibs(nullptr), fScapulae(nullptr), fSpine_cervical(nullptr),
0046   fSpine_lumbar(nullptr), fSpine_thoratic(nullptr), fSacrum(nullptr),
0047   fSternum(nullptr), fHf_upper(nullptr), fHf_lower(nullptr), fMed_lowerleg(nullptr),
0048   fMed_lowerarm(nullptr), fCartilage(nullptr), fSkin(nullptr), fBlood(nullptr),
0049   fMuscle(nullptr), fLiver(nullptr), fPancreas(nullptr), fBrain(nullptr), fHeart(nullptr), fEye(nullptr),
0050   fKidney(nullptr), fStomach(nullptr), fIntestine_sml(nullptr), fIntestine_lrg(nullptr),
0051   fSpleen(nullptr), fThyroid(nullptr), fBladder(nullptr), fOvaries_testes(nullptr), fAdrenals(nullptr),
0052   fOesophagus(nullptr), fMisc(nullptr), fUterus_prostate(nullptr), fLymph(nullptr),
0053   fBreast_glandular(nullptr), fBreast_adipose(nullptr), fGastro_content(nullptr),
0054   fUrine(nullptr)
0055 {;}
0056
0057 ICRP110PhantomMaterial_Female::~ICRP110PhantomMaterial_Female()
0058 {;}
0059
0060 void ICRP110PhantomMaterial_Female::DefineMaterials()
0061 {
0062   // Define required materials
0063
0064   G4double A;  // atomic mass
0065   G4double Z;  // atomic number
0066   G4double d;  // density
0067
0068   // General elements
0069
0070   A = 1.01*g/mole;
0071   auto elH = new G4Element ("Hydrogen","H",Z = 1.,A);
0072
0073   A = 12.011*g/mole;
0074   auto elC = new G4Element("Carbon","C",Z = 6.,A);
0075
0076   A = 14.01*g/mole;
0077   auto elN = new G4Element("Nitrogen","N",Z = 7.,A);
0078
0079   A = 16.00*g/mole;
0080   auto elO = new G4Element("Oxygen","O",Z = 8.,A);
0081
0082   A = 22.99*g/mole;
0083   auto elNa = new G4Element("Sodium","Na",Z = 11.,A);
0084
0085   A = 24.305*g/mole;
0086   auto elMg = new G4Element("Magnesium","Mg",Z = 12.,A);
0087
0088   A = 30.974*g/mole;
0089   auto elP = new G4Element("Phosphorus","P",Z = 15.,A);
0090
0091   A = 32.064*g/mole;
0092   auto elS = new G4Element("Sulfur","S",Z = 16.,A);
0093
0094   A = 35.453*g/mole;
0095   auto elCl = new G4Element("Chlorine","Cl",Z = 17.,A);
0096
0097   A = 39.098*g/mole;
0098   auto elK = new G4Element("Potassium","K",Z = 19.,A);
0099
0100   A = 40.08*g/mole;
0101   auto elCa = new G4Element("Calcium","Ca",Z = 20.,A);
0102
0103   A = 55.85*g/mole;
0104   auto elFe  = new G4Element("Iron","Fe",Z = 26.,A);
0105
0106   A = 126.90447 *g/mole;
0107   auto elI = new G4Element("Iodine","I", Z = 53.,A);
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.185 *g/cm3;
0136   fHumeri_upper = new G4Material("humeri_upper",d,11);
0137   fHumeri_upper -> AddElement(elH,0.087) ;
0138   fHumeri_upper -> AddElement(elC,0.366);
0139   fHumeri_upper -> AddElement(elN,0.025);
0140   fHumeri_upper -> AddElement(elO,0.422);
0141   fHumeri_upper -> AddElement(elNa,0.002);
0142   fHumeri_upper -> AddElement(elMg,0.001);
0143   fHumeri_upper -> AddElement(elP,0.030);
0144   fHumeri_upper -> AddElement(elS,0.003);
0145   fHumeri_upper -> AddElement(elCl,0.001);
0146   fHumeri_upper -> AddElement(elK, 0.001);
0147   fHumeri_upper -> AddElement(elCa,0.062);
0148
0149   //Humeri, lower half, spongiosa -> ID 4
0150   d = 1.117 *g/cm3;
0151   fHumeri_lower = new G4Material("humeri_lower",d,9);
0152   fHumeri_lower -> AddElement(elH,0.096);
0153   fHumeri_lower -> AddElement(elC,0.473);
0154   fHumeri_lower -> AddElement(elN,0.017);
0155   fHumeri_lower -> AddElement(elO,0.341);
0156   fHumeri_lower -> AddElement(elNa,0.002);
0157   fHumeri_lower -> AddElement(elP,0.022);
0158   fHumeri_lower -> AddElement(elS,0.002);
0159   fHumeri_lower -> AddElement(elCl,0.001);
0160   fHumeri_lower -> AddElement(elCa,0.046);
0161
0162   //Lower arm bones, spongiosa -> ID 5
0163   d = 1.117 *g/cm3;
0164   fArm_lower = new G4Material("arm_lower",d,9);
0165   fArm_lower -> AddElement(elH,0.096);
0166   fArm_lower -> AddElement(elC,0.473);
0167   fArm_lower -> AddElement(elN,0.017);
0168   fArm_lower -> AddElement(elO,0.341);
0169   fArm_lower -> AddElement(elNa,0.002);
0170   fArm_lower -> AddElement(elP,0.022);
0171   fArm_lower -> AddElement(elS,0.002);
0172   fArm_lower -> AddElement(elCl,0.001);
0173   fArm_lower -> AddElement(elCa,0.046);
0174
0175   //Hand Bones, Spongiosa ->ID 6
0176   d = 1.117 *g/cm3;
0177   fHand = new G4Material("hand",d,9);
0178   fHand -> AddElement(elH,0.096);
0179   fHand -> AddElement(elC,0.473);
0180   fHand -> AddElement(elN,0.017);
0181   fHand -> AddElement(elO,0.341);
0182   fHand -> AddElement(elNa,0.002);
0183   fHand -> AddElement(elP,0.022);
0184   fHand -> AddElement(elS,0.002);
0185   fHand -> AddElement(elCl,0.001);
0186   fHand -> AddElement(elCa,0.046);
0187
0188   //Clavicles, spongiosa -> ID 7
0189   d = 1.191 *g/cm3;
0190   fClavicle = new G4Material("clavicle",d,11);
0191   fClavicle -> AddElement(elH,0.087);
0192   fClavicle -> AddElement(elC,0.361);
0193   fClavicle -> AddElement(elN,0.025);
0194   fClavicle -> AddElement(elO,0.424);
0195   fClavicle -> AddElement(elNa,0.002);
0196   fClavicle -> AddElement(elMg,0.001);
0197   fClavicle -> AddElement(elP,0.031);
0198   fClavicle -> AddElement(elS,0.003);
0199   fClavicle -> AddElement(elCl,0.001);
0200   fClavicle -> AddElement(elK,0.001);
0201   fClavicle -> AddElement(elCa,0.064);
0202
0203   // Cranium, spongiosa -> ID 8
0204   d = 1.245 *g/cm3;
0205   fCranium = new G4Material("cranium",d,11);
0206   fCranium -> AddElement(elH,0.081);
0207   fCranium -> AddElement(elC,0.317);
0208   fCranium -> AddElement(elN,0.028);
0209   fCranium -> AddElement(elO,0.451);
0210   fCranium -> AddElement(elNa,0.002);
0211   fCranium -> AddElement(elMg, 0.001);
0212   fCranium -> AddElement(elP,0.037);
0213   fCranium -> AddElement(elS,0.003);
0214   fCranium -> AddElement(elCl,0.001);
0215   fCranium -> AddElement(elK,0.001);
0216   fCranium -> AddElement(elCa,0.078);
0217
0218   //Femora, upper half, spongiosa -> ID 9
0219   d = 1.046 *g/cm3;
0220   fFemora_upper = new G4Material("femora_upper",d,10);
0221   fFemora_upper -> AddElement(elH,0.104);
0222   fFemora_upper -> AddElement(elC,0.496);
0223   fFemora_upper -> AddElement(elN,0.018);
0224   fFemora_upper -> AddElement(elO,0.349);
0225   fFemora_upper -> AddElement(elNa,0.001);
0226   fFemora_upper -> AddElement(elP,0.009);
0227   fFemora_upper -> AddElement(elS,0.002);
0228   fFemora_upper -> AddElement(elCl,0.001);
0229   fFemora_upper -> AddElement(elK,0.001);
0230   fFemora_upper -> AddElement(elCa,0.019);
0231
0232   //Femora, lower half, spongiosa -> ID 10
0233   d = 1.117 *g/cm3;
0234   fFemora_lower = new G4Material("femora_lower",d,9);
0235   fFemora_lower -> AddElement(elH,0.096);
0236   fFemora_lower -> AddElement(elC,0.473);
0237   fFemora_lower -> AddElement(elN,0.017);
0238   fFemora_lower -> AddElement(elO,0.341);
0239   fFemora_lower -> AddElement(elNa,0.002);
0240   fFemora_lower -> AddElement(elP,0.022);
0241   fFemora_lower -> AddElement(elS,0.002);
0242   fFemora_lower -> AddElement(elCl,0.001);
0243   fFemora_lower -> AddElement(elCa,0.046);
0244
0245   //Lower leg bones, spongiosa -> ID 11
0246   d = 1.117 *g/cm3;
0247   fLeg_lower = new G4Material("leg_lower",d,9);
0248   fLeg_lower -> AddElement(elH,0.096);
0249   fLeg_lower -> AddElement(elC,0.473);
0250   fLeg_lower -> AddElement(elN,0.017);
0251   fLeg_lower -> AddElement(elO,0.341);
0252   fLeg_lower -> AddElement(elNa,0.002);
0253   fLeg_lower -> AddElement(elP,0.022);
0254   fLeg_lower -> AddElement(elS,0.002);
0255   fLeg_lower -> AddElement(elCl,0.001);
0256   fLeg_lower -> AddElement(elCa,0.046);
0257
0258   //Foot bones, spongiosa ->ID 12
0259   d = 1.117 *g/cm3;
0260   fFoot = new G4Material("foot",d,9);
0261   fFoot -> AddElement(elH,0.096);
0262   fFoot -> AddElement(elC,0.473);
0263   fFoot -> AddElement(elN,0.017);
0264   fFoot -> AddElement(elO,0.341);
0265   fFoot -> AddElement(elNa,0.002);
0266   fFoot -> AddElement(elP,0.022);
0267   fFoot -> AddElement(elS,0.002);
0268   fFoot -> AddElement(elCl,0.001);
0269   fFoot -> AddElement(elCa,0.046);
0270
0271   //Mandible, spongiosa -> ID 13
0272   d = 1.189 *g/cm3;
0273   fMandible = new G4Material("mandible",d,11);
0274   fMandible -> AddElement(elH,0.087);
0275   fMandible -> AddElement(elC,0.357);
0276   fMandible -> AddElement(elN,0.026);
0277   fMandible -> AddElement(elO,0.429);
0278   fMandible -> AddElement(elNa,0.002);
0279   fMandible -> AddElement(elMg,0.001);
0280   fMandible -> AddElement(elP,0.030);
0281   fMandible -> AddElement(elS,0.003);
0282   fMandible -> AddElement(elCl,0.001);
0283   fMandible -> AddElement(elK,0.001);
0284   fMandible -> AddElement(elCa,0.063);
0285
0286   //Pelvis, Spongiosa -> ID 14
0287   d = 1.109 *g/cm3;
0288   fPelvis = new G4Material("pelvis",d,10);
0289   fPelvis -> AddElement(elH,0.096);
0290   fPelvis -> AddElement(elC,0.406);
0291   fPelvis -> AddElement(elN,0.025);
0292   fPelvis -> AddElement(elO,0.412);
0293   fPelvis -> AddElement(elNa,0.001);
0294   fPelvis -> AddElement(elP,0.018);
0295   fPelvis -> AddElement(elS,0.002);
0296   fPelvis -> AddElement(elCl,0.001);
0297   fPelvis -> AddElement(elK,0.001);
0298   fPelvis -> AddElement(elCa,0.038);
0299
0300   //Ribs, spongiosa -> ID 15
0301   d = 1.092 *g/cm3;
0302   fRibs = new G4Material("ribs",d,11);
0303   fRibs -> AddElement(elH,0.097);
0304   fRibs -> AddElement(elC,0.381);
0305   fRibs -> AddElement(elN,0.028);
0306   fRibs -> AddElement(elO,0.445);
0307   fRibs -> AddElement(elNa,0.001);
0308   fRibs -> AddElement(elP,0.014);
0309   fRibs -> AddElement(elS,0.002);
0310   fRibs -> AddElement(elCl,0.002);
0311   fRibs -> AddElement(elK,0.001);
0312   fRibs -> AddElement(elCa,0.028);
0313   fRibs -> AddElement(elFe,0.001);
0314
0315   //Scapulae, spongiosa -> ID 16
0316   d = 1.128 *g/cm3;
0317   fScapulae = new G4Material("scapulae",d,10);
0318   fScapulae -> AddElement(elH,0.094);
0319   fScapulae -> AddElement(elC,0.406);
0320   fScapulae -> AddElement(elN,0.024);
0321   fScapulae -> AddElement(elO,0.404);
0322   fScapulae -> AddElement(elNa,0.001);
0323   fScapulae -> AddElement(elP,0.022);
0324   fScapulae -> AddElement(elS,0.002);
0325   fScapulae -> AddElement(elCl,0.001);
0326   fScapulae -> AddElement(elK,0.001);
0327   fScapulae -> AddElement(elCa,0.045);
0328
0329   //Cervical Spine, spongiosa -> ID 17
0330   d = 1.135 *g/cm3;
0331   fSpine_cervical = new G4Material("spine_cervical",d,10);
0332   fSpine_cervical -> AddElement(elH,0.092);
0333   fSpine_cervical -> AddElement(elC,0.351);
0334   fSpine_cervical -> AddElement(elN,0.029);
0335   fSpine_cervical -> AddElement(elO,0.458);
0336   fSpine_cervical -> AddElement(elNa,0.001);
0337   fSpine_cervical -> AddElement(elP,0.021);
0338   fSpine_cervical -> AddElement(elS,0.002);
0339   fSpine_cervical -> AddElement(elCl,0.002);
0340   fSpine_cervical -> AddElement(elK,0.002);
0341   fSpine_cervical -> AddElement(elCa,0.043);
0342
0343   //Thoratic Spine, spongiosa -> ID 18
0344   d = 1.084 *g/cm3;
0345   fSpine_thoratic = new G4Material("spine_thoratic",d,11);
0346   fSpine_thoratic -> AddElement(elH,0.098);
0347   fSpine_thoratic -> AddElement(elC,0.386);
0348   fSpine_thoratic -> AddElement(elN,0.028);
0349   fSpine_thoratic -> AddElement(elO,0.442);
0350   fSpine_thoratic -> AddElement(elNa,0.001);
0351   fSpine_thoratic -> AddElement(elP,0.013);
0352   fSpine_thoratic -> AddElement(elS,0.002);
0353   fSpine_thoratic -> AddElement(elCl,0.002);
0354   fSpine_thoratic -> AddElement(elK,0.001);
0355   fSpine_thoratic -> AddElement(elCa,0.026);
0356   fSpine_thoratic -> AddElement(elFe,0.001);
0357
0358   //Lumbar Spine, spongiosa -> ID 19
0359   d = 1.171 *g/cm3;
0360   fSpine_lumbar = new G4Material("spine_lumbar",d,11);
0361   fSpine_lumbar -> AddElement(elH,0.088);
0362   fSpine_lumbar -> AddElement(elC,0.329);
0363   fSpine_lumbar -> AddElement(elN,0.030);
0364   fSpine_lumbar -> AddElement(elO,0.466);
0365   fSpine_lumbar -> AddElement(elNa,0.001);
0366   fSpine_lumbar -> AddElement(elMg,0.001);
0367   fSpine_lumbar -> AddElement(elP,0.026);
0368   fSpine_lumbar -> AddElement(elS,0.003);
0369   fSpine_lumbar -> AddElement(elCl,0.001);
0370   fSpine_lumbar -> AddElement(elK,0.001);
0371   fSpine_lumbar -> AddElement(elCa,0.054);
0372
0373   //Sacrum, spongiosa -> ID 20
0374   d = 1.052 *g/cm3;
0375   fSacrum = new G4Material("sacrum",d,11);
0376   fSacrum -> AddElement(elH,0.102);
0377   fSacrum -> AddElement(elC,0.410);
0378   fSacrum -> AddElement(elN,0.027);
0379   fSacrum -> AddElement(elO,0.433);
0380   fSacrum -> AddElement(elNa,0.001);
0381   fSacrum -> AddElement(elP,0.007);
0382   fSacrum -> AddElement(elS,0.002);
0383   fSacrum -> AddElement(elCl,0.002);
0384   fSacrum -> AddElement(elK,0.001);
0385   fSacrum -> AddElement(elCa,0.014);
0386   fSacrum -> AddElement(elFe,0.001);
0387
0388   //Sternum, spongiosa -> ID 21
0389   d = 1.076 *g/cm3;
0390   fSternum = new G4Material("sternum",d,11);
0391   fSternum -> AddElement(elH,0.099);
0392   fSternum -> AddElement(elC,0.392);
0393   fSternum -> AddElement(elN,0.028);
0394   fSternum -> AddElement(elO,0.439);
0395   fSternum -> AddElement(elNa,0.001);
0396   fSternum -> AddElement(elP,0.012);
0397   fSternum -> AddElement(elS,0.002);
0398   fSternum -> AddElement(elCl,0.002);
0399   fSternum -> AddElement(elK,0.001);
0400   fSternum -> AddElement(elCa,0.023);
0401   fSternum -> AddElement(elFe,0.001);
0402
0403   //Humeri and femora, upper halves, medullary cavity -> ID 22
0404   d = 0.980 *g/cm3;
0405   fHf_upper = new G4Material("hf_upper",d,7);
0406   fHf_upper -> AddElement(elH,0.115);
0407   fHf_upper -> AddElement(elC,0.637);
0408   fHf_upper -> AddElement(elN,0.007);
0409   fHf_upper -> AddElement(elO,0.238);
0410   fHf_upper -> AddElement(elNa,0.001);
0411   fHf_upper -> AddElement(elS,0.001);
0412   fHf_upper -> AddElement(elCl,0.001);
0413
0414   //Humeri and femora, lower halves, medullary cavity -> ID 23
0415   d = 0.980 *g/cm3;
0416   fHf_lower = new G4Material("hf_lower",d,7);
0417   fHf_lower -> AddElement(elH,0.115);
0418   fHf_lower -> AddElement(elC,0.637);
0419   fHf_lower -> AddElement(elN,0.007);
0420   fHf_lower -> AddElement(elO,0.238);
0421   fHf_lower -> AddElement(elNa,0.001);
0422   fHf_lower -> AddElement(elS,0.001);
0423   fHf_lower -> AddElement(elCl,0.001);
0424
0425   //Lower arm bones, medullary cavity -> ID 24
0426   d = 0.980 *g/cm3;
0427   fMed_lowerarm = new G4Material("med_lowerarm",d,7);
0428   fMed_lowerarm -> AddElement(elH,0.115);
0429   fMed_lowerarm -> AddElement(elC,0.637);
0430   fMed_lowerarm -> AddElement(elN,0.007);
0431   fMed_lowerarm -> AddElement(elO,0.238);
0432   fMed_lowerarm -> AddElement(elNa,0.001);
0433   fMed_lowerarm -> AddElement(elS,0.001);
0434   fMed_lowerarm -> AddElement(elCl,0.001);
0435
0436   //Lower leg bones, medullary cavity -> ID 25
0437   d = 0.980 *g/cm3;
0438   fMed_lowerleg = new G4Material("med_lowerleg",d,7);
0439   fMed_lowerleg -> AddElement(elH,0.115);
0440   fMed_lowerleg -> AddElement(elC,0.637);
0441   fMed_lowerleg -> AddElement(elN,0.007);
0442   fMed_lowerleg -> AddElement(elO,0.238);
0443   fMed_lowerleg -> AddElement(elNa,0.001);
0444   fMed_lowerleg -> AddElement(elS,0.001);
0445   fMed_lowerleg -> AddElement(elCl,0.001);
0446
0447   //Cartilage -> ID 26
0448   d = 1.100 *g/cm3;
0449   fCartilage = new G4Material("cartilage",d,8);
0450   fCartilage -> AddElement(elH,0.096);
0451   fCartilage -> AddElement(elC,0.099);
0452   fCartilage -> AddElement(elN,0.022);
0453   fCartilage -> AddElement(elO,0.744);
0454   fCartilage -> AddElement(elNa,0.005);
0455   fCartilage -> AddElement(elP,0.022);
0456   fCartilage -> AddElement(elS,0.009);
0457   fCartilage -> AddElement(elCl,0.003);
0458
0459   //Skin -> Id 27
0460   d = 1.090 *g/cm3;
0461   fSkin = new G4Material("skin",d,9);
0462   fSkin -> AddElement(elH,0.100);
0463   fSkin -> AddElement(elC,0.199);
0464   fSkin -> AddElement(elN,0.042);
0465   fSkin -> AddElement(elO,0.650);
0466   fSkin -> AddElement(elNa,0.002);
0467   fSkin -> AddElement(elP,0.001);
0468   fSkin -> AddElement(elS,0.002);
0469   fSkin -> AddElement(elCl,0.003);
0470   fSkin -> AddElement(elK,0.001);
0471
0472   //Blood -> ID 28
0473   d = 1.060 *g/cm3;
0474   fBlood = new G4Material("blood",d,10);
0475   fBlood -> AddElement(elH,0.102);
0476   fBlood -> AddElement(elC,0.110);
0477   fBlood -> AddElement(elN,0.033);
0478   fBlood -> AddElement(elO,0.745);
0479   fBlood -> AddElement(elNa,0.001);
0480   fBlood -> AddElement(elP,0.001);
0481   fBlood -> AddElement(elS,0.002);
0482   fBlood -> AddElement(elCl,0.003);
0483   fBlood -> AddElement(elK,0.002);
0484   fBlood -> AddElement(elFe,0.001);
0485
0486   //Muscular Tissue -> ID 29
0487   d = 1.050 *g/cm3;
0488   fMuscle = new G4Material("muscle",d,9);
0489   fMuscle -> AddElement(elH,0.102);
0490   fMuscle -> AddElement(elC,0.142);
0491   fMuscle -> AddElement(elN,0.034);
0492   fMuscle -> AddElement(elO,0.711);
0493   fMuscle -> AddElement(elNa,0.001);
0494   fMuscle -> AddElement(elP,0.002);
0495   fMuscle -> AddElement(elS,0.003);
0496   fMuscle -> AddElement(elCl,0.001);
0497   fMuscle -> AddElement(elK,0.004);
0498
0499   //Liver -> ID 30
0500   d = 1.050 *g/cm3;
0501   fLiver = new G4Material("liver",d,9);
0502   fLiver -> AddElement(elH,0.102);
0503   fLiver -> AddElement(elC,0.131);
0504   fLiver -> AddElement(elN,0.031);
0505   fLiver -> AddElement(elO,0.724);
0506   fLiver -> AddElement(elNa,0.002);
0507   fLiver -> AddElement(elP,0.002);
0508   fLiver -> AddElement(elS,0.003);
0509   fLiver -> AddElement(elCl,0.002);
0510   fLiver -> AddElement(elK,0.003);
0511
0512   //Pancreas ->ID 31
0513   d = 1.050 *g/cm3;
0514   fPancreas = new G4Material("pancreas",d,9);
0515   fPancreas -> AddElement(elH,0.105);
0516   fPancreas -> AddElement(elC,0.157);
0517   fPancreas -> AddElement(elN,0.024);
0518   fPancreas -> AddElement(elO,0.705);
0519   fPancreas -> AddElement(elNa,0.002);
0520   fPancreas -> AddElement(elP,0.002);
0521   fPancreas -> AddElement(elS,0.001);
0522   fPancreas -> AddElement(elCl,0.002);
0523   fPancreas -> AddElement(elK,0.002);
0524
0525   //Brain -> ID 32
0526   d = 1.050 *g/cm3;
0527   fBrain = new G4Material("brain",d,9);
0528   fBrain -> AddElement(elH,0.107);
0529   fBrain -> AddElement(elC,0.144);
0530   fBrain -> AddElement(elN,0.022);
0531   fBrain -> AddElement(elO,0.713);
0532   fBrain -> AddElement(elNa,0.002);
0533   fBrain -> AddElement(elP,0.004);
0534   fBrain -> AddElement(elS,0.002);
0535   fBrain -> AddElement(elCl,0.003);
0536   fBrain -> AddElement(elK,0.003);
0537
0538   //Heart -> ID 33
0539   d = 1.050 *g/cm3;
0540   fHeart = new G4Material("heart",d,9);
0541   fHeart -> AddElement(elH,0.104);
0542   fHeart -> AddElement(elC,0.138);
0543   fHeart -> AddElement(elN,0.029);
0544   fHeart -> AddElement(elO,0.719);
0545   fHeart -> AddElement(elNa,0.001);
0546   fHeart -> AddElement(elP,0.002);
0547   fHeart -> AddElement(elS,0.002);
0548   fHeart -> AddElement(elCl,0.002);
0549   fHeart -> AddElement(elK,0.003);
0550
0551   //Eye ->ID 34
0552   d = 1.050 *g/cm3;
0553   fEye = new G4Material("eye",d,8);
0554   fEye -> AddElement(elH,0.097);
0555   fEye -> AddElement(elC,0.183);
0556   fEye -> AddElement(elN,0.054);
0557   fEye -> AddElement(elO,0.660);
0558   fEye -> AddElement(elCa,0.001);
0559   fEye -> AddElement(elP,0.001);
0560   fEye -> AddElement(elS,0.003);
0561   fEye -> AddElement(elCl,0.001);
0562
0563   //Kidneys -> ID 35
0564   d = 1.050 *g/cm3;
0565   fKidney = new G4Material("kidney",d,10);
0566   fKidney -> AddElement(elH,0.103);
0567   fKidney -> AddElement(elC,0.125);
0568   fKidney -> AddElement(elN,0.031);
0569   fKidney -> AddElement(elO,0.730);
0570   fKidney -> AddElement(elNa,0.002);
0571   fKidney -> AddElement(elP,0.002);
0572   fKidney -> AddElement(elS,0.002);
0573   fKidney -> AddElement(elCl,0.002);
0574   fKidney -> AddElement(elK,0.002);
0575   fKidney -> AddElement(elCa,0.001);
0576
0577   //Stomach ->ID 36
0578   d = 1.040 *g/cm3;
0579   fStomach = new G4Material("stomach",d,9);
0580   fStomach -> AddElement(elH,0.105);
0581   fStomach -> AddElement(elC,0.114);
0582   fStomach -> AddElement(elN,0.025);
0583   fStomach -> AddElement(elO,0.750);
0584   fStomach -> AddElement(elNa,0.001);
0585   fStomach -> AddElement(elP,0.001);
0586   fStomach -> AddElement(elS,0.001);
0587   fStomach -> AddElement(elCl,0.002);
0588   fStomach -> AddElement(elK,0.001);
0589
0590   //Small intestine ->ID 37
0591   d = 1.040 *g/cm3;
0592   fIntestine_sml = new G4Material("intestine_sml",d,9);
0593   fIntestine_sml -> AddElement(elH,0.105);
0594   fIntestine_sml -> AddElement(elC,0.114);
0595   fIntestine_sml -> AddElement(elN,0.025);
0596   fIntestine_sml -> AddElement(elO,0.750);
0597   fIntestine_sml -> AddElement(elNa,0.001);
0598   fIntestine_sml -> AddElement(elP,0.001);
0599   fIntestine_sml -> AddElement(elS,0.001);
0600   fIntestine_sml -> AddElement(elCl,0.002);
0601   fIntestine_sml -> AddElement(elK,0.001);
0602
0603   //Large intestine ->ID 38
0604   d = 1.040 *g/cm3;
0605   fIntestine_lrg = new G4Material("intestine_lrg",d,9);
0606   fIntestine_lrg -> AddElement(elH,0.105);
0607   fIntestine_lrg -> AddElement(elC,0.114);
0608   fIntestine_lrg -> AddElement(elN,0.025);
0609   fIntestine_lrg -> AddElement(elO,0.750);
0610   fIntestine_lrg -> AddElement(elNa,0.001);
0611   fIntestine_lrg -> AddElement(elP,0.001);
0612   fIntestine_lrg -> AddElement(elS,0.001);
0613   fIntestine_lrg -> AddElement(elCl,0.002);
0614   fIntestine_lrg -> AddElement(elK,0.001);
0615
0616   //Spleen -> ID 39
0617   d = 1.040 *g/cm3;
0618   fSpleen = new G4Material("spleen",d,9);
0619   fSpleen -> AddElement(elH,0.103);
0620   fSpleen -> AddElement(elC,0.112);
0621   fSpleen -> AddElement(elN,0.032);
0622   fSpleen -> AddElement(elO,0.743);
0623   fSpleen -> AddElement(elNa,0.001);
0624   fSpleen -> AddElement(elP,0.002);
0625   fSpleen -> AddElement(elS,0.002);
0626   fSpleen -> AddElement(elCl,0.002);
0627   fSpleen -> AddElement(elK,0.003);
0628
0629   //Thyroid -> ID 40
0630   d = 1.040 *g/cm3;
0631   fThyroid = new G4Material("thyroid",d,10);
0632   fThyroid -> AddElement(elH,0.104);
0633   fThyroid -> AddElement(elC,0.118);
0634   fThyroid -> AddElement(elN,0.025);
0635   fThyroid -> AddElement(elO,0.745);
0636   fThyroid -> AddElement(elNa,0.002);
0637   fThyroid -> AddElement(elP,0.001);
0638   fThyroid -> AddElement(elS,0.001);
0639   fThyroid -> AddElement(elCl,0.002);
0640   fThyroid -> AddElement(elK,0.001);
0641   fThyroid -> AddElement(elI,0.001);
0642
0643   //Urinary Bladder -> ID 41
0644   d = 1.040 *g/cm3;
0645   fBladder = new G4Material("bladder",d,9);
0646   fBladder -> AddElement(elH,0.105);
0647   fBladder -> AddElement(elC,0.096);
0648   fBladder -> AddElement(elN,0.026);
0649   fBladder -> AddElement(elO,0.761);
0650   fBladder -> AddElement(elNa,0.002);
0651   fBladder -> AddElement(elP,0.002);
0652   fBladder -> AddElement(elS,0.002);
0653   fBladder -> AddElement(elCl,0.003);
0654   fBladder -> AddElement(elK,0.003);
0655
0656   //  Ovaries (Defined as ovaries_testes to be called upon by ColourMap.dat
0657   //for visualisation purposes) -> ID 42
0658   d = 1.040 *g/cm3;
0659   fOvaries_testes = new G4Material("ovaries_testes",d,9);
0660   fOvaries_testes -> AddElement(elH,0.105);
0661   fOvaries_testes -> AddElement(elC,0.094);
0662   fOvaries_testes -> AddElement(elN,0.025);
0663   fOvaries_testes -> AddElement(elO,0.766);
0664   fOvaries_testes -> AddElement(elNa,0.002);
0665   fOvaries_testes -> AddElement(elP,0.002);
0666   fOvaries_testes -> AddElement(elS,0.002);
0667   fOvaries_testes -> AddElement(elCl,0.002);
0668   fOvaries_testes -> AddElement(elK,0.002);
0669
0670   //Adrenals -> ID 43
0671   d = 1.030 *g/cm3;
0672   fAdrenals = new G4Material("adrenals",d,9);
0673   fAdrenals -> AddElement(elH,0.104);
0674   fAdrenals -> AddElement(elC,0.228);
0675   fAdrenals -> AddElement(elN,0.028);
0676   fAdrenals -> AddElement(elO,0.630);
0677   fAdrenals -> AddElement(elNa,0.001);
0678   fAdrenals -> AddElement(elP,0.002);
0679   fAdrenals -> AddElement(elS,0.003);
0680   fAdrenals -> AddElement(elCl,0.002);
0681   fAdrenals -> AddElement(elK,0.002);
0682
0683   //Oesophagus -> ID 44
0684   d = 1.030 *g/cm3;
0685   fOesophagus = new G4Material("oesophagus",d,9);
0686   fOesophagus -> AddElement(elH,0.104);
0687   fOesophagus -> AddElement(elC,0.222);
0688   fOesophagus -> AddElement(elN,0.028);
0689   fOesophagus -> AddElement(elO,0.636);
0690   fOesophagus -> AddElement(elNa,0.001);
0691   fOesophagus -> AddElement(elP,0.002);
0692   fOesophagus -> AddElement(elS,0.003);
0693   fOesophagus -> AddElement(elCl,0.002);
0694   fOesophagus -> AddElement(elK,0.002);
0695
0696   //Miscillaneous (Gallbladder, Trachea, Thymus, Tonsils, Ureters, ...) -> ID 45
0697   d = 1.030 *g/cm3;
0698   fMisc = new G4Material("misc",d,9);
0699   fMisc -> AddElement(elH,0.105);
0700   fMisc -> AddElement(elC,0.235);
0701   fMisc -> AddElement(elN,0.028);
0702   fMisc -> AddElement(elO,0.622);
0703   fMisc -> AddElement(elNa,0.001);
0704   fMisc -> AddElement(elP,0.002);
0705   fMisc -> AddElement(elS,0.003);
0706   fMisc -> AddElement(elCl,0.002);
0707   fMisc -> AddElement(elK,0.002);
0708
0709   //Uterus (Defined as fUterus_prostate to be called upon by ColourMap.dat
0710   // for visualisation purposes) -> ID 46
0711   d = 1.030 *g/cm3;
0712   fUterus_prostate = new G4Material("uterus_prostate",d,9);
0713   fUterus_prostate -> AddElement(elH,0.105);
0714   fUterus_prostate -> AddElement(elC,0.286);
0715   fUterus_prostate -> AddElement(elN,0.025);
0716   fUterus_prostate -> AddElement(elO,0.576);
0717   fUterus_prostate -> AddElement(elNa,0.001);
0718   fUterus_prostate -> AddElement(elP,0.002);
0719   fUterus_prostate -> AddElement(elS,0.002);
0720   fUterus_prostate -> AddElement(elCl,0.001);
0721   fUterus_prostate -> AddElement(elK,0.002);
0722
0723   //Lymph -> ID 47
0724   d = 1.030 *g/cm3;
0725   fLymph = new G4Material("lymph",d,7);
0726   fLymph -> AddElement(elH,0.108);
0727   fLymph -> AddElement(elC,0.042);
0728   fLymph -> AddElement(elN,0.011);
0729   fLymph -> AddElement(elO,0.831);
0730   fLymph -> AddElement(elNa,0.003);
0731   fLymph -> AddElement(elS,0.001);
0732   fLymph -> AddElement(elCl,0.004);
0733
0734   //Breast (Mammary Gland)-> ID 48
0735   d = 1.020 *g/cm3;
0736   fBreast_glandular = new G4Material("breast_glandular",d,4);
0737   fBreast_glandular -> AddElement(elH,0.114);
0738   fBreast_glandular -> AddElement(elC,0.461);
0739   fBreast_glandular -> AddElement(elN,0.005);
0740   fBreast_glandular -> AddElement(elO,0.420);
0741
0742   //Adipose tissue (breast) -> ID 49
0743   d = 0.950 *g/cm3;
0744   fBreast_adipose = new G4Material("breast_adipose",d,7);
0745   fBreast_adipose -> AddElement(elH,0.114);
0746   fBreast_adipose -> AddElement(elC,0.589);
0747   fBreast_adipose -> AddElement(elN,0.007);
0748   fBreast_adipose -> AddElement(elO,0.287);
0749   fBreast_adipose -> AddElement(elNa,0.001);
0750   fBreast_adipose -> AddElement(elS,0.001);
0751   fBreast_adipose -> AddElement(elCl,0.001);
0752
0753   //Lung Tissue (Compressed Lung) -> ID 50
0754   d = 0.385 *g/cm3;
0755   fLung = new G4Material("lung",d,9);
0756   fLung -> AddElement(elH,0.103);
0757   fLung -> AddElement(elC,0.107);
0758   fLung -> AddElement(elN,0.032);
0759   fLung -> AddElement(elO,0.746);
0760   fLung -> AddElement(elNa,0.002);
0761   fLung -> AddElement(elP,0.002);
0762   fLung -> AddElement(elS,0.003);
0763   fLung -> AddElement(elCl,0.003);
0764   fLung -> AddElement(elK,0.002);
0765
0766   //Contents of gastro-intestinal tract -> ID 51
0767   d = 1.040 *g/cm3;
0768   fGastro_content = new G4Material("gastro_content",d,10);
0769   fGastro_content -> AddElement(elH,0.100);
0770   fGastro_content -> AddElement(elC,0.222);
0771   fGastro_content -> AddElement(elN,0.022);
0772   fGastro_content -> AddElement(elO,0.644);
0773   fGastro_content -> AddElement(elNa,0.001);
0774   fGastro_content -> AddElement(elP,0.002);
0775   fGastro_content -> AddElement(elS,0.003);
0776   fGastro_content -> AddElement(elCl,0.001);
0777   fGastro_content -> AddElement(elK,0.004);
0778   fGastro_content -> AddElement(elCa,0.001);
0779
0780   //Urine -> ID 52
0781   d = 1.040 *g/cm3;
0782   fUrine = new G4Material("urine",d,7);
0783   fUrine -> AddElement(elH,0.107);
0784   fUrine -> AddElement(elC,0.003);
0785   fUrine -> AddElement(elN,0.010);
0786   fUrine -> AddElement(elO,0.873);
0787   fUrine -> AddElement(elNa,0.004);
0788   fUrine -> AddElement(elP,0.001);
0789   fUrine -> AddElement(elK,0.002);
0790
0791 }
0792
0793 G4Material* ICRP110PhantomMaterial_Female::GetMaterial(G4String material)
0794 {
0795   // Returns a material
0796   G4Material* pttoMaterial = G4Material::GetMaterial(material);
0797   if (!pttoMaterial) G4cout << "WARNING: material '" << material << "' is not defined!" << G4endl;
0798   return pttoMaterial;
0799 }