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. ***** Table : Nb of materials = 8 ***** Material: Iron density: 7.870 g/cm3 RadL: 1.758 cm Nucl.Int.Length: 16.999 cm Imean: 286.000 eV temperature: 293.15 K pressure: 1.00 atm ---> Element: Fe (Fe) Z = 26.0 N = 56 A = 55.845 g/mole ---> Isotope: Fe54 Z = 26 N = 54 A = 53.94 g/mole abundance: 5.845 % ---> Isotope: Fe56 Z = 26 N = 56 A = 55.93 g/mole abundance: 91.754 % ---> Isotope: Fe57 Z = 26 N = 57 A = 56.94 g/mole abundance: 2.119 % ---> Isotope: Fe58 Z = 26 N = 58 A = 57.93 g/mole abundance: 0.282 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: H2O density: 1.000 g/cm3 RadL: 36.092 cm Nucl.Int.Length: 75.356 cm Imean: 68.998 eV temperature: 293.15 K pressure: 1.00 atm ---> Element: Hydrogen (H) Z = 1.0 N = 1 A = 1.010 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: 11.21 % ElmAbundance 66.67 % ---> Element: Oxygen (O) Z = 8.0 N = 16 A = 16.000 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: 88.79 % ElmAbundance 33.33 % Material: Air density: 1.290 mg/cm3 RadL: 285.161 m Nucl.Int.Length: 662.680 m Imean: 85.703 eV temperature: 293.15 K pressure: 1.00 atm ---> Element: Nitrogen (N) Z = 7.0 N = 14 A = 14.010 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: 70.00 % ElmAbundance 72.71 % ---> Element: Oxygen (O) Z = 8.0 N = 16 A = 16.000 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: 30.00 % ElmAbundance 27.29 % Material: vacuum density: 0.010 mg/cm3 RadL: 36.786 km Nucl.Int.Length: 85.486 km Imean: 85.703 eV temperature: 273.15 K pressure: 0.02 atm ---> Element: Nitrogen (N) Z = 7.0 N = 14 A = 14.010 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: 70.00 % ElmAbundance 72.71 % ---> Element: Oxygen (O) Z = 8.0 N = 16 A = 16.000 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: 30.00 % ElmAbundance 27.29 % Material: Argon density: 1.784 mg/cm3 RadL: 109.604 m Nucl.Int.Length: 670.814 m Imean: 188.000 eV temperature: 293.15 K pressure: 1.00 atm ---> Element: Argon (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: 100.00 % ElmAbundance 100.00 % Material: N2 density: 1.251 mg/cm3 RadL: 303.846 m Nucl.Int.Length: 674.714 m Imean: 82.000 eV temperature: 293.15 K pressure: 1.00 atm ---> Element: Nitrogen (N) Z = 7.0 N = 14 A = 14.010 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: 100.00 % ElmAbundance 100.00 % Material: O2 density: 1.429 mg/cm3 RadL: 239.621 m Nucl.Int.Length: 617.219 m Imean: 95.000 eV temperature: 293.15 K pressure: 1.00 atm ---> Element: Oxygen (O) Z = 8.0 N = 16 A = 16.000 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: 100.00 % ElmAbundance 100.00 % Material: LaboratoryVacuum density: 0.000 mg/cm3 RadL: 28329262.634 km Nucl.Int.Length: 66178558.148 km Imean: 85.664 eV temperature: 273.15 K pressure: 0.00 atm ---> Element: Nitrogen (N) Z = 7.0 N = 14 A = 14.010 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.57 % ElmAbundance 78.48 % ---> Element: Oxygen (O) Z = 8.0 N = 16 A = 16.000 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.15 % ElmAbundance 21.05 % ---> Element: Argon (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 % end material ----------------------------------------------------------- Geometry and materials ----------------------------------------------------------- ---> World: ---> LaboratoryVacuum in World ---> WorldSizeXY: 3 m ---> WorldSizeZ: 3 m ----------------------------------------------------------- ---> Purging Magnet: ---> Gap made of Iron ---> GapSizeY1: 10 cm ---> GapSizeY2: 10 cm ---> GapSizeX1: 10 cm ---> GapSizeX2: 18.37 cm ----------------------------------------------------------- ---> Measurement Volume: ---> LaboratoryVacuum in Measurement volume ---> MeasureVolumeXY: 2.8 m ---> MeasureVolumeZ: 1 cm ---> At SSD = -50.5 cm ----------------------------------------------------------- ----------------------------------------------------------- Magnetic field ----------------------------------------------------------- ---> Reading the field grid from PurgMag3D.TABLE ... [ Number of values x,y,z: 6 6 10 ] ---> ... done reading ---> assumed the order: x, y, z, Bx, By, Bz ---> Min values x,y,z: -5 -5 -26 cm ---> Max values x,y,z: 5 17 10 cm ---> The field will be offset by 0.44 cm After reordering if neccesary ---> Min values x,y,z: -5 -5 -26 cm ---> Max values x,y,z: 5 17 10 cm ---> Dif values x,y,z (range): 10 22 36 cm in z ----------------------------------------------------------- DeltaStep 0.01mm G4ChordFinder: stepperDriverId: 2 PurgMagPhysicsList::SetCuts:CutLength : 1 um -------- WWWW ------- G4Exception-START -------- WWWW ------- *** G4Exception : ProcCuts110 issued by : G4ProductionCuts::SetProductionCut Setting cuts for particles other than photon, e-, e+ or proton has no effect. *** This is just a warning message. *** -------- WWWW -------- G4Exception-END --------- WWWW ------- /run/verbose 0 /event/verbose 0 /run/initialize /run/beamOn 100000 ======================================================================= ====== 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 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 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 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 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 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 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 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 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 : LaboratoryVacuum Range cuts : gamma 1 um e- 1 um e+ 1 um proton 1 um Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 100 eV Region(s) which use this couple : DefaultRegionForTheWorld Index : 1 used in the geometry : Yes Material : Iron Range cuts : gamma 1 um e- 1 um e+ 1 um proton 1 um Energy thresholds : gamma 990 eV e- 5.58698 keV e+ 5.47819 keV proton 100 eV Region(s) which use this couple : DefaultRegionForTheWorld ================================================================== Ntuple-1 created Ntuple-2 created Ntuple-3 created ---> Run 0 (master) start. ------- MixMaxRng engine status ------- Current state vector is: mixmax state, file version 1.0 N=17 V[N]={906770732717044781, 629165745432651234, 1235682547346241386, 68420008233404568, 2262190991329242458, 2266470399991071809, 1976726662926872232, 245458862506414172, 1955974201201518530, 2155248512522080758, 604170912935414061, 1116171330120743511, 1861018313684488333, 1296715403254578286, 1549011045957234151, 370819759640195970, 2230139271784837643} counter= 17sumtot= 1977567618660788324 --------------------------------------- ---> Begin Of Event: 0 ---> Begin Of Event: 10000 ---> Begin Of Event: 20000 ---> Begin Of Event: 30000 ---> Begin Of Event: 40000 ---> Begin Of Event: 50000 ---> Begin Of Event: 60000 ---> Begin Of Event: 70000 ---> Begin Of Event: 80000 ---> Begin Of Event: 90000 Total number of event = 100000 ------- MixMaxRng engine status ------- Current state vector is: mixmax state, file version 1.0 N=17 V[N]={77748014430414613, 1916023684659400533, 1170895291302529645, 1308717968961542533, 617576989286624977, 4159304169773412, 1773860259953165952, 1170460085642752549, 1346723549001032558, 2164775407857314692, 1400227920416410176, 1742052209028641158, 51329690927869269, 1772437326104352635, 406661536190202989, 1018565785029837114, 1261036022485472146} counter= 10sumtot= 756506971737785343 --------------------------------------- Graphics systems deleted. Visualization Manager deleting...