Gabriele Cosmo
696 lines
34 KiB
Plaintext
696 lines
34 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/
|
|
**************************************************************
|
|
|
|
PhysicsList:
|
|
/run/verbose 2
|
|
/tracking/verbose 0
|
|
#
|
|
/testhadr/TargetMat G4_W
|
|
/testhadr/TargetRadius 1 cm
|
|
/testhadr/TargetLength 10 cm
|
|
/testhadr/NumberDivZ 100
|
|
/testhadr/PrintModulo 10
|
|
#
|
|
/run/setCut 1 mm
|
|
#/testhadr/Physics QBBC
|
|
#
|
|
/physics_lists/factory/addNeutrino
|
|
/physics_lists/factory/addChargeExchange
|
|
#
|
|
/run/initialize
|
|
userDetector->Construct() start.
|
|
### Target consist of 100 disks of G4_W with R(mm)= 10 Width(mm)= 1 Total Length(mm)= 100 ###
|
|
|
|
***** Table : Nb of materials = 4 *****
|
|
|
|
Material: G4_Al density: 2.699 g/cm3 RadL: 8.896 cm Nucl.Int.Length: 38.894 cm
|
|
Imean: 166.000 eV temperature: 293.15 K pressure: 1.00 atm
|
|
|
|
---> Element: Al (Al) Z = 13.0 N = 27 A = 26.982 g/mole
|
|
---> Isotope: Al27 Z = 13 N = 27 A = 26.98 g/mole abundance: 100.000 %
|
|
ElmMassFraction: 100.00 % ElmAbundance 100.00 %
|
|
|
|
|
|
Material: G4_Galactic density: 0.000 mg/cm3 RadL: 204310098.490 pc Nucl.Int.Length: 3240.779 pc
|
|
Imean: 21.800 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: Battery density: 2.165 g/cm3 RadL: 10.553 cm Nucl.Int.Length: 46.079 cm
|
|
Imean: 144.880 eV temperature: 293.15 K pressure: 1.00 atm
|
|
|
|
---> Element: C (C) Z = 6.0 N = 12 A = 12.011 g/mole
|
|
---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.930 %
|
|
---> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.070 %
|
|
ElmMassFraction: 19.52 % ElmAbundance 30.66 %
|
|
|
|
---> Element: Al (Al) Z = 13.0 N = 27 A = 26.982 g/mole
|
|
---> Isotope: Al27 Z = 13 N = 27 A = 26.98 g/mole abundance: 100.000 %
|
|
ElmMassFraction: 39.80 % ElmAbundance 27.83 %
|
|
|
|
---> Element: Ti (Ti) Z = 22.0 N = 48 A = 47.867 g/mole
|
|
---> Isotope: Ti46 Z = 22 N = 46 A = 45.95 g/mole abundance: 8.250 %
|
|
---> Isotope: Ti47 Z = 22 N = 47 A = 46.95 g/mole abundance: 7.440 %
|
|
---> Isotope: Ti48 Z = 22 N = 48 A = 47.95 g/mole abundance: 73.720 %
|
|
---> Isotope: Ti49 Z = 22 N = 49 A = 48.95 g/mole abundance: 5.410 %
|
|
---> Isotope: Ti50 Z = 22 N = 50 A = 49.94 g/mole abundance: 5.180 %
|
|
ElmMassFraction: 2.00 % ElmAbundance 0.79 %
|
|
|
|
---> Element: Cu (Cu) Z = 29.0 N = 64 A = 63.546 g/mole
|
|
---> Isotope: Cu63 Z = 29 N = 63 A = 62.93 g/mole abundance: 69.170 %
|
|
---> Isotope: Cu65 Z = 29 N = 65 A = 64.93 g/mole abundance: 30.830 %
|
|
ElmMassFraction: 8.40 % ElmAbundance 2.49 %
|
|
|
|
---> Element: Li (Li) Z = 3.0 N = 7 A = 6.940 g/mole
|
|
---> Isotope: Li6 Z = 3 N = 6 A = 6.02 g/mole abundance: 7.590 %
|
|
---> Isotope: Li7 Z = 3 N = 7 A = 7.02 g/mole abundance: 92.410 %
|
|
ElmMassFraction: 1.70 % ElmAbundance 4.62 %
|
|
|
|
---> Element: Co (Co) Z = 27.0 N = 59 A = 58.933 g/mole
|
|
---> Isotope: Co59 Z = 27 N = 59 A = 58.93 g/mole abundance: 100.000 %
|
|
ElmMassFraction: 14.46 % ElmAbundance 4.63 %
|
|
|
|
---> Element: O (O) Z = 8.0 N = 16 A = 15.999 g/mole
|
|
---> Isotope: O16 Z = 8 N = 16 A = 15.99 g/mole abundance: 99.757 %
|
|
---> Isotope: O17 Z = 8 N = 17 A = 17.00 g/mole abundance: 0.038 %
|
|
---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.205 %
|
|
ElmMassFraction: 13.42 % ElmAbundance 15.82 %
|
|
|
|
---> 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: 0.70 % ElmAbundance 13.16 %
|
|
|
|
|
|
Material: G4_W density: 19.300 g/cm3 RadL: 3.504 mm Nucl.Int.Length: 10.312 cm
|
|
Imean: 727.000 eV temperature: 293.15 K pressure: 1.00 atm
|
|
|
|
---> Element: W (W) Z = 74.0 N = 184 A = 183.842 g/mole
|
|
---> Isotope: W180 Z = 74 N = 180 A = 179.95 g/mole abundance: 0.120 %
|
|
---> Isotope: W182 Z = 74 N = 182 A = 181.95 g/mole abundance: 26.500 %
|
|
---> Isotope: W183 Z = 74 N = 183 A = 182.95 g/mole abundance: 14.310 %
|
|
---> Isotope: W184 Z = 74 N = 184 A = 183.95 g/mole abundance: 30.640 %
|
|
---> Isotope: W186 Z = 74 N = 186 A = 185.95 g/mole abundance: 28.430 %
|
|
ElmMassFraction: 100.00 % ElmAbundance 100.00 %
|
|
|
|
|
|
|
|
World is registered to the default region.
|
|
physicsList->Construct() start.
|
|
physicsList->CheckParticleList() start.
|
|
physicsList->setCut() start.
|
|
#
|
|
/gun/particle proton
|
|
/gun/energy 3. 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 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
|
|
|
|
Region <DefaultRegionForTheWorld> -- -- appears in <World> world volume
|
|
This region is in the mass world.
|
|
Root logical volume(s) : World
|
|
Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
|
|
Materials : G4_Galactic G4_W
|
|
Production cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
|
|
|
|
Region <DefaultRegionForParallelWorld> -- -- is not associated to any world.
|
|
Root logical volume(s) :
|
|
Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
|
|
Materials :
|
|
Production cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
|
|
|
|
========= Table of registered couples ============================
|
|
|
|
Index : 0 used in the geometry : Yes
|
|
Material : G4_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_W
|
|
Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
|
|
Energy thresholds : gamma 106.43 keV e- 2.31004 MeV e+ 2.14587 MeV proton 100 keV
|
|
Region(s) which use this couple :
|
|
DefaultRegionForTheWorld
|
|
|
|
==================================================================
|
|
|
|
Start closing geometry.
|
|
--------------------------------------------------------------------------------
|
|
G4VoxelisationHelper::ReportVoxelStats -- Voxel Statistics
|
|
|
|
Total memory consumed for geometry optimisation: 10 kByte
|
|
Total CPU time elapsed for geometry optimisation: 0 seconds
|
|
|
|
Voxelisation: top CPU users:
|
|
Percent Total CPU System CPU Memory Volume
|
|
------- ---------- ---------- -------- ----------
|
|
0.00 0.00 0.00 11k Check
|
|
|
|
Voxelisation: top memory users:
|
|
Percent Memory Heads Nodes Pointers Total CPU Volume
|
|
------- -------- ------ ------ -------- ---------- ----------
|
|
100.00 10k 1 197 200 0.00 Check
|
|
--------------------------------------------------------------------------------
|
|
|
|
### Run 0 starts.
|
|
### Run 0 start
|
|
Initialisation time: User=0.280000s Real=0.764216s Sys=0.000000s
|
|
EventAction: Event #0 started
|
|
EventAction: Event #10 started
|
|
EventAction: Event #20 started
|
|
EventAction: Event #30 started
|
|
EventAction: Event #40 started
|
|
EventAction: Event #50 started
|
|
EventAction: Event #60 started
|
|
EventAction: Event #70 started
|
|
EventAction: Event #80 started
|
|
EventAction: Event #90 started
|
|
Run terminated.
|
|
Run Summary
|
|
Number of events processed : 100
|
|
User=0.020000s Real=0.041764s Sys=0.000000s
|
|
RunAction::EndOfRunAction: User=0.020000s Real=0.041806s Sys=0.000000s
|
|
RunAction: End of run actions are started
|
|
HistoManager: End of run actions are started
|
|
========================================================
|
|
Beam particle proton
|
|
Beam Energy(MeV) 3000
|
|
Number of events 100
|
|
Average energy deposit (MeV) 219.7 RMS(MeV) 8.754
|
|
Average number of steps 174.3
|
|
Average number of gamma 4.02
|
|
Average number of e- 8.35
|
|
Average number of e+ 0.05
|
|
Average number of neutrons 0
|
|
Average number of protons 0
|
|
Average number of antiprotons 0
|
|
Average number of pi+ & pi- 0
|
|
Average number of pi0 0
|
|
Average number of kaons 0
|
|
Average number of muons 0
|
|
Average number of deuterons+tritons 0
|
|
Average number of He3+alpha 0
|
|
Average number of ions 0.02
|
|
Average number of forward neutrons 0
|
|
Average number of reflected neutrons 0
|
|
Average number of leaked neutrons 0
|
|
Average number of proton leak 0
|
|
Average number of pion leak 0
|
|
========================================================
|
|
|
|
#
|
|
G4 kernel has come to Quit state.
|
|
Deleting G4Run (id:0)
|
|
UserDetectorConstruction deleted 0xb6e830
|
|
UserPhysicsList deleted 0xb91a88
|
|
UserActionInitialization deleted 0
|
|
UserWorkerInitialization deleted 0
|
|
UserWorkerThreadInitialization deleted 0
|
|
UserRunAction deleted.
|
|
UserPrimaryGenerator deleted.
|
|
RunManager is deleting RunManagerKernel.
|
|
G4SDManager deleted.
|
|
EventManager deleted.
|
|
Units table cleared.
|
|
TransportationManager deleted.
|
|
Total navigation history collections cleaned: 16
|
|
G4RNGHelper object is deleted.
|
|
================== Deleting memory pools ===================
|
|
Pool ID '20G4NavigationLevelRep', size : 0.0231 MB
|
|
Pool ID '24G4ReferenceCountedHandleIvE', size : 0.000961 MB
|
|
Pool ID '17G4DynamicParticle', size : 0.00288 MB
|
|
Pool ID '7G4Event', size : 0.000961 MB
|
|
Pool ID '15G4PrimaryVertex', size : 0.000961 MB
|
|
Pool ID '17G4PrimaryParticle', size : 0.000961 MB
|
|
Pool ID '15G4HCofThisEvent', size : 0.000961 MB
|
|
Pool ID '7G4Track', size : 0.00481 MB
|
|
Pool ID '18G4TouchableHistory', size : 0.00192 MB
|
|
Pool ID '15G4CountedObjectIvE', size : 0.000961 MB
|
|
Number of memory pools allocated: 10 of which, static: 0
|
|
Dynamic pools deleted: 10 / Total memory freed: 0.038 MB
|
|
============================================================
|
|
G4Allocator objects are deleted.
|
|
UImanager deleted.
|
|
StateManager deleted.
|
|
RunManagerKernel is deleted. Good bye :)
|
|
RunManager is deleted.
|