Environment variable "G4FORCE_RUN_MANAGER_TYPE" enabled with value == Serial. Forcing G4RunManager type... ############################################ !!! WARNING - FPE detection is activated !!! ############################################ ################################ !!! G4Backtrace is activated !!! ################################ ************************************************************** Geant4 version Name: geant4-11-04-patch-01 (13-March-2026) Copyright : Geant4 Collaboration References : NIM A 506 (2003), 250-303 : IEEE-TNS 53 (2006), 270-278 : NIM A 835 (2016), 186-225 WWW : http://geant4.org/ ************************************************************** Visualization Manager instantiating with verbosity "warnings (3)"... Visualization Manager initialising... Registering graphics systems... You have successfully registered the following graphics systems. Registered graphics systems are: ASCIITree (ATree) DAWNFILE (DAWNFILE) RayTracer (RT) VRML2FILE (VRML2FILE) gMocrenFile (gMocrenFile) TOOLSSG_OFFSCREEN (TSG_OFFSCREEN, TSG_FILE) OpenGLImmediateQt (OGLIQt, OGLI) OpenGLStoredQt (OGLSQt, OGLS) OpenGLImmediateXm (OGLIXm, OGLIQt_FALLBACK) OpenGLStoredXm (OGLSXm, OGLSQt_FALLBACK) OpenGLImmediateX (OGLIX, OGLIQt_FALLBACK, OGLIXm_FALLBACK) OpenGLStoredX (OGLSX, OGLSQt_FALLBACK, OGLSXm_FALLBACK) RayTracerX (RTX) RayTracerQt (RTQt) TOOLSSG_X11_GLES (TSG_X11_GLES, TSGX11, TSG_XT_GLES_FALLBACK) TOOLSSG_X11_ZB (TSG_X11_ZB, TSGX11ZB) TOOLSSG_XT_GLES (TSG_XT_GLES, TSGXt, TSG_QT_GLES_FALLBACK) TOOLSSG_XT_ZB (TSG_XT_ZB, TSGXtZB) TOOLSSG_QT_GLES (TSG_QT_GLES, TSGQt, TSG, OGL) TOOLSSG_QT_ZB (TSG_QT_ZB, TSGQtZB, TSGZB) You may choose a graphics system (driver) with a parameter of the command "/vis/open" or "/vis/sceneHandler/create", or you may omit the driver parameter and choose at run time: - by argument in the construction of G4VisExecutive - by environment variable "G4VIS_DEFAULT_DRIVER" - by entry in "~/.g4session" - by build flags. - Note: This feature is not allowed in batch mode. For further information see "examples/basic/B1/exampleB1.cc" and "vis.mac". Registering model factories... You have successfully registered the following model factories. Registered model factories: generic drawByAttribute drawByCharge drawByOriginVolume drawByParticleID drawByEncounteredVolume Registered models: None Registered filter factories: attributeFilter chargeFilter originVolumeFilter particleFilter encounteredVolumeFilter Registered filters: None You have successfully registered the following user vis actions. Run Duration User Vis Actions: none End of Event User Vis Actions: none End of Run User Vis Actions: none Some /vis commands (optionally) take a string to specify colour. "/vis/list" to see available colours. # # # # Initialize New Phantom # # Choose model: ORNLFemale, ORNLMale, MIRD, MIRDHead, ORNLHead # /phantom/setPhantomModel MIRD >> Phantom MIRD will be built. # # Choose Sex of Phantom: Male or Female # /phantom/setPhantomSex Male >> Male Phantom will be built. # Construct your Phantom # /phantom/buildNewPhantom ****************>>>> NEW PHANTOM CONSTRUCTION <<<<***************** /run/initialize MotherVolume: physicalWorld sensitivity : 0 Construct Head with mother physicalWorld Checking overlaps for volume physicalHead:0 (G4UnionSolid) ... OK! Volume of Head = 4654.79 cm^3 Material of Head = soft_tissue Density of Material = 0.9869 g/cm^3 Mass of Head = 4593.82 g MotherVolume: physicalHead sensitivity : 0 Construct Skull with mother volume physicalHead Checking overlaps for volume physicalSkull:0 (G4SubtractionSolid) ... OK! Skull created !!!!!! Volume of Skull = 847.724 cm^3 Material of Skull = skeleton Density of Material = 1.4862 g/cm^3 Mass of Skull = 1259.89 g Construct Brain with mother physicalHead Checking overlaps for volume physicalBrain:0 (G4Ellipsoid) ... OK! Volume of Brain = 1470.27 cm^3 Material of Brain = soft_tissue Density of Material = 0.9869 g/cm^3 Mass of Brain = 1451 g MotherVolume: physicalWorld sensitivity : 0 Construct Trunk with mother volume physicalWorld Checking overlaps for volume physicalTrunk:0 (G4EllipticalTube) ... OK! Trunk created !!!!!! Volume of Trunk = 43982.3 cm^3 Material of Trunk = soft_tissue Density of Material = 0.9869 g/cm^3 Mass of Trunk = 43406.1 g MotherVolume: physicalWorld sensitivity : 0 Construct LeftLeg with mother volume physicalWorld Checking overlaps for volume physicalLeftLeg:0 (G4Cons) ... OK! LeftLeg created !!!!!! Volume of LeftLeg = 10388.2 cm^3 Material of LeftLeg = soft_tissue Density of Material = 0.9869 g/cm^3 Mass of LeftLeg = 10252.1 g MotherVolume: physicalWorld sensitivity : 0 Construct RightLeg with mother volume physicalWorld Checking overlaps for volume physicalRightLeg:0 (G4Cons) ... OK! RightLeg created !!!!!! Volume of RightLeg = 10388.2 cm^3 Material of RightLeg = soft_tissue Density of Material = 0.9869 g/cm^3 Mass of RightLeg = 10252.1 g MotherVolume: physicalTrunk sensitivity : 0 Construct LeftArmBone with mother physicalTrunk Checking overlaps for volume physicalLeftArmBone:0 (G4EllipticalTube) ... OK! LeftArmBone created !!!!!! Volume of LeftArmBone = 819.39 cm^3 Material of LeftArmBone = skeleton Density of Material = 1.4862 g/cm^3 Mass of LeftArmBone = 1217.78 g MotherVolume: physicalTrunk sensitivity : 0 Construct RightArmBone with mother volume physicalTrunk Checking overlaps for volume physicalRightArmBone:0 (G4EllipticalTube) ... OK! RightArmBone created !!!!!! Volume of RightArmBone = 819.39 cm^3 Material of RightArmBone = skeleton Density of Material = 1.4862 g/cm^3 Mass of RightArmBone = 1217.78 g MotherVolume: physicalLeftLeg sensitivity : 0 Construct LeftLegBone with mother volume physicalLeftLeg Checking overlaps for volume physicalLeftLegBone:0 (G4Cons) ... OK! LeftLegBone created !!!!!! Volume of LeftLegBone = 1399.74 cm^3 Material of LeftLegBone = skeleton Density of Material = 1.4862 g/cm^3 Mass of LeftLegBone = 2080.29 g MotherVolume: physicalRightLeg sensitivity : 0 Construct RightLegBone with mother volume physicalRightLeg Checking overlaps for volume physicalRightLegBone:0 (G4Cons) ... OK! RightLegBone created !!!!!! Volume of RightLegBone = 1399.74 cm^3 Material of RightLegBone = skeleton Density of Material = 1.4862 g/cm^3 Mass of RightLegBone = 2080.29 g MotherVolume: physicalHead sensitivity : 0 Construct UpperSpine with mother volume physicalHead Checking overlaps for volume physicalUpperSpine:0 (G4SubtractionSolid) ... OK! UpperSpine created !!!!!! Volume of UpperSpine = 126.465 cm^3 Material of UpperSpine = skeleton Density of Material = 1.4862 g/cm^3 Mass of UpperSpine = 187.952 g MotherVolume: physicalTrunk sensitivity : 0 Construct LeftScapula with mother volume physicalTrunk Checking overlaps for volume physicalLeftScapula:0 (G4SubtractionSolid) ... OK! LeftScapula created !!!!!! Volume of LeftScapula = 104.415 cm^3 Material of LeftScapula = skeleton Density of Material = 1.4862 g/cm^3 Mass of LeftScapula = 155.181 g MotherVolume: physicalTrunk sensitivity : 0 Construct RightScapula with mother volume physicalTrunk Checking overlaps for volume physicalRightScapula:0 (G4SubtractionSolid) ... OK! RightScapula created !!!!!! Volume of RightScapula = 106.924 cm^3 Material of RightScapula = skeleton Density of Material = 1.4862 g/cm^3 Mass of RightScapula = 158.91 g Construct LeftAdrenal with mother physicalTrunk Checking overlaps for volume physicalLeftAdrenal:0 (G4Ellipsoid) ... OK! Left LeftAdrenal created !!!!!! Volume of LeftAdrenal = 7.85398 cm^3 Material of LeftAdrenal = soft_tissue Density of Material = 0.9869 g/cm^3 Mass of LeftAdrenal = 7.75109 g Construct RightAdrenal with mother volume physicalTrunk Checking overlaps for volume physicalRightAdrenal:0 (G4Ellipsoid) ... OK! Right RightAdrenal created !!!!!! Volume of RightAdrenal = 7.85398 cm^3 Material of RightAdrenal = soft_tissue Density of Material = 0.9869 g/cm^3 Mass of RightAdrenal = 7.75109 g Construct Thymus with mother volume physicalTrunk Checking overlaps for volume physicalThymus:0 (G4Ellipsoid) ... OK! Thymus created !!!!!! Volume of Thymus = 25.1327 cm^3 Material of Thymus = soft_tissue Density of Material = 0.9869 g/cm^3 Mass of Thymus = 24.8035 g MotherVolume: physicalTrunk sensitivity : 0 Construct LeftClavicle with mother volume physicalTrunk Checking overlaps for volume physicalLeftClavicle:0 (G4Torus) ... OK! LeftClavicle created !!!!!! Volume of LeftClavicle = 13.6657 cm^3 Material of LeftClavicle = skeleton Density of Material = 1.4862 g/cm^3 Mass of LeftClavicle = 20.3099 g MotherVolume: physicalTrunk sensitivity : 0 Construct RightClavicle with mother volume physicalTrunk Checking overlaps for volume physicalRightClavicle:0 (G4Torus) ... OK! RightClavicle created !!!!!! Volume of RightClavicle = 13.6657 cm^3 Material of RightClavicle = skeleton Density of Material = 1.4862 g/cm^3 Mass of RightClavicle = 20.3099 g Construct SmallIntestine with mother volume physicalTrunk Checking overlaps for volume physicalSmallIntestine:0 (G4SubtractionSolid) ... OK! SmallIntestine created !!!!!! Volume of SmallIntestine = 1019.24 cm^3 Material of SmallIntestine = soft_tissue Density of Material = 0.9869 g/cm^3 Mass of SmallIntestine = 1005.89 g Construct RibCage with mother volume physicalTrunk Checking overlaps for volume physicalRibCage:0 (G4SubtractionSolid) ... OK! Checking overlaps for volume physicalRib:0 (G4SubtractionSolid) ... OK! Checking overlaps for volume physicalRib:0 (G4SubtractionSolid) ... OK! Checking overlaps for volume physicalRib:0 (G4SubtractionSolid) ... OK! Checking overlaps for volume physicalRib:0 (G4SubtractionSolid) ... OK! Checking overlaps for volume physicalRib:0 (G4SubtractionSolid) ... OK! Checking overlaps for volume physicalRib:0 (G4SubtractionSolid) ... OK! Checking overlaps for volume physicalRib:0 (G4SubtractionSolid) ... OK! Checking overlaps for volume physicalRib:0 (G4SubtractionSolid) ... OK! Checking overlaps for volume physicalRib:0 (G4SubtractionSolid) ... OK! Checking overlaps for volume physicalRib:0 (G4SubtractionSolid) ... OK! Checking overlaps for volume physicalRib:0 (G4SubtractionSolid) ... OK! Checking overlaps for volume physicalRib:0 (G4SubtractionSolid) ... OK! RibCage created !!!!!! Volume of RibCage = 697.036 cm^3 Material of RibCage = skeleton Density of Material = 1.4862 g/cm^3 Mass of RibCage = 1035.94 g MotherVolume: physicalTrunk sensitivity : 0 Construct MiddleLowerSpine with mother volume physicalTrunk Checking overlaps for volume physicalMiddleLowerSpine:0 (G4EllipticalTube) ... OK! MiddleLowerSpine created !!!!!! Volume of MiddleLowerSpine = 753.982 cm^3 Material of MiddleLowerSpine = skeleton Density of Material = 1.4862 g/cm^3 Mass of MiddleLowerSpine = 1120.57 g Construct Pelvis with mother volume physicalTrunk Checking overlaps for volume physicalPelvis:0 (G4SubtractionSolid) ... OK! Pelvis created !!!!!! Volume of Pelvis = 610.442 cm^3 Material of Pelvis = skeleton Density of Material = 1.4862 g/cm^3 Mass of Pelvis = 907.239 g Construct Stomach with mother volume physicalTrunk Checking overlaps for volume physicalStomach:0 (G4Ellipsoid) ... OK! Stomach created !!!!!! Volume of Stomach = 402.124 cm^3 Material of Stomach = soft_tissue Density of Material = 0.9869 g/cm^3 Mass of Stomach = 396.856 g Construct UpperLargeIntestine with mother volume physicalTrunk Checking overlaps for volume physicalUpperLargeIntestine:0 (G4UnionSolid) ... OK! UpperLargeIntestine created !!!!!! Volume of UpperLargeIntestine = 434.914 cm^3 Material of UpperLargeIntestine = soft_tissue Density of Material = 0.9869 g/cm^3 Mass of UpperLargeIntestine = 429.217 g Construct LowerLargeIntestine with mother volume physicalTrunk Checking overlaps for volume physicalLowerLargeIntestine:0 (G4UnionSolid) ... OK! LowerLargeIntestine created !!!!!! Volume of LowerLargeIntestine = 344.316 cm^3 Material of LowerLargeIntestine = soft_tissue Density of Material = 0.9869 g/cm^3 Mass of LowerLargeIntestine = 339.805 g Construct Spleen with mother volume physicalTrunk Checking overlaps for volume physicalSpleen:0 (G4Ellipsoid) ... OK! Spleen created !!!!!! Volume of Spleen = 175.929 cm^3 Material of Spleen = soft_tissue Density of Material = 0.9869 g/cm^3 Mass of Spleen = 173.625 g Construct Pancreas with mother volume physicalTrunk Checking overlaps for volume physicalPancreas:0 (G4SubtractionSolid) ... OK! Pancreas created !!!!!! Volume of Pancreas = 61.0986 cm^3 Material of Pancreas = soft_tissue Density of Material = 0.9869 g/cm^3 Mass of Pancreas = 60.2982 g Construct LeftKidney with mother volume physicalTrunk Checking overlaps for volume physicalLeftKidney:0 (G4SubtractionSolid) ... OK! Left LeftKidney created !!!!!! Volume of LeftKidney = 143.998 cm^3 Material of LeftKidney = soft_tissue Density of Material = 0.9869 g/cm^3 Mass of LeftKidney = 142.112 g Construct RightKidney with mother volume physicalTrunk Checking overlaps for volume physicalRightKidney:0 (G4SubtractionSolid) ... OK! RightKidney created !!!!!! Volume of RightKidney = 143.996 cm^3 Material of RightKidney = soft_tissue Density of Material = 0.9869 g/cm^3 Mass of RightKidney = 142.109 g Construct UrinaryBladder with mother volume physicalTrunk Checking overlaps for volume physicalUrinaryBladder:0 (G4SubtractionSolid) ... OK! UrinaryBladder created !!!!!! Volume of UrinaryBladder = 45.8796 cm^3 Material of UrinaryBladder = soft_tissue Density of Material = 0.9869 g/cm^3 Mass of UrinaryBladder = 45.2785 g MotherVolume: physicalWorld sensitivity : 0 Construct MaleGenitalia with mother volume physicalWorld Checking overlaps for volume physicalMaleGenitalia:0 (G4SubtractionSolid) ... OK! MaleGenitalia created !!!!!! Volume of MaleGenitalia = 228.823 cm^3 Material of MaleGenitalia = soft_tissue Density of Material = 0.9869 g/cm^3 Mass of MaleGenitalia = 225.826 g MotherVolume: physicalWorld sensitivity : 0 Construct LeftTeste with mother volume physicalMaleGenitalia Checking overlaps for volume physicalLeftTeste:0 (G4Ellipsoid) ... OK! LeftTeste created !!!!!! Volume of LeftTeste = 18.7867 cm^3 Material of LeftTeste = soft_tissue Density of Material = 0.9869 g/cm^3 Mass of LeftTeste = 18.5406 g MotherVolume: physicalWorld sensitivity : 0 Construct RightTeste with mother volume physicalMaleGenitalia Checking overlaps for volume physicalRightTeste:0 (G4Ellipsoid) ... OK! RightTeste created !!!!!! Volume of RightTeste = 18.7867 cm^3 Material of RightTeste = soft_tissue Density of Material = 0.9869 g/cm^3 Mass of RightTeste = 18.5406 g # or ... #/control/execute adultMIRDFemale.mac # or ... #/control/execute adultORNLMale.mac # or... #/control/execute adultORNLFemale.mac # # Primary radiation field /control/execute primary.mac # ### Radiation field defined ### by means of GPS /gps/verbose 0 /gps/particle e- /gps/energy 1 MeV # /gps/pos/type Volume # /gps/pos/shape Para /gps/pos/halfx 20 cm /gps/pos/halfy 1 mm /gps/pos/halfz 100 cm /gps/pos/centre 0. 30. 0. cm /gps/direction 0 -1 0 # /run/beamOn 100 ======================================================================= ====== Electromagnetic Physics Parameters ======== ======================================================================= LPM effect enabled 1 Enable creation and use of sampling tables 0 Apply cuts on all EM processes 0 Use combined TransportationWithMsc Disabled Use general process 1 Enable linear polarisation for gamma 0 Enable photoeffect sampling below K-shell 1 Enable sampling of quantum entanglement 0 X-section factor for integral approach 0.8 Min kinetic energy for tables 100 eV Max kinetic energy for tables 100 TeV Number of bins per decade of a table 20 Verbose level 1 Verbose level for worker thread 0 Bremsstrahlung energy threshold above which primary e+- is added to the list of secondary 100 TeV Bremsstrahlung energy threshold above which primary muon/hadron is added to the list of secondary 100 TeV Positron annihilation at rest model AllisonPositronium Enable 3 gamma annihilation on fly 1 Lowest triplet kinetic energy 1 MeV Enable sampling of gamma linear polarisation 0 5D gamma conversion model type 0 5D gamma conversion model on isolated ion 0 Use RiGe 5D e+e- pair production model by muons 0 Livermore data directory epics_2017 ======================================================================= ====== Ionisation Parameters ======== ======================================================================= Step function for e+- (0.2, 0.01 mm) Step function for muons/hadrons (0.1, 0.05 mm) Step function for light ions (0.1, 0.02 mm) Step function for general ions (0.1, 0.001 mm) Lowest e+e- kinetic energy 100 eV Lowest muon/hadron kinetic energy 1 keV Use ICRU90 data 1 Fluctuations of dE/dx are enabled 1 Type of fluctuation model for leptons and hadrons Urban Use built-in Birks saturation 0 Build CSDA range enabled 0 Use cut as a final range enabled 0 Enable angular generator interface 1 Max kinetic energy for CSDA tables 1 GeV Max kinetic energy for NIEL computation 1 MeV Linear loss limit 0.01 Read data from file for e+e- pair production by mu 0 ======================================================================= ====== Multiple Scattering Parameters ======== ======================================================================= Type of msc step limit algorithm for e+- 2 Type of msc step limit algorithm for muons/hadrons 0 Msc lateral displacement for e+- enabled 1 Msc lateral displacement for muons and hadrons 1 Urban msc model lateral displacement alg96 1 Range factor for msc step limit for e+- 0.08 Range factor for msc step limit for muons/hadrons 0.2 Geometry factor for msc step limitation of e+- 2.5 Safety factor for msc step limit for e+- 0.6 Skin parameter for msc step limitation of e+- 3 Lambda limit for msc step limit for e+- 1 mm Use Mott correction for e- scattering 1 Factor used for dynamic computation of angular limit between single and multiple scattering 1 Fixed angular limit between single and multiple scattering 3.1416 rad Upper energy limit for e+- multiple scattering 100 MeV Type of electron single scattering model 0 Type of nuclear form-factor 1 Screening factor 1 ======================================================================= ====== Atomic Deexcitation Parameters ======== ======================================================================= Fluorescence enabled 1 Directory in G4LEDATA for fluorescence data files fluor Auger electron cascade enabled 0 PIXE atomic de-excitation enabled 0 De-excitation module ignores cuts 0 Type of PIXE cross section for hadrons Empirical Type of PIXE cross section for e+- Livermore ======================================================================= ### === Deexcitation model UAtomDeexcitation is activated for 1 region: DefaultRegionForTheWorld 1 0 0 ### === Ignore cuts flag: 0 phot: for gamma SubType=12 BuildTable=0 LambdaPrime table from 200 keV to 100 TeV in 174 bins ===== EM models for the G4Region DefaultRegionForTheWorld ====== LivermorePhElectric : Emin= 0 eV Emax= 100 TeV SauterGavrila Fluo compt: for gamma SubType=13 BuildTable=1 Lambda table from 100 eV to 1 MeV, 20 bins/decade, spline: 1 LambdaPrime table from 1 MeV to 100 TeV in 160 bins ===== EM models for the G4Region DefaultRegionForTheWorld ====== LowEPComptonModel : Emin= 0 eV Emax= 20 MeV Fluo KleinNishina : Emin= 20 MeV Emax= 100 TeV Fluo conv: for gamma SubType=14 BuildTable=1 Lambda table from 1.022 MeV to 100 TeV, 20 bins/decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== BetheHeitler5D : Emin= 0 eV Emax= 100 TeV ModifiedTsai Rayl: for gamma SubType=11 BuildTable=1 Lambda table from 100 eV to 150 keV, 20 bins/decade, spline: 0 LambdaPrime table from 150 keV to 100 TeV in 176 bins ===== EM models for the G4Region DefaultRegionForTheWorld ====== LivermoreRayleigh : Emin= 0 eV Emax= 100 TeV CullenGenerator msc: for e- SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== GoudsmitSaunderson : Emin= 0 eV Emax= 100 MeV Nbins=120 100 eV - 100 MeV StepLim=SafetyPlus Rfact=0.08 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm WentzelVIUni : Emin= 100 MeV Emax= 100 TeV Nbins=120 100 MeV - 100 TeV StepLim=SafetyPlus Rfact=0.08 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm eIoni: for e- XStype:3 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 StepFunction=(0.2, 0.01 mm), integ: 3, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== PenIoni : Emin= 0 eV Emax= 100 keV MollerBhabha : Emin= 100 keV Emax= 100 TeV deltaVI eBrem: for e- XStype:4 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 LPM flag: 1 for E > 1 GeV, VertexHighEnergyTh(GeV)= 100000 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eBremSB : Emin= 0 eV Emax= 1 GeV AngularGen2BS eBremLPM : Emin= 1 GeV Emax= 100 TeV AngularGen2BS ePairProd: for e- XStype:1 SubType=4 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 0 Sampling table 25x1001 from 0.1 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== ePairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi CoulombScat: for e- XStype:1 SubType=1 BuildTable=1 Lambda table from 100 MeV to 100 TeV, 20 bins/decade, spline: 0 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 100 MeV Emax= 100 TeV msc: for e+ SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== GoudsmitSaunderson : Emin= 0 eV Emax= 100 MeV Nbins=120 100 eV - 100 MeV StepLim=SafetyPlus Rfact=0.08 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm WentzelVIUni : Emin= 100 MeV Emax= 100 TeV Nbins=120 100 MeV - 100 TeV StepLim=SafetyPlus Rfact=0.08 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm eIoni: for e+ XStype:3 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 StepFunction=(0.2, 0.01 mm), integ: 3, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== PenIoni : Emin= 0 eV Emax= 100 keV MollerBhabha : Emin= 100 keV Emax= 100 TeV deltaVI eBrem: for e+ XStype:4 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 LPM flag: 1 for E > 1 GeV, VertexHighEnergyTh(GeV)= 100000 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eBremSB : Emin= 0 eV Emax= 1 GeV AngularGen2BS eBremLPM : Emin= 1 GeV Emax= 100 TeV AngularGen2BS ePairProd: for e+ XStype:1 SubType=4 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 0 Sampling table 25x1001 from 0.1 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== ePairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi annihil: for e+ XStype:2 SubType=5 AtRestModel:Allison BuildTable=0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eplusTo2or3gamma : Emin= 0 eV Emax= 100 TeV CoulombScat: for e+ XStype:1 SubType=1 BuildTable=1 Lambda table from 100 MeV to 100 TeV, 20 bins/decade, spline: 0 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 100 MeV Emax= 100 TeV msc: for proton SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm hIoni: for proton XStype:3 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 2 MeV deltaVI BetheBloch : Emin= 2 MeV Emax= 100 TeV deltaVI hBrems: for proton XStype:1 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi ===== Limit on energy threshold has been applied hPairProd: for proton XStype:1 SubType=4 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 Sampling table 17x1001, from 7.50618 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi CoulombScat: for proton XStype:1 SubType=1 BuildTable=1 Lambda table from threshold to 100 TeV, 20 bins/decade, spline: 0 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 100 TeV nuclearStopping: for proton SubType=8 BuildTable=0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU49NucStopping : Emin= 0 eV Emax= 1 MeV msc: for GenericIon SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm ionIoni: for GenericIon XStype:3 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 StepFunction=(0.1, 0.001 mm), integ: 3, fluct: 1, linLossLim= 0.02 ===== EM models for the G4Region DefaultRegionForTheWorld ====== LindhardSorensen : Emin= 0 eV Emax= 100 TeV deltaVI nuclearStopping: for GenericIon SubType=8 BuildTable=0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU49NucStopping : Emin= 0 eV Emax= 1 MeV msc: for alpha SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm ionIoni: for alpha XStype:3 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 StepFunction=(0.1, 0.02 mm), integ: 3, fluct: 1, linLossLim= 0.02 ===== EM models for the G4Region DefaultRegionForTheWorld ====== BraggIon : Emin= 0 eV Emax=7.9452 MeV deltaVI BetheBloch : Emin=7.9452 MeV Emax= 100 TeV deltaVI nuclearStopping: for alpha SubType=8 BuildTable=0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU49NucStopping : Emin= 0 eV Emax= 1 MeV msc: for anti_proton SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm hIoni: for anti_proton XStype:3 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax= 2 MeV deltaVI BetheBloch : Emin= 2 MeV Emax= 100 TeV deltaVI hBrems: for anti_proton XStype:1 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi ===== Limit on energy threshold has been applied hPairProd: for anti_proton XStype:1 SubType=4 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 Sampling table 17x1001, from 7.50618 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi CoulombScat: for anti_proton XStype:1 SubType=1 BuildTable=1 Lambda table from threshold to 100 TeV, 20 bins/decade, spline: 0 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 100 TeV msc: for kaon+ SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm hIoni: for kaon+ XStype:3 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax=1.05231 MeV deltaVI BetheBloch : Emin=1.05231 MeV Emax= 100 TeV deltaVI hBrems: for kaon+ XStype:1 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi ===== Limit on energy threshold has been applied hPairProd: for kaon+ XStype:1 SubType=4 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 Sampling table 18x1001, from 3.94942 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi CoulombScat: for kaon+ XStype:1 SubType=1 BuildTable=1 Lambda table from threshold to 100 TeV, 20 bins/decade, spline: 0 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 100 TeV msc: for kaon- SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm hIoni: for kaon- XStype:3 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax=1.05231 MeV deltaVI BetheBloch : Emin=1.05231 MeV Emax= 100 TeV deltaVI hBrems: for kaon- XStype:1 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi ===== Limit on energy threshold has been applied hPairProd: for kaon- XStype:1 SubType=4 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 Sampling table 18x1001, from 3.94942 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi CoulombScat: for kaon- XStype:1 SubType=1 BuildTable=1 Used Lambda table of kaon+ ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 100 TeV msc: for mu+ SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm muIoni: for mu+ XStype:3 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 200 keV deltaVI MuBetheBloch : Emin= 200 keV Emax= 100 TeV deltaVI muBrems: for mu+ XStype:1 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MuBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi ===== Limit on energy threshold has been applied muPairProd: for mu+ XStype:1 SubType=4 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 Sampling table 21x1001, from 0.85 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== muPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi CoulombScat: for mu+ XStype:1 SubType=1 BuildTable=1 Lambda table from threshold to 100 TeV, 20 bins/decade, spline: 0 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 100 TeV msc: for mu- SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm muIoni: for mu- XStype:3 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax= 200 keV deltaVI MuBetheBloch : Emin= 200 keV Emax= 100 TeV deltaVI muBrems: for mu- XStype:1 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MuBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi ===== Limit on energy threshold has been applied muPairProd: for mu- XStype:1 SubType=4 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 Sampling table 21x1001, from 0.85 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== muPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi CoulombScat: for mu- XStype:1 SubType=1 BuildTable=1 Used Lambda table of mu+ ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 100 TeV msc: for pi+ SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm hIoni: for pi+ XStype:3 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax=297.505 keV deltaVI BetheBloch : Emin=297.505 keV Emax= 100 TeV deltaVI hBrems: for pi+ XStype:1 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi ===== Limit on energy threshold has been applied hPairProd: for pi+ XStype:1 SubType=4 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 Sampling table 20x1001, from 1.11656 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi CoulombScat: for pi+ XStype:1 SubType=1 BuildTable=1 Lambda table from threshold to 100 TeV, 20 bins/decade, spline: 0 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 100 TeV msc: for pi- SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm hIoni: for pi- XStype:3 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax=297.505 keV deltaVI BetheBloch : Emin=297.505 keV Emax= 100 TeV deltaVI hBrems: for pi- XStype:1 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi ===== Limit on energy threshold has been applied hPairProd: for pi- XStype:1 SubType=4 dE/dx and range tables from 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 Sampling table 20x1001, from 1.11656 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi CoulombScat: for pi- XStype:1 SubType=1 BuildTable=1 Used Lambda table of pi+ ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 100 TeV ========= Table of registered couples ============================ Index : 0 used in the geometry : Yes Material : Air Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm Energy thresholds : gamma 250 eV e- 329.56 eV e+ 325.472 eV proton 100 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 1 used in the geometry : Yes Material : soft_tissue Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm Energy thresholds : gamma 2.80001 keV e- 347.402 keV e+ 337.896 keV proton 100 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 2 used in the geometry : Yes Material : skeleton Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm Energy thresholds : gamma 4.00973 keV e- 441.639 keV e+ 427.712 keV proton 100 keV Region(s) which use this couple : DefaultRegionForTheWorld ================================================================== ### Run 0 start. ... set ntuple merging mode : G4NtupleMergeMode::kNone - done ... set ntuple merging mode : G4NtupleMergeMode::kNone - done ... set ntuple merging row mode : row-wise - done ... create directory for ntuples : human_phantom_ntuple - done ... create file : human_phantom.root - done ... open analysis file : human_phantom.root - done ... open analysis file : human_phantom.root - done ... create ntuple booking : 1 ntupleId 1 - done ... create ntuple T column : organID ntupleId 1 - done ... create ntuple T column : edep ntupleId 1 - done Number of events = 100 Energy Total in Run:logicalBrain, ID: 0, Energy Deposition (MeV): 0 Energy Total in Run:logicalHead, ID: 1, Energy Deposition (MeV): 4.30138 Analysis: organ 1 edep: 4.30138 Energy Total in Run:logicalHeart, ID: 2, Energy Deposition (MeV): 0 Energy Total in Run:logicalLeftAdrenal, ID: 3, Energy Deposition (MeV): 0 Energy Total in Run:logicalLeftArmBone, ID: 4, Energy Deposition (MeV): 0 Energy Total in Run:logicalLeftBreast, ID: 5, Energy Deposition (MeV): 0 Energy Total in Run:logicalLeftClavicle, ID: 6, Energy Deposition (MeV): 0 Energy Total in Run:logicalLeftKidney, ID: 7, Energy Deposition (MeV): 0 Energy Total in Run:logicalLeftLeg, ID: 8, Energy Deposition (MeV): 11.6454 Analysis: organ 8 edep: 11.6454 Energy Total in Run:logicalLeftLegBone, ID: 9, Energy Deposition (MeV): 0 Energy Total in Run:logicalLeftLung, ID: 10, Energy Deposition (MeV): 0 Energy Total in Run:logicalLeftOvary, ID: 11, Energy Deposition (MeV): 0 Energy Total in Run:logicalLeftScapula, ID: 12, Energy Deposition (MeV): 0 Energy Total in Run:logicalLeftTeste, ID: 13, Energy Deposition (MeV): 0 Energy Total in Run:logicalLowerLargeIntestine, ID: 14, Energy Deposition (MeV): 0 Energy Total in Run:logicalMaleGenitalia, ID: 15, Energy Deposition (MeV): 2.08718 Analysis: organ 15 edep: 2.08718 Energy Total in Run:logicalMiddleLowerSpine, ID: 16, Energy Deposition (MeV): 0 Energy Total in Run:logicalPancreas, ID: 17, Energy Deposition (MeV): 0 Energy Total in Run:logicalPelvis, ID: 18, Energy Deposition (MeV): 0 Energy Total in Run:logicalRibCage, ID: 19, Energy Deposition (MeV): 0.935689 Analysis: organ 19 edep: 0.935689 Energy Total in Run:logicalRightAdrenal, ID: 20, Energy Deposition (MeV): 0 Energy Total in Run:logicalRightArmBone, ID: 21, Energy Deposition (MeV): 0 Energy Total in Run:logicalRightBreast, ID: 22, Energy Deposition (MeV): 0 Energy Total in Run:logicalRightClavicle, ID: 23, Energy Deposition (MeV): 0 Energy Total in Run:logicalRightKidney, ID: 24, Energy Deposition (MeV): 0 Energy Total in Run:logicalRightLeg, ID: 25, Energy Deposition (MeV): 7.5709 Analysis: organ 25 edep: 7.5709 Energy Total in Run:logicalRightLegBone, ID: 26, Energy Deposition (MeV): 0 Energy Total in Run:logicalRightLung, ID: 27, Energy Deposition (MeV): 0 Energy Total in Run:logicalRightOvary, ID: 28, Energy Deposition (MeV): 0 Energy Total in Run:logicalRightScapula, ID: 29, Energy Deposition (MeV): 0 Energy Total in Run:logicalRightTeste, ID: 30, Energy Deposition (MeV): 0 Energy Total in Run:logicalSkull, ID: 31, Energy Deposition (MeV): 0 Energy Total in Run:logicalSmallIntestine, ID: 32, Energy Deposition (MeV): 0 Energy Total in Run:logicalSpleen, ID: 33, Energy Deposition (MeV): 0 Energy Total in Run:logicalStomach, ID: 34, Energy Deposition (MeV): 0 Energy Total in Run:logicalThymus, ID: 35, Energy Deposition (MeV): 0 Energy Total in Run:logicalThyroid, ID: 36, Energy Deposition (MeV): 0 Energy Total in Run:logicalTrunk, ID: 37, Energy Deposition (MeV): 28.8373 Analysis: organ 37 edep: 28.8373 Energy Total in Run:logicalUpperLargeIntestine, ID: 38, Energy Deposition (MeV): 0 Energy Total in Run:logicalUpperSpine, ID: 39, Energy Deposition (MeV): 0 Energy Total in Run:logicalUrinaryBladder, ID: 40, Energy Deposition (MeV): 0 Energy Total in Run:logicalUterus, ID: 41, Energy Deposition (MeV): 0 Total Energy deposit in the body is: 55.3779 MeV ... write file : human_phantom.root - done ... close file : human_phantom.root - done ... clear files - done ... clear h1 - done ... clear h2 - done ... clear h3 - done ... clear p1 - done ... clear p2 - done ... clear ntupleBookings - done ... clear ntuples - done ... clear files - done ... clear all data - done Graphics systems deleted. Visualization Manager deleting...