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-02-ref-06 (28-June-2024)
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: 113427284.261 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 Lowest triplet kinetic energy 1 MeV
0092 Enable sampling of gamma linear polarisation 0
0093 5D gamma conversion model type 0
0094 5D gamma conversion model on isolated ion 0
0095 Livermore data directory epics_2017
0096 =======================================================================
0097 ====== Ionisation Parameters ========
0098 =======================================================================
0099 Step function for e+- (0.2, 0.1 mm)
0100 Step function for muons/hadrons (0.2, 0.05 mm)
0101 Step function for light ions (0.1, 0.02 mm)
0102 Step function for general ions (0.1, 0.001 mm)
0103 Lowest e+e- kinetic energy 100 eV
0104 Lowest muon/hadron kinetic energy 1 keV
0105 Use ICRU90 data 1
0106 Fluctuations of dE/dx are enabled 1
0107 Type of fluctuation model for leptons and hadrons Urban
0108 Use built-in Birks satuaration 0
0109 Build CSDA range enabled 0
0110 Use cut as a final range enabled 0
0111 Enable angular generator interface 1
0112 Max kinetic energy for CSDA tables 1 GeV
0113 Max kinetic energy for NIEL computation 1 MeV
0114 Linear loss limit 0.01
0115 Read data from file for e+e- pair production by mu 0
0116 =======================================================================
0117 ====== Multiple Scattering Parameters ========
0118 =======================================================================
0119 Type of msc step limit algorithm for e+- 3
0120 Type of msc step limit algorithm for muons/hadrons 0
0121 Msc lateral displacement for e+- enabled 1
0122 Msc lateral displacement for muons and hadrons 1
0123 Urban msc model lateral displacement alg96 1
0124 Range factor for msc step limit for e+- 0.04
0125 Range factor for msc step limit for muons/hadrons 0.2
0126 Geometry factor for msc step limitation of e+- 2.5
0127 Safety factor for msc step limit for e+- 0.6
0128 Skin parameter for msc step limitation of e+- 1
0129 Lambda limit for msc step limit for e+- 1 mm
0130 Use Mott correction for e- scattering 1
0131 Factor used for dynamic computation of angular
0132 limit between single and multiple scattering 1
0133 Fixed angular limit between single
0134 and multiple scattering 3.1416 rad
0135 Upper energy limit for e+- multiple scattering 100 MeV
0136 Type of electron single scattering model 0
0137 Type of nuclear form-factor 1
0138 Screening factor 1
0139 =======================================================================
0140 ====== Atomic Deexcitation Parameters ========
0141 =======================================================================
0142 Fluorescence enabled 1
0143 Directory in G4LEDATA for fluorescence data files fluor
0144 Auger electron cascade enabled 0
0145 PIXE atomic de-excitation enabled 0
0146 De-excitation module ignores cuts 0
0147 Type of PIXE cross section for hadrons Empirical
0148 Type of PIXE cross section for e+- Livermore
0149 =======================================================================
0150
0151 ### === Deexcitation model UAtomDeexcitation is activated for 1 region:
0152 DefaultRegionForTheWorld 1 0 0
0153 ### === Ignore cuts flag: 0
0154
0155 phot: for gamma SubType=12 BuildTable=0
0156 LambdaPrime table from 200 keV to 100 TeV in 174 bins
0157 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0158 LivermorePhElectric : Emin= 0 eV Emax= 100 TeV SauterGavrila Fluo
0159
0160 compt: for gamma SubType=13 BuildTable=1
0161 Lambda table from 10 eV to 1 MeV, 20 bins/decade, spline: 1
0162 LambdaPrime table from 1 MeV to 100 TeV in 160 bins
0163 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0164 KleinNishina : Emin= 0 eV Emax= 100 TeV Fluo
0165
0166 conv: for gamma SubType=14 BuildTable=1
0167 Lambda table from 1.022 MeV to 100 TeV, 20 bins/decade, spline: 1
0168 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0169 BetheHeitlerLPM : Emin= 0 eV Emax= 100 TeV ModifiedTsai
0170
0171 Rayl: for gamma SubType=11 BuildTable=1
0172 Lambda table from 10 eV to 150 keV, 20 bins/decade, spline: 0
0173 LambdaPrime table from 150 keV to 100 TeV in 176 bins
0174 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0175 LivermoreRayleigh : Emin= 0 eV Emax= 100 TeV CullenGenerator
0176
0177 msc: for e- SubType= 10
0178 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0179 UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV
0180 StepLim=DistanceToBoundary Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0181
0182 eIoni: for e- XStype:3 SubType=2
0183 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0184 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0185 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0186 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0187 MollerBhabha : Emin= 0 eV Emax= 100 TeV deltaVI
0188
0189 eBrem: for e- XStype:4 SubType=3
0190 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0191 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0192 LPM flag: 1 for E > 1 GeV, VertexHighEnergyTh(GeV)= 100000
0193 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0194 eBremSB : Emin= 0 eV Emax= 1 GeV AngularGen2BS
0195 eBremLPM : Emin= 1 GeV Emax= 100 TeV AngularGen2BS
0196
0197 ePairProd: for e- XStype:1 SubType=4
0198 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0199 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 0
0200 Sampling table 25x1001 from 0.1 GeV to 100 TeV
0201 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0202 ePairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0203
0204 msc: for e+ SubType= 10
0205 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0206 UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV
0207 StepLim=DistanceToBoundary Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0208
0209 eIoni: for e+ XStype:3 SubType=2
0210 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0211 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0212 StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
0213 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0214 MollerBhabha : Emin= 0 eV Emax= 100 TeV deltaVI
0215
0216 eBrem: for e+ XStype:4 SubType=3
0217 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0218 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0219 LPM flag: 1 for E > 1 GeV, VertexHighEnergyTh(GeV)= 100000
0220 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0221 eBremSB : Emin= 0 eV Emax= 1 GeV AngularGen2BS
0222 eBremLPM : Emin= 1 GeV Emax= 100 TeV AngularGen2BS
0223
0224 ePairProd: for e+ XStype:1 SubType=4
0225 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0226 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 0
0227 Sampling table 25x1001 from 0.1 GeV to 100 TeV
0228 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0229 ePairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0230
0231 annihil: for e+ XStype:2 SubType=5 AtRestModel:Allison BuildTable=0
0232 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0233 eplus2gg : Emin= 0 eV Emax= 100 TeV
0234
0235 msc: for proton SubType= 10
0236 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0237 UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV
0238 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0239
0240 hIoni: for proton XStype:3 SubType=2
0241 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0242 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0243 StepFunction=(0.2, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
0244 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0245 Bragg : Emin= 0 eV Emax= 2 MeV deltaVI
0246 BetheBloch : Emin= 2 MeV Emax= 100 TeV deltaVI
0247
0248 hBrems: for proton XStype:1 SubType=3
0249 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0250 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0251 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0252 hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0253
0254 hPairProd: for proton XStype:1 SubType=4
0255 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0256 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0257 Sampling table 17x1001 from 7.50618 GeV to 100 TeV
0258 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0259 hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0260
0261 nuclearStopping: for proton SubType=8 BuildTable=0
0262 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0263 ICRU49NucStopping : Emin= 0 eV Emax= 1 MeV
0264
0265 msc: for GenericIon SubType= 10
0266 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0267 UrbanMsc : Emin= 0 eV Emax= 100 TeV
0268 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0269
0270 ionIoni: for GenericIon XStype:3 SubType=2
0271 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0272 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0273 StepFunction=(0.1, 0.001 mm), integ: 3, fluct: 1, linLossLim= 0.02
0274 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0275 LindhardSorensen : Emin= 0 eV Emax= 100 TeV deltaVI
0276
0277 nuclearStopping: for GenericIon SubType=8 BuildTable=0
0278 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0279 ICRU49NucStopping : Emin= 0 eV Emax= 1 MeV
0280
0281 msc: for alpha SubType= 10
0282 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0283 UrbanMsc : Emin= 0 eV Emax= 100 TeV
0284 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0285
0286 ionIoni: for alpha XStype:3 SubType=2
0287 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0288 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0289 StepFunction=(0.1, 0.02 mm), integ: 3, fluct: 1, linLossLim= 0.02
0290 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0291 BraggIon : Emin= 0 eV Emax=7.9452 MeV deltaVI
0292 BetheBloch : Emin=7.9452 MeV Emax= 100 TeV deltaVI
0293
0294 nuclearStopping: for alpha SubType=8 BuildTable=0
0295 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0296 ICRU49NucStopping : Emin= 0 eV Emax= 1 MeV
0297
0298 msc: for anti_proton SubType= 10
0299 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0300 UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV
0301 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0302
0303 hIoni: for anti_proton XStype:3 SubType=2
0304 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0305 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0306 StepFunction=(0.2, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
0307 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0308 ICRU73QO : Emin= 0 eV Emax= 2 MeV deltaVI
0309 BetheBloch : Emin= 2 MeV Emax= 100 TeV deltaVI
0310
0311 hBrems: for anti_proton XStype:1 SubType=3
0312 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0313 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0314 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0315 hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0316
0317 hPairProd: for anti_proton XStype:1 SubType=4
0318 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0319 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0320 Sampling table 17x1001 from 7.50618 GeV to 100 TeV
0321 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0322 hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0323
0324 msc: for kaon+ SubType= 10
0325 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0326 UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV
0327 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0328
0329 hIoni: for kaon+ XStype:3 SubType=2
0330 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0331 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0332 StepFunction=(0.2, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
0333 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0334 Bragg : Emin= 0 eV Emax=1.05231 MeV deltaVI
0335 BetheBloch : Emin=1.05231 MeV Emax= 100 TeV deltaVI
0336
0337 hBrems: for kaon+ XStype:1 SubType=3
0338 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0339 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0340 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0341 hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0342
0343 hPairProd: for kaon+ XStype:1 SubType=4
0344 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0345 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0346 Sampling table 18x1001 from 3.94942 GeV to 100 TeV
0347 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0348 hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0349
0350 msc: for kaon- SubType= 10
0351 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0352 UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV
0353 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0354
0355 hIoni: for kaon- XStype:3 SubType=2
0356 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0357 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0358 StepFunction=(0.2, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
0359 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0360 ICRU73QO : Emin= 0 eV Emax=1.05231 MeV deltaVI
0361 BetheBloch : Emin=1.05231 MeV Emax= 100 TeV deltaVI
0362
0363 hBrems: for kaon- XStype:1 SubType=3
0364 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0365 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0366 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0367 hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0368
0369 hPairProd: for kaon- XStype:1 SubType=4
0370 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0371 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0372 Sampling table 18x1001 from 3.94942 GeV to 100 TeV
0373 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0374 hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0375
0376 msc: for mu+ SubType= 10
0377 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0378 UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV
0379 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0380
0381 muIoni: for mu+ XStype:3 SubType=2
0382 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0383 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0384 StepFunction=(0.2, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
0385 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0386 Bragg : Emin= 0 eV Emax= 200 keV deltaVI
0387 MuBetheBloch : Emin= 200 keV Emax= 100 TeV deltaVI
0388
0389 muBrems: for mu+ XStype:1 SubType=3
0390 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0391 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0392 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0393 MuBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0394
0395 muPairProd: for mu+ XStype:1 SubType=4
0396 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0397 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0398 Sampling table 21x1001 from 0.85 GeV to 100 TeV
0399 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0400 muPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0401
0402 msc: for mu- SubType= 10
0403 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0404 UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV
0405 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0406
0407 muIoni: for mu- XStype:3 SubType=2
0408 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0409 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0410 StepFunction=(0.2, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
0411 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0412 ICRU73QO : Emin= 0 eV Emax= 200 keV deltaVI
0413 MuBetheBloch : Emin= 200 keV Emax= 100 TeV deltaVI
0414
0415 muBrems: for mu- XStype:1 SubType=3
0416 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0417 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0418 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0419 MuBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0420
0421 muPairProd: for mu- XStype:1 SubType=4
0422 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0423 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0424 Sampling table 21x1001 from 0.85 GeV to 100 TeV
0425 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0426 muPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0427
0428 msc: for pi+ SubType= 10
0429 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0430 UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV
0431 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0432
0433 hIoni: for pi+ XStype:3 SubType=2
0434 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0435 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0436 StepFunction=(0.2, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
0437 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0438 Bragg : Emin= 0 eV Emax=297.505 keV deltaVI
0439 BetheBloch : Emin=297.505 keV Emax= 100 TeV deltaVI
0440
0441 hBrems: for pi+ XStype:1 SubType=3
0442 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0443 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0444 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0445 hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0446
0447 hPairProd: for pi+ XStype:1 SubType=4
0448 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0449 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0450 Sampling table 20x1001 from 1.11656 GeV to 100 TeV
0451 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0452 hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0453
0454 msc: for pi- SubType= 10
0455 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0456 UrbanMsc : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV
0457 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
0458
0459 hIoni: for pi- XStype:3 SubType=2
0460 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0461 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0462 StepFunction=(0.2, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
0463 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0464 ICRU73QO : Emin= 0 eV Emax=297.505 keV deltaVI
0465 BetheBloch : Emin=297.505 keV Emax= 100 TeV deltaVI
0466
0467 hBrems: for pi- XStype:1 SubType=3
0468 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0469 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0470 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0471 hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0472
0473 hPairProd: for pi- XStype:1 SubType=4
0474 dE/dx and range tables from 10 eV to 100 TeV in 260 bins
0475 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0476 Sampling table 20x1001 from 1.11656 GeV to 100 TeV
0477 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0478 hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0479
0480 ========= Table of registered couples ============================
0481
0482 Index : 0 used in the geometry : Yes
0483 Material : G4_Galactic
0484 Range cuts : gamma 7 um e- 7 um e+ 7 um proton 7 um
0485 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 700 eV
0486 Region(s) which use this couple :
0487 DefaultRegionForTheWorld
0488
0489 Index : 1 used in the geometry : Yes
0490 Material : Silicon
0491 Range cuts : gamma 7 um e- 7 um e+ 7 um proton 7 um
0492 Energy thresholds : gamma 990 eV e- 21.0638 keV e+ 20.5957 keV proton 700 eV
0493 Region(s) which use this couple :
0494 DefaultRegionForTheWorld
0495
0496 ==================================================================
0497
0498 ### Run 0 starts.
0499
0500 ------- MixMaxRng engine status -------
0501 Current state vector is:
0502 mixmax state, file version 1.0
0503 N=17 V[N]={906770732717044781, 629165745432651234, 1235682547346241386, 68420008233404568, 2262190991329242458, 2266470399991071809, 1976726662926872232, 245458862506414172, 1955974201201518530, 2155248512522080758, 604170912935414061, 1116171330120743511, 1861018313684488333, 1296715403254578286, 1549011045957234151, 370819759640195970, 2230139271784837643} counter= 17sumtot= 1977567618660788324
0504 ---------------------------------------
0505 --> Event 0 starts.
0506 --> Event 10000 starts.
0507 --> Event 20000 starts.
0508 --> Event 30000 starts.
0509 --> Event 40000 starts.
0510 --> Event 50000 starts.
0511 --> Event 60000 starts.
0512 --> Event 70000 starts.
0513 --> Event 80000 starts.
0514 --> Event 90000 starts.
0515 Run terminated.
0516 Run Summary
0517 Number of events processed : 100000
0518 User=0.890000s Real=0.938960s Sys=0.000000s
0519
0520 ======================== run summary ======================
0521
0522 The run was 100000 pi+ of 5 GeV through 20 um of Silicon (density: 2.33 g/cm3 )
0523
0524 Total energy deposit in absorber per event = 6.346 keV +- 16.23 eV
0525
0526 -----> Mean dE/dx = 3.173 MeV/cm (1.362 MeV*cm2/g)
0527
0528 From formulas :
0529 restricted dEdx = 2.96 MeV/cm (1.27 MeV*cm2/g)
0530 full dEdx = 4.727 MeV/cm (2.029 MeV*cm2/g)
0531
0532 Leakage : primary = 5 GeV +- 1.301 keV secondaries = 3.621 keV +- 1.301 keV
0533 Energy balance : edep + eleak = 5 GeV
0534
0535 Total track length (charged) in absorber per event = 20.26 um +- 7.733 nm
0536 Total track length (neutral) in absorber per event = 1.273 Ang +- 7.208e+04 fm
0537
0538 Number of steps (charged) in absorber per event = 1.129 +- 0.003973
0539 Number of steps (neutral) in absorber per event = 4e-05 +- 2e-05
0540
0541 Number of secondaries per event : Gammas = 4e-05 electrons = 0.01731 positrons = 0
0542
0543 Number of events with the primary particle transmitted = 100 %
0544 Number of events with at least 1 particle transmitted (same charge as primary) = 100 %
0545
0546 Number of events with the primary particle reflected = 0 %
0547 Number of events with at least 1 particle reflected (same charge as primary) = 0 %
0548
0549 MultipleScattering:
0550 rms proj angle of transmit primary particle = 0.02511 mrad (central part only)
0551 computed theta0 (Highland formula) = 0.02627 mrad
0552 central part defined as +- 0.07881 mrad Tail ratio = 1.945 %
0553
0554 Gamma process counts:
0555 Photoeffect 3
0556 Compton 0
0557 Conversion 0
0558 Rayleigh 0
0559
0560 ------- MixMaxRng engine status -------
0561 Current state vector is:
0562 mixmax state, file version 1.0
0563 N=17 V[N]={821612629383806282, 1355057195926553987, 1269729022158354381, 44614384821324673, 1303746237006252945, 287979836684896906, 2187769835849609370, 723285766230658542, 843848951070026816, 411140667534285683, 1214416906255647570, 627651330024082620, 2171256320394340744, 456036330617358175, 101021189349451762, 395972390457516759, 480661350402717899} counter= 4sumtot= 860742288884721408
0564 ---------------------------------------
0565 #
0566 /testem/det/setAbsMat Gold
0567 /testem/det/setAbsThick 9.658 um
0568 /testem/det/setAbsYZ 1 mm
0569 /testem/gun/setDefault
0570 # Field Should be issued after /run/initialize
0571 #/testem/det/setField 1 tesla
0572 /gun/particle e-
0573 /gun/energy 15.7 MeV
0574 /run/beamOn 50000
0575
0576 ### === Deexcitation model UAtomDeexcitation is activated for 1 region:
0577 DefaultRegionForTheWorld 1 0 0
0578 ### === Ignore cuts flag: 0
0579
0580 ========= Table of registered couples ============================
0581
0582 Index : 0 used in the geometry : Yes
0583 Material : G4_Galactic
0584 Range cuts : gamma 7 um e- 7 um e+ 7 um proton 7 um
0585 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 700 eV
0586 Region(s) which use this couple :
0587 DefaultRegionForTheWorld
0588
0589 Index : 1 used in the geometry : No
0590 Material : Silicon
0591 Range cuts : gamma 7 um e- 7 um e+ 7 um proton 7 um
0592 Energy thresholds : gamma 990 eV e- 21.0638 keV e+ 20.5957 keV proton 700 eV
0593
0594 Index : 2 used in the geometry : Yes
0595 Material : Gold
0596 Range cuts : gamma 7 um e- 7 um e+ 7 um proton 7 um
0597 Energy thresholds : gamma 6.52544 keV e- 65.1638 keV e+ 63.9139 keV proton 700 eV
0598 Region(s) which use this couple :
0599 DefaultRegionForTheWorld
0600
0601 ==================================================================
0602
0603 ### Run 1 starts.
0604
0605 ------- MixMaxRng engine status -------
0606 Current state vector is:
0607 mixmax state, file version 1.0
0608 N=17 V[N]={821612629383806282, 1355057195926553987, 1269729022158354381, 44614384821324673, 1303746237006252945, 287979836684896906, 2187769835849609370, 723285766230658542, 843848951070026816, 411140667534285683, 1214416906255647570, 627651330024082620, 2171256320394340744, 456036330617358175, 101021189349451762, 395972390457516759, 480661350402717899} counter= 4sumtot= 860742288884721408
0609 ---------------------------------------
0610 --> Event 0 starts.
0611 --> Event 10000 starts.
0612 --> Event 20000 starts.
0613 --> Event 30000 starts.
0614 --> Event 40000 starts.
0615 Run terminated.
0616 Run Summary
0617 Number of events processed : 50000
0618 User=0.650000s Real=0.669512s Sys=0.000000s
0619
0620 ======================== run summary ======================
0621
0622 The run was 50000 e- of 15.7 MeV through 9.66 um of Gold (density: 19.3 g/cm3 )
0623
0624 Total energy deposit in absorber per event = 18.5 keV +- 65.14 eV
0625
0626 -----> Mean dE/dx = 19.16 MeV/cm (0.9915 MeV*cm2/g)
0627
0628 From formulas :
0629 restricted dEdx = 17.91 MeV/cm (0.927 MeV*cm2/g)
0630 full dEdx = 61.41 MeV/cm (3.179 MeV*cm2/g)
0631
0632 Leakage : primary = 15.64 MeV +- 2.184 keV secondaries = 40.21 keV +- 2.179 keV
0633 Energy balance : edep + eleak = 15.7 MeV
0634
0635 Total track length (charged) in absorber per event = 9.932 um +- 12.05 nm
0636 Total track length (neutral) in absorber per event = 130.4 nm +- 8.558 nm
0637
0638 Number of steps (charged) in absorber per event = 2.708 +- 0.00813
0639 Number of steps (neutral) in absorber per event = 0.02784 +- 0.0007708
0640
0641 Number of secondaries per event : Gammas = 0.02774 electrons = 0.02262 positrons = 0
0642
0643 Number of events with the primary particle transmitted = 99.99 %
0644 Number of events with at least 1 particle transmitted (same charge as primary) = 99.99 %
0645
0646 Number of events with the primary particle reflected = 0.008 %
0647 Number of events with at least 1 particle reflected (same charge as primary) = 0.322 %
0648
0649 MultipleScattering:
0650 rms proj angle of transmit primary particle = 34.83 mrad (central part only)
0651 computed theta0 (Highland formula) = 35.1 mrad
0652 central part defined as +- 105.3 mrad Tail ratio = 2.914 %
0653
0654 Gamma process counts:
0655 Photoeffect 233
0656 Compton 0
0657 Conversion 0
0658 Rayleigh 5
0659
0660 ------- MixMaxRng engine status -------
0661 Current state vector is:
0662 mixmax state, file version 1.0
0663 N=17 V[N]={1350763663012815612, 1294348897496478787, 1497574536371361083, 1311753551330282605, 1604272835801072024, 478382643570230237, 174370242934593580, 1118323367948153675, 183190413119826769, 6592242593151900, 954709532099726697, 1955739663614221939, 708984632719044540, 2058128678670554479, 968441555857833481, 1086034042937008217, 1747224124524366416} counter= 14sumtot= 52090550891170433
0664 ---------------------------------------
0665 #
0666 /run/setCut 10 um
0667 #
0668 /testem/det/setAbsMat G4_POLYETHYLENE
0669 /testem/det/setAbsThick 4.74 mm
0670 /testem/det/setAbsYZ 50 mm
0671 /testem/gun/setDefault
0672 /gun/particle mu+
0673 /gun/energy 96.2 MeV
0674 /run/beamOn 50000
0675
0676 ### === Deexcitation model UAtomDeexcitation is activated for 1 region:
0677 DefaultRegionForTheWorld 1 0 0
0678 ### === Ignore cuts flag: 0
0679
0680 ========= Table of registered couples ============================
0681
0682 Index : 0 used in the geometry : Yes
0683 Material : G4_Galactic
0684 Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um
0685 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 1 keV
0686 Region(s) which use this couple :
0687 DefaultRegionForTheWorld
0688
0689 Index : 1 used in the geometry : No
0690 Material : Silicon
0691 Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um
0692 Energy thresholds : gamma 990 eV e- 31.6994 keV e+ 31.3501 keV proton 1 keV
0693
0694 Index : 2 used in the geometry : No
0695 Material : Gold
0696 Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um
0697 Energy thresholds : gamma 8.32988 keV e- 80.7452 keV e+ 79.0938 keV proton 1 keV
0698
0699 Index : 3 used in the geometry : Yes
0700 Material : G4_POLYETHYLENE
0701 Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um
0702 Energy thresholds : gamma 990 eV e- 13.7224 keV e+ 13.4607 keV proton 1 keV
0703 Region(s) which use this couple :
0704 DefaultRegionForTheWorld
0705
0706 ==================================================================
0707
0708 ### Run 2 starts.
0709
0710 ------- MixMaxRng engine status -------
0711 Current state vector is:
0712 mixmax state, file version 1.0
0713 N=17 V[N]={1350763663012815612, 1294348897496478787, 1497574536371361083, 1311753551330282605, 1604272835801072024, 478382643570230237, 174370242934593580, 1118323367948153675, 183190413119826769, 6592242593151900, 954709532099726697, 1955739663614221939, 708984632719044540, 2058128678670554479, 968441555857833481, 1086034042937008217, 1747224124524366416} counter= 14sumtot= 52090550891170433
0714 ---------------------------------------
0715 --> Event 0 starts.
0716 --> Event 10000 starts.
0717 --> Event 20000 starts.
0718 --> Event 30000 starts.
0719 --> Event 40000 starts.
0720 Run terminated.
0721 Run Summary
0722 Number of events processed : 50000
0723 User=1.820000s Real=1.980485s Sys=0.000000s
0724
0725 ======================== run summary ======================
0726
0727 The run was 50000 mu+ of 96.2 MeV through 4.74 mm of G4_POLYETHYLENE (density: 940 mg/cm3)
0728
0729 Total energy deposit in absorber per event = 1.053 MeV +- 987.7 eV
0730
0731 -----> Mean dE/dx = 2.221 MeV/cm (2.363 MeV*cm2/g)
0732
0733 From formulas :
0734 restricted dEdx = 1.761 MeV/cm (1.873 MeV*cm2/g)
0735 full dEdx = 2.28 MeV/cm (2.425 MeV*cm2/g)
0736
0737 Leakage : primary = 95.12 MeV +- 1.349 keV secondaries = 27.37 keV +- 729.9 eV
0738 Energy balance : edep + eleak = 96.2 MeV
0739
0740 Total track length (charged) in absorber per event = 5.168 mm +- 3.713 um
0741 Total track length (neutral) in absorber per event = 18.36 um +- 2.016 um
0742
0743 Number of steps (charged) in absorber per event = 11.46 +- 0.03112
0744 Number of steps (neutral) in absorber per event = 0.00928 +- 0.0004498
0745
0746 Number of secondaries per event : Gammas = 0.0089 electrons = 4.187 positrons = 0
0747
0748 Number of events with the primary particle transmitted = 100 %
0749 Number of events with at least 1 particle transmitted (same charge as primary) = 100 %
0750
0751 Number of events with the primary particle reflected = 0.004 %
0752 Number of events with at least 1 particle reflected (same charge as primary) = 0.004 %
0753
0754 MultipleScattering:
0755 rms proj angle of transmit primary particle = 7.321 mrad (central part only)
0756 computed theta0 (Highland formula) = 7.636 mrad
0757 central part defined as +- 22.91 mrad Tail ratio = 1.233 %
0758
0759 Gamma process counts:
0760 Photoeffect 290
0761 Compton 16
0762 Conversion 0
0763 Rayleigh 3
0764
0765 ------- MixMaxRng engine status -------
0766 Current state vector is:
0767 mixmax state, file version 1.0
0768 N=17 V[N]={259808161633608381, 1699288445904995528, 1090904378916437161, 706721855983804764, 650736876350838486, 2058478167238721866, 1715703218711355282, 390147856004027483, 1510430911949678032, 2292085793239906216, 1954290827584041406, 1081447968060984640, 598848276166094160, 570192178062982954, 971153001893022519, 484014818966045887, 1454776790956935882} counter= 3sumtot= 1042285453913929039
0769 ---------------------------------------
0770 #
0771 /run/setCut 1 mm
0772 #
0773 /testem/stack/killSecondaries
0774 #
0775 /testem/det/setAbsMat Aluminium
0776 /testem/det/setAbsThick 8.004 mm
0777 /testem/gun/setDefault
0778 /gun/particle proton
0779 /gun/energy 174.664 GeV
0780 /run/beamOn 50000
0781
0782 ### === Deexcitation model UAtomDeexcitation is activated for 1 region:
0783 DefaultRegionForTheWorld 1 0 0
0784 ### === Ignore cuts flag: 0
0785
0786 ========= Table of registered couples ============================
0787
0788 Index : 0 used in the geometry : Yes
0789 Material : G4_Galactic
0790 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
0791 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 100 keV
0792 Region(s) which use this couple :
0793 DefaultRegionForTheWorld
0794
0795 Index : 1 used in the geometry : No
0796 Material : Silicon
0797 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
0798 Energy thresholds : gamma 6.93672 keV e- 546.311 keV e+ 525.593 keV proton 100 keV
0799
0800 Index : 2 used in the geometry : No
0801 Material : Gold
0802 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
0803 Energy thresholds : gamma 112.405 keV e- 2.28837 MeV e+ 2.12466 MeV proton 100 keV
0804
0805 Index : 3 used in the geometry : No
0806 Material : G4_POLYETHYLENE
0807 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
0808 Energy thresholds : gamma 2.25978 keV e- 351.127 keV e+ 341.776 keV proton 100 keV
0809
0810 Index : 4 used in the geometry : Yes
0811 Material : Aluminium
0812 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
0813 Energy thresholds : gamma 6.91691 keV e- 597.603 keV e+ 574.29 keV proton 100 keV
0814 Region(s) which use this couple :
0815 DefaultRegionForTheWorld
0816
0817 ==================================================================
0818
0819 ### Run 3 starts.
0820
0821 ------- MixMaxRng engine status -------
0822 Current state vector is:
0823 mixmax state, file version 1.0
0824 N=17 V[N]={259808161633608381, 1699288445904995528, 1090904378916437161, 706721855983804764, 650736876350838486, 2058478167238721866, 1715703218711355282, 390147856004027483, 1510430911949678032, 2292085793239906216, 1954290827584041406, 1081447968060984640, 598848276166094160, 570192178062982954, 971153001893022519, 484014818966045887, 1454776790956935882} counter= 3sumtot= 1042285453913929039
0825 ---------------------------------------
0826 --> Event 0 starts.
0827 --> Event 10000 starts.
0828 --> Event 20000 starts.
0829 --> Event 30000 starts.
0830 --> Event 40000 starts.
0831 Run terminated.
0832 Run Summary
0833 Number of events processed : 50000
0834 User=0.480000s Real=0.511424s Sys=0.000000s
0835
0836 ======================== run summary ======================
0837
0838 The run was 50000 proton of 175 GeV through 8 mm of Aluminium (density: 2.7 g/cm3 )
0839
0840 Total energy deposit in absorber per event = 5.161 MeV +- 331.3 keV
0841
0842 -----> Mean dE/dx = 6.449 MeV/cm (2.388 MeV*cm2/g)
0843
0844 From formulas :
0845 restricted dEdx = 4.024 MeV/cm (1.49 MeV*cm2/g)
0846 full dEdx = 6.049 MeV/cm (2.24 MeV*cm2/g)
0847
0848 Leakage : primary = 174.7 GeV +- 331.4 keV secondaries = 0 eV +- 0 eV
0849 Energy balance : edep + eleak = 174.7 GeV
0850
0851 Total track length (charged) in absorber per event = 8.004 mm +- 0 fm
0852 Total track length (neutral) in absorber per event = 0 fm +- 0 fm
0853
0854 Number of steps (charged) in absorber per event = 1.264 +- 0.002314
0855 Number of steps (neutral) in absorber per event = 0 +- 0
0856
0857 Number of secondaries per event : Gammas = 0 electrons = 0.2644 positrons = 0.00072
0858
0859 Number of events with the primary particle transmitted = 100 %
0860 Number of events with at least 1 particle transmitted (same charge as primary) = 100 %
0861
0862 Number of events with the primary particle reflected = 0 %
0863 Number of events with at least 1 particle reflected (same charge as primary) = 0 %
0864
0865 MultipleScattering:
0866 rms proj angle of transmit primary particle = 0.02057 mrad (central part only)
0867 computed theta0 (Highland formula) = 0.02111 mrad
0868 central part defined as +- 0.06333 mrad Tail ratio = 1.352 %
0869
0870 Gamma process counts:
0871 Photoeffect 0
0872 Compton 0
0873 Conversion 0
0874 Rayleigh 0
0875
0876 ------- MixMaxRng engine status -------
0877 Current state vector is:
0878 mixmax state, file version 1.0
0879 N=17 V[N]={208646914259207529, 44936455847492391, 617853924538367756, 122840575292410615, 554625140930138278, 425420266314921049, 1100409250128904939, 88037367507246226, 886917685562055761, 484189764832677705, 1494479114852379083, 856245190962181511, 63808615407432736, 1760858753733414007, 1397402278650583195, 1407199650688570504, 212537499099074009} counter= 10sumtot= 197193402538587539
0880 ---------------------------------------
0881 #
0882 /testem/det/setAbsMat Iron
0883 /testem/det/setAbsThick 1 m
0884 /testem/gun/setDefault
0885 /gun/particle mu+
0886 /gun/energy 100 GeV
0887 /run/beamOn 50000
0888
0889 ### === Deexcitation model UAtomDeexcitation is activated for 1 region:
0890 DefaultRegionForTheWorld 1 0 0
0891 ### === Ignore cuts flag: 0
0892
0893 ========= Table of registered couples ============================
0894
0895 Index : 0 used in the geometry : Yes
0896 Material : G4_Galactic
0897 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
0898 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 100 keV
0899 Region(s) which use this couple :
0900 DefaultRegionForTheWorld
0901
0902 Index : 1 used in the geometry : No
0903 Material : Silicon
0904 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
0905 Energy thresholds : gamma 6.93672 keV e- 546.311 keV e+ 525.593 keV proton 100 keV
0906
0907 Index : 2 used in the geometry : No
0908 Material : Gold
0909 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
0910 Energy thresholds : gamma 112.405 keV e- 2.28837 MeV e+ 2.12466 MeV proton 100 keV
0911
0912 Index : 3 used in the geometry : No
0913 Material : G4_POLYETHYLENE
0914 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
0915 Energy thresholds : gamma 2.25978 keV e- 351.127 keV e+ 341.776 keV proton 100 keV
0916
0917 Index : 4 used in the geometry : No
0918 Material : Aluminium
0919 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
0920 Energy thresholds : gamma 6.91691 keV e- 597.603 keV e+ 574.29 keV proton 100 keV
0921
0922 Index : 5 used in the geometry : Yes
0923 Material : Iron
0924 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
0925 Energy thresholds : gamma 20.7519 keV e- 1.28828 MeV e+ 1.21724 MeV proton 100 keV
0926 Region(s) which use this couple :
0927 DefaultRegionForTheWorld
0928
0929 ==================================================================
0930
0931 ### Run 4 starts.
0932
0933 ------- MixMaxRng engine status -------
0934 Current state vector is:
0935 mixmax state, file version 1.0
0936 N=17 V[N]={208646914259207529, 44936455847492391, 617853924538367756, 122840575292410615, 554625140930138278, 425420266314921049, 1100409250128904939, 88037367507246226, 886917685562055761, 484189764832677705, 1494479114852379083, 856245190962181511, 63808615407432736, 1760858753733414007, 1397402278650583195, 1407199650688570504, 212537499099074009} counter= 10sumtot= 197193402538587539
0937 ---------------------------------------
0938 --> Event 0 starts.
0939 --> Event 10000 starts.
0940 --> Event 20000 starts.
0941 --> Event 30000 starts.
0942 --> Event 40000 starts.
0943 Run terminated.
0944 Run Summary
0945 Number of events processed : 50000
0946 User=7.540000s Real=7.847703s Sys=0.010000s
0947
0948 ======================== run summary ======================
0949
0950 The run was 50000 mu+ of 100 GeV through 1 m of Iron (density: 7.87 g/cm3 )
0951
0952 Total energy deposit in absorber per event = 2.166 GeV +- 14.98 MeV
0953
0954 -----> Mean dE/dx = 21.66 MeV/cm (2.752 MeV*cm2/g)
0955
0956 From formulas :
0957 restricted dEdx = 11.21 MeV/cm (1.425 MeV*cm2/g)
0958 full dEdx = 21.61 MeV/cm (2.746 MeV*cm2/g)
0959
0960 Leakage : primary = 97.83 GeV +- 14.99 MeV secondaries = 0 eV +- 0 eV
0961 Energy balance : edep + eleak = 100 GeV
0962
0963 Total track length (charged) in absorber per event = 100 cm +- 7.917 um
0964 Total track length (neutral) in absorber per event = 0 fm +- 0 fm
0965
0966 Number of steps (charged) in absorber per event = 47.58 +- 0.03057
0967 Number of steps (neutral) in absorber per event = 0 +- 0
0968
0969 Number of secondaries per event : Gammas = 0.0488 electrons = 46.53 positrons = 1.27
0970
0971 Number of events with the primary particle transmitted = 100 %
0972 Number of events with at least 1 particle transmitted (same charge as primary) = 100 %
0973
0974 Number of events with the primary particle reflected = 0 %
0975 Number of events with at least 1 particle reflected (same charge as primary) = 0 %
0976
0977 MultipleScattering:
0978 rms proj angle of transmit primary particle = 1.173 mrad (central part only)
0979 computed theta0 (Highland formula) = 1.182 mrad
0980 central part defined as +- 3.546 mrad Tail ratio = 0.562 %
0981
0982 Gamma process counts:
0983 Photoeffect 0
0984 Compton 0
0985 Conversion 0
0986 Rayleigh 0
0987
0988 ------- MixMaxRng engine status -------
0989 Current state vector is:
0990 mixmax state, file version 1.0
0991 N=17 V[N]={1677103876344921596, 2145489090854365033, 481303451215822613, 1616489373288118465, 384234704660471306, 839699839614129357, 1492764307969312906, 821359852248163654, 784933776173016720, 1581977921863174715, 1177477298979774340, 2171462053635775432, 513278174562618836, 299019185349657766, 561025169451506329, 147941292274911013, 1054042942077394166} counter= 14sumtot= 1608701246067276590
0992 ---------------------------------------
0993 #
0994 ================== Deleting memory pools ===================
0995 Number of memory pools allocated: 9 of which, static: 0
0996 Dynamic pools deleted: 9 / Total memory freed: 0.053 MB
0997 ============================================================
