Warning, /geant4/examples/extended/electromagnetic/TestEm5/TestEm5.out is written in an unsupported language. File is not indexed.
0001 Environment variable "G4FORCE_RUN_MANAGER_TYPE" enabled with value == Serial. Forcing G4RunManager type...
0002
0003 ############################################
0004 !!! WARNING - FPE detection is activated !!!
0005 ############################################
0006
0007
0008 ################################
0009 !!! G4Backtrace is activated !!!
0010 ################################
0011
0012
0013 **************************************************************
0014 Geant4 version Name: geant4-11-04-patch-01 (13-March-2026)
0015 Copyright : Geant4 Collaboration
0016 References : NIM A 506 (2003), 250-303
0017 : IEEE-TNS 53 (2006), 270-278
0018 : NIM A 835 (2016), 186-225
0019 WWW : http://geant4.org/
0020 **************************************************************
0021
0022 /control/cout/ignoreThreadsExcept 0
0023 /run/numberOfThreads 2
0024 *** /run/numberOfThreads command is issued in sequential mode.
0025 Command is ignored.
0026 /run/verbose 1
0027 #
0028 /testem/det/setAbsMat Silicon
0029 /testem/det/setAbsThick 20 um
0030 /testem/det/setAbsYZ 50 um
0031 #
0032 /testem/phys/addPhysics emstandard_opt3
0033 PhysicsList::AddPhysicsList: <emstandard_opt3>
0034 /process/em/verbose 1
0035 /process/em/workerVerbose 0
0036 #
0037 /run/setCut 7 um
0038 #
0039 /run/initialize
0040
0041 Material: G4_Galactic density: 0.000 mg/cm3 RadL: 204310098.490 pc Nucl.Int.Length: 3240.779 pc
0042 Imean: 21.800 eV temperature: 2.73 K pressure: 0.00 atm
0043
0044 ---> Element: H (H) Z = 1.0 N = 1 A = 1.008 g/mole
0045 ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.989 %
0046 ---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.011 %
0047 ElmMassFraction: 100.00 % ElmAbundance 100.00 %
0048
0049
0050 Material: Silicon density: 2.330 g/cm3 RadL: 9.366 cm Nucl.Int.Length: 45.660 cm
0051 Imean: 173.000 eV temperature: 293.15 K pressure: 1.00 atm
0052
0053 ---> Element: Si (Si) Z = 14.0 N = 28 A = 28.085 g/mole
0054 ---> Isotope: Si28 Z = 14 N = 28 A = 27.98 g/mole abundance: 92.230 %
0055 ---> Isotope: Si29 Z = 14 N = 29 A = 28.98 g/mole abundance: 4.683 %
0056 ---> Isotope: Si30 Z = 14 N = 30 A = 29.97 g/mole abundance: 3.087 %
0057 ElmMassFraction: 100.00 % ElmAbundance 100.00 %
0058
0059
0060 The WORLD is made of 24 um of G4_Galactic. The transverse size (YZ) of the world is 60 um
0061 The ABSORBER is made of 20 um of Silicon. The transverse size (YZ) is 50 um
0062 X position of the middle of the absorber 0 fm
0063 #
0064 /testem/gun/setDefault
0065 /gun/particle pi+
0066 /gun/energy 5 GeV
0067 #
0068 /run/printProgress 10000
0069 /run/beamOn 100000
0070 =======================================================================
0071 ====== Electromagnetic Physics Parameters ========
0072 =======================================================================
0073 LPM effect enabled 1
0074 Enable creation and use of sampling tables 0
0075 Apply cuts on all EM processes 0
0076 Use combined TransportationWithMsc Disabled
0077 Use general process 1
0078 Enable linear polarisation for gamma 0
0079 Enable photoeffect sampling below K-shell 1
0080 Enable sampling of quantum entanglement 0
0081 X-section factor for integral approach 0.8
0082 Min kinetic energy for tables 10 eV
0083 Max kinetic energy for tables 100 TeV
0084 Number of bins per decade of a table 20
0085 Verbose level 1
0086 Verbose level for worker thread 0
0087 Bremsstrahlung energy threshold above which
0088 primary e+- is added to the list of secondary 100 TeV
0089 Bremsstrahlung energy threshold above which primary
0090 muon/hadron is added to the list of secondary 100 TeV
0091 Positron annihilation at rest model AllisonPositronium
0092 Enable 3 gamma annihilation on fly 0
0093 Lowest triplet kinetic energy 1 MeV
0094 Enable sampling of gamma linear polarisation 0
0095 5D gamma conversion model type 0
0096 5D gamma conversion model on isolated ion 0
0097 Use RiGe 5D e+e- pair production model by muons 0
0098 Livermore data directory epics_2017
0099 =======================================================================
0100 ====== Ionisation Parameters ========
0101 =======================================================================
0102 Step function for e+- (0.2, 0.1 mm)
0103 Step function for muons/hadrons (0.2, 0.05 mm)
0104 Step function for light ions (0.1, 0.02 mm)
0105 Step function for general ions (0.1, 0.001 mm)
0106 Lowest e+e- kinetic energy 100 eV
0107 Lowest muon/hadron kinetic energy 1 keV
0108 Use ICRU90 data 1
0109 Fluctuations of dE/dx are enabled 1
0110 Type of fluctuation model for leptons and hadrons Urban
0111 Use built-in Birks saturation 0
0112 Build CSDA range enabled 0
0113 Use cut as a final range enabled 0
0114 Enable angular generator interface 1
0115 Max kinetic energy for CSDA tables 1 GeV
0116 Max kinetic energy for NIEL computation 1 MeV
0117 Linear loss limit 0.01
0118 Read data from file for e+e- pair production by mu 0
0119 =======================================================================
0120 ====== Multiple Scattering Parameters ========
0121 =======================================================================
0122 Type of msc step limit algorithm for e+- 3
0123 Type of msc step limit algorithm for muons/hadrons 0
0124 Msc lateral displacement for e+- enabled 1
0125 Msc lateral displacement for muons and hadrons 1
0126 Urban msc model lateral displacement alg96 1
0127 Range factor for msc step limit for e+- 0.04
0128 Range factor for msc step limit for muons/hadrons 0.2
0129 Geometry factor for msc step limitation of e+- 2.5
0130 Safety factor for msc step limit for e+- 0.6
0131 Skin parameter for msc step limitation of e+- 1
0132 Lambda limit for msc step limit for e+- 1 mm
0133 Use Mott correction for e- scattering 1
0134 Factor used for dynamic computation of angular
0135 limit between single and multiple scattering 1
0136 Fixed angular limit between single
0137 and multiple scattering 3.1416 rad
0138 Upper energy limit for e+- multiple scattering 100 MeV
0139 Type of electron single scattering model 0
0140 Type of nuclear form-factor 1
0141 Screening factor 1
0142 =======================================================================
0143 ====== Atomic Deexcitation Parameters ========
0144 =======================================================================
0145 Fluorescence enabled 1
0146 Directory in G4LEDATA for fluorescence data files fluor
0147 Auger electron cascade enabled 0
0148 PIXE atomic de-excitation enabled 0
0149 De-excitation module ignores cuts 0
0150 Type of PIXE cross section for hadrons Empirical
0151 Type of PIXE cross section for e+- Livermore
0152 =======================================================================
0153
0154 ### === Deexcitation model UAtomDeexcitation is activated for 1 region:
0155 DefaultRegionForTheWorld 1 0 0
0156 ### === Ignore cuts flag: 0
0157
0158 phot: for gamma SubType=12 BuildTable=0
0159 LambdaPrime table from 200 keV to 100 TeV in 174 bins
0160 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0161 LivermorePhElectric : Emin= 0 eV Emax= 100 TeV SauterGavrila Fluo
0162
0163 compt: for gamma SubType=13 BuildTable=1
0164 Lambda table from 10 eV to 1 MeV, 20 bins/decade, spline: 1
0165 LambdaPrime table from 1 MeV to 100 TeV in 160 bins
0166 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0167 KleinNishina : Emin= 0 eV Emax= 100 TeV Fluo
0168
0169 conv: for gamma SubType=14 BuildTable=1
0170 Lambda table from 1.022 MeV to 100 TeV, 20 bins/decade, spline: 1
0171 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0172 BetheHeitlerLPM : Emin= 0 eV Emax= 100 TeV ModifiedTsai
0173
0174 Rayl: for gamma SubType=11 BuildTable=1
0175 Lambda table from 10 eV to 150 keV, 20 bins/decade, spline: 0
0176 LambdaPrime table from 150 keV to 100 TeV in 176 bins
0177 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0178 LivermoreRayleigh : Emin= 0 eV Emax= 100 TeV CullenGenerator
0179
0180 msc: for e- SubType= 10
0181 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0182 UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV
0183 StepLim=DistanceToBoundary Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0184
0185 eIoni: for e- XStype:3 SubType=2
0186 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0187 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0188 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0189 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0190 MollerBhabha : Emin= 0 eV Emax= 100 TeV deltaVI
0191
0192 eBrem: for e- XStype:4 SubType=3
0193 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0194 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0195 LPM flag: 1 for E > 1 GeV, VertexHighEnergyTh(GeV)= 100000
0196 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0197 eBremSB : Emin= 0 eV Emax= 1 GeV AngularGen2BS
0198 eBremLPM : Emin= 1 GeV Emax= 100 TeV AngularGen2BS
0199
0200 ePairProd: for e- XStype:1 SubType=4
0201 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0202 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 0
0203 Sampling table 25x1001 from 0.1 GeV to 100 TeV
0204 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0205 ePairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0206
0207 msc: for e+ SubType= 10
0208 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0209 UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV
0210 StepLim=DistanceToBoundary Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0211
0212 eIoni: for e+ XStype:3 SubType=2
0213 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0214 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0215 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0216 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0217 MollerBhabha : Emin= 0 eV Emax= 100 TeV deltaVI
0218
0219 eBrem: for e+ XStype:4 SubType=3
0220 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0221 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0222 LPM flag: 1 for E > 1 GeV, VertexHighEnergyTh(GeV)= 100000
0223 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0224 eBremSB : Emin= 0 eV Emax= 1 GeV AngularGen2BS
0225 eBremLPM : Emin= 1 GeV Emax= 100 TeV AngularGen2BS
0226
0227 ePairProd: for e+ XStype:1 SubType=4
0228 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0229 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 0
0230 Sampling table 25x1001 from 0.1 GeV to 100 TeV
0231 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0232 ePairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0233
0234 annihil: for e+ XStype:2 SubType=5 AtRestModel:Allison BuildTable=0
0235 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0236 eplus2gg : Emin= 0 eV Emax= 100 TeV
0237
0238 msc: for proton SubType= 10
0239 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0240 UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV
0241 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0242
0243 hIoni: for proton XStype:3 SubType=2
0244 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0245 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0246 StepFunction=(0.2, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
0247 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0248 Bragg : Emin= 0 eV Emax= 2 MeV deltaVI
0249 BetheBloch : Emin= 2 MeV Emax= 100 TeV deltaVI
0250
0251 hBrems: for proton XStype:1 SubType=3
0252 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0253 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0254 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0255 hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0256
0257 hPairProd: for proton XStype:1 SubType=4
0258 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0259 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0260 Sampling table 17x1001, from 7.50618 GeV to 100 TeV
0261 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0262 hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0263
0264 nuclearStopping: for proton SubType=8 BuildTable=0
0265 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0266 ICRU49NucStopping : Emin= 0 eV Emax= 1 MeV
0267
0268 msc: for GenericIon SubType= 10
0269 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0270 UrbanMsc : Emin= 0 eV Emax= 100 TeV
0271 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0272
0273 ionIoni: for GenericIon XStype:3 SubType=2
0274 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0275 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0276 StepFunction=(0.1, 0.001 mm), integ: 3, fluct: 1, linLossLim= 0.02
0277 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0278 LindhardSorensen : Emin= 0 eV Emax= 100 TeV deltaVI
0279
0280 nuclearStopping: for GenericIon SubType=8 BuildTable=0
0281 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0282 ICRU49NucStopping : Emin= 0 eV Emax= 1 MeV
0283
0284 msc: for alpha SubType= 10
0285 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0286 UrbanMsc : Emin= 0 eV Emax= 100 TeV
0287 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0288
0289 ionIoni: for alpha XStype:3 SubType=2
0290 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0291 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0292 StepFunction=(0.1, 0.02 mm), integ: 3, fluct: 1, linLossLim= 0.02
0293 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0294 BraggIon : Emin= 0 eV Emax=7.9452 MeV deltaVI
0295 BetheBloch : Emin=7.9452 MeV Emax= 100 TeV deltaVI
0296
0297 nuclearStopping: for alpha SubType=8 BuildTable=0
0298 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0299 ICRU49NucStopping : Emin= 0 eV Emax= 1 MeV
0300
0301 msc: for anti_proton SubType= 10
0302 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0303 UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV
0304 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0305
0306 hIoni: for anti_proton XStype:3 SubType=2
0307 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0308 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0309 StepFunction=(0.2, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
0310 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0311 ICRU73QO : Emin= 0 eV Emax= 2 MeV deltaVI
0312 BetheBloch : Emin= 2 MeV Emax= 100 TeV deltaVI
0313
0314 hBrems: for anti_proton XStype:1 SubType=3
0315 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0316 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0317 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0318 hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0319
0320 hPairProd: for anti_proton XStype:1 SubType=4
0321 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0322 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0323 Sampling table 17x1001, from 7.50618 GeV to 100 TeV
0324 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0325 hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0326
0327 msc: for kaon+ SubType= 10
0328 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0329 UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV
0330 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0331
0332 hIoni: for kaon+ XStype:3 SubType=2
0333 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0334 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0335 StepFunction=(0.2, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
0336 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0337 Bragg : Emin= 0 eV Emax=1.05231 MeV deltaVI
0338 BetheBloch : Emin=1.05231 MeV Emax= 100 TeV deltaVI
0339
0340 hBrems: for kaon+ XStype:1 SubType=3
0341 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0342 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0343 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0344 hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0345
0346 hPairProd: for kaon+ XStype:1 SubType=4
0347 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0348 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0349 Sampling table 18x1001, from 3.94942 GeV to 100 TeV
0350 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0351 hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0352
0353 msc: for kaon- SubType= 10
0354 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0355 UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV
0356 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0357
0358 hIoni: for kaon- XStype:3 SubType=2
0359 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0360 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0361 StepFunction=(0.2, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
0362 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0363 ICRU73QO : Emin= 0 eV Emax=1.05231 MeV deltaVI
0364 BetheBloch : Emin=1.05231 MeV Emax= 100 TeV deltaVI
0365
0366 hBrems: for kaon- XStype:1 SubType=3
0367 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0368 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0369 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0370 hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0371
0372 hPairProd: for kaon- XStype:1 SubType=4
0373 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0374 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0375 Sampling table 18x1001, from 3.94942 GeV to 100 TeV
0376 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0377 hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0378
0379 msc: for mu+ SubType= 10
0380 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0381 UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV
0382 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0383
0384 muIoni: for mu+ XStype:3 SubType=2
0385 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0386 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0387 StepFunction=(0.2, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
0388 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0389 Bragg : Emin= 0 eV Emax= 200 keV deltaVI
0390 MuBetheBloch : Emin= 200 keV Emax= 100 TeV deltaVI
0391
0392 muBrems: for mu+ XStype:1 SubType=3
0393 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0394 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0395 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0396 MuBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0397
0398 muPairProd: for mu+ XStype:1 SubType=4
0399 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0400 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0401 Sampling table 21x1001, from 0.85 GeV to 100 TeV
0402 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0403 muPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0404
0405 msc: for mu- SubType= 10
0406 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0407 UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV
0408 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0409
0410 muIoni: for mu- XStype:3 SubType=2
0411 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0412 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0413 StepFunction=(0.2, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
0414 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0415 ICRU73QO : Emin= 0 eV Emax= 200 keV deltaVI
0416 MuBetheBloch : Emin= 200 keV Emax= 100 TeV deltaVI
0417
0418 muBrems: for mu- XStype:1 SubType=3
0419 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0420 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0421 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0422 MuBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0423
0424 muPairProd: for mu- XStype:1 SubType=4
0425 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0426 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0427 Sampling table 21x1001, from 0.85 GeV to 100 TeV
0428 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0429 muPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0430
0431 msc: for pi+ SubType= 10
0432 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0433 UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV
0434 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0435
0436 hIoni: for pi+ XStype:3 SubType=2
0437 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0438 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0439 StepFunction=(0.2, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
0440 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0441 Bragg : Emin= 0 eV Emax=297.505 keV deltaVI
0442 BetheBloch : Emin=297.505 keV Emax= 100 TeV deltaVI
0443
0444 hBrems: for pi+ XStype:1 SubType=3
0445 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0446 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0447 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0448 hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0449
0450 hPairProd: for pi+ XStype:1 SubType=4
0451 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0452 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0453 Sampling table 20x1001, from 1.11656 GeV to 100 TeV
0454 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0455 hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0456
0457 msc: for pi- SubType= 10
0458 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0459 UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV
0460 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0461
0462 hIoni: for pi- XStype:3 SubType=2
0463 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0464 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0465 StepFunction=(0.2, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
0466 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0467 ICRU73QO : Emin= 0 eV Emax=297.505 keV deltaVI
0468 BetheBloch : Emin=297.505 keV Emax= 100 TeV deltaVI
0469
0470 hBrems: for pi- XStype:1 SubType=3
0471 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0472 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0473 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0474 hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0475
0476 hPairProd: for pi- XStype:1 SubType=4
0477 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0478 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0479 Sampling table 20x1001, from 1.11656 GeV to 100 TeV
0480 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0481 hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0482
0483 ========= Table of registered couples ============================
0484
0485 Index : 0 used in the geometry : Yes
0486 Material : G4_Galactic
0487 Range cuts : gamma 7 um e- 7 um e+ 7 um proton 7 um
0488 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 700 eV
0489 Region(s) which use this couple :
0490 DefaultRegionForTheWorld
0491
0492 Index : 1 used in the geometry : Yes
0493 Material : Silicon
0494 Range cuts : gamma 7 um e- 7 um e+ 7 um proton 7 um
0495 Energy thresholds : gamma 990 eV e- 21.0638 keV e+ 20.5957 keV proton 700 eV
0496 Region(s) which use this couple :
0497 DefaultRegionForTheWorld
0498
0499 ==================================================================
0500
0501 ### Run 0 starts.
0502
0503 ------- MixMaxRng engine status -------
0504 Current state vector is:
0505 mixmax state, file version 1.0
0506 N=17 V[N]={906770732717044781, 629165745432651234, 1235682547346241386, 68420008233404568, 2262190991329242458, 2266470399991071809, 1976726662926872232, 245458862506414172, 1955974201201518530, 2155248512522080758, 604170912935414061, 1116171330120743511, 1861018313684488333, 1296715403254578286, 1549011045957234151, 370819759640195970, 2230139271784837643} counter= 17sumtot= 1977567618660788324
0507 ---------------------------------------
0508 --> Event 0 starts.
0509 --> Event 10000 starts.
0510 --> Event 20000 starts.
0511 --> Event 30000 starts.
0512 --> Event 40000 starts.
0513 --> Event 50000 starts.
0514 --> Event 60000 starts.
0515 --> Event 70000 starts.
0516 --> Event 80000 starts.
0517 --> Event 90000 starts.
0518 Run terminated.
0519 Run Summary
0520 Number of events processed : 100000
0521 User=0.430000s Real=0.434872s Sys=0.000000s
0522
0523 ======================== run summary ======================
0524
0525 The run was 100000 pi+ of 5 GeV through 20 um of Silicon (density: 2.33 g/cm3 )
0526
0527 Total energy deposit in absorber per event = 6.353 keV +- 16.34 eV
0528
0529 -----> Mean dE/dx = 3.177 MeV/cm (1.363 MeV*cm2/g)
0530
0531 From formulas :
0532 restricted dEdx = 2.96 MeV/cm (1.27 MeV*cm2/g)
0533 full dEdx = 4.726 MeV/cm (2.028 MeV*cm2/g)
0534
0535 Leakage : primary = 5 GeV +- 725.3 eV secondaries = 2.672 keV +- 724.6 eV
0536 Energy balance : edep + eleak = 5 GeV
0537
0538 Total track length (charged) in absorber per event = 20.26 um +- 7.771 nm
0539 Total track length (neutral) in absorber per event = 7.453e+04 fm +- 7.453e+04 fm
0540
0541 Number of steps (charged) in absorber per event = 1.179 +- 0.005528
0542 Number of steps (neutral) in absorber per event = 1e-05 +- 1e-05
0543
0544 Number of secondaries per event : Gammas = 1e-05 electrons = 0.017 positrons = 0
0545
0546 Number of events with the primary particle transmitted = 100 %
0547 Number of events with at least 1 particle transmitted (same charge as primary) = 100 %
0548
0549 Number of events with the primary particle reflected = 0 %
0550 Number of events with at least 1 particle reflected (same charge as primary) = 0 %
0551
0552 MultipleScattering:
0553 rms proj angle of transmit primary particle = 0.02504 mrad (central part only)
0554 computed theta0 (Highland formula) = 0.02627 mrad
0555 central part defined as +- 0.07881 mrad Tail ratio = 1.951 %
0556
0557 Gamma process counts:
0558 Photoeffect 0
0559 Compton 0
0560 Conversion 0
0561 Rayleigh 0
0562
0563 ------- MixMaxRng engine status -------
0564 Current state vector is:
0565 mixmax state, file version 1.0
0566 N=17 V[N]={1378456096184970916, 645125664873386088, 289313246081823319, 692181631145386720, 1554795592127445348, 429523111405020466, 1767267687826005316, 887948550155982521, 745800083929026403, 309730682002155907, 2175523263042675480, 1439287206468397773, 1375105445691351298, 782374144723864041, 1715829858399311388, 655882193793759711, 1581270085403564578} counter= 8sumtot= 2284513478758269616
0567 ---------------------------------------
0568 #
0569 /testem/det/setAbsMat Gold
0570 /testem/det/setAbsThick 9.658 um
0571 /testem/det/setAbsYZ 1 mm
0572 /testem/gun/setDefault
0573 # Field Should be issued after /run/initialize
0574 #/testem/det/setField 1 tesla
0575 /gun/particle e-
0576 /gun/energy 15.7 MeV
0577 /run/beamOn 50000
0578
0579 ### === Deexcitation model UAtomDeexcitation is activated for 1 region:
0580 DefaultRegionForTheWorld 1 0 0
0581 ### === Ignore cuts flag: 0
0582
0583 ========= Table of registered couples ============================
0584
0585 Index : 0 used in the geometry : Yes
0586 Material : G4_Galactic
0587 Range cuts : gamma 7 um e- 7 um e+ 7 um proton 7 um
0588 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 700 eV
0589 Region(s) which use this couple :
0590 DefaultRegionForTheWorld
0591
0592 Index : 1 used in the geometry : No
0593 Material : Silicon
0594 Range cuts : gamma 7 um e- 7 um e+ 7 um proton 7 um
0595 Energy thresholds : gamma 990 eV e- 21.0638 keV e+ 20.5957 keV proton 700 eV
0596
0597 Index : 2 used in the geometry : Yes
0598 Material : Gold
0599 Range cuts : gamma 7 um e- 7 um e+ 7 um proton 7 um
0600 Energy thresholds : gamma 6.52544 keV e- 65.1638 keV e+ 63.9139 keV proton 700 eV
0601 Region(s) which use this couple :
0602 DefaultRegionForTheWorld
0603
0604 ==================================================================
0605
0606 ### Run 1 starts.
0607
0608 ------- MixMaxRng engine status -------
0609 Current state vector is:
0610 mixmax state, file version 1.0
0611 N=17 V[N]={1378456096184970916, 645125664873386088, 289313246081823319, 692181631145386720, 1554795592127445348, 429523111405020466, 1767267687826005316, 887948550155982521, 745800083929026403, 309730682002155907, 2175523263042675480, 1439287206468397773, 1375105445691351298, 782374144723864041, 1715829858399311388, 655882193793759711, 1581270085403564578} counter= 8sumtot= 2284513478758269616
0612 ---------------------------------------
0613 --> Event 0 starts.
0614 --> Event 10000 starts.
0615 --> Event 20000 starts.
0616 --> Event 30000 starts.
0617 --> Event 40000 starts.
0618 Run terminated.
0619 Run Summary
0620 Number of events processed : 50000
0621 User=0.300000s Real=0.298352s Sys=0.000000s
0622
0623 ======================== run summary ======================
0624
0625 The run was 50000 e- of 15.7 MeV through 9.66 um of Gold (density: 19.3 g/cm3 )
0626
0627 Total energy deposit in absorber per event = 18.38 keV +- 64.44 eV
0628
0629 -----> Mean dE/dx = 19.03 MeV/cm (0.9848 MeV*cm2/g)
0630
0631 From formulas :
0632 restricted dEdx = 17.91 MeV/cm (0.927 MeV*cm2/g)
0633 full dEdx = 61.41 MeV/cm (3.179 MeV*cm2/g)
0634
0635 Leakage : primary = 15.64 MeV +- 2.11 keV secondaries = 37.48 keV +- 2.102 keV
0636 Energy balance : edep + eleak = 15.7 MeV
0637
0638 Total track length (charged) in absorber per event = 9.911 um +- 11.24 nm
0639 Total track length (neutral) in absorber per event = 141.4 nm +- 13.24 nm
0640
0641 Number of steps (charged) in absorber per event = 2.602 +- 0.009024
0642 Number of steps (neutral) in absorber per event = 0.02698 +- 0.0007528
0643
0644 Number of secondaries per event : Gammas = 0.0269 electrons = 0.02112 positrons = 0
0645
0646 Number of events with the primary particle transmitted = 99.98 %
0647 Number of events with at least 1 particle transmitted (same charge as primary) = 99.98 %
0648
0649 Number of events with the primary particle reflected = 0.02 %
0650 Number of events with at least 1 particle reflected (same charge as primary) = 0.31 %
0651
0652 MultipleScattering:
0653 rms proj angle of transmit primary particle = 34.71 mrad (central part only)
0654 computed theta0 (Highland formula) = 35.1 mrad
0655 central part defined as +- 105.3 mrad Tail ratio = 2.948 %
0656
0657 Gamma process counts:
0658 Photoeffect 198
0659 Compton 0
0660 Conversion 0
0661 Rayleigh 4
0662
0663 ------- MixMaxRng engine status -------
0664 Current state vector is:
0665 mixmax state, file version 1.0
0666 N=17 V[N]={187465526861526162, 1919470187343084000, 1200600146555496649, 2254698630955541145, 1010695742965213971, 386425232931658995, 1773838893809632618, 1545813557997812675, 1661416122218187049, 1838452107617613130, 852398828744869876, 2030305144316100, 2300250181027420193, 126900021668951897, 823556986218677264, 1235812847093715199, 1130805818374825746} counter= 11sumtot= 1803887063818991061
0667 ---------------------------------------
0668 #
0669 /run/setCut 10 um
0670 #
0671 /testem/det/setAbsMat G4_POLYETHYLENE
0672 /testem/det/setAbsThick 4.74 mm
0673 /testem/det/setAbsYZ 50 mm
0674 /testem/gun/setDefault
0675 /gun/particle mu+
0676 /gun/energy 96.2 MeV
0677 /run/beamOn 50000
0678
0679 ### === Deexcitation model UAtomDeexcitation is activated for 1 region:
0680 DefaultRegionForTheWorld 1 0 0
0681 ### === Ignore cuts flag: 0
0682
0683 ========= Table of registered couples ============================
0684
0685 Index : 0 used in the geometry : Yes
0686 Material : G4_Galactic
0687 Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um
0688 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 1 keV
0689 Region(s) which use this couple :
0690 DefaultRegionForTheWorld
0691
0692 Index : 1 used in the geometry : No
0693 Material : Silicon
0694 Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um
0695 Energy thresholds : gamma 990 eV e- 31.6994 keV e+ 31.3501 keV proton 1 keV
0696
0697 Index : 2 used in the geometry : No
0698 Material : Gold
0699 Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um
0700 Energy thresholds : gamma 8.32988 keV e- 80.7452 keV e+ 79.0938 keV proton 1 keV
0701
0702 Index : 3 used in the geometry : Yes
0703 Material : G4_POLYETHYLENE
0704 Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um
0705 Energy thresholds : gamma 990 eV e- 13.7224 keV e+ 13.4607 keV proton 1 keV
0706 Region(s) which use this couple :
0707 DefaultRegionForTheWorld
0708
0709 ==================================================================
0710
0711 ### Run 2 starts.
0712
0713 ------- MixMaxRng engine status -------
0714 Current state vector is:
0715 mixmax state, file version 1.0
0716 N=17 V[N]={187465526861526162, 1919470187343084000, 1200600146555496649, 2254698630955541145, 1010695742965213971, 386425232931658995, 1773838893809632618, 1545813557997812675, 1661416122218187049, 1838452107617613130, 852398828744869876, 2030305144316100, 2300250181027420193, 126900021668951897, 823556986218677264, 1235812847093715199, 1130805818374825746} counter= 11sumtot= 1803887063818991061
0717 ---------------------------------------
0718 --> Event 0 starts.
0719 --> Event 10000 starts.
0720 --> Event 20000 starts.
0721 --> Event 30000 starts.
0722 --> Event 40000 starts.
0723 Run terminated.
0724 Run Summary
0725 Number of events processed : 50000
0726 User=0.910000s Real=0.936989s Sys=0.000000s
0727
0728 ======================== run summary ======================
0729
0730 The run was 50000 mu+ of 96.2 MeV through 4.74 mm of G4_POLYETHYLENE (density: 940 mg/cm3)
0731
0732 Total energy deposit in absorber per event = 1.048 MeV +- 993.7 eV
0733
0734 -----> Mean dE/dx = 2.21 MeV/cm (2.351 MeV*cm2/g)
0735
0736 From formulas :
0737 restricted dEdx = 1.749 MeV/cm (1.861 MeV*cm2/g)
0738 full dEdx = 2.265 MeV/cm (2.409 MeV*cm2/g)
0739
0740 Leakage : primary = 95.13 MeV +- 1.332 keV secondaries = 26.01 keV +- 701.3 eV
0741 Energy balance : edep + eleak = 96.2 MeV
0742
0743 Total track length (charged) in absorber per event = 5.17 mm +- 3.775 um
0744 Total track length (neutral) in absorber per event = 22.35 um +- 2.273 um
0745
0746 Number of steps (charged) in absorber per event = 12.17 +- 0.03721
0747 Number of steps (neutral) in absorber per event = 0.00988 +- 0.0004846
0748
0749 Number of secondaries per event : Gammas = 0.00914 electrons = 4.181 positrons = 0
0750
0751 Number of events with the primary particle transmitted = 100 %
0752 Number of events with at least 1 particle transmitted (same charge as primary) = 100 %
0753
0754 Number of events with the primary particle reflected = 0 %
0755 Number of events with at least 1 particle reflected (same charge as primary) = 0 %
0756
0757 MultipleScattering:
0758 rms proj angle of transmit primary particle = 7.373 mrad (central part only)
0759 computed theta0 (Highland formula) = 7.636 mrad
0760 central part defined as +- 22.91 mrad Tail ratio = 1.214 %
0761
0762 Gamma process counts:
0763 Photoeffect 305
0764 Compton 23
0765 Conversion 0
0766 Rayleigh 14
0767
0768 ------- MixMaxRng engine status -------
0769 Current state vector is:
0770 mixmax state, file version 1.0
0771 N=17 V[N]={2083361880407413973, 558014575596605443, 155743980189950201, 2182957000861726697, 1838655335096136537, 1946843552728816280, 1681406852924851701, 474130547566224162, 1018240738768540513, 1872432096900372193, 472993835442224762, 10222201860268061, 1779133025925577121, 2263978871037324573, 1473228915816656057, 35068858945786700, 1174357358231345229} counter= 2sumtot= 268182545376574644
0772 ---------------------------------------
0773 #
0774 /run/setCut 1 mm
0775 #
0776 /testem/stack/killSecondaries
0777 #
0778 /testem/det/setAbsMat Aluminium
0779 /testem/det/setAbsThick 8.004 mm
0780 /testem/gun/setDefault
0781 /gun/particle proton
0782 /gun/energy 174.664 GeV
0783 /run/beamOn 50000
0784
0785 ### === Deexcitation model UAtomDeexcitation is activated for 1 region:
0786 DefaultRegionForTheWorld 1 0 0
0787 ### === Ignore cuts flag: 0
0788
0789 ========= Table of registered couples ============================
0790
0791 Index : 0 used in the geometry : Yes
0792 Material : G4_Galactic
0793 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
0794 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 100 keV
0795 Region(s) which use this couple :
0796 DefaultRegionForTheWorld
0797
0798 Index : 1 used in the geometry : No
0799 Material : Silicon
0800 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
0801 Energy thresholds : gamma 6.93672 keV e- 546.311 keV e+ 525.593 keV proton 100 keV
0802
0803 Index : 2 used in the geometry : No
0804 Material : Gold
0805 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
0806 Energy thresholds : gamma 112.405 keV e- 2.28837 MeV e+ 2.12466 MeV proton 100 keV
0807
0808 Index : 3 used in the geometry : No
0809 Material : G4_POLYETHYLENE
0810 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
0811 Energy thresholds : gamma 2.25978 keV e- 351.127 keV e+ 341.776 keV proton 100 keV
0812
0813 Index : 4 used in the geometry : Yes
0814 Material : Aluminium
0815 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
0816 Energy thresholds : gamma 6.91691 keV e- 597.603 keV e+ 574.29 keV proton 100 keV
0817 Region(s) which use this couple :
0818 DefaultRegionForTheWorld
0819
0820 ==================================================================
0821
0822 ### Run 3 starts.
0823
0824 ------- MixMaxRng engine status -------
0825 Current state vector is:
0826 mixmax state, file version 1.0
0827 N=17 V[N]={2083361880407413973, 558014575596605443, 155743980189950201, 2182957000861726697, 1838655335096136537, 1946843552728816280, 1681406852924851701, 474130547566224162, 1018240738768540513, 1872432096900372193, 472993835442224762, 10222201860268061, 1779133025925577121, 2263978871037324573, 1473228915816656057, 35068858945786700, 1174357358231345229} counter= 2sumtot= 268182545376574644
0828 ---------------------------------------
0829 --> Event 0 starts.
0830 --> Event 10000 starts.
0831 --> Event 20000 starts.
0832 --> Event 30000 starts.
0833 --> Event 40000 starts.
0834 Run terminated.
0835 Run Summary
0836 Number of events processed : 50000
0837 User=0.220000s Real=0.217912s Sys=0.000000s
0838
0839 ======================== run summary ======================
0840
0841 The run was 50000 proton of 175 GeV through 8 mm of Aluminium (density: 2.7 g/cm3 )
0842
0843 Total energy deposit in absorber per event = 5.362 MeV +- 381.9 keV
0844
0845 -----> Mean dE/dx = 6.699 MeV/cm (2.481 MeV*cm2/g)
0846
0847 From formulas :
0848 restricted dEdx = 4.024 MeV/cm (1.49 MeV*cm2/g)
0849 full dEdx = 6.049 MeV/cm (2.24 MeV*cm2/g)
0850
0851 Leakage : primary = 174.7 GeV +- 381.9 keV secondaries = 0 eV +- 0 eV
0852 Energy balance : edep + eleak = 174.7 GeV
0853
0854 Total track length (charged) in absorber per event = 8.004 mm +- 0 fm
0855 Total track length (neutral) in absorber per event = 0 fm +- 0 fm
0856
0857 Number of steps (charged) in absorber per event = 1.268 +- 0.002316
0858 Number of steps (neutral) in absorber per event = 0 +- 0
0859
0860 Number of secondaries per event : Gammas = 0 electrons = 0.268 positrons = 0.00088
0861
0862 Number of events with the primary particle transmitted = 100 %
0863 Number of events with at least 1 particle transmitted (same charge as primary) = 100 %
0864
0865 Number of events with the primary particle reflected = 0 %
0866 Number of events with at least 1 particle reflected (same charge as primary) = 0 %
0867
0868 MultipleScattering:
0869 rms proj angle of transmit primary particle = 0.02062 mrad (central part only)
0870 computed theta0 (Highland formula) = 0.02111 mrad
0871 central part defined as +- 0.06333 mrad Tail ratio = 1.327 %
0872
0873 Gamma process counts:
0874 Photoeffect 0
0875 Compton 0
0876 Conversion 0
0877 Rayleigh 0
0878
0879 ------- MixMaxRng engine status -------
0880 Current state vector is:
0881 mixmax state, file version 1.0
0882 N=17 V[N]={864120598200464168, 189420282748189000, 2125101903705486721, 1863466081725669111, 1383960382039097673, 550892631529147885, 2244069534527202395, 933979936320874605, 791647544461653037, 286412522794523854, 1616800613386878212, 478416970401061220, 1548587089976077386, 917812281992815632, 2284924600166281660, 94168853074452663, 1674934093723585259} counter= 15sumtot= 1401971847063908873
0883 ---------------------------------------
0884 #
0885 /testem/det/setAbsMat Iron
0886 /testem/det/setAbsThick 1 m
0887 /testem/gun/setDefault
0888 /gun/particle mu+
0889 /gun/energy 100 GeV
0890 /run/beamOn 50000
0891
0892 ### === Deexcitation model UAtomDeexcitation is activated for 1 region:
0893 DefaultRegionForTheWorld 1 0 0
0894 ### === Ignore cuts flag: 0
0895
0896 ========= Table of registered couples ============================
0897
0898 Index : 0 used in the geometry : Yes
0899 Material : G4_Galactic
0900 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
0901 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 100 keV
0902 Region(s) which use this couple :
0903 DefaultRegionForTheWorld
0904
0905 Index : 1 used in the geometry : No
0906 Material : Silicon
0907 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
0908 Energy thresholds : gamma 6.93672 keV e- 546.311 keV e+ 525.593 keV proton 100 keV
0909
0910 Index : 2 used in the geometry : No
0911 Material : Gold
0912 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
0913 Energy thresholds : gamma 112.405 keV e- 2.28837 MeV e+ 2.12466 MeV proton 100 keV
0914
0915 Index : 3 used in the geometry : No
0916 Material : G4_POLYETHYLENE
0917 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
0918 Energy thresholds : gamma 2.25978 keV e- 351.127 keV e+ 341.776 keV proton 100 keV
0919
0920 Index : 4 used in the geometry : No
0921 Material : Aluminium
0922 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
0923 Energy thresholds : gamma 6.91691 keV e- 597.603 keV e+ 574.29 keV proton 100 keV
0924
0925 Index : 5 used in the geometry : Yes
0926 Material : Iron
0927 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
0928 Energy thresholds : gamma 20.7519 keV e- 1.28828 MeV e+ 1.21724 MeV proton 100 keV
0929 Region(s) which use this couple :
0930 DefaultRegionForTheWorld
0931
0932 ==================================================================
0933
0934 ### Run 4 starts.
0935
0936 ------- MixMaxRng engine status -------
0937 Current state vector is:
0938 mixmax state, file version 1.0
0939 N=17 V[N]={864120598200464168, 189420282748189000, 2125101903705486721, 1863466081725669111, 1383960382039097673, 550892631529147885, 2244069534527202395, 933979936320874605, 791647544461653037, 286412522794523854, 1616800613386878212, 478416970401061220, 1548587089976077386, 917812281992815632, 2284924600166281660, 94168853074452663, 1674934093723585259} counter= 15sumtot= 1401971847063908873
0940 ---------------------------------------
0941 --> Event 0 starts.
0942 --> Event 10000 starts.
0943 --> Event 20000 starts.
0944 --> Event 30000 starts.
0945 --> Event 40000 starts.
0946 Run terminated.
0947 Run Summary
0948 Number of events processed : 50000
0949 User=3.840000s Real=3.908897s Sys=0.000000s
0950
0951 ======================== run summary ======================
0952
0953 The run was 50000 mu+ of 100 GeV through 1 m of Iron (density: 7.87 g/cm3 )
0954
0955 Total energy deposit in absorber per event = 2.163 GeV +- 15.06 MeV
0956
0957 -----> Mean dE/dx = 21.63 MeV/cm (2.748 MeV*cm2/g)
0958
0959 From formulas :
0960 restricted dEdx = 11.21 MeV/cm (1.425 MeV*cm2/g)
0961 full dEdx = 21.61 MeV/cm (2.746 MeV*cm2/g)
0962
0963 Leakage : primary = 97.84 GeV +- 15.06 MeV secondaries = 0 eV +- 0 eV
0964 Energy balance : edep + eleak = 100 GeV
0965
0966 Total track length (charged) in absorber per event = 1 m +- 7.06 nm
0967 Total track length (neutral) in absorber per event = 0 fm +- 0 fm
0968
0969 Number of steps (charged) in absorber per event = 47.58 +- 0.03044
0970 Number of steps (neutral) in absorber per event = 0 +- 0
0971
0972 Number of secondaries per event : Gammas = 0.05044 electrons = 46.52 positrons = 1.254
0973
0974 Number of events with the primary particle transmitted = 100 %
0975 Number of events with at least 1 particle transmitted (same charge as primary) = 100 %
0976
0977 Number of events with the primary particle reflected = 0 %
0978 Number of events with at least 1 particle reflected (same charge as primary) = 0 %
0979
0980 MultipleScattering:
0981 rms proj angle of transmit primary particle = 1.178 mrad (central part only)
0982 computed theta0 (Highland formula) = 1.182 mrad
0983 central part defined as +- 3.546 mrad Tail ratio = 0.579 %
0984
0985 Gamma process counts:
0986 Photoeffect 0
0987 Compton 0
0988 Conversion 0
0989 Rayleigh 0
0990
0991 ------- MixMaxRng engine status -------
0992 Current state vector is:
0993 mixmax state, file version 1.0
0994 N=17 V[N]={1046585659322464193, 1700724626221482459, 1710146587566228662, 2274062821915192324, 1093598942862882934, 859975426035031537, 1845450543346184012, 1347610440888355163, 1747669745252223770, 1589127766525861706, 1326010817860615296, 367647672722605028, 58019342554185087, 2022010207749596709, 2304966030133393681, 135872583825542285, 522422170598193944} counter= 6sumtot= 1199314302456793231
0995 ---------------------------------------
0996 #
0997 ================== Deleting memory pools ===================
0998 Number of memory pools allocated: 9 of which, static: 0
0999 Dynamic pools deleted: 9 / Total memory freed: 0.056 MB
1000 ============================================================
