/control/execute exercise5a.in /gamos/setParam GmReadPhantomGeometry:FileName ../exercise3/world.geom /gamos/setParam GmSimplePhantomGeometry:MaterialNames G4_WATER /gamos/setParam GmSimplePhantomGeometry:MaterialZVoxels 100 /gamos/setParam GmSimplePhantomGeometry:MaterialDensities 1. /gamos/setParam GmSimplePhantomGeometry:NVoxels 100 100 100 /gamos/setParam GmSimplePhantomGeometry:PhantomDims -50 50 -50 50 -50 50 /gamos/geometry GmSimplePhantomGeometry /gamos/physicsList GmEMExtendedPhysics /gamos/generator GmGenerator GmGenUtils::FileInPath ./isotopes.dat in . GmGenUtils::FileInPath /home/arce/gamosnew/GAMOS.5.0.0/data/isotopes.dat in /home/arce/gamosnew/GAMOS.5.0.0/data /run/initialize GmGenUtils::FileInPath ./../exercise3/world.geom in . GmEMPhysics PhysicsList:DEFAULT = gamma-lowener is registered GmEMPhysics PhysicsList: DEFAULT = electron-lowener is registered GmEMPhysics PhysicsList: DEFAULT = positron-standard is registered GmPhysicsGammaLowEner::ConstructProcess /gamos/physics/addPhysics radioactiveDecay /gamos/filter nu_eF GmParticleFilter nu_e /gamos/userAction GmKillAtStackingActionUA nu_eF /gamos/generator/addSingleParticleSource I125_source I125[0.0] 0.1*keV /gamos/scoring/createMFDetector PDDdet phantom /gamos/scoring/addScorer2MFD PDDscorer GmG4PSDoseDeposit PDDdet GmVPrimitiveScorer::SetParameters: no parameters /gamos/scoring/addPrinter2Scorer RTPSPDoseHistos PDDscorer phantomContainer GmRegularParamUtils::GetPhantomMotherTranslation (0,0,0) global (0,0,0) phantomContainer GmRegularParamUtils::GetPhantomMotherRotation [ ( 1 0 0) ( 0 1 0) ( 0 0 1) ] global [ ( 1 0 0) ( 0 1 0) ( 0 0 1) ] RTPSPDoseHistos nvoxel 100 100 100 RTPSPDoseHistos half dim 0.5 0.5 0.5 RTPSPDoseHistos translation (0,0,0) RTPSPDoseHistos rotation [ ( 1 0 0) ( 0 1 0) ( 0 0 1) ] CreateHisto1D RTPSPDoseHistos: Dose = 36201 CreateHisto1D RTPSPDoseHistos: Dose-volume = 36202 /gamos/userAction GmCountProcessesUA /run/beamOn 10000 phot: for gamma SubType= 12 LambdaPrime table from 200 keV to 10 TeV in 54 bins ===== EM models for the G4Region DefaultRegionForTheWorld ====== LivermorePhElectric : Emin= 0 eV Emax= 1 GeV AngularGenSauterGavrila FluoActive PhotoElectric : Emin= 1 GeV Emax= 10 TeV AngularGenSauterGavrila FluoActive compt: for gamma SubType= 13 Lambda table from 100 eV to 1 MeV in 28 bins, spline: 1 LambdaPrime table from 1 MeV to 10 TeV in 49 bins ===== EM models for the G4Region DefaultRegionForTheWorld ====== LivermoreCompton : Emin= 0 eV Emax= 1 GeV FluoActive Klein-Nishina : Emin= 1 GeV Emax= 10 TeV conv: for gamma SubType= 14 Lambda table from 1.022 MeV to 10 TeV in 49 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== LivermoreConversion : Emin= 0 eV Emax= 1 GeV BetheHeitler : Emin= 1 GeV Emax= 80 GeV BetheHeitlerLPM : Emin= 80 GeV Emax= 10 TeV Rayl: for gamma SubType= 11 Lambda table from 100 eV to 100 keV in 21 bins, spline: 0 LambdaPrime table from 100 keV to 10 TeV in 56 bins ===== EM models for the G4Region DefaultRegionForTheWorld ====== LivermoreRayleigh : Emin= 0 eV Emax= 1 GeV CullenGenerator LivermoreRayleigh : Emin= 1 GeV Emax= 10 TeV CullenGenerator msc: for e- SubType= 10 RangeFactor= 0.04, step limit type: 2, lateralDisplacement: 1, skin= 1, geomFactor= 2.5 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 77 bins Emin= 100 eV Emax= 10 TeV ### === Deexcitation model UAtomDeexcitation is activated for 1 region: DefaultRegionForTheWorld eIoni: for e- SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== LowEnergyIoni : Emin= 0 eV Emax= 1 GeV MollerBhabha : Emin= 1 GeV Emax= 10 TeV eBrem: for e- SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 LPM flag: 1 for E > 1 GeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== LowEnBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban eBremLPM : Emin= 1 GeV Emax= 10 TeV DipBustGen msc: for e+ SubType= 10 RangeFactor= 0.04, step limit type: 1, lateralDisplacement: 1, skin= 1, geomFactor= 2.5 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 77 bins Emin= 100 eV Emax= 10 TeV eIoni: for e+ SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MollerBhabha : Emin= 0 eV Emax= 10 TeV eBrem: for e+ SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 LPM flag: 1 for E > 1 GeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== eBremSB : Emin= 0 eV Emax= 1 GeV DipBustGen eBremLPM : Emin= 1 GeV Emax= 10 TeV DipBustGen annihil: for e+ SubType= 5 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eplus2gg : Emin= 0 eV Emax= 10 TeV msc: for proton SubType= 10 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 77 bins Emin= 100 eV Emax= 10 TeV hIoni: for proton SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 2 MeV BetheBloch : Emin= 2 MeV Emax= 10 TeV hBrems: for proton SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 10 TeV hPairProd: for proton SubType= 4 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 10 TeV ### === Deexcitation model UAtomDeexcitation is activated for 1 region: DefaultRegionForTheWorld msc: for GenericIon SubType= 10 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc95 : Emin= 0 eV Emax= 10 TeV ionIoni: for GenericIon SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 finalRange(mm)= 0.01, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02 Stopping Power data for 17 ion/material pairs ===== EM models for the G4Region DefaultRegionForTheWorld ====== ParamICRU73 : Emin= 0 eV Emax= 10 TeV nuclearStopping: for GenericIon SubType= 8 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV msc: for alpha SubType= 10 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 77 bins Emin= 100 eV Emax= 10 TeV ionIoni: for alpha SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 finalRange(mm)= 0.02, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02 ===== EM models for the G4Region DefaultRegionForTheWorld ====== BraggIon : Emin= 0 eV Emax= 7.9452 MeV BetheBloch : Emin= 7.9452 MeV Emax= 10 TeV nuclearStopping: for alpha SubType= 8 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV msc: for anti_proton SubType= 10 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 77 bins Emin= 100 eV Emax= 10 TeV hIoni: for anti_proton SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax= 2 MeV BetheBloch : Emin= 2 MeV Emax= 10 TeV msc: for kaon+ SubType= 10 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 77 bins Emin= 100 eV Emax= 10 TeV hIoni: for kaon+ SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 1.05231 MeV BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV hBrems: for kaon+ SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 10 TeV hPairProd: for kaon+ SubType= 4 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 10 TeV msc: for kaon- SubType= 10 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 77 bins Emin= 100 eV Emax= 10 TeV hIoni: for kaon- SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax= 1.05231 MeV BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV hBrems: for kaon- SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 10 TeV hPairProd: for kaon- SubType= 4 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 10 TeV msc: for mu+ SubType= 10 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 10 TeV Table with 77 bins Emin= 100 eV Emax= 10 TeV muIoni: for mu+ SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 200 keV BetheBloch : Emin= 200 keV Emax= 1 GeV MuBetheBloch : Emin= 1 GeV Emax= 10 TeV muBrems: for mu+ SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MuBrem : Emin= 0 eV Emax= 10 TeV muPairProd: for mu+ SubType= 4 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== muPairProd : Emin= 0 eV Emax= 10 TeV CoulombScat: for mu+ SubType= 1 Lambda table from 100 eV to 10 TeV in 77 bins, spline: 1 0 < Theta(degree) < 180; pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 10 TeV msc: for mu- SubType= 10 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 10 TeV Table with 77 bins Emin= 100 eV Emax= 10 TeV muIoni: for mu- SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax= 200 keV BetheBloch : Emin= 200 keV Emax= 1 GeV MuBetheBloch : Emin= 1 GeV Emax= 10 TeV muBrems: for mu- SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MuBrem : Emin= 0 eV Emax= 10 TeV muPairProd: for mu- SubType= 4 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== muPairProd : Emin= 0 eV Emax= 10 TeV CoulombScat: for mu- SubType= 1 Lambda table from 100 eV to 10 TeV in 77 bins, spline: 1 0 < Theta(degree) < 180; pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 10 TeV msc: for pi+ SubType= 10 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 77 bins Emin= 100 eV Emax= 10 TeV hIoni: for pi+ SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 297.505 keV BetheBloch : Emin= 297.505 keV Emax= 10 TeV hBrems: for pi+ SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 10 TeV hPairProd: for pi+ SubType= 4 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 10 TeV msc: for pi- SubType= 10 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 77 bins Emin= 100 eV Emax= 10 TeV hIoni: for pi- SubType= 2 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax= 297.505 keV BetheBloch : Emin= 297.505 keV Emax= 10 TeV hBrems: for pi- SubType= 3 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 10 TeV hPairProd: for pi- SubType= 4 dE/dx and range tables from 100 eV to 10 TeV in 77 bins Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 10 TeV ========= Table of registered couples ============================== Index : 0 used in the geometry : Yes recalculation needed : No Material : G4_AIR Range cuts : gamma 100 um e- 100 um e+ 100 um proton 100 um Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 10 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 1 used in the geometry : Yes recalculation needed : No Material : G4_WATER1 Range cuts : gamma 100 um e- 100 um e+ 100 um proton 100 um Energy thresholds : gamma 1.10981 keV e- 84.6589 keV e+ 83.5287 keV proton 10 keV Region(s) which use this couple : DefaultRegionForTheWorld ==================================================================== G4RunManagerKernel -- G4ScoreSplittingProcess is appended to all particles. @@@@@ GmCountProcessesUA GmCountProcessesUA Start Run processing. Run terminated. Run Summary Number of events processed : 10000 User=113.66s Real=114.77s Sys=0.12s RTPSPDoseHistos MINIMUM DOSE 2.19801e-16 Gy RTPSPDoseHistos MAXIMUM DOSE 3.96138e-10 Gy RTPSPDoseHistos AVERAGE ERROR 20% = 0.0223639 RTPSPDoseHistos AVERAGE ERROR 50% = 0.0223639 RTPSPDoseHistos AVERAGE ERROR 90% = 0.0223639 CreateHisto1D RTPSPDoseHistos: Dose Profile X_merged = 66301 CreateHisto1D RTPSPDoseHistos: Dose Profile Y_merged = 66302 CreateHisto1D RTPSPDoseHistos: Dose Profile Z_merged = 66303 CreateHisto2D RTPSPDoseHistos: Dose XY_merged = 66304 CreateHisto2D RTPSPDoseHistos: Dose XZ_merged = 66305 CreateHisto2D RTPSPDoseHistos: Dose YZ_merged = 66306 @@@@@ GmCountProcessesUA PART_COUNT: I125[0.0] = 10000 PART_COUNT: Te125[0.0] = 10000 PART_COUNT: Te125[35.5] = 10000 PART_COUNT: e- = 74676 PART_COUNT: gamma = 16223 PROC_COUNT I125[0.0] : RadioactiveDecay = 10000 PROC_COUNT I125[0.0] : ionIoni_NoSeco = 10000 PROC_COUNT I125[0.0] : ALL = 20000 PROC_COUNT Te125[0.0] : RadioactiveDecay = 10000 PROC_COUNT Te125[0.0] : ALL = 10000 PROC_COUNT Te125[35.5] : RadioactiveDecay = 10000 PROC_COUNT Te125[35.5] : Transportation = 222 PROC_COUNT Te125[35.5] : ionIoni_NoSeco = 9776 PROC_COUNT Te125[35.5] : ALL = 19998 PROC_COUNT e- : Transportation = 70273 PROC_COUNT e- : eBrem = 21 PROC_COUNT e- : eBrem_NoSeco = 2 PROC_COUNT e- : eIoni = 89 PROC_COUNT e- : eIoni_NoSeco = 70880 PROC_COUNT e- : msc = 251337 PROC_COUNT e- : ALL = 392602 PROC_COUNT gamma : Rayl = 3221 PROC_COUNT gamma : Transportation = 882245 PROC_COUNT gamma : compt = 10657 PROC_COUNT gamma : phot = 13050 PROC_COUNT gamma : ALL = 909173 PROC_CREATOR_COUNT I125[0.0] : Primary = 10000 PROC_CREATOR_COUNT Te125[0.0] : RadioactiveDecay = 10000 PROC_CREATOR_COUNT Te125[35.5] : RadioactiveDecay = 10000 PROC_CREATOR_COUNT e- : RadioactiveDecay = 50872 PROC_CREATOR_COUNT e- : compt = 10657 PROC_CREATOR_COUNT e- : eIoni = 89 PROC_CREATOR_COUNT e- : phot = 13058 PROC_CREATOR_COUNT gamma : RadioactiveDecay = 16201 PROC_CREATOR_COUNT gamma : eBrem = 21 PROC_CREATOR_COUNT gamma : phot = 1 Graphics systems deleted. Visualization Manager deleting... === saving histograms in file === dose_RTPSPDoseHistos.root