Gabriele Cosmo
661 lines
32 KiB
Plaintext
661 lines
32 KiB
Plaintext
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/
|
|
**************************************************************
|
|
|
|
/control/cout/ignoreThreadsExcept 0
|
|
/run/verbose 1
|
|
#
|
|
/testem/phys/addPhysics emstandard_opt0
|
|
PhysicsList::AddPhysicsList: <emstandard_opt0>
|
|
#
|
|
/process/em/verbose 1
|
|
/process/msc/Skin 1
|
|
/run/physicsModified
|
|
/run/initialize
|
|
|
|
-------------------------------------------------------------
|
|
---> The calorimeter is 50 layers of:
|
|
G4_Pb: 2.3 mm
|
|
G4_lAr: 5.7 mm
|
|
-------------------------------------------------------------
|
|
|
|
Material: Galactic density: 0.000 mg/cm3 RadL: 204310098.490 pc Nucl.Int.Length: 3240.779 pc
|
|
Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm
|
|
|
|
---> Element: H (H) Z = 1.0 N = 1 A = 1.008 g/mole
|
|
---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.989 %
|
|
---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.011 %
|
|
ElmMassFraction: 100.00 % ElmAbundance 100.00 %
|
|
|
|
|
|
Material: G4_Pb density: 11.350 g/cm3 RadL: 5.613 mm Nucl.Int.Length: 18.248 cm
|
|
Imean: 823.000 eV temperature: 293.15 K pressure: 1.00 atm
|
|
|
|
---> Element: Pb (Pb) Z = 82.0 N = 207 A = 207.217 g/mole
|
|
---> Isotope: Pb204 Z = 82 N = 204 A = 203.97 g/mole abundance: 1.400 %
|
|
---> Isotope: Pb206 Z = 82 N = 206 A = 205.97 g/mole abundance: 24.100 %
|
|
---> Isotope: Pb207 Z = 82 N = 207 A = 206.98 g/mole abundance: 22.100 %
|
|
---> Isotope: Pb208 Z = 82 N = 208 A = 207.98 g/mole abundance: 52.400 %
|
|
ElmMassFraction: 100.00 % ElmAbundance 100.00 %
|
|
|
|
|
|
Material: G4_lAr density: 1.396 g/cm3 RadL: 14.003 cm Nucl.Int.Length: 85.706 cm
|
|
Imean: 188.000 eV temperature: 293.15 K pressure: 1.00 atm
|
|
|
|
---> Element: Ar (Ar) Z = 18.0 N = 40 A = 39.948 g/mole
|
|
---> Isotope: Ar36 Z = 18 N = 36 A = 35.97 g/mole abundance: 0.337 %
|
|
---> Isotope: Ar38 Z = 18 N = 38 A = 37.96 g/mole abundance: 0.063 %
|
|
---> Isotope: Ar40 Z = 18 N = 40 A = 39.96 g/mole abundance: 99.600 %
|
|
ElmMassFraction: 100.00 % ElmAbundance 100.00 %
|
|
|
|
|
|
-------------------------------------------------------------
|
|
#
|
|
/gun/particle e-
|
|
/gun/energy 1 GeV
|
|
#
|
|
###/analysis/setFileName testem3
|
|
###/analysis/h1/set 1 100 700. 1000. MeV #edep in absor1
|
|
###/analysis/h1/set 2 100 0. 300. MeV #edep in absor2
|
|
###/analysis/h1/set 11 52 0. 52. none #long. profile in absor1
|
|
###/analysis/h1/set 12 52 0. 52. none #long. profile in absor2
|
|
###/analysis/h1/set 21 102 0. 102. none #energy flow
|
|
###/analysis/h1/set 22 102 0. 102. none #lateral energy leakage
|
|
#
|
|
/run/printProgress 50
|
|
/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 7
|
|
Verbose level 1
|
|
Verbose level for worker thread 0
|
|
Bremsstrahlung energy threshold above which
|
|
primary e+- is added to the list of secondary 100 TeV
|
|
Bremsstrahlung energy threshold above which primary
|
|
muon/hadron is added to the list of secondary 100 TeV
|
|
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.2, 1 mm)
|
|
Step function for muons/hadrons (0.2, 0.1 mm)
|
|
Step function for light ions (0.2, 0.1 mm)
|
|
Step function for general ions (0.2, 0.1 mm)
|
|
Lowest e+e- kinetic energy 1 keV
|
|
Lowest muon/hadron kinetic energy 1 keV
|
|
Use ICRU90 data 0
|
|
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 0
|
|
Max kinetic energy for CSDA tables 1 GeV
|
|
Max kinetic energy for NIEL computation 0 eV
|
|
Linear loss limit 0.01
|
|
Read data from file for e+e- pair production by mu 0
|
|
=======================================================================
|
|
====== Multiple Scattering Parameters ========
|
|
=======================================================================
|
|
Type of msc step limit algorithm for e+- 1
|
|
Type of msc step limit algorithm for muons/hadrons 0
|
|
Msc lateral displacement for e+- enabled 1
|
|
Msc lateral displacement for muons and hadrons 0
|
|
Urban msc model lateral displacement alg96 1
|
|
Range factor for msc step limit for e+- 0.04
|
|
Range factor for msc step limit for muons/hadrons 0.2
|
|
Geometry factor for msc step limitation of e+- 2.5
|
|
Safety factor for msc step limit for e+- 0.6
|
|
Skin parameter for msc step limitation of e+- 1
|
|
Lambda limit for msc step limit for e+- 1 mm
|
|
Use Mott correction for e- scattering 0
|
|
Factor used for dynamic computation of angular
|
|
limit between single and multiple scattering 1
|
|
Fixed angular limit between single
|
|
and multiple scattering 3.1416 rad
|
|
Upper energy limit for e+- multiple scattering 100 MeV
|
|
Type of electron single scattering model 0
|
|
Type of nuclear form-factor 1
|
|
Screening factor 1
|
|
=======================================================================
|
|
|
|
phot: for gamma SubType=12 BuildTable=0
|
|
LambdaPrime table from 200 keV to 100 TeV in 61 bins
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
LivermorePhElectric : Emin= 0 eV Emax= 100 TeV SauterGavrila Fluo
|
|
|
|
compt: for gamma SubType=13 BuildTable=1
|
|
Lambda table from 100 eV to 1 MeV, 7 bins/decade, spline: 1
|
|
LambdaPrime table from 1 MeV to 100 TeV in 56 bins
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
Klein-Nishina : Emin= 0 eV Emax= 100 TeV
|
|
|
|
conv: for gamma SubType=14 BuildTable=1
|
|
Lambda table from 1.022 MeV to 100 TeV, 18 bins/decade, spline: 1
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
BetheHeitlerLPM : Emin= 0 eV Emax= 100 TeV ModifiedTsai
|
|
|
|
Rayl: for gamma SubType=11 BuildTable=1
|
|
Lambda table from 100 eV to 150 keV, 7 bins/decade, spline: 0
|
|
LambdaPrime table from 150 keV to 100 TeV in 62 bins
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
LivermoreRayleigh : Emin= 0 eV Emax= 100 TeV CullenGenerator
|
|
|
|
msc: for e- SubType= 10
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
UrbanMsc : Emin= 0 eV Emax= 100 MeV Nbins=42 100 eV - 100 MeV
|
|
StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
|
|
WentzelVIUni : Emin= 100 MeV Emax= 100 TeV Nbins=42 100 MeV - 100 TeV
|
|
StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
|
|
|
|
eIoni: for e- XStype:3 SubType=2
|
|
dE/dx and range tables from 100 eV to 100 TeV in 84 bins
|
|
Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
|
|
StepFunction=(0.2, 1 mm), integ: 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 84 bins
|
|
Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
|
|
LPM flag: 1 for E > 1 GeV, VertexHighEnergyTh(GeV)= 100000
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
eBremSB : Emin= 0 eV Emax= 1 GeV ModifiedTsai
|
|
eBremLPM : Emin= 1 GeV Emax= 100 TeV ModifiedTsai
|
|
|
|
CoulombScat: for e- XStype:1 SubType=1 BuildTable=1
|
|
Lambda table from 100 MeV to 100 TeV, 7 bins/decade, spline: 0
|
|
ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
eCoulombScattering : Emin= 100 MeV Emax= 100 TeV
|
|
|
|
msc: for e+ SubType= 10
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
UrbanMsc : Emin= 0 eV Emax= 100 MeV Nbins=42 100 eV - 100 MeV
|
|
StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
|
|
WentzelVIUni : Emin= 100 MeV Emax= 100 TeV Nbins=42 100 MeV - 100 TeV
|
|
StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
|
|
|
|
eIoni: for e+ XStype:3 SubType=2
|
|
dE/dx and range tables from 100 eV to 100 TeV in 84 bins
|
|
Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
|
|
StepFunction=(0.2, 1 mm), integ: 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 84 bins
|
|
Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
|
|
LPM flag: 1 for E > 1 GeV, VertexHighEnergyTh(GeV)= 100000
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
eBremSB : Emin= 0 eV Emax= 1 GeV ModifiedTsai
|
|
eBremLPM : Emin= 1 GeV Emax= 100 TeV ModifiedTsai
|
|
|
|
annihil: for e+ XStype:2 SubType=5 AtRestModel:Simple BuildTable=0
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
eplus2gg : Emin= 0 eV Emax= 100 TeV
|
|
|
|
CoulombScat: for e+ XStype:1 SubType=1 BuildTable=1
|
|
Lambda table from 100 MeV to 100 TeV, 7 bins/decade, spline: 0
|
|
ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
eCoulombScattering : Emin= 100 MeV Emax= 100 TeV
|
|
|
|
msc: for proton SubType= 10
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV
|
|
StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
|
|
|
|
hIoni: for proton XStype:3 SubType=2
|
|
dE/dx and range tables from 100 eV to 100 TeV in 84 bins
|
|
Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
|
|
StepFunction=(0.2, 0.1 mm), integ: 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 84 bins
|
|
Lambda tables from threshold to 100 TeV, 7 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 84 bins
|
|
Lambda tables from threshold to 100 TeV, 7 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, 7 bins/decade, spline: 0
|
|
ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
eCoulombScattering : Emin= 0 eV Emax= 100 TeV
|
|
|
|
msc: for GenericIon SubType= 10
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
UrbanMsc : Emin= 0 eV Emax= 100 TeV
|
|
StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
|
|
|
|
ionIoni: for GenericIon XStype:3 SubType=2
|
|
dE/dx and range tables from 100 eV to 100 TeV in 84 bins
|
|
Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
|
|
StepFunction=(0.2, 0.1 mm), integ: 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 84 bins
|
|
Lambda tables from threshold to 100 TeV, 7 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=84 100 eV - 100 TeV
|
|
StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
|
|
|
|
hIoni: for anti_proton XStype:3 SubType=2
|
|
dE/dx and range tables from 100 eV to 100 TeV in 84 bins
|
|
Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
|
|
StepFunction=(0.2, 0.1 mm), integ: 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 84 bins
|
|
Lambda tables from threshold to 100 TeV, 7 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 84 bins
|
|
Lambda tables from threshold to 100 TeV, 7 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, 7 bins/decade, spline: 0
|
|
ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
eCoulombScattering : Emin= 0 eV Emax= 100 TeV
|
|
|
|
msc: for kaon+ SubType= 10
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV
|
|
StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
|
|
|
|
hIoni: for kaon+ XStype:3 SubType=2
|
|
dE/dx and range tables from 100 eV to 100 TeV in 84 bins
|
|
Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
|
|
StepFunction=(0.2, 0.1 mm), integ: 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 84 bins
|
|
Lambda tables from threshold to 100 TeV, 7 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 84 bins
|
|
Lambda tables from threshold to 100 TeV, 7 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, 7 bins/decade, spline: 0
|
|
ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
eCoulombScattering : Emin= 0 eV Emax= 100 TeV
|
|
|
|
msc: for kaon- SubType= 10
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV
|
|
StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
|
|
|
|
hIoni: for kaon- XStype:3 SubType=2
|
|
dE/dx and range tables from 100 eV to 100 TeV in 84 bins
|
|
Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
|
|
StepFunction=(0.2, 0.1 mm), integ: 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 84 bins
|
|
Lambda tables from threshold to 100 TeV, 7 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 84 bins
|
|
Lambda tables from threshold to 100 TeV, 7 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=84 100 eV - 100 TeV
|
|
StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
|
|
|
|
muIoni: for mu+ XStype:3 SubType=2
|
|
dE/dx and range tables from 100 eV to 100 TeV in 84 bins
|
|
Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
|
|
StepFunction=(0.2, 0.1 mm), integ: 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 84 bins
|
|
Lambda tables from threshold to 100 TeV, 7 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 84 bins
|
|
Lambda tables from threshold to 100 TeV, 7 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, 7 bins/decade, spline: 0
|
|
ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
eCoulombScattering : Emin= 0 eV Emax= 100 TeV
|
|
|
|
msc: for mu- SubType= 10
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV
|
|
StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
|
|
|
|
muIoni: for mu- XStype:3 SubType=2
|
|
dE/dx and range tables from 100 eV to 100 TeV in 84 bins
|
|
Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
|
|
StepFunction=(0.2, 0.1 mm), integ: 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 84 bins
|
|
Lambda tables from threshold to 100 TeV, 7 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 84 bins
|
|
Lambda tables from threshold to 100 TeV, 7 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=84 100 eV - 100 TeV
|
|
StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
|
|
|
|
hIoni: for pi+ XStype:3 SubType=2
|
|
dE/dx and range tables from 100 eV to 100 TeV in 84 bins
|
|
Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
|
|
StepFunction=(0.2, 0.1 mm), integ: 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 84 bins
|
|
Lambda tables from threshold to 100 TeV, 7 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 84 bins
|
|
Lambda tables from threshold to 100 TeV, 7 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, 7 bins/decade, spline: 0
|
|
ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
eCoulombScattering : Emin= 0 eV Emax= 100 TeV
|
|
|
|
msc: for pi- SubType= 10
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=84 100 eV - 100 TeV
|
|
StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm
|
|
|
|
hIoni: for pi- XStype:3 SubType=2
|
|
dE/dx and range tables from 100 eV to 100 TeV in 84 bins
|
|
Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
|
|
StepFunction=(0.2, 0.1 mm), integ: 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 84 bins
|
|
Lambda tables from threshold to 100 TeV, 7 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 84 bins
|
|
Lambda tables from threshold to 100 TeV, 7 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 : Galactic
|
|
Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
|
|
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 100 keV
|
|
Region(s) which use this couple :
|
|
DefaultRegionForTheWorld
|
|
|
|
Index : 1 used in the geometry : Yes
|
|
Material : G4_Pb
|
|
Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
|
|
Energy thresholds : gamma 101.797 keV e- 1.37158 MeV e+ 1.2826 MeV proton 100 keV
|
|
Region(s) which use this couple :
|
|
DefaultRegionForTheWorld
|
|
|
|
Index : 2 used in the geometry : Yes
|
|
Material : G4_lAr
|
|
Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
|
|
Energy thresholds : gamma 6.20202 keV e- 347.457 keV e+ 336.377 keV proton 100 keV
|
|
Region(s) which use this couple :
|
|
DefaultRegionForTheWorld
|
|
|
|
==================================================================
|
|
|
|
### Run 0 starts.
|
|
--> Event 0 starts.
|
|
--> Event 50 starts.
|
|
Run terminated.
|
|
Run Summary
|
|
Number of events processed : 100
|
|
User=0.710000s Real=0.703590s Sys=0.000000s
|
|
|
|
------------------------------------------------------------
|
|
material Edep RMS sqrt(E0(GeV))*rmsE/Emean total tracklen
|
|
|
|
G4_Pb: 774.03 MeV : 14.58 MeV 1.883 +- 0.1883 % 54.6 cm +- 1.17 cm
|
|
G4_lAr: 215.74 MeV : 13.37 MeV 6.199 +- 0.6199 % 1.06 m +- 6.84 cm
|
|
|
|
Total energy deposited = 989.8 MeV +- 6.905 MeV
|
|
Energy leakage = 10.22 MeV +- 6.905 MeV
|
|
Total energy : Edep + Eleak = 1 GeV +- 61.04 eV
|
|
------------------------------------------------------------
|
|
Beam particle e- E = 1 GeV
|
|
Mean number of gamma 508.8
|
|
Mean number of e- 870
|
|
Mean number of e+ 53.89
|
|
Mean number of charged steps 3710.83
|
|
Mean number of neutral steps 3715.22
|
|
------------------------------------------------------------
|
|
|
|
Energy deposition from Energy flow balance :
|
|
material Total Edep
|
|
|
|
G4_Pb: 774.038 MeV
|
|
G4_lAr: 215.736 MeV
|
|
|
|
------------------------------------------------------------
|
|
|
|
#
|
|
/process/em/setSecBiasing eIoni World 0.3 2 MeV
|
|
/run/physicsModified
|
|
/process/em/setSecBiasing eBrem World 0.2 1 MeV
|
|
/run/physicsModified
|
|
/run/physicsModified
|
|
/run/beamOn 100
|
|
Secondary biasing is activated for e+ and eIoni inside G4Regions:
|
|
DefaultRegionForTheWorld BiasingWeight= 3.33333
|
|
Secondary biasing is activated for e+ and eBrem inside G4Regions:
|
|
DefaultRegionForTheWorld BiasingWeight= 5
|
|
Secondary biasing is activated for e- and eIoni inside G4Regions:
|
|
DefaultRegionForTheWorld BiasingWeight= 3.33333
|
|
Secondary biasing is activated for e- and eBrem inside G4Regions:
|
|
DefaultRegionForTheWorld BiasingWeight= 5
|
|
|
|
========= Table of registered couples ============================
|
|
|
|
Index : 0 used in the geometry : Yes
|
|
Material : Galactic
|
|
Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
|
|
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 100 keV
|
|
Region(s) which use this couple :
|
|
DefaultRegionForTheWorld
|
|
|
|
Index : 1 used in the geometry : Yes
|
|
Material : G4_Pb
|
|
Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
|
|
Energy thresholds : gamma 101.797 keV e- 1.37158 MeV e+ 1.2826 MeV proton 100 keV
|
|
Region(s) which use this couple :
|
|
DefaultRegionForTheWorld
|
|
|
|
Index : 2 used in the geometry : Yes
|
|
Material : G4_lAr
|
|
Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
|
|
Energy thresholds : gamma 6.20202 keV e- 347.457 keV e+ 336.377 keV proton 100 keV
|
|
Region(s) which use this couple :
|
|
DefaultRegionForTheWorld
|
|
|
|
==================================================================
|
|
|
|
### Run 1 starts.
|
|
--> Event 0 starts.
|
|
--> Event 50 starts.
|
|
Run terminated.
|
|
Run Summary
|
|
Number of events processed : 100
|
|
User=0.580000s Real=0.589049s Sys=0.000000s
|
|
|
|
------------------------------------------------------------
|
|
material Edep RMS sqrt(E0(GeV))*rmsE/Emean total tracklen
|
|
|
|
G4_Pb: 778.62 MeV : 23.49 MeV 3.017 +- 0.3017 % 52.3 cm +- 1.35 cm
|
|
G4_lAr: 214.09 MeV : 18.28 MeV 8.54 +- 0.854 % 1.01 m +- 8.3 cm
|
|
|
|
Total energy deposited = 992.7 MeV +- 16.58 MeV
|
|
Energy leakage = 8.722 MeV +- 5.448 MeV
|
|
Total energy : Edep + Eleak = 1.001 GeV +- 15.52 MeV
|
|
------------------------------------------------------------
|
|
Beam particle e- E = 1 GeV
|
|
Mean number of gamma 292.4
|
|
Mean number of e- 557.7
|
|
Mean number of e+ 53.54
|
|
Mean number of charged steps 3140.74
|
|
Mean number of neutral steps 2920.78
|
|
------------------------------------------------------------
|
|
|
|
Energy deposition from Energy flow balance :
|
|
material Total Edep
|
|
|
|
G4_Pb: 775.556 MeV
|
|
G4_lAr: 215.711 MeV
|
|
|
|
------------------------------------------------------------
|
|
|
|
#
|
|
================== Deleting memory pools ===================
|
|
Number of memory pools allocated: 9 of which, static: 0
|
|
Dynamic pools deleted: 9 / Total memory freed: 0.088 MB
|
|
============================================================
|