################################ !!! G4Backtrace is activated !!! ################################ ************************************************************** Geant4 version Name: geant4-11-00-patch-02 [MT] (25-May-2022) << in Multi-threaded mode >> 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) G4HepRepFile (HepRepFile) RayTracer (RayTracer) VRML2FILE (VRML2FILE) gMocrenFile (gMocrenFile) OpenGLImmediateQt (OGLIQt, OGLI) OpenGLStoredQt (OGLSQt, OGL, OGLS) OpenGLImmediateX (OGLIX, OGLIQt_FALLBACK) OpenGLStoredX (OGLSX, OGLSQt_FALLBACK) Qt3D (Qt3D) 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. Checking overlaps for volume crystal:0 (G4Box) ... OK! Checking overlaps for volume crystal:1 (G4Box) ... OK! Checking overlaps for volume crystal:2 (G4Box) ... OK! Checking overlaps for volume crystal:3 (G4Box) ... OK! Checking overlaps for volume crystal:4 (G4Box) ... OK! Checking overlaps for volume crystal:5 (G4Box) ... OK! Checking overlaps for volume crystal:6 (G4Box) ... OK! Checking overlaps for volume crystal:7 (G4Box) ... OK! Checking overlaps for volume crystal:8 (G4Box) ... OK! Checking overlaps for volume crystal:9 (G4Box) ... OK! Checking overlaps for volume crystal:10 (G4Box) ... OK! Checking overlaps for volume crystal:11 (G4Box) ... OK! Checking overlaps for volume crystal:12 (G4Box) ... OK! Checking overlaps for volume crystal:13 (G4Box) ... OK! Checking overlaps for volume crystal:14 (G4Box) ... OK! Checking overlaps for volume crystal:15 (G4Box) ... OK! Checking overlaps for volume crystal:16 (G4Box) ... OK! Checking overlaps for volume crystal:17 (G4Box) ... OK! Checking overlaps for volume crystal:18 (G4Box) ... OK! Checking overlaps for volume crystal:19 (G4Box) ... OK! Checking overlaps for volume crystal:20 (G4Box) ... OK! Checking overlaps for volume crystal:21 (G4Box) ... OK! Checking overlaps for volume crystal:22 (G4Box) ... OK! Checking overlaps for volume crystal:23 (G4Box) ... OK! Checking overlaps for volume crystal:24 (G4Box) ... OK! Checking overlaps for volume crystal:25 (G4Box) ... OK! Checking overlaps for volume crystal:26 (G4Box) ... OK! Checking overlaps for volume crystal:27 (G4Box) ... OK! Checking overlaps for volume crystal:28 (G4Box) ... OK! Checking overlaps for volume crystal:29 (G4Box) ... OK! Checking overlaps for volume crystal:30 (G4Box) ... OK! Checking overlaps for volume crystal:31 (G4Box) ... OK! Checking overlaps for volume ring:0 (G4Tubs) ... OK! Checking overlaps for volume ring:1 (G4Tubs) ... OK! Checking overlaps for volume ring:2 (G4Tubs) ... OK! Checking overlaps for volume ring:3 (G4Tubs) ... OK! Checking overlaps for volume ring:4 (G4Tubs) ... OK! Checking overlaps for volume ring:5 (G4Tubs) ... OK! Checking overlaps for volume ring:6 (G4Tubs) ... OK! Checking overlaps for volume ring:7 (G4Tubs) ... OK! Checking overlaps for volume ring:8 (G4Tubs) ... OK! Checking overlaps for volume Detector:0 (G4Tubs) ... OK! Checking overlaps for volume Patient:0 (G4Tubs) ... OK! ***** Table : Nb of materials = 3 ***** Material: Lu2SiO5 density: 7.400 g/cm3 RadL: 1.143 cm Nucl.Int.Length: 20.920 cm Imean: 418.807 eV temperature: 293.15 K pressure: 1.00 atm ---> Element: Lu (Lu) Z = 71.0 N = 175 A = 174.967 g/mole ---> Isotope: Lu175 Z = 71 N = 175 A = 174.94 g/mole abundance: 97.410 % ---> Isotope: Lu176 Z = 71 N = 176 A = 175.94 g/mole abundance: 2.590 % ElmMassFraction: 76.40 % ElmAbundance 25.00 % ---> Element: Si (Si) Z = 14.0 N = 28 A = 28.085 g/mole ---> Isotope: Si28 Z = 14 N = 28 A = 27.98 g/mole abundance: 92.230 % ---> Isotope: Si29 Z = 14 N = 29 A = 28.98 g/mole abundance: 4.683 % ---> Isotope: Si30 Z = 14 N = 30 A = 29.97 g/mole abundance: 3.087 % ElmMassFraction: 6.13 % ElmAbundance 12.50 % ---> Element: O (O) Z = 8.0 N = 16 A = 15.999 g/mole ---> Isotope: O16 Z = 8 N = 16 A = 15.99 g/mole abundance: 99.757 % ---> Isotope: O17 Z = 8 N = 17 A = 17.00 g/mole abundance: 0.038 % ---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.205 % ElmMassFraction: 17.47 % ElmAbundance 62.50 % Material: G4_AIR density: 1.205 mg/cm3 RadL: 303.921 m Nucl.Int.Length: 710.095 m Imean: 85.700 eV temperature: 293.15 K pressure: 1.00 atm ---> Element: C (C) Z = 6.0 N = 12 A = 12.011 g/mole ---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.930 % ---> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.070 % ElmMassFraction: 0.01 % ElmAbundance 0.02 % ---> Element: N (N) Z = 7.0 N = 14 A = 14.007 g/mole ---> Isotope: N14 Z = 7 N = 14 A = 14.00 g/mole abundance: 99.632 % ---> Isotope: N15 Z = 7 N = 15 A = 15.00 g/mole abundance: 0.368 % ElmMassFraction: 75.53 % ElmAbundance 78.44 % ---> Element: O (O) Z = 8.0 N = 16 A = 15.999 g/mole ---> Isotope: O16 Z = 8 N = 16 A = 15.99 g/mole abundance: 99.757 % ---> Isotope: O17 Z = 8 N = 17 A = 17.00 g/mole abundance: 0.038 % ---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.205 % ElmMassFraction: 23.18 % ElmAbundance 21.07 % ---> Element: Ar (Ar) Z = 18.0 N = 40 A = 39.948 g/mole ---> Isotope: Ar36 Z = 18 N = 36 A = 35.97 g/mole abundance: 0.337 % ---> Isotope: Ar38 Z = 18 N = 38 A = 37.96 g/mole abundance: 0.063 % ---> Isotope: Ar40 Z = 18 N = 40 A = 39.96 g/mole abundance: 99.600 % ElmMassFraction: 1.28 % ElmAbundance 0.47 % Material: G4_BRAIN_ICRP density: 1.040 g/cm3 RadL: 35.403 cm Nucl.Int.Length: 72.156 cm Imean: 73.300 eV temperature: 293.15 K pressure: 1.00 atm ---> Element: H (H) Z = 1.0 N = 1 A = 1.008 g/mole ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.989 % ---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.011 % ElmMassFraction: 10.70 % ElmAbundance 64.44 % ---> Element: C (C) Z = 6.0 N = 12 A = 12.011 g/mole ---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.930 % ---> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.070 % ElmMassFraction: 14.50 % ElmAbundance 7.33 % ---> Element: N (N) Z = 7.0 N = 14 A = 14.007 g/mole ---> Isotope: N14 Z = 7 N = 14 A = 14.00 g/mole abundance: 99.632 % ---> Isotope: N15 Z = 7 N = 15 A = 15.00 g/mole abundance: 0.368 % ElmMassFraction: 2.20 % ElmAbundance 0.95 % ---> Element: O (O) Z = 8.0 N = 16 A = 15.999 g/mole ---> Isotope: O16 Z = 8 N = 16 A = 15.99 g/mole abundance: 99.757 % ---> Isotope: O17 Z = 8 N = 17 A = 17.00 g/mole abundance: 0.038 % ---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.205 % ElmMassFraction: 71.20 % ElmAbundance 27.01 % ---> Element: Na (Na) Z = 11.0 N = 23 A = 22.990 g/mole ---> Isotope: Na23 Z = 11 N = 23 A = 22.99 g/mole abundance: 100.000 % ElmMassFraction: 0.20 % ElmAbundance 0.05 % ---> Element: P (P) Z = 15.0 N = 31 A = 30.974 g/mole ---> Isotope: P31 Z = 15 N = 31 A = 30.97 g/mole abundance: 100.000 % ElmMassFraction: 0.40 % ElmAbundance 0.08 % ---> Element: S (S) Z = 16.0 N = 32 A = 32.066 g/mole ---> Isotope: S32 Z = 16 N = 32 A = 31.97 g/mole abundance: 94.930 % ---> Isotope: S33 Z = 16 N = 33 A = 32.97 g/mole abundance: 0.760 % ---> Isotope: S34 Z = 16 N = 34 A = 33.97 g/mole abundance: 4.290 % ---> Isotope: S36 Z = 16 N = 36 A = 35.97 g/mole abundance: 0.020 % ElmMassFraction: 0.20 % ElmAbundance 0.04 % ---> Element: Cl (Cl) Z = 17.0 N = 35 A = 35.453 g/mole ---> Isotope: Cl35 Z = 17 N = 35 A = 34.97 g/mole abundance: 75.780 % ---> Isotope: Cl37 Z = 17 N = 37 A = 36.97 g/mole abundance: 24.220 % ElmMassFraction: 0.30 % ElmAbundance 0.05 % ---> Element: K (K) Z = 19.0 N = 39 A = 39.098 g/mole ---> Isotope: K39 Z = 19 N = 39 A = 38.96 g/mole abundance: 93.258 % ---> Isotope: K40 Z = 19 N = 40 A = 39.96 g/mole abundance: 0.012 % ---> Isotope: K41 Z = 19 N = 41 A = 40.96 g/mole abundance: 6.730 % ElmMassFraction: 0.30 % ElmAbundance 0.05 % G4SDManager::AddNewCollection : the collection is registered at 1 G4SDManager::AddNewCollection : the collection is registered at 2 ======================================================================= ====== Electromagnetic Physics Parameters ======== ======================================================================= LPM effect enabled 1 Enable creation and use of sampling tables 0 Apply cuts on all EM processes 0 Use general process 0 Enable linear polarisation for gamma 0 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 7 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 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 Livermore data directory livermore ======================================================================= ====== Ionisation Parameters ======== ======================================================================= Step function for e+- (0.2, 1 mm) Step function for muons/hadrons (0.2, 0.1 mm) Step function for light ions (0.2, 0.1 mm) Step function for general ions (0.2, 0.1 mm) Lowest e+e- kinetic energy 1 keV Lowest muon/hadron kinetic energy 1 keV Fluctuations of dE/dx are enabled 1 Use ICRU90 data 0 Use built-in Birks satuaration 0 Build CSDA range enabled 0 Use cut as a final range enabled 0 Enable angular generator interface 0 Max kinetic energy for CSDA tables 1 GeV Max kinetic energy for NIEL computation 0 eV 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+- 1 Type of msc step limit algorithm for muons/hadrons 0 Msc lateral displacement for e+- enabled 1 Msc lateral displacement for muons and hadrons 0 Urban msc model lateral displacement alg96 1 Range factor for msc step limit for e+- 0.04 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+- 1 Lambda limit for msc step limit for e+- 1 mm Use Mott correction for e- scattering 0 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 Fluorescence Bearden data files enabled 0 Fluorescence ANSTO data files enabled 0 Auger electron cascade enabled 1 PIXE atomic de-excitation enabled 0 De-excitation module ignores cuts 1 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 1 0 ### === Auger flag: 1 ### === Ignore cuts flag: 1 phot: for gamma SubType=12 BuildTable=0 LambdaPrime table from 200 keV to 100 TeV in 61 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, 7 bins/decade, spline: 1 LambdaPrime table from 1 MeV to 100 TeV in 56 bins ===== EM models for the G4Region DefaultRegionForTheWorld ====== Klein-Nishina : Emin= 0 eV Emax= 100 TeV conv: for gamma SubType=14 BuildTable=1 Lambda table from 1.022 MeV to 100 TeV, 18 bins/decade, spline: 1 ===== EM models for the G4Region DefaultRegionForTheWorld ====== BetheHeitlerLPM : Emin= 0 eV Emax= 100 TeV ModifiedTsai Rayl: for gamma SubType=11 BuildTable=1 Lambda table from 100 eV to 100 keV, 7 bins/decade, spline: 0 LambdaPrime table from 100 keV to 100 TeV in 63 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 ====== UrbanMsc : Emin= 0 eV Emax= 100 MeV Nbins=42 100 eV - 100 MeV StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm WentzelVIUni : Emin= 100 MeV Emax= 100 TeV Nbins=42 100 MeV - 100 TeV StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm eIoni: for e- XStype:1 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 1 mm), integ: 1, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MollerBhabha : Emin= 0 eV Emax= 100 TeV eBrem: for e- XStype:4 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 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 ModifiedTsai eBremLPM : Emin= 1 GeV Emax= 100 TeV ModifiedTsai CoulombScat: for e- XStype:1 SubType=1 BuildTable=1 Lambda table from 100 MeV to 100 TeV, 7 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 ====== UrbanMsc : Emin= 0 eV Emax= 100 MeV Nbins=42 100 eV - 100 MeV StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm WentzelVIUni : Emin= 100 MeV Emax= 100 TeV Nbins=42 100 MeV - 100 TeV StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm eIoni: for e+ XStype:1 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 1 mm), integ: 1, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MollerBhabha : Emin= 0 eV Emax= 100 TeV eBrem: for e+ XStype:4 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 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 ModifiedTsai eBremLPM : Emin= 1 GeV Emax= 100 TeV ModifiedTsai annihil: for e+ XStype:2 SubType=5 BuildTable=0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eplus2gg : Emin= 0 eV Emax= 100 TeV CoulombScat: for e+ XStype:1 SubType=1 BuildTable=1 Lambda table from 100 MeV to 100 TeV, 7 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=84 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm hIoni: for proton XStype:1 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax= 2 MeV BetheBloch : Emin= 2 MeV Emax= 100 TeV hBrems: for proton XStype:1 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi hPairProd: for proton XStype:1 SubType=4 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 0 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, 7 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 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:0 Skin=1 Llim=1 mm ionIoni: for GenericIon XStype:1 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.02 Stopping Power data for 17 ion/material pairs ===== EM models for the G4Region DefaultRegionForTheWorld ====== BraggIon : Emin= 0 eV Emax= 2 MeV BetheBloch : Emin= 2 MeV Emax= 100 TeV ====================================================================== ====== Radioactive Decay Physics Parameters ======= ====================================================================== Max life time 1.4427e+06 ps Internal e- conversion flag 1 Stored internal conversion coefficients 1 Enable correlated gamma emission 0 Max 2J for sampling of angular correlations 10 Atomic de-excitation enabled 1 Auger electron emission enabled 1 Check EM cuts disabled for atomic de-excitation 1 Use Bearden atomic level energies 0 Use ANSTO fluorescence model 0 Threshold for very long decay time at rest 1e+27 ns ====================================================================== 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:0 Skin=1 Llim=1 mm ionIoni: for alpha XStype:1 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.02 ===== EM models for the G4Region DefaultRegionForTheWorld ====== BraggIon : Emin= 0 eV Emax=7.9452 MeV BetheBloch : Emin=7.9452 MeV Emax= 100 TeV msc: for anti_proton SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm hIoni: for anti_proton XStype:1 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax= 2 MeV BetheBloch : Emin= 2 MeV Emax= 100 TeV hBrems: for anti_proton XStype:1 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi hPairProd: for anti_proton XStype:1 SubType=4 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 0 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, 7 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=84 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm hIoni: for kaon+ XStype:1 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax=1.05231 MeV BetheBloch : Emin=1.05231 MeV Emax= 100 TeV hBrems: for kaon+ XStype:1 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi hPairProd: for kaon+ XStype:1 SubType=4 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 0 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, 7 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=84 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm hIoni: for kaon- XStype:1 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax=1.05231 MeV BetheBloch : Emin=1.05231 MeV Emax= 100 TeV hBrems: for kaon- XStype:1 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi hPairProd: for kaon- XStype:1 SubType=4 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 0 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=84 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm muIoni: for mu+ XStype:1 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 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= 100 TeV muBrems: for mu+ XStype:1 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MuBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi muPairProd: for mu+ XStype:1 SubType=4 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 0 Sampling table 21x1001; from 1 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, 7 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=84 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm muIoni: for mu- XStype:1 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 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= 100 TeV muBrems: for mu- XStype:1 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== MuBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi muPairProd: for mu- XStype:1 SubType=4 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 0 Sampling table 21x1001; from 1 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=84 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm hIoni: for pi+ XStype:1 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== Bragg : Emin= 0 eV Emax=297.505 keV BetheBloch : Emin=297.505 keV Emax= 100 TeV hBrems: for pi+ XStype:1 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi hPairProd: for pi+ XStype:1 SubType=4 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 0 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, 7 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=84 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm hIoni: for pi- XStype:1 SubType=2 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1 StepFunction=(0.2, 0.1 mm), integ: 1, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU73QO : Emin= 0 eV Emax=297.505 keV BetheBloch : Emin=297.505 keV Emax= 100 TeV hBrems: for pi- XStype:1 SubType=3 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi hPairProd: for pi- XStype:1 SubType=4 dE/dx and range tables from 100 eV to 100 TeV in 84 bins Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 0 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 ========================================================================= --> G4TaskRunManager::CreateAndStartWorkers() --> Initializing workers... ========================================================================= # /run/printProgress 1000 /run/beamOn 10000 G4VisManager: Using G4TrajectoryDrawByCharge as fallback trajectory model. See commands in /vis/modeling/trajectories/ for other options. ### Run 0 starts. ### Run 0 start. ... set ntuple merging row mode : row-wise - done ... create file : scoring.root - done ... open analysis file : scoring.root - done ... open analysis file : scoring.root - done G4WT7 > /run/printProgress 1000 G4WT6 > /run/printProgress 1000 G4WT0 > /run/printProgress 1000 G4WT5 > /run/printProgress 1000 G4WT2 > /run/printProgress 1000 G4WT4 > /run/printProgress 1000 G4WT3 > /run/printProgress 1000 G4WT1 > /run/printProgress 1000 ============================================================================================ --> G4TaskRunManager::CreateAndStartWorkers() --> Creating 100 tasks with 100 events/task... ============================================================================================ G4WT7 > ### Run 0 starts on worker thread 7. G4WT1 > ### Run 0 starts on worker thread 1. G4WT1 > ### Run 0 start. G4WT6 > ### Run 0 starts on worker thread 6. G4WT7 > ### Run 0 start. G4WT2 > ### Run 0 starts on worker thread 2. G4WT4 > ### Run 0 starts on worker thread 4. G4WT1 > ... set ntuple merging row mode : row-wise - done G4WT5 > ### Run 0 starts on worker thread 5. G4WT3 > ### Run 0 starts on worker thread 3. G4WT6 > ### Run 0 start. G4WT7 > ... set ntuple merging row mode : row-wise - done G4WT6 > ... set ntuple merging row mode : row-wise - done G4WT2 > ### Run 0 start. G4WT0 > ### Run 0 starts on worker thread 0. G4WT4 > ### Run 0 start. G4WT5 > ### Run 0 start. G4WT3 > ### Run 0 start. G4WT2 > ... set ntuple merging row mode : row-wise - done G4WT3 > ... set ntuple merging row mode : row-wise - done G4WT0 > ### Run 0 start. G4WT5 > ... set ntuple merging row mode : row-wise - done G4WT0 > ... set ntuple merging row mode : row-wise - done G4WT4 > ... set ntuple merging row mode : row-wise - done G4WT1 > ... create file : scoring_t1.root - done G4WT1 > ... open analysis file : scoring.root - done G4WT1 > ... open analysis file : scoring.root - done G4WT1 > --> Event 0 starts with initial seeds (13049039,61775110). G4WT6 > ... create file : scoring_t6.root - done G4WT5 > ... create file : scoring_t5.root - done G4WT0 > ... create file : scoring_t0.root - done G4WT5 > ... open analysis file : scoring.root - done G4WT4 > ... create file : scoring_t4.root - done G4WT6 > ... open analysis file : scoring.root - done G4WT4 > ... open analysis file : scoring.root - done G4WT7 > ... create file : scoring_t7.root - done G4WT7 > ... open analysis file : scoring.root - done G4WT3 > ... create file : scoring_t3.root - done G4WT3 > ... open analysis file : scoring.root - done G4WT0 > ... open analysis file : scoring.root - done G4WT2 > ... create file : scoring_t2.root - done G4WT2 > ... open analysis file : scoring.root - done G4WT5 > ... open analysis file : scoring.root - done G4WT2 > ... open analysis file : scoring.root - done G4WT4 > ... open analysis file : scoring.root - done G4WT0 > ... open analysis file : scoring.root - done G4WT6 > ... open analysis file : scoring.root - done G4WT7 > ... open analysis file : scoring.root - done G4WT3 > ... open analysis file : scoring.root - done G4WT5 > --> Event 1000 starts with initial seeds (32766986,29233990). G4WT7 > --> Event 2000 starts with initial seeds (49764252,54850316). G4WT6 > --> Event 3000 starts with initial seeds (74171463,44032866). G4WT7 > --> Event 4000 starts with initial seeds (66846944,33828183). G4WT6 > --> Event 5000 starts with initial seeds (8383287,34842104). G4WT7 > --> Event 6000 starts with initial seeds (99382317,80379040). G4WT4 > --> Event 7000 starts with initial seeds (79437796,16414614). G4WT6 > --> Event 8000 starts with initial seeds (15511632,88960611). G4WT2 > --> Event 9000 starts with initial seeds (38426545,76991870). G4WT7 > G4WT7 > --------------------End of Local Run------------------------ G4WT7 > The run was 1100 F18; Nb of 'good' e+ annihilations: 132 G4WT7 > Total dose in patient : 32.882 picoGy G4WT7 > ------------------------------------------------------------ G4WT2 > G4WT1 > G4WT6 > G4WT6 > --------------------End of Local Run------------------------ G4WT6 > The run was 1300 F18; Nb of 'good' e+ annihilations: 161 G4WT6 > Total dose in patient : 39.8074 picoGy G4WT4 > G4WT3 > G4WT3 > --------------------End of Local Run------------------------ G4WT3 > The run was 1100 F18; Nb of 'good' e+ annihilations: 126 G4WT3 > Total dose in patient : 34.0866 picoGy G4WT5 > G4WT5 > --------------------End of Local Run------------------------ G4WT5 > The run was 1300 F18; Nb of 'good' e+ annihilations: 164 G4WT0 > G4WT0 > --------------------End of Local Run------------------------ G4WT0 > The run was 1300 F18; Nb of 'good' e+ annihilations: 155 G4WT2 > --------------------End of Local Run------------------------ G4WT2 > The run was 1300 F18; Nb of 'good' e+ annihilations: 175 G4WT2 > Total dose in patient : 39.4576 picoGy G4WT2 > ------------------------------------------------------------ G4WT2 > G4WT1 > --------------------End of Local Run------------------------ G4WT2 > ... merge all H1 - done G4WT2 > ... merge all H2 - done G4WT2 > ... merge all H3 - done G4WT4 > --------------------End of Local Run------------------------ G4WT5 > Total dose in patient : 40.3098 picoGy G4WT5 > ------------------------------------------------------------ G4WT5 > G4WT0 > Total dose in patient : 40.4257 picoGy G4WT0 > ------------------------------------------------------------ G4WT0 > G4WT6 > ------------------------------------------------------------ G4WT1 > The run was 1400 F18; Nb of 'good' e+ annihilations: 177 G4WT1 > Total dose in patient : 43.1911 picoGy G4WT1 > ------------------------------------------------------------ G4WT1 > G4WT3 > ------------------------------------------------------------ G4WT3 > G4WT2 > ... merge all P1 - done G4WT2 > ... merge all P2 - done G4WT7 > G4WT4 > The run was 1200 F18; Nb of 'good' e+ annihilations: 142 G4WT4 > Total dose in patient : 36.505 picoGy G4WT4 > ------------------------------------------------------------ G4WT4 > G4WT5 > ... merge all H1 - done G4WT5 > ... merge all H2 - done G4WT5 > ... merge all H3 - done G4WT5 > ... merge all P1 - done G4WT6 > G4WT5 > ... merge all P2 - done G4WT3 > ... merge all H1 - done G4WT3 > ... merge all H2 - done G4WT3 > ... merge all H3 - done G4WT3 > ... merge all P1 - done G4WT3 > ... merge all P2 - done G4WT0 > ... merge all H1 - done G4WT0 > ... merge all H2 - done G4WT0 > ... merge all H3 - done G4WT0 > ... merge all P1 - done G4WT0 > ... merge all P2 - done G4WT1 > ... merge all H1 - done G4WT1 > ... merge all H2 - done G4WT1 > ... merge all H3 - done G4WT1 > ... merge all P1 - done G4WT1 > ... merge all P2 - done G4WT6 > ... merge all H1 - done G4WT6 > ... merge all H2 - done G4WT6 > ... merge all H3 - done G4WT4 > ... merge all H1 - done G4WT4 > ... merge all H2 - done G4WT4 > ... merge all H3 - done G4WT4 > ... merge all P1 - done G4WT4 > ... merge all P2 - done G4WT6 > ... merge all P1 - done G4WT6 > ... merge all P2 - done G4WT7 > ... merge all H1 - done G4WT7 > ... merge all H2 - done G4WT7 > ... merge all H3 - done G4WT7 > ... merge all P1 - done G4WT7 > ... merge all P2 - done G4WT3 > ... write file : scoring_t3.root - done G4WT2 > ... write file : scoring_t2.root - done G4WT5 > ... write file : scoring_t5.root - done G4WT1 > ... write file : scoring_t1.root - done G4WT4 > ... write file : scoring_t4.root - done G4WT0 > ... write file : scoring_t0.root - done G4WT6 > ... write file : scoring_t6.root - done G4WT7 > ... write file : scoring_t7.root - done G4WT3 > ... close file : scoring_t3.root - done G4WT1 > ... close file : scoring_t1.root - done G4WT5 > ... close file : scoring_t5.root - done G4WT0 > ... close file : scoring_t0.root - done G4WT2 > ... close file : scoring_t2.root - done G4WT7 > ... close file : scoring_t7.root - done G4WT6 > ... close file : scoring_t6.root - done G4WT4 > ... close file : scoring_t4.root - done --------------------End of Global Run----------------------- The run was 10000 events ; Nb of 'good' e+ annihilations: 1232 Total dose in patient : 306.665 picoGy ------------------------------------------------------------ ... write file : scoring.root - done ... close file : scoring.root - done ... delete empty file : scoring.root - done 0 events have been kept for refreshing and/or reviewing. "/vis/reviewKeptEvents" to review them one by one. "/vis/enable", then "/vis/viewer/flush" or "/vis/viewer/rebuild" to see them accumulated. Graphics systems deleted. Visualization Manager deleting... Deleting G4ThreadLocalSingletons for type G4GeometryTolerance ... Deleting G4ThreadLocalSingletons for type G4ProcessTable ... Deleting G4ThreadLocalSingletons for type G4PhysicsListHelper ... Deleting G4ThreadLocalSingletons for type G4AccumulableManager ... Deleting G4ThreadLocalSingletons for type G4LossTableManager ... Deleting G4ThreadLocalSingletons for type G4GenericAnalysisManager ... Deleting G4ThreadLocalSingletons for type G4String ... Deleting G4ThreadLocalSingletons for type std::__cxx11::basic_ostringstream, std::allocator > ... ThreadPool destroyed