Files
2026-03-19 16:51:22 +01:00

988 lines
45 KiB
Plaintext

############################################
!!! 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/
**************************************************************
/run/verbose 1
#/tracking/verbose 1
#
# select EM PhysicsList
#
/testem/phys/addPhysics emstandard_opt1
PhysicsList::Set emstandard_opt1 EM physics
#
/process/em/AddEmRegion VertexDetector G4EmStandard_opt4
/process/em/AddEmRegion MuonDetector G4EmStandard
#
/run/initialize
Ecal is 36 cm of G4_CESIUM_IODIDE
Absorber is 2 mm of G4_Al
Vertex is 3 cm of 3 layers of Si of 100 um npads= 601
### New geometry is constructed
#
/run/setCut 0.3 mm
/run/setCutForRegion VertexDetector 0.01 mm
/run/setCutForRegion MuonDetector 10 mm
/testem/stepMax 100 mm
#
/testem/det/acceptance1 0.8271 0.02974 2.5
/testem/det/acceptance9 0.9570 0.01286 2.5
/testem/det/acceptance25 0.9793 0.00884 2.5
#
/gun/particle e-
/gun/energy 1 GeV
/run/beamOn 100
=======================================================================
====== Electromagnetic Physics Parameters ========
=======================================================================
LPM effect enabled 1
Enable creation and use of sampling tables 0
Apply cuts on all EM processes 1
Use combined TransportationWithMsc Enabled
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 SimplePositronium
Enable 3 gamma annihilation on fly 0
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.8, 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
Use ICRU90 data 0
Fluctuations of dE/dx are enabled 1
Type of fluctuation model for leptons and hadrons Universal
Use built-in Birks saturation 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+- 0
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.2
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 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 3 regions:
DefaultRegionForTheWorld 0 0 0
VertexDetector 1 0 0
MuonDetector 0 0 0
### === Ignore cuts flag: 0
phot: for gamma applyCuts:1 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 applyCuts:1 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 ======
Klein-Nishina : Emin= 0 eV Emax= 100 TeV
===== EM models for the G4Region VertexDetector ======
LowEPComptonModel : Emin= 0 eV Emax= 20 MeV Fluo
KleinNishina : Emin= 20 MeV Emax= 100 TeV Fluo
conv: for gamma applyCuts:1 SubType=14 BuildTable=1
Lambda table from 1.022 MeV to 100 TeV, 20 bins/decade, spline: 1
===== EM models for the G4Region DefaultRegionForTheWorld ======
BetheHeitlerLPM : Emin= 0 eV Emax= 100 TeV ModifiedTsai
===== EM models for the G4Region VertexDetector ======
BetheHeitler5D : Emin= 0 eV Emax= 100 TeV ModifiedTsai
Rayl: for gamma applyCuts:1 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 ======
DummyModel : Emin= 0 eV Emax= 100 TeV
StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
===== EM models for the G4Region VertexDetector ======
LivermoreRayleigh : Emin= 0 eV Emax= 100 TeV CullenGenerator
TransportationWithMsc: for e-
===== EM models for the G4Region DefaultRegionForTheWorld ======
UrbanMsc : Emin= 0 eV Emax= 100 MeV
StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
WentzelVIUni : Emin= 100 MeV Emax= 100 TeV
StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
===== EM models for the G4Region VertexDetector ======
GoudsmitSaunderson : Emin= 0 eV Emax= 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
StepLim=SafetyPlus Rfact=0.08 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
===== EM models for the G4Region MuonDetector ======
UrbanMsc : Emin= 0 eV Emax= 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
StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 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.8, 1 mm), integ: 3, 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 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 ModifiedTsai
eBremLPM : Emin= 1 GeV Emax= 100 TeV ModifiedTsai
CoulombScat: for e- XStype:1 applyCuts: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
===== EM models for the G4Region VertexDetector ======
eCoulombScattering : Emin= 100 MeV Emax= 100 TeV
===== EM models for the G4Region MuonDetector ======
eCoulombScattering : Emin= 100 MeV Emax= 100 TeV
TransportationWithMsc: for e+
===== EM models for the G4Region DefaultRegionForTheWorld ======
UrbanMsc : Emin= 0 eV Emax= 100 MeV
StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
WentzelVIUni : Emin= 100 MeV Emax= 100 TeV
StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
===== EM models for the G4Region VertexDetector ======
GoudsmitSaunderson : Emin= 0 eV Emax= 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
StepLim=SafetyPlus Rfact=0.08 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
===== EM models for the G4Region MuonDetector ======
UrbanMsc : Emin= 0 eV Emax= 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
StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 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.8, 1 mm), integ: 3, 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 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 ModifiedTsai
eBremLPM : Emin= 1 GeV Emax= 100 TeV ModifiedTsai
annihil: for e+ XStype:2 applyCuts:1 SubType=5 AtRestModel:Simple BuildTable=0
===== EM models for the G4Region DefaultRegionForTheWorld ======
eplus2gg : Emin= 0 eV Emax= 100 TeV
CoulombScat: for e+ XStype:1 applyCuts: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
===== EM models for the G4Region VertexDetector ======
eCoulombScattering : Emin= 100 MeV Emax= 100 TeV
===== EM models for the G4Region MuonDetector ======
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:0 Skin=1 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.2, 0.1 mm), integ: 3, 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 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 applyCuts: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 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: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.1 mm), integ: 3, fluct: 1, linLossLim= 0.02
===== EM models for the G4Region DefaultRegionForTheWorld ======
Bragg : Emin= 0 eV Emax= 2 MeV
BetheBloch : Emin= 2 MeV Emax= 100 TeV
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: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.1 mm), integ: 3, 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=240 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: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.1 mm), integ: 3, 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 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 applyCuts: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:0 Skin=1 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.2, 0.1 mm), integ: 3, 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 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 applyCuts: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:0 Skin=1 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.2, 0.1 mm), integ: 3, 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 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 applyCuts: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:0 Skin=1 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.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
===== EM models for the G4Region DefaultRegionForTheWorld ======
Bragg : Emin= 0 eV Emax= 200 keV
MuBetheBloch : Emin= 200 keV Emax= 100 TeV
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 applyCuts: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:0 Skin=1 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.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
===== EM models for the G4Region DefaultRegionForTheWorld ======
ICRU73QO : Emin= 0 eV Emax= 200 keV
MuBetheBloch : Emin= 200 keV Emax= 100 TeV
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 applyCuts: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:0 Skin=1 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.2, 0.1 mm), integ: 3, 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 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 applyCuts: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:0 Skin=1 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.2, 0.1 mm), integ: 3, 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 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 applyCuts: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 : G4_AIR
Range cuts : gamma 300 um e- 300 um e+ 300 um proton 300 um
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 30 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 1 used in the geometry : Yes
Material : G4_CESIUM_IODIDE
Range cuts : gamma 300 um e- 300 um e+ 300 um proton 300 um
Energy thresholds : gamma 18.2699 keV e- 291.571 keV e+ 281.583 keV proton 30 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 2 used in the geometry : Yes
Material : G4_Al
Range cuts : gamma 300 um e- 300 um e+ 300 um proton 300 um
Energy thresholds : gamma 3.9261 keV e- 260.03 keV e+ 253.129 keV proton 30 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 3 used in the geometry : Yes
Material : G4_AIR
Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 1 keV
Region(s) which use this couple :
VertexDetector
Index : 4 used in the geometry : Yes
Material : G4_Al
Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um
Energy thresholds : gamma 990 eV e- 34.1058 keV e+ 33.7225 keV proton 1 keV
Region(s) which use this couple :
VertexDetector
Index : 5 used in the geometry : Yes
Material : G4_Si
Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um
Energy thresholds : gamma 990 eV e- 31.6994 keV e+ 31.3501 keV proton 1 keV
Region(s) which use this couple :
VertexDetector
Index : 6 used in the geometry : Yes
Material : G4_Fe
Range cuts : gamma 1 cm e- 1 cm e+ 1 cm proton 1 cm
Energy thresholds : gamma 60.6431 keV e- 13.0433 MeV e+ 12.1452 MeV proton 1 MeV
Region(s) which use this couple :
MuonDetector
Index : 7 used in the geometry : Yes
Material : G4_Al
Range cuts : gamma 1 cm e- 1 cm e+ 1 cm proton 1 cm
Energy thresholds : gamma 20.4222 keV e- 4.56335 MeV e+ 4.27483 MeV proton 1 MeV
Region(s) which use this couple :
MuonDetector
==================================================================
### Run 0 starts.
### Run 0 start
### PrimaryGeneratorAction::GeneratePrimaries ##########
### e- E(MeV)= 1000
########################################################
Run terminated.
Run Summary
Number of events processed : 100
User=0.310000s Real=0.311598s Sys=0.000000s
RunAction: End of run actions are started
HistoManager: End of run actions are started RunID# 0
=================================================================
Number of events 100
Average number of e- 403.7
Average number of gamma 635.5
Average number of e+ 47.22
Average number of steps 2946
Edep 1x1 = 0.8329 +- 0.003391 res= 4.071 % 100
Edep 3x3 = 0.959 +- 0.001318 res= 1.374 % 100
Edep 5x5 = 0.9796 +- 0.0009483 res= 0.9681 % 100
=========== Ratios without trancating ===========================
E1/E9 = 0.8685 +- 0.003149
E1/E25 = 0.8503 +- 0.003311
E9/E25 = 0.979 +- 0.0008637
Beam Energy 1 GeV
==================================================================
<<<<<ACCEPTANCE>>>>> 100 events for Crystal Calorimeter
Edep1x1: 0.8329 delEdep1x1= -0.1671 nrms= -1.671
Erms1x1: 0.03391 delErms1x1= -0.9661 nrms= -9.661
Edep3x3: 0.959 delEdep3x3= -0.04099 nrms= -0.4099
Erms3x3: 0.01318 delErms3x3= -0.9868 nrms= -9.868
Edep5x5: 0.9796 delEdep5x5= -0.02042 nrms= -0.2042
Erms5x5: 0.009483 delErms5x5= -0.9905 nrms= -9.905
<<<<<END>>>>> IS ACCEPTED
Z bremsstrahlung photoeffect compton conversion
13 0 0 160 0
53 261 42 0 44
55 278 48 0 48
#
/testem/det/acceptance1 0.0136 0.0001 500
/testem/det/acceptance9 0.0139 0.0001 500
/testem/det/acceptance25 0.0139 0.0001 500
#
/gun/particle mu-
/gun/energy 20 GeV
/run/beamOn 1000
========= Table of registered couples ============================
Index : 0 used in the geometry : Yes
Material : G4_AIR
Range cuts : gamma 300 um e- 300 um e+ 300 um proton 300 um
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 30 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 1 used in the geometry : Yes
Material : G4_CESIUM_IODIDE
Range cuts : gamma 300 um e- 300 um e+ 300 um proton 300 um
Energy thresholds : gamma 18.27 keV e- 291.6 keV e+ 281.6 keV proton 30 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 2 used in the geometry : Yes
Material : G4_Al
Range cuts : gamma 300 um e- 300 um e+ 300 um proton 300 um
Energy thresholds : gamma 3.926 keV e- 260 keV e+ 253.1 keV proton 30 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 3 used in the geometry : Yes
Material : G4_AIR
Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 1 keV
Region(s) which use this couple :
VertexDetector
Index : 4 used in the geometry : Yes
Material : G4_Al
Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um
Energy thresholds : gamma 990 eV e- 34.11 keV e+ 33.72 keV proton 1 keV
Region(s) which use this couple :
VertexDetector
Index : 5 used in the geometry : Yes
Material : G4_Si
Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um
Energy thresholds : gamma 990 eV e- 31.7 keV e+ 31.35 keV proton 1 keV
Region(s) which use this couple :
VertexDetector
Index : 6 used in the geometry : Yes
Material : G4_Fe
Range cuts : gamma 1 cm e- 1 cm e+ 1 cm proton 1 cm
Energy thresholds : gamma 60.64 keV e- 13.04 MeV e+ 12.15 MeV proton 1 MeV
Region(s) which use this couple :
MuonDetector
Index : 7 used in the geometry : Yes
Material : G4_Al
Range cuts : gamma 1 cm e- 1 cm e+ 1 cm proton 1 cm
Energy thresholds : gamma 20.42 keV e- 4.563 MeV e+ 4.275 MeV proton 1 MeV
Region(s) which use this couple :
MuonDetector
==================================================================
### Run 1 starts.
### Run 1 start
### PrimaryGeneratorAction::GeneratePrimaries ##########
### mu- E(MeV)= 2e+04
########################################################
Run terminated.
Run Summary
Number of events processed : 1000
User=0.470000s Real=0.472948s Sys=0.000000s
RunAction: End of run actions are started
HistoManager: End of run actions are started RunID# 1
=================================================================
Number of events 1000
Average number of e- 69.02
Average number of gamma 61.77
Average number of e+ 3.389
Average number of steps 421.8
Edep 1x1 = 0.0141 +- 0.0003603
Edep 3x3 = 0.01445 +- 0.0004103
Edep 5x5 = 0.01451 +- 0.0004186
=========== Ratios without trancating ===========================
E1/E9 = 0.984 +- 0.0006927
E1/E25 = 0.9811 +- 0.0007849
E9/E25 = 0.997 +- 0.0001965
Beam Energy 20 GeV
==================================================================
#
/process/em/setSecBiasing eBrem World 0.1 2 MeV
/run/physicsModified
/process/em/setSecBiasing eIoni World 0.0 1 MeV
/run/physicsModified
/process/em/setSecBiasing compt World 0.0 1 MeV
/run/physicsModified
/process/em/setSecBiasing phot World 0.0 1 MeV
/run/physicsModified
/process/em/setSecBiasing conv World 0.0 1 MeV
/run/physicsModified
/run/physicsModified
/gun/particle e-
/gun/energy 1 GeV
/run/beamOn 100
### === Deexcitation model UAtomDeexcitation is activated for 3 regions:
DefaultRegionForTheWorld 0 0 0
VertexDetector 1 0 0
MuonDetector 0 0 0
### === Ignore cuts flag: 0
Secondary biasing is activated for e+ and eIoni inside G4Regions:
DefaultRegionForTheWorld BiasingWeight= 0
Secondary biasing is activated for e+ and eBrem inside G4Regions:
DefaultRegionForTheWorld BiasingWeight= 10
Secondary biasing is activated for e- and eIoni inside G4Regions:
DefaultRegionForTheWorld BiasingWeight= 0
Secondary biasing is activated for e- and eBrem inside G4Regions:
DefaultRegionForTheWorld BiasingWeight= 10
Secondary biasing is activated for gamma and phot inside G4Regions:
DefaultRegionForTheWorld BiasingWeight= 0
Secondary biasing is activated for gamma and compt inside G4Regions:
DefaultRegionForTheWorld BiasingWeight= 0
Secondary biasing is activated for gamma and conv inside G4Regions:
DefaultRegionForTheWorld BiasingWeight= 0
========= Table of registered couples ============================
Index : 0 used in the geometry : Yes
Material : G4_AIR
Range cuts : gamma 300 um e- 300 um e+ 300 um proton 300 um
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 30 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 1 used in the geometry : Yes
Material : G4_CESIUM_IODIDE
Range cuts : gamma 300 um e- 300 um e+ 300 um proton 300 um
Energy thresholds : gamma 18.27 keV e- 291.6 keV e+ 281.6 keV proton 30 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 2 used in the geometry : Yes
Material : G4_Al
Range cuts : gamma 300 um e- 300 um e+ 300 um proton 300 um
Energy thresholds : gamma 3.926 keV e- 260 keV e+ 253.1 keV proton 30 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 3 used in the geometry : Yes
Material : G4_AIR
Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 1 keV
Region(s) which use this couple :
VertexDetector
Index : 4 used in the geometry : Yes
Material : G4_Al
Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um
Energy thresholds : gamma 990 eV e- 34.11 keV e+ 33.72 keV proton 1 keV
Region(s) which use this couple :
VertexDetector
Index : 5 used in the geometry : Yes
Material : G4_Si
Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um
Energy thresholds : gamma 990 eV e- 31.7 keV e+ 31.35 keV proton 1 keV
Region(s) which use this couple :
VertexDetector
Index : 6 used in the geometry : Yes
Material : G4_Fe
Range cuts : gamma 1 cm e- 1 cm e+ 1 cm proton 1 cm
Energy thresholds : gamma 60.64 keV e- 13.04 MeV e+ 12.15 MeV proton 1 MeV
Region(s) which use this couple :
MuonDetector
Index : 7 used in the geometry : Yes
Material : G4_Al
Range cuts : gamma 1 cm e- 1 cm e+ 1 cm proton 1 cm
Energy thresholds : gamma 20.42 keV e- 4.563 MeV e+ 4.275 MeV proton 1 MeV
Region(s) which use this couple :
MuonDetector
==================================================================
### Run 2 starts.
### Run 2 start
### PrimaryGeneratorAction::GeneratePrimaries ##########
### e- E(MeV)= 1000
########################################################
Run terminated.
Run Summary
Number of events processed : 100
User=0.220000s Real=0.217922s Sys=0.000000s
RunAction: End of run actions are started
HistoManager: End of run actions are started RunID# 2
=================================================================
Number of events 100
Average number of e- 173.2
Average number of gamma 223
Average number of e+ 48.24
Average number of steps 1838
Edep 1x1 = 0.8322 +- 0.004077 res= 4.898 % 100
Edep 3x3 = 0.9578 +- 0.003436 res= 3.587 % 100
Edep 5x5 = 0.9806 +- 0.003509 res= 3.579 % 100
=========== Ratios without trancating ===========================
E1/E9 = 0.869 +- 0.003095
E1/E25 = 0.8488 +- 0.003257
E9/E25 = 0.9768 +- 0.001139
Beam Energy 1 GeV
==================================================================
#
/process/em/deexcitation VertexDetector true false true
/run/physicsModified
/process/em/fluo true
/run/physicsModified
/process/em/pixe true
/run/physicsModified
#/process/em/auger false
/run/physicsModified
#
/run/beamOn 100
### === Deexcitation model UAtomDeexcitation is activated for 3 regions:
DefaultRegionForTheWorld 0 0 0
VertexDetector 1 0 1
MuonDetector 0 0 0
### === G4UAtomicDeexcitation::InitialiseForNewRun()
### === Ignore cuts flag: 0
### === PIXE model for hadrons: Empirical
### === PIXE model for e+-: Livermore
Secondary biasing is activated for e+ and eIoni inside G4Regions:
DefaultRegionForTheWorld BiasingWeight= 0
Secondary biasing is activated for e+ and eBrem inside G4Regions:
DefaultRegionForTheWorld BiasingWeight= 10
Secondary biasing is activated for e- and eIoni inside G4Regions:
DefaultRegionForTheWorld BiasingWeight= 0
Secondary biasing is activated for e- and eBrem inside G4Regions:
DefaultRegionForTheWorld BiasingWeight= 10
Secondary biasing is activated for gamma and phot inside G4Regions:
DefaultRegionForTheWorld BiasingWeight= 0
Secondary biasing is activated for gamma and compt inside G4Regions:
DefaultRegionForTheWorld BiasingWeight= 0
Secondary biasing is activated for gamma and conv inside G4Regions:
DefaultRegionForTheWorld BiasingWeight= 0
========= Table of registered couples ============================
Index : 0 used in the geometry : Yes
Material : G4_AIR
Range cuts : gamma 300 um e- 300 um e+ 300 um proton 300 um
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 30 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 1 used in the geometry : Yes
Material : G4_CESIUM_IODIDE
Range cuts : gamma 300 um e- 300 um e+ 300 um proton 300 um
Energy thresholds : gamma 18.27 keV e- 291.6 keV e+ 281.6 keV proton 30 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 2 used in the geometry : Yes
Material : G4_Al
Range cuts : gamma 300 um e- 300 um e+ 300 um proton 300 um
Energy thresholds : gamma 3.926 keV e- 260 keV e+ 253.1 keV proton 30 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 3 used in the geometry : Yes
Material : G4_AIR
Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 1 keV
Region(s) which use this couple :
VertexDetector
Index : 4 used in the geometry : Yes
Material : G4_Al
Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um
Energy thresholds : gamma 990 eV e- 34.11 keV e+ 33.72 keV proton 1 keV
Region(s) which use this couple :
VertexDetector
Index : 5 used in the geometry : Yes
Material : G4_Si
Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um
Energy thresholds : gamma 990 eV e- 31.7 keV e+ 31.35 keV proton 1 keV
Region(s) which use this couple :
VertexDetector
Index : 6 used in the geometry : Yes
Material : G4_Fe
Range cuts : gamma 1 cm e- 1 cm e+ 1 cm proton 1 cm
Energy thresholds : gamma 60.64 keV e- 13.04 MeV e+ 12.15 MeV proton 1 MeV
Region(s) which use this couple :
MuonDetector
Index : 7 used in the geometry : Yes
Material : G4_Al
Range cuts : gamma 1 cm e- 1 cm e+ 1 cm proton 1 cm
Energy thresholds : gamma 20.42 keV e- 4.563 MeV e+ 4.275 MeV proton 1 MeV
Region(s) which use this couple :
MuonDetector
==================================================================
### Run 3 starts.
### Run 3 start
### PrimaryGeneratorAction::GeneratePrimaries ##########
### e- E(MeV)= 1000
########################################################
Run terminated.
Run Summary
Number of events processed : 100
User=0.220000s Real=0.223399s Sys=0.000000s
RunAction: End of run actions are started
HistoManager: End of run actions are started RunID# 3
=================================================================
Number of events 100
Average number of e- 172.7
Average number of gamma 223
Average number of e+ 47.33
Average number of steps 1842
Edep 1x1 = 0.8241 +- 0.003772 res= 4.577 % 100
Edep 3x3 = 0.9537 +- 0.003438 res= 3.605 % 100
Edep 5x5 = 0.9756 +- 0.003368 res= 3.452 % 100
=========== Ratios without trancating ===========================
E1/E9 = 0.8642 +- 0.00281
E1/E25 = 0.8449 +- 0.003046
E9/E25 = 0.9775 +- 0.001322
Beam Energy 1 GeV
==================================================================
#
================== Deleting memory pools ===================
Number of memory pools allocated: 9 of which, static: 0
Dynamic pools deleted: 9 / Total memory freed: 0.09 MB
============================================================