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/ ************************************************************** <<< Geant4 Physics List simulation engine: FTFP_BERT G4VModularPhysicsList::ReplacePhysics: G4EmStandard with type : 2 is replaced with G4EmStandard_opt4 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. /tracking/verbose 0 /run/verbose 1 /process/optical/verbose 1 /run/physicsModified /control/cout/ignoreThreadsExcept 0 /opnovice2/boxMaterial G4_BGO Tank material set to G4_BGO /opnovice2/worldMaterial G4_AIR /opnovice2/boxProperty ABSLENGTH 0.000002 1 0.000005 2 0.000008 3 The MPT for the box is now: 16: ABSLENGTH 2e-06 1 5e-06 2 8e-06 3 ............. /opnovice2/boxProperty RAYLEIGH 0.000002 1 0.000008 1 The MPT for the box is now: 11: RAYLEIGH 2e-06 1 8e-06 1 16: ABSLENGTH 2e-06 1 5e-06 2 8e-06 3 ............. /opnovice2/boxProperty RINDEX 0.000002 1.3 0.000008 1.4 The MPT for the box is now: 0: RINDEX 2e-06 1.3 8e-06 1.4 9: GROUPVEL 2e-06 218.486 8e-06 203.645 11: RAYLEIGH 2e-06 1 8e-06 1 16: ABSLENGTH 2e-06 1 5e-06 2 8e-06 3 ............. /opnovice2/boxProperty SCINTILLATIONCOMPONENT1 0.000002 1.0 0.000005 1.1 0.000008 1.3 The MPT for the box is now: 0: RINDEX 2e-06 1.3 8e-06 1.4 9: GROUPVEL 2e-06 218.486 8e-06 203.645 11: RAYLEIGH 2e-06 1 8e-06 1 16: ABSLENGTH 2e-06 1 5e-06 2 8e-06 3 23: SCINTILLATIONCOMPONENT1 2e-06 1 5e-06 1.1 8e-06 1.3 ............. /opnovice2/boxProperty SCINTILLATIONCOMPONENT2 0.000002 0.1 0.000003 0.2 0.000004 0.4 0.000005 0.6 0.000006 0.8 0.000007 0.9 .000008 1.0 The MPT for the box is now: 0: RINDEX 2e-06 1.3 8e-06 1.4 9: GROUPVEL 2e-06 218.486 8e-06 203.645 11: RAYLEIGH 2e-06 1 8e-06 1 16: ABSLENGTH 2e-06 1 5e-06 2 8e-06 3 23: SCINTILLATIONCOMPONENT1 2e-06 1 5e-06 1.1 8e-06 1.3 24: SCINTILLATIONCOMPONENT2 2e-06 0.1 3e-06 0.2 4e-06 0.4 5e-06 0.6 6e-06 0.8 7e-06 0.9 8e-06 1 ............. /opnovice2/boxProperty SCINTILLATIONCOMPONENT3 0.000002 0.2 0.000005 0.1 0.000008 0.05 The MPT for the box is now: 0: RINDEX 2e-06 1.3 8e-06 1.4 9: GROUPVEL 2e-06 218.486 8e-06 203.645 11: RAYLEIGH 2e-06 1 8e-06 1 16: ABSLENGTH 2e-06 1 5e-06 2 8e-06 3 23: SCINTILLATIONCOMPONENT1 2e-06 1 5e-06 1.1 8e-06 1.3 24: SCINTILLATIONCOMPONENT2 2e-06 0.1 3e-06 0.2 4e-06 0.4 5e-06 0.6 6e-06 0.8 7e-06 0.9 8e-06 1 25: SCINTILLATIONCOMPONENT3 2e-06 0.2 5e-06 0.1 8e-06 0.05 ............. /opnovice2/boxConstProperty SCINTILLATIONTIMECONSTANT1 20 The MPT for the box is now: 0: RINDEX 2e-06 1.3 8e-06 1.4 9: GROUPVEL 2e-06 218.486 8e-06 203.645 11: RAYLEIGH 2e-06 1 8e-06 1 16: ABSLENGTH 2e-06 1 5e-06 2 8e-06 3 23: SCINTILLATIONCOMPONENT1 2e-06 1 5e-06 1.1 8e-06 1.3 24: SCINTILLATIONCOMPONENT2 2e-06 0.1 3e-06 0.2 4e-06 0.4 5e-06 0.6 6e-06 0.8 7e-06 0.9 8e-06 1 25: SCINTILLATIONCOMPONENT3 2e-06 0.2 5e-06 0.1 8e-06 0.05 30: SCINTILLATIONTIMECONSTANT1 20 ............. /opnovice2/boxConstProperty SCINTILLATIONTIMECONSTANT2 100 The MPT for the box is now: 0: RINDEX 2e-06 1.3 8e-06 1.4 9: GROUPVEL 2e-06 218.486 8e-06 203.645 11: RAYLEIGH 2e-06 1 8e-06 1 16: ABSLENGTH 2e-06 1 5e-06 2 8e-06 3 23: SCINTILLATIONCOMPONENT1 2e-06 1 5e-06 1.1 8e-06 1.3 24: SCINTILLATIONCOMPONENT2 2e-06 0.1 3e-06 0.2 4e-06 0.4 5e-06 0.6 6e-06 0.8 7e-06 0.9 8e-06 1 25: SCINTILLATIONCOMPONENT3 2e-06 0.2 5e-06 0.1 8e-06 0.05 30: SCINTILLATIONTIMECONSTANT1 20 31: SCINTILLATIONTIMECONSTANT2 100 ............. /opnovice2/boxConstProperty SCINTILLATIONTIMECONSTANT3 200 The MPT for the box is now: 0: RINDEX 2e-06 1.3 8e-06 1.4 9: GROUPVEL 2e-06 218.486 8e-06 203.645 11: RAYLEIGH 2e-06 1 8e-06 1 16: ABSLENGTH 2e-06 1 5e-06 2 8e-06 3 23: SCINTILLATIONCOMPONENT1 2e-06 1 5e-06 1.1 8e-06 1.3 24: SCINTILLATIONCOMPONENT2 2e-06 0.1 3e-06 0.2 4e-06 0.4 5e-06 0.6 6e-06 0.8 7e-06 0.9 8e-06 1 25: SCINTILLATIONCOMPONENT3 2e-06 0.2 5e-06 0.1 8e-06 0.05 30: SCINTILLATIONTIMECONSTANT1 20 31: SCINTILLATIONTIMECONSTANT2 100 32: SCINTILLATIONTIMECONSTANT3 200 ............. /opnovice2/boxConstProperty SCINTILLATIONYIELD 5000.0 The MPT for the box is now: 0: RINDEX 2e-06 1.3 8e-06 1.4 9: GROUPVEL 2e-06 218.486 8e-06 203.645 11: RAYLEIGH 2e-06 1 8e-06 1 16: ABSLENGTH 2e-06 1 5e-06 2 8e-06 3 23: SCINTILLATIONCOMPONENT1 2e-06 1 5e-06 1.1 8e-06 1.3 24: SCINTILLATIONCOMPONENT2 2e-06 0.1 3e-06 0.2 4e-06 0.4 5e-06 0.6 6e-06 0.8 7e-06 0.9 8e-06 1 25: SCINTILLATIONCOMPONENT3 2e-06 0.2 5e-06 0.1 8e-06 0.05 10: SCINTILLATIONYIELD 5000 30: SCINTILLATIONTIMECONSTANT1 20 31: SCINTILLATIONTIMECONSTANT2 100 32: SCINTILLATIONTIMECONSTANT3 200 ............. /opnovice2/boxConstProperty SCINTILLATIONYIELD1 1.0 The MPT for the box is now: 0: RINDEX 2e-06 1.3 8e-06 1.4 9: GROUPVEL 2e-06 218.486 8e-06 203.645 11: RAYLEIGH 2e-06 1 8e-06 1 16: ABSLENGTH 2e-06 1 5e-06 2 8e-06 3 23: SCINTILLATIONCOMPONENT1 2e-06 1 5e-06 1.1 8e-06 1.3 24: SCINTILLATIONCOMPONENT2 2e-06 0.1 3e-06 0.2 4e-06 0.4 5e-06 0.6 6e-06 0.8 7e-06 0.9 8e-06 1 25: SCINTILLATIONCOMPONENT3 2e-06 0.2 5e-06 0.1 8e-06 0.05 10: SCINTILLATIONYIELD 5000 30: SCINTILLATIONTIMECONSTANT1 20 31: SCINTILLATIONTIMECONSTANT2 100 32: SCINTILLATIONTIMECONSTANT3 200 36: SCINTILLATIONYIELD1 1 ............. /opnovice2/boxConstProperty SCINTILLATIONYIELD2 1.0 The MPT for the box is now: 0: RINDEX 2e-06 1.3 8e-06 1.4 9: GROUPVEL 2e-06 218.486 8e-06 203.645 11: RAYLEIGH 2e-06 1 8e-06 1 16: ABSLENGTH 2e-06 1 5e-06 2 8e-06 3 23: SCINTILLATIONCOMPONENT1 2e-06 1 5e-06 1.1 8e-06 1.3 24: SCINTILLATIONCOMPONENT2 2e-06 0.1 3e-06 0.2 4e-06 0.4 5e-06 0.6 6e-06 0.8 7e-06 0.9 8e-06 1 25: SCINTILLATIONCOMPONENT3 2e-06 0.2 5e-06 0.1 8e-06 0.05 10: SCINTILLATIONYIELD 5000 30: SCINTILLATIONTIMECONSTANT1 20 31: SCINTILLATIONTIMECONSTANT2 100 32: SCINTILLATIONTIMECONSTANT3 200 36: SCINTILLATIONYIELD1 1 37: SCINTILLATIONYIELD2 1 ............. /opnovice2/boxConstProperty SCINTILLATIONYIELD3 0.1 The MPT for the box is now: 0: RINDEX 2e-06 1.3 8e-06 1.4 9: GROUPVEL 2e-06 218.486 8e-06 203.645 11: RAYLEIGH 2e-06 1 8e-06 1 16: ABSLENGTH 2e-06 1 5e-06 2 8e-06 3 23: SCINTILLATIONCOMPONENT1 2e-06 1 5e-06 1.1 8e-06 1.3 24: SCINTILLATIONCOMPONENT2 2e-06 0.1 3e-06 0.2 4e-06 0.4 5e-06 0.6 6e-06 0.8 7e-06 0.9 8e-06 1 25: SCINTILLATIONCOMPONENT3 2e-06 0.2 5e-06 0.1 8e-06 0.05 10: SCINTILLATIONYIELD 5000 30: SCINTILLATIONTIMECONSTANT1 20 31: SCINTILLATIONTIMECONSTANT2 100 32: SCINTILLATIONTIMECONSTANT3 200 36: SCINTILLATIONYIELD1 1 37: SCINTILLATIONYIELD2 1 38: SCINTILLATIONYIELD3 0.1 ............. /opnovice2/boxConstProperty RESOLUTIONSCALE 1 The MPT for the box is now: 0: RINDEX 2e-06 1.3 8e-06 1.4 9: GROUPVEL 2e-06 218.486 8e-06 203.645 11: RAYLEIGH 2e-06 1 8e-06 1 16: ABSLENGTH 2e-06 1 5e-06 2 8e-06 3 23: SCINTILLATIONCOMPONENT1 2e-06 1 5e-06 1.1 8e-06 1.3 24: SCINTILLATIONCOMPONENT2 2e-06 0.1 3e-06 0.2 4e-06 0.4 5e-06 0.6 6e-06 0.8 7e-06 0.9 8e-06 1 25: SCINTILLATIONCOMPONENT3 2e-06 0.2 5e-06 0.1 8e-06 0.05 10: SCINTILLATIONYIELD 5000 11: RESOLUTIONSCALE 1 30: SCINTILLATIONTIMECONSTANT1 20 31: SCINTILLATIONTIMECONSTANT2 100 32: SCINTILLATIONTIMECONSTANT3 200 36: SCINTILLATIONYIELD1 1 37: SCINTILLATIONYIELD2 1 38: SCINTILLATIONYIELD3 0.1 ............. /opnovice2/boxConstProperty SCINTILLATIONRISETIME1 3 The MPT for the box is now: 0: RINDEX 2e-06 1.3 8e-06 1.4 9: GROUPVEL 2e-06 218.486 8e-06 203.645 11: RAYLEIGH 2e-06 1 8e-06 1 16: ABSLENGTH 2e-06 1 5e-06 2 8e-06 3 23: SCINTILLATIONCOMPONENT1 2e-06 1 5e-06 1.1 8e-06 1.3 24: SCINTILLATIONCOMPONENT2 2e-06 0.1 3e-06 0.2 4e-06 0.4 5e-06 0.6 6e-06 0.8 7e-06 0.9 8e-06 1 25: SCINTILLATIONCOMPONENT3 2e-06 0.2 5e-06 0.1 8e-06 0.05 10: SCINTILLATIONYIELD 5000 11: RESOLUTIONSCALE 1 30: SCINTILLATIONTIMECONSTANT1 20 31: SCINTILLATIONTIMECONSTANT2 100 32: SCINTILLATIONTIMECONSTANT3 200 33: SCINTILLATIONRISETIME1 3 36: SCINTILLATIONYIELD1 1 37: SCINTILLATIONYIELD2 1 38: SCINTILLATIONYIELD3 0.1 ............. /opnovice2/boxConstProperty SCINTILLATIONRISETIME2 10 The MPT for the box is now: 0: RINDEX 2e-06 1.3 8e-06 1.4 9: GROUPVEL 2e-06 218.486 8e-06 203.645 11: RAYLEIGH 2e-06 1 8e-06 1 16: ABSLENGTH 2e-06 1 5e-06 2 8e-06 3 23: SCINTILLATIONCOMPONENT1 2e-06 1 5e-06 1.1 8e-06 1.3 24: SCINTILLATIONCOMPONENT2 2e-06 0.1 3e-06 0.2 4e-06 0.4 5e-06 0.6 6e-06 0.8 7e-06 0.9 8e-06 1 25: SCINTILLATIONCOMPONENT3 2e-06 0.2 5e-06 0.1 8e-06 0.05 10: SCINTILLATIONYIELD 5000 11: RESOLUTIONSCALE 1 30: SCINTILLATIONTIMECONSTANT1 20 31: SCINTILLATIONTIMECONSTANT2 100 32: SCINTILLATIONTIMECONSTANT3 200 33: SCINTILLATIONRISETIME1 3 34: SCINTILLATIONRISETIME2 10 36: SCINTILLATIONYIELD1 1 37: SCINTILLATIONYIELD2 1 38: SCINTILLATIONYIELD3 0.1 ............. /opnovice2/boxConstProperty SCINTILLATIONRISETIME3 20 The MPT for the box is now: 0: RINDEX 2e-06 1.3 8e-06 1.4 9: GROUPVEL 2e-06 218.486 8e-06 203.645 11: RAYLEIGH 2e-06 1 8e-06 1 16: ABSLENGTH 2e-06 1 5e-06 2 8e-06 3 23: SCINTILLATIONCOMPONENT1 2e-06 1 5e-06 1.1 8e-06 1.3 24: SCINTILLATIONCOMPONENT2 2e-06 0.1 3e-06 0.2 4e-06 0.4 5e-06 0.6 6e-06 0.8 7e-06 0.9 8e-06 1 25: SCINTILLATIONCOMPONENT3 2e-06 0.2 5e-06 0.1 8e-06 0.05 10: SCINTILLATIONYIELD 5000 11: RESOLUTIONSCALE 1 30: SCINTILLATIONTIMECONSTANT1 20 31: SCINTILLATIONTIMECONSTANT2 100 32: SCINTILLATIONTIMECONSTANT3 200 33: SCINTILLATIONRISETIME1 3 34: SCINTILLATIONRISETIME2 10 35: SCINTILLATIONRISETIME3 20 36: SCINTILLATIONYIELD1 1 37: SCINTILLATIONYIELD2 1 38: SCINTILLATIONYIELD3 0.1 ............. /opnovice2/worldProperty RINDEX 0.000002 1.01 0.000008 1.01 The MPT for the world is now: 0: RINDEX 2e-06 1.01 8e-06 1.01 9: GROUPVEL 2e-06 296.824 8e-06 296.824 ............. /opnovice2/worldProperty ABSLENGTH 0.000002 100 0.000005 100 0.000008 100 The MPT for the world is now: 0: RINDEX 2e-06 1.01 8e-06 1.01 9: GROUPVEL 2e-06 296.824 8e-06 296.824 16: ABSLENGTH 2e-06 100 5e-06 100 8e-06 100 ............. /opnovice2/surfaceFinish ground /run/geometryModified /opnovice2/surfaceSigmaAlpha 0.2 /run/geometryModified Surface sigma alpha set to: 0.2 /opnovice2/surfaceProperty SPECULARLOBECONSTANT 0.000002 0.1 0.000008 0.1 SPECULARLOBECONSTANT 0.000002 0.1 0.000008 0.1 The MPT for the surface is now: 6: SPECULARLOBECONSTANT 2e-06 0.1 8e-06 0.1 ............. /opnovice2/surfaceProperty SPECULARSPIKECONSTANT 0.000002 0.1 0.000008 0.1 SPECULARSPIKECONSTANT 0.000002 0.1 0.000008 0.1 The MPT for the surface is now: 6: SPECULARLOBECONSTANT 2e-06 0.1 8e-06 0.1 7: SPECULARSPIKECONSTANT 2e-06 0.1 8e-06 0.1 ............. /opnovice2/surfaceProperty BACKSCATTERCONSTANT 0.000002 0.1 0.000008 0.1 BACKSCATTERCONSTANT 0.000002 0.1 0.000008 0.1 The MPT for the surface is now: 6: SPECULARLOBECONSTANT 2e-06 0.1 8e-06 0.1 7: SPECULARSPIKECONSTANT 2e-06 0.1 8e-06 0.1 8: BACKSCATTERCONSTANT 2e-06 0.1 8e-06 0.1 ............. /opnovice2/surfaceProperty TRANSMITTANCE 0.000002 0.1 0.000008 0.1 TRANSMITTANCE 0.000002 0.1 0.000008 0.1 The MPT for the surface is now: 5: TRANSMITTANCE 2e-06 0.1 8e-06 0.1 6: SPECULARLOBECONSTANT 2e-06 0.1 8e-06 0.1 7: SPECULARSPIKECONSTANT 2e-06 0.1 8e-06 0.1 8: BACKSCATTERCONSTANT 2e-06 0.1 8e-06 0.1 ............. /opnovice2/surfaceProperty REFLECTIVITY 0.000002 0.8 0.000008 0.8 REFLECTIVITY 0.000002 0.8 0.000008 0.8 The MPT for the surface is now: 1: REFLECTIVITY 2e-06 0.8 8e-06 0.8 5: TRANSMITTANCE 2e-06 0.1 8e-06 0.1 6: SPECULARLOBECONSTANT 2e-06 0.1 8e-06 0.1 7: SPECULARSPIKECONSTANT 2e-06 0.1 8e-06 0.1 8: BACKSCATTERCONSTANT 2e-06 0.1 8e-06 0.1 ............. /opnovice2/surfaceProperty EFFICIENCY 0.000002 0.1 0.000008 0.1 EFFICIENCY 0.000002 0.1 0.000008 0.1 The MPT for the surface is now: 1: REFLECTIVITY 2e-06 0.8 8e-06 0.8 4: EFFICIENCY 2e-06 0.1 8e-06 0.1 5: TRANSMITTANCE 2e-06 0.1 8e-06 0.1 6: SPECULARLOBECONSTANT 2e-06 0.1 8e-06 0.1 7: SPECULARSPIKECONSTANT 2e-06 0.1 8e-06 0.1 8: BACKSCATTERCONSTANT 2e-06 0.1 8e-06 0.1 ............. /process/optical/cerenkov/setMaxPhotons 3 Cerenkov max photons: 3 /run/physicsModified /process/optical/cerenkov/setMaxBetaChange 10 /run/physicsModified /process/optical/scintillation/setByParticleType false /run/physicsModified /process/optical/scintillation/setTrackInfo false /run/physicsModified /process/optical/scintillation/setFiniteRiseTime true /run/physicsModified /process/optical/scintillation/setStackPhotons true /run/physicsModified /run/initialize ****** opticalSurface->DumpInfo: Surface type = 1 Surface finish = 3 Surface model = 1 Surface parameter ----------------- sigma_alpha: 0.2 ****** end of opticalSurface->DumpInfo hInelastic FTFP_BERT : threshold between BERT and FTFP is over the interval for pions : 3 to 6 GeV for kaons : 3 to 6 GeV for proton : 3 to 6 GeV for neutron : 3 to 6 GeV ### Adding tracking cuts for neutron TimeCut(ns)= 10000 KinEnergyCut(MeV)= 0 OpAbsorption is created OpRayleigh is created OpMieHG is created OpBoundary is created OpWLS is created OpWLS2 is created Cerenkov is created. ### Birks coefficients used in run time G4_BGO 0.126 mm/MeV 0.089838 g/cm^2/MeV massFactor= 40.0239 effCharge= 1652.42 /analysis/setFileName electron /analysis/h1/set 1 100 0 10 /analysis/h1/set 2 100 0 10 /analysis/h1/set 3 400 0 200 /gun/particle e- /gun/energy 1 MeV /gun/position -1 0 0 m /gun/direction 1 0 0 # /run/printProgress 100 /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 1 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 ==================================================================== HADRONIC PROCESSES SUMMARY (verbose level 1) ----------------------------------------------------------------------- Hadronic Processes for neutron Process: hadElastic Model: hElasticCHIPS: 0 eV ---> 100 TeV Cr_sctns: G4NeutronElasticXS: 0 eV ---> 100 TeV Process: neutronInelastic Model: FTFP: 3 GeV ---> 100 TeV Model: BertiniCascade: 0 eV ---> 6 GeV Cr_sctns: G4NeutronInelasticXS: 0 eV ---> 100 TeV Process: nCapture Model: nRadCapture: 0 eV ---> 100 TeV Cr_sctns: G4NeutronCaptureXS: 0 eV ---> 100 TeV Process: nKiller ----------------------------------------------------------------------- Hadronic Processes for B- Process: hadElastic Model: hElasticLHEP: 0 eV ---> 100 TeV Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV Process: B-Inelastic Model: FTFP: 0 eV ---> 100 TeV Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV ----------------------------------------------------------------------- Hadronic Processes for D- Process: hadElastic Model: hElasticLHEP: 0 eV ---> 100 TeV Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV Process: D-Inelastic Model: FTFP: 0 eV ---> 100 TeV Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV ----------------------------------------------------------------------- Hadronic Processes for GenericIon Process: ionInelastic Model: Binary Light Ion Cascade: 0 eV /n ---> 6 GeV/n Model: FTFP: 3 GeV/n ---> 100 TeV/n Cr_sctns: Glauber-Gribov Nucl-nucl: 0 eV ---> 25.6 PeV ----------------------------------------------------------------------- Hadronic Processes for He3 Process: hadElastic Model: hElasticLHEP: 0 eV /n ---> 100 TeV/n Cr_sctns: Glauber-Gribov Nucl-nucl: 0 eV ---> 25.6 PeV Process: He3Inelastic Model: Binary Light Ion Cascade: 0 eV /n ---> 6 GeV/n Model: FTFP: 3 GeV/n ---> 100 TeV/n Cr_sctns: Glauber-Gribov Nucl-nucl: 0 eV ---> 25.6 PeV ----------------------------------------------------------------------- Hadronic Processes for alpha Process: hadElastic Model: hElasticLHEP: 0 eV /n ---> 100 TeV/n Cr_sctns: Glauber-Gribov Nucl-nucl: 0 eV ---> 25.6 PeV Process: alphaInelastic Model: Binary Light Ion Cascade: 0 eV /n ---> 6 GeV/n Model: FTFP: 3 GeV/n ---> 100 TeV/n Cr_sctns: Glauber-Gribov Nucl-nucl: 0 eV ---> 25.6 PeV ----------------------------------------------------------------------- Hadronic Processes for anti_He3 Process: hadElastic Model: hElasticLHEP: 0 eV /n ---> 100.1 MeV/n Model: AntiAElastic: 100 MeV/n ---> 100 TeV/n Cr_sctns: AntiAGlauber: 0 eV ---> 25.6 PeV Process: anti_He3Inelastic Model: FTFP: 0 eV /n ---> 100 TeV/n Cr_sctns: AntiAGlauber: 0 eV ---> 25.6 PeV Process: hFritiofCaptureAtRest ----------------------------------------------------------------------- Hadronic Processes for anti_alpha Process: hadElastic Model: hElasticLHEP: 0 eV /n ---> 100.1 MeV/n Model: AntiAElastic: 100 MeV/n ---> 100 TeV/n Cr_sctns: AntiAGlauber: 0 eV ---> 25.6 PeV Process: anti_alphaInelastic Model: FTFP: 0 eV /n ---> 100 TeV/n Cr_sctns: AntiAGlauber: 0 eV ---> 25.6 PeV Process: hFritiofCaptureAtRest ----------------------------------------------------------------------- Hadronic Processes for anti_deuteron Process: hadElastic Model: hElasticLHEP: 0 eV /n ---> 100.1 MeV/n Model: AntiAElastic: 100 MeV/n ---> 100 TeV/n Cr_sctns: AntiAGlauber: 0 eV ---> 25.6 PeV Process: anti_deuteronInelastic Model: FTFP: 0 eV /n ---> 100 TeV/n Cr_sctns: AntiAGlauber: 0 eV ---> 25.6 PeV Process: hFritiofCaptureAtRest ------------------------------------------------------------------------- Hadronic Processes for anti_hypertriton Process: hFritiofCaptureAtRest ----------------------------------------------------------------------- Hadronic Processes for anti_lambda Process: hadElastic Model: hElasticLHEP: 0 eV ---> 100 TeV Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV Process: anti_lambdaInelastic Model: FTFP: 0 eV ---> 100 TeV Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV Process: hFritiofCaptureAtRest ----------------------------------------------------------------------- Hadronic Processes for anti_neutron Process: hadElastic Model: hElasticLHEP: 0 eV ---> 100.1 MeV Model: AntiAElastic: 100 MeV ---> 100 TeV Cr_sctns: AntiAGlauber: 0 eV ---> 25.6 PeV Process: anti_neutronInelastic Model: FTFP: 0 eV ---> 100 TeV Cr_sctns: AntiAGlauber: 0 eV ---> 25.6 PeV Process: hFritiofCaptureAtRest ----------------------------------------------------------------------- Hadronic Processes for anti_proton Process: hadElastic Model: hElasticLHEP: 0 eV ---> 100.1 MeV Model: AntiAElastic: 100 MeV ---> 100 TeV Cr_sctns: AntiAGlauber: 0 eV ---> 25.6 PeV Process: anti_protonInelastic Model: FTFP: 0 eV ---> 100 TeV Cr_sctns: AntiAGlauber: 0 eV ---> 25.6 PeV Process: hFritiofCaptureAtRest ----------------------------------------------------------------------- Hadronic Processes for anti_triton Process: hadElastic Model: hElasticLHEP: 0 eV /n ---> 100.1 MeV/n Model: AntiAElastic: 100 MeV/n ---> 100 TeV/n Cr_sctns: AntiAGlauber: 0 eV ---> 25.6 PeV Process: anti_tritonInelastic Model: FTFP: 0 eV /n ---> 100 TeV/n Cr_sctns: AntiAGlauber: 0 eV ---> 25.6 PeV Process: hFritiofCaptureAtRest ----------------------------------------------------------------------- Hadronic Processes for deuteron Process: hadElastic Model: hElasticLHEP: 0 eV /n ---> 100 TeV/n Cr_sctns: Glauber-Gribov Nucl-nucl: 0 eV ---> 25.6 PeV Process: dInelastic Model: Binary Light Ion Cascade: 0 eV /n ---> 6 GeV/n Model: FTFP: 3 GeV/n ---> 100 TeV/n Cr_sctns: Glauber-Gribov Nucl-nucl: 0 eV ---> 25.6 PeV ----------------------------------------------------------------------- Hadronic Processes for e+ Process: positronNuclear Model: G4ElectroVDNuclearModel: 0 eV ---> 1 PeV Cr_sctns: ElectroNuclearXS: 0 eV ---> 100 TeV ----------------------------------------------------------------------- Hadronic Processes for e- Process: electronNuclear Model: G4ElectroVDNuclearModel: 0 eV ---> 1 PeV Cr_sctns: ElectroNuclearXS: 0 eV ---> 100 TeV ----------------------------------------------------------------------- Hadronic Processes for gamma Process: photonNuclear Model: GammaNPreco: 0 eV ---> 200 MeV Model: BertiniCascade: 199 MeV ---> 6 GeV Model: TheoFSGenerator: 3 GeV ---> 100 TeV Cr_sctns: GammaNuclearXS: 0 eV ---> 100 TeV ----------------------------------------------------------------------- Hadronic Processes for kaon+ Process: hadElastic Model: hElasticLHEP: 0 eV ---> 100 TeV Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV Process: kaon+Inelastic Model: FTFP: 3 GeV ---> 100 TeV Model: BertiniCascade: 0 eV ---> 6 GeV Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV ----------------------------------------------------------------------- Hadronic Processes for kaon- Process: hadElastic Model: hElasticLHEP: 0 eV ---> 100 TeV Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV Process: kaon-Inelastic Model: FTFP: 3 GeV ---> 100 TeV Model: BertiniCascade: 0 eV ---> 6 GeV Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV Process: hBertiniCaptureAtRest ----------------------------------------------------------------------- Hadronic Processes for lambda Process: hadElastic Model: hElasticLHEP: 0 eV ---> 100 TeV Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV Process: lambdaInelastic Model: FTFP: 3 GeV ---> 100 TeV Model: BertiniCascade: 0 eV ---> 6 GeV Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV ----------------------------------------------------------------------- Hadronic Processes for mu+ Process: muonNuclear Model: G4MuonVDNuclearModel: 0 eV ---> 1 PeV Cr_sctns: KokoulinMuonNuclearXS: 0 eV ---> 100 TeV ----------------------------------------------------------------------- Hadronic Processes for mu- Process: muonNuclear Model: G4MuonVDNuclearModel: 0 eV ---> 1 PeV Cr_sctns: KokoulinMuonNuclearXS: 0 eV ---> 100 TeV Process: muMinusCaptureAtRest ----------------------------------------------------------------------- Hadronic Processes for pi+ Process: hadElastic Model: hElasticGlauber: 0 eV ---> 100 TeV Cr_sctns: BarashenkovGlauberGribov: 0 eV ---> 100 TeV Process: pi+Inelastic Model: FTFP: 3 GeV ---> 100 TeV Model: BertiniCascade: 0 eV ---> 6 GeV Cr_sctns: BarashenkovGlauberGribov: 0 eV ---> 100 TeV ----------------------------------------------------------------------- Hadronic Processes for pi- Process: hadElastic Model: hElasticGlauber: 0 eV ---> 100 TeV Cr_sctns: BarashenkovGlauberGribov: 0 eV ---> 100 TeV Process: pi-Inelastic Model: FTFP: 3 GeV ---> 100 TeV Model: BertiniCascade: 0 eV ---> 6 GeV Cr_sctns: BarashenkovGlauberGribov: 0 eV ---> 100 TeV Process: hBertiniCaptureAtRest ----------------------------------------------------------------------- Hadronic Processes for proton Process: hadElastic Model: hElasticCHIPS: 0 eV ---> 100 TeV Cr_sctns: BarashenkovGlauberGribov: 0 eV ---> 100 TeV Process: protonInelastic Model: FTFP: 3 GeV ---> 100 TeV Model: BertiniCascade: 0 eV ---> 6 GeV Cr_sctns: BarashenkovGlauberGribov: 0 eV ---> 100 TeV ----------------------------------------------------------------------- Hadronic Processes for sigma- Process: hadElastic Model: hElasticLHEP: 0 eV ---> 100 TeV Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV Process: sigma-Inelastic Model: FTFP: 3 GeV ---> 100 TeV Model: BertiniCascade: 0 eV ---> 6 GeV Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV Process: hBertiniCaptureAtRest ----------------------------------------------------------------------- Hadronic Processes for triton Process: hadElastic Model: hElasticLHEP: 0 eV /n ---> 100 TeV/n Cr_sctns: Glauber-Gribov Nucl-nucl: 0 eV ---> 25.6 PeV Process: tInelastic Model: Binary Light Ion Cascade: 0 eV /n ---> 6 GeV/n Model: FTFP: 3 GeV/n ---> 100 TeV/n Cr_sctns: Glauber-Gribov Nucl-nucl: 0 eV ---> 25.6 PeV ======================================================================= ====== Geant4 Native Pre-compound Model Parameters ======== ======================================================================= Type of pre-compound model 0 Type of pre-compound inverse x-section 1 Pre-compound model active 1 Pre-compound excitation low energy 0.1 MeV Pre-compound excitation high energy 15 MeV Angular generator for pre-compound model 1 Use NeverGoBack option for pre-compound model 0 Use SoftCutOff option for pre-compound model 0 Use CEM transitions for pre-compound model 1 Use GNASH transitions for pre-compound model 0 Use HETC submodel for pre-compound model 0 ======================================================================= ====== Nuclear De-excitation Module Parameters ======== ======================================================================= Type of de-excitation inverse x-section 3 Type of de-excitation factory Evaporation+GEM Number of de-excitation channels 68 Type of Fermi BreakUp model ModelVI Min excitation energy 0.01 keV Min energy per nucleon for multifragmentation 2e+05 MeV Level density (1/MeV) 0.075 Use simple level density model 1 Use discrete excitation energy of the residual 0 Time limit for long lived isomeres 1 ns Isomer production flag 1 Internal e- conversion flag 1 Store e- internal conversion data 1 Correlated gamma emission flag 0 Max 2J for sampling of angular correlations 10 ======================================================================= ========= Table of registered couples ============================ Index : 0 used in the geometry : Yes Material : G4_AIR Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 70 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 1 used in the geometry : Yes Material : G4_BGO Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um Energy thresholds : gamma 64.4521 keV e- 770.195 keV e+ 731.433 keV proton 70 keV Region(s) which use this couple : DefaultRegionForTheWorld ================================================================== ### Run 0 starts. ... set ntuple merging row mode : row-wise - done ... create file : electron.root - done ... open analysis file : electron.root - done ... open analysis file : electron.root - done --> Event 0 starts. Run terminated. Run Summary Number of events processed : 100 User=1.320000s Real=1.319979s Sys=0.000000s Run Summary --------------------------------- Primary particle was: e- with energy 1 MeV. Number of events: 100 Material of world: G4_AIR Material of tank: G4_BGO Average energy of Cerenkov photons created per event: 2.3e+02 eV. Average number of Cerenkov photons created per event: 45 Average energy per photon: 5.2 eV. Average energy of scintillation photons created per event: 1.7e+04 eV. Average number of scintillation photons created per event: 3.1e+03 Average energy per photon: 5.4 eV. Average number of photons absorbed by WLS per event: 0 Average number of photons created by WLS per event: 0 Average energy of WLS photons created per event: 0 eV. Average number of photons absorbed by WLS2 per event: 0 Average number of photons created by WLS2 per event: 0 Average energy of WLS2 photons created per event: 0 eV. Average number of OpRayleigh per event: 4.8e+03 Average number of OpAbsorption per event: 3e+03 Surface events (on +X surface, maximum one per photon) this run: # of primary particles: 1e+02 OpAbsorption before surface: 168802 Total # of surface events: 147495 Surface events by process: Transmission: 14801 Fresnel refraction: 48242 Lambertian reflection: 49321 Lobe reflection: 6148 Spike reflection: 7111 Backscattering: 7078 Absorption: 13376 Detection: 1418 Sum: 147495 Unaccounted for: 0 --------------------------------- Histogram statistics for the entire run: Cerenkov spectrum: mean = 5.19026 eV rms = 1.73685 eV. Scintillation spectrum: mean = 5.4313 eV rms = 1.66662 eV. Scintillation time: mean = 48.407 ns rms = 46.9642 ns. ... write file : electron.root - done ... close file : electron.root - done Graphics systems deleted. Visualization Manager deleting... ================== Deleting memory pools =================== Number of memory pools allocated: 10 of which, static: 0 Dynamic pools deleted: 10 / Total memory freed: 0.46 MB ============================================================