source: trunk/examples/extended/electromagnetic/TestEm5/TestEm5.out@ 1346

Last change on this file since 1346 was 1342, checked in by garnier, 15 years ago

update ti head

File size: 95.8 KB
RevLine 
[807]1
[1337]2 ############################################
3 !!! WARNING - FPE detection is activated !!!
4 ############################################
5
[807]6*************************************************************
[1342]7 Geant4 version Name: geant4-09-03-ref-09 (25-June-2010)
[807]8 Copyright : Geant4 Collaboration
9 Reference : NIM A 506 (2003), 250-303
10 WWW : http://cern.ch/geant4
11*************************************************************
12
13/run/verbose 2
14#
15/testem/det/setAbsMat Silicon
16/testem/det/setAbsThick 20 um
17/testem/det/setAbsYZ 50 um
18#
[1337]19/testem/phys/addPhysics emstandard_opt3
20PhysicsList::AddPhysicsList: <emstandard_opt3>
[807]21#
[1337]22/testem/phys/setCuts 7 um
[807]23#
24/run/initialize
25userDetector->Construct() start.
26
[1342]27 Material: Galactic density: 0.000 kg/m3 RadL: 204727512.315 pc Nucl.Int.Length: 114561548.020 pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm
[807]28 ---> Element: Galactic ( ) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
29
[1342]30 Material: Silicon density: 2.330 g/cm3 RadL: 9.368 cm Nucl.Int.Length: 45.761 cm Imean: 173.000 eV
[807]31 ---> Element: Silicon ( ) Z = 14.0 N = 28.1 A = 28.09 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
32
33 The WORLD is made of 30 um of Galactic. The transverse size (YZ) of the world is 60 um
34 The ABSORBER is made of 20 um of Silicon. The transverse size (YZ) is 50 um
35 X position of the middle of the absorber 0 fm
36World is registered to the default region.
37physicsList->Construct() start.
38physicsList->Construct() start.
[1337]39physicsList->CheckParticleList() start.
[807]40physicsList->setCut() start.
41PhysicsList::SetCuts:CutLength : 1 mm
42#
43/testem/gun/setDefault
44/gun/particle pi+
45/gun/energy 5 GeV
46#
47/testem/histo/setFileName testem5
[1337]48/testem/histo/setFileType root
[807]49/testem/histo/setHisto 13 100 -0.25 0.25 mrad
50----> SetHisto 13: (transmit, charged) : projected angle at exit (mrad); 100 bins from -0.25 mrad to 0.25 mrad
51#
52/run/beamOn 100000
53
[1337]54phot: for gamma SubType= 12
55 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]56 PhotoElectric : Emin= 0 eV Emax= 10 TeV
[807]57
[1337]58compt: for gamma SubType= 13
59 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
60 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]61 Klein-Nishina : Emin= 0 eV Emax= 10 TeV
[807]62
[1337]63conv: for gamma SubType= 14
64 Lambda tables from 1.022 MeV to 10 TeV in 220 bins, spline: 1
65 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]66 BetheHeitler : Emin= 0 eV Emax= 10 TeV
[807]67
[1337]68Rayl: for gamma SubType= 11
69 Lambda tables from 100 eV to 10 TeV in 200 bins, spline: 1
70 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]71 LivermoreRayleigh : Emin= 0 eV Emax= 100 GeV
[807]72
[1337]73msc: for e- SubType= 10
74 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
75 RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5
76 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]77 UrbanMsc93 : Emin= 0 eV Emax= 10 TeV
[807]78
[1337]79eIoni: for e- SubType= 2
80 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
81 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
82 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
83 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]84 MollerBhabha : Emin= 0 eV Emax= 10 TeV
[807]85
[1337]86eBrem: for e- SubType= 3
87 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
88 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
89 LPM flag: 1 for E > 1 GeV
90 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]91 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban
92 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban
[807]93
[1337]94eIoni: for e+ SubType= 2
95 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
96 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
97 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
98 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]99 MollerBhabha : Emin= 0 eV Emax= 10 TeV
[807]100
[1337]101eBrem: for e+ SubType= 3
102 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
103 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
104 LPM flag: 1 for E > 1 GeV
105 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]106 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban
107 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban
[807]108
[1337]109annihil: for e+ SubType= 5
110 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
111 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]112 eplus2gg : Emin= 0 eV Emax= 10 TeV
[807]113
[1337]114msc: for proton SubType= 10
115 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
116 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
117 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]118 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV
[807]119
[1337]120hIoni: for proton SubType= 2
121 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
122 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
123 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
124 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]125 Bragg : Emin= 0 eV Emax= 2 MeV
126 BetheBloch : Emin= 2 MeV Emax= 10 TeV
[807]127
[1337]128hBrems: for proton SubType= 3
129 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
130 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
131 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]132 hBrem : Emin= 0 eV Emax= 10 TeV
[807]133
[1337]134hPairProd: for proton SubType= 4
135 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
136 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
137 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]138 hPairProd : Emin= 0 eV Emax= 10 TeV
[807]139
[1337]140msc: for GenericIon SubType= 10
141 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
142 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]143 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV
[807]144
[1337]145ionIoni: for GenericIon SubType= 2
146 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
147 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
148 finalRange(mm)= 0.01, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
[1342]149 Stopping Power data for 17 ion/material pairs
[1337]150 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]151 ParamICRU73 : Emin= 0 eV Emax= 10 TeV
[807]152
[1337]153nuclearStopping: for GenericIon SubType= 1
154 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]155 ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV
[807]156
[1337]157nuclearStopping: for alpha SubType= 1
158 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]159 ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV
[807]160
[1337]161hIoni: for anti_proton SubType= 2
162 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
163 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
164 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
165 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]166 ICRU73QO : Emin= 0 eV Emax= 2 MeV
167 BetheBloch : Emin= 2 MeV Emax= 10 TeV
[807]168
[1342]169msc: for kaon+ SubType= 10
170 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
171 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
172 ===== EM models for the G4Region DefaultRegionForTheWorld ======
173 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV
174
[1337]175hIoni: for kaon+ SubType= 2
176 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
177 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
178 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
179 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]180 Bragg : Emin= 0 eV Emax= 1.05231 MeV
181 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV
[807]182
[1337]183hBrems: for kaon+ SubType= 3
184 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
185 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
186 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]187 hBrem : Emin= 0 eV Emax= 10 TeV
[807]188
[1337]189hPairProd: for kaon+ SubType= 4
190 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
191 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
192 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]193 hPairProd : Emin= 0 eV Emax= 10 TeV
[807]194
[1337]195hIoni: for kaon- SubType= 2
196 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
197 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
198 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
199 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]200 ICRU73QO : Emin= 0 eV Emax= 1.05231 MeV
201 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV
[807]202
[1337]203hBrems: for kaon- SubType= 3
204 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
205 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
206 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]207 hBrem : Emin= 0 eV Emax= 10 TeV
[807]208
[1337]209hPairProd: for kaon- SubType= 4
210 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
211 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
212 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]213 hPairProd : Emin= 0 eV Emax= 10 TeV
[807]214
[1337]215muMsc: for mu+ SubType= 10
216 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
217 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.4592
218 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]219 WentzelVIUni : Emin= 0 eV Emax= 10 TeV
[1337]220
221muIoni: for mu+ SubType= 2
222 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
223 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
224 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
225 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]226 Bragg : Emin= 0 eV Emax= 200 keV
227 BetheBloch : Emin= 200 keV Emax= 1 GeV
228 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV
[1337]229
230muBrems: for mu+ SubType= 3
231 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
232 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
233 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]234 MuBrem : Emin= 0 eV Emax= 10 TeV
[1337]235
236muPairProd: for mu+ SubType= 4
237 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
238 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
239 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]240 muPairProd : Emin= 0 eV Emax= 10 TeV
[1337]241
242CoulombScat: for mu+ SubType= 1
243 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
244 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531
245 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]246 eCoulombScattering : Emin= 0 eV Emax= 10 TeV
[1337]247
248muIoni: for mu- SubType= 2
249 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
250 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
251 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
252 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]253 ICRU73QO : Emin= 0 eV Emax= 200 keV
254 BetheBloch : Emin= 200 keV Emax= 1 GeV
255 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV
[1337]256
257muBrems: for mu- SubType= 3
258 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
259 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
260 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]261 MuBrem : Emin= 0 eV Emax= 10 TeV
[1337]262
263muPairProd: for mu- SubType= 4
264 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
265 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
266 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]267 muPairProd : Emin= 0 eV Emax= 10 TeV
[1337]268
269CoulombScat: for mu- SubType= 1
270 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
271 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531
272 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]273 eCoulombScattering : Emin= 0 eV Emax= 10 TeV
[1337]274
275hIoni: for pi+ SubType= 2
276 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
277 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
278 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
279 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]280 Bragg : Emin= 0 eV Emax= 297.505 keV
281 BetheBloch : Emin= 297.505 keV Emax= 10 TeV
[1337]282
283hBrems: for pi+ SubType= 3
284 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
285 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
286 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]287 hBrem : Emin= 0 eV Emax= 10 TeV
[1337]288
289hPairProd: for pi+ SubType= 4
290 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
291 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
292 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]293 hPairProd : Emin= 0 eV Emax= 10 TeV
[1337]294
295msc: for pi- SubType= 10
296 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
297 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
298 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]299 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV
[1337]300
301hIoni: for pi- SubType= 2
302 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
303 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
304 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
305 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]306 ICRU73QO : Emin= 0 eV Emax= 297.505 keV
307 BetheBloch : Emin= 297.505 keV Emax= 10 TeV
[1337]308
309hBrems: for pi- SubType= 3
310 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
311 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
312 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]313 hBrem : Emin= 0 eV Emax= 10 TeV
[1337]314
315hPairProd: for pi- SubType= 4
316 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
317 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
318 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]319 hPairProd : Emin= 0 eV Emax= 10 TeV
[1337]320
[807]321Region <DefaultRegionForTheWorld> -- appears in <World> world volume
322 Root logical volume(s) : World
323 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
324 Materials : Galactic Silicon
[1337]325 Production cuts : gamma 7 um e- 7 um e+ 7 um proton 0 fm
[807]326
327========= Table of registered couples ==============================
328
329Index : 0 used in the geometry : Yes recalculation needed : No
330 Material : Galactic
[1337]331 Range cuts : gamma 7 um e- 7 um e+ 7 um proton 0 fm
332 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 0 eV
[807]333 Region(s) which use this couple :
334 DefaultRegionForTheWorld
335
336Index : 1 used in the geometry : Yes recalculation needed : No
337 Material : Silicon
[1337]338 Range cuts : gamma 7 um e- 7 um e+ 7 um proton 0 fm
339 Energy thresholds : gamma 990 eV e- 21.4886 keV e+ 20.8954 keV proton 0 eV
[807]340 Region(s) which use this couple :
341 DefaultRegionForTheWorld
342
343====================================================================
344
345Start closing geometry.
346G4GeometryManager::ReportVoxelStats -- Voxel Statistics
347
348 Total memory consumed for geometry optimisation: 0 kByte
349 Total CPU time elapsed for geometry optimisation: 0 seconds
350### Run 0 start.
351
352--------- Ranecu engine status ---------
353 Initial seed (index) = 0
354 Current couple of seeds = 9876, 54321
355----------------------------------------
356Start Run processing.
357
358---> Begin of Event: 0
359
360---> Begin of Event: 10000
361
362---> Begin of Event: 20000
363
364---> Begin of Event: 30000
365
366---> Begin of Event: 40000
367
368---> Begin of Event: 50000
369
370---> Begin of Event: 60000
371
372---> Begin of Event: 70000
373
374---> Begin of Event: 80000
375
376---> Begin of Event: 90000
377Run terminated.
378Run Summary
379 Number of events processed : 100000
[1342]380 User=4.2s Real=6.41s Sys=2.21s
[807]381
382 ======================== run summary ======================
383
384 The run was 100000 pi+ of 5 GeV through 20 um of Silicon (density: 2.33 g/cm3 )
385
[1342]386 Total energy deposit in absorber per event = 6.355 keV +- 16.41 eV
[807]387
[1342]388 -----> Mean dE/dx = 3.177 MeV/cm (1.364 MeV*cm2/g)
[807]389
390 From formulas :
[1337]391 restricted dEdx = 2.963 MeV/cm (1.272 MeV*cm2/g)
392 full dEdx = 4.726 MeV/cm (2.028 MeV*cm2/g)
[807]393
[1342]394 Leakage : primary = 5 GeV +- 1.274 keV secondaries = 3.64 keV +- 1.273 keV
[807]395 Energy balance : edep + eleak = 5 GeV
396
[1342]397 Total track length (charged) in absorber per event = 20.26 um +- 7.821 nm
398 Total track length (neutral) in absorber per event = 5.528 Ang +- 3.51 Ang
[807]399
[1342]400 Number of steps (charged) in absorber per event = 1.159 +- 8.843e-06
401 Number of steps (neutral) in absorber per event = 3e-05 +- 8.843e-06
[807]402
[1342]403 Number of secondaries per event : Gammas = 3e-05; electrons = 0.01661; positrons = 0
[807]404
405 Number of events with the primary particle transmitted = 100 %
406 Number of events with at least 1 particle transmitted (same charge as primary) = 100 %
407
408 Number of events with the primary particle reflected = 0 %
409 Number of events with at least 1 particle reflected (same charge as primary) = 0 %
410
411 MultipleScattering:
[1342]412 rms proj angle of transmit primary particle = 0.02446 mrad (central part only)
[807]413 computed theta0 (Highland formula) = 0.02627 mrad
[1342]414 central part defined as +- 0.0788 mrad; Tail ratio = 1.553 %
[807]415
416--------- Ranecu engine status ---------
417 Initial seed (index) = 0
[1342]418 Current couple of seeds = 1096967505, 715451002
[807]419----------------------------------------
[1337]420#
421/testem/det/setAbsMat Gold
422/testem/det/setAbsThick 9.658 um
423/testem/det/setAbsYZ 1 mm
424/testem/det/update
425
[1342]426 Material: Galactic density: 0.000 kg/m3 RadL: 204727512.315 pc Nucl.Int.Length: 114561548.020 pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm
[1337]427 ---> Element: Galactic ( ) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
428
[1342]429 Material: Gold density: 19.320 g/cm3 RadL: 3.344 mm Nucl.Int.Length: 10.539 cm Imean: 790.000 eV
[1337]430 ---> Element: Gold ( ) Z = 79.0 N = 197.0 A = 196.97 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
431
432 The WORLD is made of 14.487 um of Galactic. The transverse size (YZ) of the world is 1.2 mm
433 The ABSORBER is made of 9.658 um of Gold. The transverse size (YZ) is 1 mm
434 X position of the middle of the absorber 0 fm
435Obsolete world logical volume is removed from the default region.
436World is registered to the default region.
437/testem/gun/setDefault
438/gun/particle e-
439/gun/energy 15.7 MeV
440/run/beamOn 50000
441
442phot: for gamma SubType= 12
443 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]444 PhotoElectric : Emin= 0 eV Emax= 10 TeV
[1337]445
446compt: for gamma SubType= 13
447 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
448 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]449 Klein-Nishina : Emin= 0 eV Emax= 10 TeV
[1337]450
451conv: for gamma SubType= 14
[1342]452 Lambda tables from 1.022 MeV to 10 TeV in 220 bins, spline: 1
[1337]453 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]454 BetheHeitler : Emin= 0 eV Emax= 10 TeV
[1337]455
456Rayl: for gamma SubType= 11
457 Lambda tables from 100 eV to 10 TeV in 200 bins, spline: 1
458 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]459 LivermoreRayleigh : Emin= 0 eV Emax= 100 GeV
[1337]460
461msc: for e- SubType= 10
462 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
463 RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5
464 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]465 UrbanMsc93 : Emin= 0 eV Emax= 10 TeV
[1337]466
467eIoni: for e- SubType= 2
468 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
469 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
470 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
471 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]472 MollerBhabha : Emin= 0 eV Emax= 10 TeV
[1337]473
474eBrem: for e- SubType= 3
475 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
476 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
477 LPM flag: 1 for E > 1 GeV
478 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]479 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban
480 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban
[1337]481
482eIoni: for e+ SubType= 2
483 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
484 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
485 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
486 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]487 MollerBhabha : Emin= 0 eV Emax= 10 TeV
[1337]488
489eBrem: for e+ SubType= 3
490 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
491 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
492 LPM flag: 1 for E > 1 GeV
493 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]494 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban
495 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban
[1337]496
497annihil: for e+ SubType= 5
498 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
499 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]500 eplus2gg : Emin= 0 eV Emax= 10 TeV
[1337]501
502msc: for proton SubType= 10
503 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
504 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
505 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]506 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV
[1337]507
508hIoni: for proton SubType= 2
509 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
510 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
511 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
512 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]513 Bragg : Emin= 0 eV Emax= 2 MeV
514 BetheBloch : Emin= 2 MeV Emax= 10 TeV
[1337]515
516hBrems: for proton SubType= 3
517 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
518 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
519 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]520 hBrem : Emin= 0 eV Emax= 10 TeV
[1337]521
522hPairProd: for proton SubType= 4
523 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
524 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
525 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]526 hPairProd : Emin= 0 eV Emax= 10 TeV
[1337]527
528msc: for GenericIon SubType= 10
529 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
530 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]531 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV
[1337]532
533ionIoni: for GenericIon SubType= 2
534 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
535 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
536 finalRange(mm)= 0.01, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
[1342]537 Stopping Power data for 17 ion/material pairs
[1337]538 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]539 ParamICRU73 : Emin= 0 eV Emax= 10 TeV
[1337]540
541nuclearStopping: for GenericIon SubType= 1
542 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]543 ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV
[1337]544
545nuclearStopping: for alpha SubType= 1
546 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]547 ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV
[1337]548
549hIoni: for anti_proton SubType= 2
550 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
551 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
552 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
553 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]554 ICRU73QO : Emin= 0 eV Emax= 2 MeV
555 BetheBloch : Emin= 2 MeV Emax= 10 TeV
[1337]556
[1342]557msc: for kaon+ SubType= 10
558 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
559 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
560 ===== EM models for the G4Region DefaultRegionForTheWorld ======
561 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV
562
[1337]563hIoni: for kaon+ SubType= 2
564 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
565 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
566 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
567 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]568 Bragg : Emin= 0 eV Emax= 1.05231 MeV
569 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV
[1337]570
571hBrems: for kaon+ SubType= 3
572 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
573 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
574 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]575 hBrem : Emin= 0 eV Emax= 10 TeV
[1337]576
577hPairProd: for kaon+ SubType= 4
578 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
579 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
580 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]581 hPairProd : Emin= 0 eV Emax= 10 TeV
[1337]582
583hIoni: for kaon- SubType= 2
584 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
585 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
586 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
587 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]588 ICRU73QO : Emin= 0 eV Emax= 1.05231 MeV
589 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV
[1337]590
591hBrems: for kaon- SubType= 3
592 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
593 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
594 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]595 hBrem : Emin= 0 eV Emax= 10 TeV
[1337]596
597hPairProd: for kaon- SubType= 4
598 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
599 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
600 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]601 hPairProd : Emin= 0 eV Emax= 10 TeV
[1337]602
603muMsc: for mu+ SubType= 10
604 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
[1342]605 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.4592
[1337]606 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]607 WentzelVIUni : Emin= 0 eV Emax= 10 TeV
[1337]608
609muIoni: for mu+ SubType= 2
610 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
611 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
612 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
613 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]614 Bragg : Emin= 0 eV Emax= 200 keV
615 BetheBloch : Emin= 200 keV Emax= 1 GeV
616 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV
[1337]617
618muBrems: for mu+ SubType= 3
619 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
620 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
621 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]622 MuBrem : Emin= 0 eV Emax= 10 TeV
[1337]623
624muPairProd: for mu+ SubType= 4
625 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
626 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
627 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]628 muPairProd : Emin= 0 eV Emax= 10 TeV
[1337]629
630CoulombScat: for mu+ SubType= 1
631 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
[1342]632 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531
[1337]633 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]634 eCoulombScattering : Emin= 0 eV Emax= 10 TeV
[1337]635
636muIoni: for mu- SubType= 2
637 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
638 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
639 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
640 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]641 ICRU73QO : Emin= 0 eV Emax= 200 keV
642 BetheBloch : Emin= 200 keV Emax= 1 GeV
643 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV
[1337]644
645muBrems: for mu- SubType= 3
646 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
647 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
648 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]649 MuBrem : Emin= 0 eV Emax= 10 TeV
[1337]650
651muPairProd: for mu- SubType= 4
652 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
653 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
654 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]655 muPairProd : Emin= 0 eV Emax= 10 TeV
[1337]656
657CoulombScat: for mu- SubType= 1
658 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
[1342]659 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531
[1337]660 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]661 eCoulombScattering : Emin= 0 eV Emax= 10 TeV
[1337]662
663hIoni: for pi+ SubType= 2
664 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
665 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
666 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
667 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]668 Bragg : Emin= 0 eV Emax= 297.505 keV
669 BetheBloch : Emin= 297.505 keV Emax= 10 TeV
[1337]670
671hBrems: for pi+ SubType= 3
672 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
673 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
674 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]675 hBrem : Emin= 0 eV Emax= 10 TeV
[1337]676
677hPairProd: for pi+ SubType= 4
678 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
679 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
680 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]681 hPairProd : Emin= 0 eV Emax= 10 TeV
[1337]682
683msc: for pi- SubType= 10
684 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
685 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
686 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]687 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV
[1337]688
689hIoni: for pi- SubType= 2
690 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
691 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
692 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
693 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]694 ICRU73QO : Emin= 0 eV Emax= 297.505 keV
695 BetheBloch : Emin= 297.505 keV Emax= 10 TeV
[1337]696
697hBrems: for pi- SubType= 3
698 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
699 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
700 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]701 hBrem : Emin= 0 eV Emax= 10 TeV
[1337]702
703hPairProd: for pi- SubType= 4
704 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
705 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
706 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]707 hPairProd : Emin= 0 eV Emax= 10 TeV
[1337]708
709Region <DefaultRegionForTheWorld> -- appears in <World> world volume
710 Root logical volume(s) : World
711 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
712 Materials : Galactic Gold
713 Production cuts : gamma 7 um e- 7 um e+ 7 um proton 0 fm
714
715========= Table of registered couples ==============================
716
717Index : 0 used in the geometry : Yes recalculation needed : No
718 Material : Galactic
719 Range cuts : gamma 7 um e- 7 um e+ 7 um proton 0 fm
720 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 0 eV
721 Region(s) which use this couple :
722 DefaultRegionForTheWorld
723
724Index : 1 used in the geometry : No recalculation needed : No
725 Material : Silicon
726 Range cuts : gamma 7 um e- 7 um e+ 7 um proton 0 fm
[1342]727 Energy thresholds : gamma 990 eV e- 21.4886 keV e+ 20.8954 keV proton 0 eV
[1337]728
729Index : 2 used in the geometry : Yes recalculation needed : No
730 Material : Gold
731 Range cuts : gamma 7 um e- 7 um e+ 7 um proton 0 fm
[1342]732 Energy thresholds : gamma 6.54229 keV e- 65.1406 keV e+ 63.8405 keV proton 0 eV
[1337]733 Region(s) which use this couple :
734 DefaultRegionForTheWorld
735
736====================================================================
737
738Start closing geometry.
739G4GeometryManager::ReportVoxelStats -- Voxel Statistics
740
741 Total memory consumed for geometry optimisation: 0 kByte
742 Total CPU time elapsed for geometry optimisation: 0 seconds
743### Run 1 start.
744
745--------- Ranecu engine status ---------
746 Initial seed (index) = 0
[1342]747 Current couple of seeds = 1096967505, 715451002
[1337]748----------------------------------------
749Start Run processing.
750
751---> Begin of Event: 0
752
753---> Begin of Event: 10000
754
755---> Begin of Event: 20000
756
757---> Begin of Event: 30000
758
759---> Begin of Event: 40000
760Run terminated.
761Run Summary
762 Number of events processed : 50000
[1342]763 User=3.35s Real=4.56s Sys=1.21s
[1337]764
765 ======================== run summary ======================
766
767 The run was 50000 e- of 15.7 MeV through 9.66 um of Gold (density: 19.3 g/cm3 )
768
[1342]769 Total energy deposit in absorber per event = 18.54 keV +- 65.36 eV
[1337]770
[1342]771 -----> Mean dE/dx = 19.2 MeV/cm (0.9938 MeV*cm2/g)
[1337]772
773 From formulas :
774 restricted dEdx = 17.9 MeV/cm (0.9266 MeV*cm2/g)
775 full dEdx = 59.7 MeV/cm (3.09 MeV*cm2/g)
776
[1342]777 Leakage : primary = 15.63 MeV +- 2.365 keV secondaries = 47.06 keV +- 2.361 keV
[1337]778 Energy balance : edep + eleak = 15.7 MeV
779
[1342]780 Total track length (charged) in absorber per event = 9.938 um +- 12.01 nm
781 Total track length (neutral) in absorber per event = 131.3 nm +- 4.242 nm
[1337]782
[1342]783 Number of steps (charged) in absorber per event = 5.069 +- 1.55e-05
784 Number of steps (neutral) in absorber per event = 0.02996 +- 1.55e-05
[1337]785
[1342]786 Number of secondaries per event : Gammas = 0.02986; electrons = 0.023; positrons = 0
[1337]787
[1342]788 Number of events with the primary particle transmitted = 99.98 %
789 Number of events with at least 1 particle transmitted (same charge as primary) = 99.98 %
[1337]790
[1342]791 Number of events with the primary particle reflected = 0.016 %
792 Number of events with at least 1 particle reflected (same charge as primary) = 0.294 %
[1337]793
794 MultipleScattering:
[1342]795 rms proj angle of transmit primary particle = 32.82 mrad (central part only)
[1337]796 computed theta0 (Highland formula) = 35.1 mrad
[1342]797 central part defined as +- 105.3 mrad; Tail ratio = 2.774 %
[1337]798
799--------- Ranecu engine status ---------
800 Initial seed (index) = 0
[1342]801 Current couple of seeds = 1228101883, 1478657312
[1337]802----------------------------------------
803#
804/testem/phys/setCuts 10 um
805#
806/testem/det/setAbsMat G4_POLYETHYLENE
807/testem/det/setAbsThick 4.74 mm
808/testem/det/setAbsYZ 50 mm
809/testem/det/update
810
[1342]811 Material: Galactic density: 0.000 kg/m3 RadL: 204727512.315 pc Nucl.Int.Length: 114561548.020 pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm
[1337]812 ---> Element: Galactic ( ) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
813
[1342]814 Material: G4_POLYETHYLENE (C_2H_4)_N-Polyethylene density: 940.000 mg/cm3 RadL: 47.632 cm Nucl.Int.Length: 72.122 cm Imean: 57.400 eV
[1337]815 ---> Element: H (H) Z = 1.0 N = 1.0 A = 1.01 g/mole
816 ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.99 %
817 ---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.01 % ElmMassFraction: 14.37 % ElmAbundance 66.67 %
818 ---> Element: C (C) Z = 6.0 N = 12.0 A = 12.01 g/mole
819 ---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.93 %
820 ---> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.07 % ElmMassFraction: 85.63 % ElmAbundance 33.33 %
821
822 The WORLD is made of 7.11 mm of Galactic. The transverse size (YZ) of the world is 6 cm
823 The ABSORBER is made of 4.74 mm of G4_POLYETHYLENE. The transverse size (YZ) is 5 cm
824 X position of the middle of the absorber 0 fm
825Obsolete world logical volume is removed from the default region.
826World is registered to the default region.
827/testem/gun/setDefault
828/gun/particle mu+
829/gun/energy 96.2 MeV
830/run/beamOn 50000
831
832phot: for gamma SubType= 12
833 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]834 PhotoElectric : Emin= 0 eV Emax= 10 TeV
[1337]835
836compt: for gamma SubType= 13
837 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
838 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]839 Klein-Nishina : Emin= 0 eV Emax= 10 TeV
[1337]840
841conv: for gamma SubType= 14
[1342]842 Lambda tables from 1.022 MeV to 10 TeV in 220 bins, spline: 1
[1337]843 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]844 BetheHeitler : Emin= 0 eV Emax= 10 TeV
[1337]845
846Rayl: for gamma SubType= 11
847 Lambda tables from 100 eV to 10 TeV in 200 bins, spline: 1
848 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]849 LivermoreRayleigh : Emin= 0 eV Emax= 100 GeV
[1337]850
851msc: for e- SubType= 10
852 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
853 RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5
854 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]855 UrbanMsc93 : Emin= 0 eV Emax= 10 TeV
[1337]856
857eIoni: for e- SubType= 2
858 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
859 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
860 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
861 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]862 MollerBhabha : Emin= 0 eV Emax= 10 TeV
[1337]863
864eBrem: for e- SubType= 3
865 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
866 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
867 LPM flag: 1 for E > 1 GeV
868 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]869 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban
870 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban
[1337]871
872eIoni: for e+ SubType= 2
873 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
874 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
875 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
876 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]877 MollerBhabha : Emin= 0 eV Emax= 10 TeV
[1337]878
879eBrem: for e+ SubType= 3
880 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
881 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
882 LPM flag: 1 for E > 1 GeV
883 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]884 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban
885 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban
[1337]886
887annihil: for e+ SubType= 5
888 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
889 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]890 eplus2gg : Emin= 0 eV Emax= 10 TeV
[1337]891
892msc: for proton SubType= 10
893 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
894 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
895 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]896 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV
[1337]897
898hIoni: for proton SubType= 2
899 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
900 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
901 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
902 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]903 Bragg : Emin= 0 eV Emax= 2 MeV
904 BetheBloch : Emin= 2 MeV Emax= 10 TeV
[1337]905
906hBrems: for proton SubType= 3
907 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
908 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
909 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]910 hBrem : Emin= 0 eV Emax= 10 TeV
[1337]911
912hPairProd: for proton SubType= 4
913 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
914 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
915 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]916 hPairProd : Emin= 0 eV Emax= 10 TeV
[1337]917
918msc: for GenericIon SubType= 10
919 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
920 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]921 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV
[1337]922
923ionIoni: for GenericIon SubType= 2
924 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
925 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
926 finalRange(mm)= 0.01, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
[1342]927 Stopping Power data for 17 ion/material pairs
[1337]928 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]929 ParamICRU73 : Emin= 0 eV Emax= 10 TeV
[1337]930
931nuclearStopping: for GenericIon SubType= 1
932 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]933 ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV
[1337]934
935nuclearStopping: for alpha SubType= 1
936 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]937 ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV
[1337]938
939hIoni: for anti_proton SubType= 2
940 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
941 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
942 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
943 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]944 ICRU73QO : Emin= 0 eV Emax= 2 MeV
945 BetheBloch : Emin= 2 MeV Emax= 10 TeV
[1337]946
[1342]947msc: for kaon+ SubType= 10
948 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
949 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
950 ===== EM models for the G4Region DefaultRegionForTheWorld ======
951 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV
952
[1337]953hIoni: for kaon+ SubType= 2
954 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
955 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
956 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
957 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]958 Bragg : Emin= 0 eV Emax= 1.05231 MeV
959 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV
[1337]960
961hBrems: for kaon+ SubType= 3
962 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
963 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
964 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]965 hBrem : Emin= 0 eV Emax= 10 TeV
[1337]966
967hPairProd: for kaon+ SubType= 4
968 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
969 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
970 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]971 hPairProd : Emin= 0 eV Emax= 10 TeV
[1337]972
973hIoni: for kaon- SubType= 2
974 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
975 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
976 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
977 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]978 ICRU73QO : Emin= 0 eV Emax= 1.05231 MeV
979 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV
[1337]980
981hBrems: for kaon- SubType= 3
982 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
983 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
984 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]985 hBrem : Emin= 0 eV Emax= 10 TeV
[1337]986
987hPairProd: for kaon- SubType= 4
988 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
989 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
990 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]991 hPairProd : Emin= 0 eV Emax= 10 TeV
[1337]992
993muMsc: for mu+ SubType= 10
994 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
[1342]995 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.4592
[1337]996 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]997 WentzelVIUni : Emin= 0 eV Emax= 10 TeV
[1337]998
999muIoni: for mu+ SubType= 2
1000 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1001 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1002 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1003 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1004 Bragg : Emin= 0 eV Emax= 200 keV
1005 BetheBloch : Emin= 200 keV Emax= 1 GeV
1006 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV
[1337]1007
1008muBrems: for mu+ SubType= 3
1009 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1010 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1011 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1012 MuBrem : Emin= 0 eV Emax= 10 TeV
[1337]1013
1014muPairProd: for mu+ SubType= 4
1015 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1016 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1017 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1018 muPairProd : Emin= 0 eV Emax= 10 TeV
[1337]1019
1020CoulombScat: for mu+ SubType= 1
1021 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
[1342]1022 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531
[1337]1023 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1024 eCoulombScattering : Emin= 0 eV Emax= 10 TeV
[1337]1025
1026muIoni: for mu- SubType= 2
1027 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1028 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1029 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1030 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1031 ICRU73QO : Emin= 0 eV Emax= 200 keV
1032 BetheBloch : Emin= 200 keV Emax= 1 GeV
1033 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV
[1337]1034
1035muBrems: for mu- SubType= 3
1036 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1037 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1038 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1039 MuBrem : Emin= 0 eV Emax= 10 TeV
[1337]1040
1041muPairProd: for mu- SubType= 4
1042 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1043 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1044 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1045 muPairProd : Emin= 0 eV Emax= 10 TeV
[1337]1046
1047CoulombScat: for mu- SubType= 1
1048 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
[1342]1049 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531
[1337]1050 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1051 eCoulombScattering : Emin= 0 eV Emax= 10 TeV
[1337]1052
1053hIoni: for pi+ SubType= 2
1054 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1055 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1056 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1057 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1058 Bragg : Emin= 0 eV Emax= 297.505 keV
1059 BetheBloch : Emin= 297.505 keV Emax= 10 TeV
[1337]1060
1061hBrems: for pi+ SubType= 3
1062 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1063 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1064 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1065 hBrem : Emin= 0 eV Emax= 10 TeV
[1337]1066
1067hPairProd: for pi+ SubType= 4
1068 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1069 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1070 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1071 hPairProd : Emin= 0 eV Emax= 10 TeV
[1337]1072
1073msc: for pi- SubType= 10
1074 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
1075 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
1076 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1077 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV
[1337]1078
1079hIoni: for pi- SubType= 2
1080 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1081 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1082 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1083 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1084 ICRU73QO : Emin= 0 eV Emax= 297.505 keV
1085 BetheBloch : Emin= 297.505 keV Emax= 10 TeV
[1337]1086
1087hBrems: for pi- SubType= 3
1088 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1089 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1090 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1091 hBrem : Emin= 0 eV Emax= 10 TeV
[1337]1092
1093hPairProd: for pi- SubType= 4
1094 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1095 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1096 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1097 hPairProd : Emin= 0 eV Emax= 10 TeV
[1337]1098
1099Region <DefaultRegionForTheWorld> -- appears in <World> world volume
1100 Root logical volume(s) : World
1101 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
1102 Materials : Galactic G4_POLYETHYLENE
1103 Production cuts : gamma 10 um e- 10 um e+ 10 um proton 0 fm
1104
1105========= Table of registered couples ==============================
1106
1107Index : 0 used in the geometry : Yes recalculation needed : No
1108 Material : Galactic
1109 Range cuts : gamma 10 um e- 10 um e+ 10 um proton 0 fm
1110 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 0 eV
1111 Region(s) which use this couple :
1112 DefaultRegionForTheWorld
1113
1114Index : 1 used in the geometry : No recalculation needed : No
1115 Material : Silicon
1116 Range cuts : gamma 10 um e- 10 um e+ 10 um proton 0 fm
[1342]1117 Energy thresholds : gamma 990 eV e- 31.9516 keV e+ 31.7376 keV proton 0 eV
[1337]1118
1119Index : 2 used in the geometry : No recalculation needed : No
1120 Material : Gold
1121 Range cuts : gamma 10 um e- 10 um e+ 10 um proton 0 fm
[1342]1122 Energy thresholds : gamma 8.27707 keV e- 80.7687 keV e+ 79.1567 keV proton 0 eV
[1337]1123
1124Index : 3 used in the geometry : Yes recalculation needed : No
1125 Material : G4_POLYETHYLENE
1126 Range cuts : gamma 10 um e- 10 um e+ 10 um proton 0 fm
[1342]1127 Energy thresholds : gamma 990 eV e- 14.1748 keV e+ 13.7548 keV proton 0 eV
[1337]1128 Region(s) which use this couple :
1129 DefaultRegionForTheWorld
1130
1131====================================================================
1132
1133Start closing geometry.
1134G4GeometryManager::ReportVoxelStats -- Voxel Statistics
1135
1136 Total memory consumed for geometry optimisation: 0 kByte
1137 Total CPU time elapsed for geometry optimisation: 0 seconds
1138### Run 2 start.
1139
1140--------- Ranecu engine status ---------
1141 Initial seed (index) = 0
[1342]1142 Current couple of seeds = 1228101883, 1478657312
[1337]1143----------------------------------------
1144Start Run processing.
1145
1146---> Begin of Event: 0
1147
1148---> Begin of Event: 10000
1149
1150---> Begin of Event: 20000
1151
1152---> Begin of Event: 30000
1153
1154---> Begin of Event: 40000
1155Run terminated.
1156Run Summary
1157 Number of events processed : 50000
[1342]1158 User=6.36s Real=7.53s Sys=1.16s
[1337]1159
1160 ======================== run summary ======================
1161
1162 The run was 50000 mu+ of 96.2 MeV through 4.74 mm of G4_POLYETHYLENE (density: 940 mg/cm3)
1163
[1342]1164 Total energy deposit in absorber per event = 1.031 MeV +- 1.011 keV
[1337]1165
[1342]1166 -----> Mean dE/dx = 2.175 MeV/cm (2.314 MeV*cm2/g)
[1337]1167
1168 From formulas :
[1342]1169 restricted dEdx = 1.747 MeV/cm (1.859 MeV*cm2/g)
1170 full dEdx = 2.259 MeV/cm (2.403 MeV*cm2/g)
[1337]1171
[1342]1172 Leakage : primary = 95.14 MeV +- 1.348 keV secondaries = 26.07 keV +- 707.1 eV
[1337]1173 Energy balance : edep + eleak = 96.2 MeV
1174
[1342]1175 Total track length (charged) in absorber per event = 5.102 mm +- 3.895 um
1176 Total track length (neutral) in absorber per event = 26.38 um +- 2.481 um
[1337]1177
[1342]1178 Number of steps (charged) in absorber per event = 11.27 +- 0.0002791
1179 Number of steps (neutral) in absorber per event = 0.01202 +- 0.0002791
[1337]1180
[1342]1181 Number of secondaries per event : Gammas = 0.01128; electrons = 3.989; positrons = 0
[1337]1182
1183 Number of events with the primary particle transmitted = 100 %
1184 Number of events with at least 1 particle transmitted (same charge as primary) = 100 %
1185
1186 Number of events with the primary particle reflected = 0 %
1187 Number of events with at least 1 particle reflected (same charge as primary) = 0 %
1188
1189 MultipleScattering:
[1342]1190 rms proj angle of transmit primary particle = 7.275 mrad (central part only)
[1337]1191 computed theta0 (Highland formula) = 7.636 mrad
[1342]1192 central part defined as +- 22.91 mrad; Tail ratio = 1.259 %
[1337]1193
1194--------- Ranecu engine status ---------
1195 Initial seed (index) = 0
[1342]1196 Current couple of seeds = 1091561076, 776987475
[1337]1197----------------------------------------
1198#
1199/testem/phys/setCuts 1 mm
1200#
1201/testem/stack/killSecondaries
1202#
1203/testem/det/setAbsMat Aluminium
1204/testem/det/setAbsThick 8.004 mm
1205/testem/det/update
1206
[1342]1207 Material: Galactic density: 0.000 kg/m3 RadL: 204727512.315 pc Nucl.Int.Length: 114561548.020 pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm
[1337]1208 ---> Element: Galactic ( ) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
1209
[1342]1210 Material: Aluminium density: 2.700 g/cm3 RadL: 8.893 cm Nucl.Int.Length: 38.860 cm Imean: 166.000 eV
[1337]1211 ---> Element: Aluminium ( ) Z = 13.0 N = 27.0 A = 26.98 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
1212
1213 The WORLD is made of 1.2006 cm of Galactic. The transverse size (YZ) of the world is 6 cm
1214 The ABSORBER is made of 8.004 mm of Aluminium. The transverse size (YZ) is 5 cm
1215 X position of the middle of the absorber 0 fm
1216Obsolete world logical volume is removed from the default region.
1217World is registered to the default region.
1218/testem/gun/setDefault
1219/gun/particle proton
1220/gun/energy 174.664 GeV
1221/run/beamOn 50000
1222
1223phot: for gamma SubType= 12
1224 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1225 PhotoElectric : Emin= 0 eV Emax= 10 TeV
[1337]1226
1227compt: for gamma SubType= 13
1228 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
1229 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1230 Klein-Nishina : Emin= 0 eV Emax= 10 TeV
[1337]1231
1232conv: for gamma SubType= 14
[1342]1233 Lambda tables from 1.022 MeV to 10 TeV in 220 bins, spline: 1
[1337]1234 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1235 BetheHeitler : Emin= 0 eV Emax= 10 TeV
[1337]1236
1237Rayl: for gamma SubType= 11
1238 Lambda tables from 100 eV to 10 TeV in 200 bins, spline: 1
1239 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1240 LivermoreRayleigh : Emin= 0 eV Emax= 100 GeV
[1337]1241
1242msc: for e- SubType= 10
1243 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
1244 RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5
1245 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1246 UrbanMsc93 : Emin= 0 eV Emax= 10 TeV
[1337]1247
1248eIoni: for e- SubType= 2
1249 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1250 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1251 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1252 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1253 MollerBhabha : Emin= 0 eV Emax= 10 TeV
[1337]1254
1255eBrem: for e- SubType= 3
1256 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1257 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1258 LPM flag: 1 for E > 1 GeV
1259 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1260 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban
1261 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban
[1337]1262
1263eIoni: for e+ SubType= 2
1264 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1265 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1266 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1267 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1268 MollerBhabha : Emin= 0 eV Emax= 10 TeV
[1337]1269
1270eBrem: for e+ SubType= 3
1271 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1272 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1273 LPM flag: 1 for E > 1 GeV
1274 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1275 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban
1276 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban
[1337]1277
1278annihil: for e+ SubType= 5
1279 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
1280 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1281 eplus2gg : Emin= 0 eV Emax= 10 TeV
[1337]1282
1283msc: for proton SubType= 10
1284 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
1285 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
1286 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1287 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV
[1337]1288
1289hIoni: for proton SubType= 2
1290 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1291 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1292 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1293 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1294 Bragg : Emin= 0 eV Emax= 2 MeV
1295 BetheBloch : Emin= 2 MeV Emax= 10 TeV
[1337]1296
1297hBrems: for proton SubType= 3
1298 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1299 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1300 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1301 hBrem : Emin= 0 eV Emax= 10 TeV
[1337]1302
1303hPairProd: for proton SubType= 4
1304 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1305 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1306 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1307 hPairProd : Emin= 0 eV Emax= 10 TeV
[1337]1308
1309msc: for GenericIon SubType= 10
1310 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
1311 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1312 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV
[1337]1313
1314ionIoni: for GenericIon SubType= 2
1315 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1316 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1317 finalRange(mm)= 0.01, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
[1342]1318 Stopping Power data for 17 ion/material pairs
[1337]1319 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1320 ParamICRU73 : Emin= 0 eV Emax= 10 TeV
[1337]1321
1322nuclearStopping: for GenericIon SubType= 1
1323 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1324 ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV
[1337]1325
1326nuclearStopping: for alpha SubType= 1
1327 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1328 ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV
[1337]1329
1330hIoni: for anti_proton SubType= 2
1331 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1332 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1333 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1334 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1335 ICRU73QO : Emin= 0 eV Emax= 2 MeV
1336 BetheBloch : Emin= 2 MeV Emax= 10 TeV
[1337]1337
[1342]1338msc: for kaon+ SubType= 10
1339 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
1340 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
1341 ===== EM models for the G4Region DefaultRegionForTheWorld ======
1342 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV
1343
[1337]1344hIoni: for kaon+ SubType= 2
1345 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1346 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1347 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1348 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1349 Bragg : Emin= 0 eV Emax= 1.05231 MeV
1350 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV
[1337]1351
1352hBrems: for kaon+ SubType= 3
1353 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1354 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1355 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1356 hBrem : Emin= 0 eV Emax= 10 TeV
[1337]1357
1358hPairProd: for kaon+ SubType= 4
1359 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1360 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1361 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1362 hPairProd : Emin= 0 eV Emax= 10 TeV
[1337]1363
1364hIoni: for kaon- SubType= 2
1365 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1366 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1367 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1368 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1369 ICRU73QO : Emin= 0 eV Emax= 1.05231 MeV
1370 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV
[1337]1371
1372hBrems: for kaon- SubType= 3
1373 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1374 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1375 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1376 hBrem : Emin= 0 eV Emax= 10 TeV
[1337]1377
1378hPairProd: for kaon- SubType= 4
1379 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1380 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1381 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1382 hPairProd : Emin= 0 eV Emax= 10 TeV
[1337]1383
1384muMsc: for mu+ SubType= 10
1385 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
[1342]1386 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.4592
[1337]1387 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1388 WentzelVIUni : Emin= 0 eV Emax= 10 TeV
[1337]1389
1390muIoni: for mu+ SubType= 2
1391 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1392 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1393 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1394 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1395 Bragg : Emin= 0 eV Emax= 200 keV
1396 BetheBloch : Emin= 200 keV Emax= 1 GeV
1397 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV
[1337]1398
1399muBrems: for mu+ SubType= 3
1400 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1401 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1402 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1403 MuBrem : Emin= 0 eV Emax= 10 TeV
[1337]1404
1405muPairProd: for mu+ SubType= 4
1406 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1407 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1408 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1409 muPairProd : Emin= 0 eV Emax= 10 TeV
[1337]1410
1411CoulombScat: for mu+ SubType= 1
1412 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
[1342]1413 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531
[1337]1414 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1415 eCoulombScattering : Emin= 0 eV Emax= 10 TeV
[1337]1416
1417muIoni: for mu- SubType= 2
1418 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1419 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1420 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1421 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1422 ICRU73QO : Emin= 0 eV Emax= 200 keV
1423 BetheBloch : Emin= 200 keV Emax= 1 GeV
1424 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV
[1337]1425
1426muBrems: for mu- SubType= 3
1427 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1428 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1429 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1430 MuBrem : Emin= 0 eV Emax= 10 TeV
[1337]1431
1432muPairProd: for mu- SubType= 4
1433 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1434 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1435 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1436 muPairProd : Emin= 0 eV Emax= 10 TeV
[1337]1437
1438CoulombScat: for mu- SubType= 1
1439 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
[1342]1440 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531
[1337]1441 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1442 eCoulombScattering : Emin= 0 eV Emax= 10 TeV
[1337]1443
1444hIoni: for pi+ SubType= 2
1445 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1446 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1447 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1448 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1449 Bragg : Emin= 0 eV Emax= 297.505 keV
1450 BetheBloch : Emin= 297.505 keV Emax= 10 TeV
[1337]1451
1452hBrems: for pi+ SubType= 3
1453 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1454 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1455 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1456 hBrem : Emin= 0 eV Emax= 10 TeV
[1337]1457
1458hPairProd: for pi+ SubType= 4
1459 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1460 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1461 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1462 hPairProd : Emin= 0 eV Emax= 10 TeV
[1337]1463
1464msc: for pi- SubType= 10
1465 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
1466 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
1467 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1468 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV
[1337]1469
1470hIoni: for pi- SubType= 2
1471 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1472 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1473 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1474 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1475 ICRU73QO : Emin= 0 eV Emax= 297.505 keV
1476 BetheBloch : Emin= 297.505 keV Emax= 10 TeV
[1337]1477
1478hBrems: for pi- SubType= 3
1479 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1480 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1481 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1482 hBrem : Emin= 0 eV Emax= 10 TeV
[1337]1483
1484hPairProd: for pi- SubType= 4
1485 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1486 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1487 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1488 hPairProd : Emin= 0 eV Emax= 10 TeV
[1337]1489
1490Region <DefaultRegionForTheWorld> -- appears in <World> world volume
1491 Root logical volume(s) : World
1492 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
1493 Materials : Galactic Aluminium
1494 Production cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm
1495
1496========= Table of registered couples ==============================
1497
1498Index : 0 used in the geometry : Yes recalculation needed : No
1499 Material : Galactic
1500 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm
1501 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 0 eV
1502 Region(s) which use this couple :
1503 DefaultRegionForTheWorld
1504
1505Index : 1 used in the geometry : No recalculation needed : No
1506 Material : Silicon
1507 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm
[1342]1508 Energy thresholds : gamma 6.95018 keV e- 548.291 keV e+ 526.624 keV proton 0 eV
[1337]1509
1510Index : 2 used in the geometry : No recalculation needed : No
1511 Material : Gold
1512 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm
[1342]1513 Energy thresholds : gamma 112.644 keV e- 2.30976 MeV e+ 2.13081 MeV proton 0 eV
[1337]1514
1515Index : 3 used in the geometry : No recalculation needed : No
1516 Material : G4_POLYETHYLENE
1517 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm
[1342]1518 Energy thresholds : gamma 2.26261 keV e- 351.877 keV e+ 342.545 keV proton 0 eV
[1337]1519
1520Index : 4 used in the geometry : Yes recalculation needed : No
1521 Material : Aluminium
1522 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm
[1342]1523 Energy thresholds : gamma 6.90363 keV e- 598.345 keV e+ 570.85 keV proton 0 eV
[1337]1524 Region(s) which use this couple :
1525 DefaultRegionForTheWorld
1526
1527====================================================================
1528
1529Start closing geometry.
1530G4GeometryManager::ReportVoxelStats -- Voxel Statistics
1531
1532 Total memory consumed for geometry optimisation: 0 kByte
1533 Total CPU time elapsed for geometry optimisation: 0 seconds
1534### Run 3 start.
1535
1536--------- Ranecu engine status ---------
1537 Initial seed (index) = 0
[1342]1538 Current couple of seeds = 1091561076, 776987475
[1337]1539----------------------------------------
1540Start Run processing.
1541
1542---> Begin of Event: 0
1543
1544---> Begin of Event: 10000
1545
1546---> Begin of Event: 20000
1547
1548---> Begin of Event: 30000
1549
1550---> Begin of Event: 40000
1551Run terminated.
1552Run Summary
1553 Number of events processed : 50000
[1342]1554 User=2.09s Real=3.2s Sys=1.12s
[1337]1555
1556 ======================== run summary ======================
1557
1558 The run was 50000 proton of 175 GeV through 8 mm of Aluminium (density: 2.7 g/cm3 )
1559
[1342]1560 Total energy deposit in absorber per event = 4.881 MeV +- 246.9 keV
[1337]1561
[1342]1562 -----> Mean dE/dx = 6.098 MeV/cm (2.259 MeV*cm2/g)
[1337]1563
1564 From formulas :
[1342]1565 restricted dEdx = 4.025 MeV/cm (1.491 MeV*cm2/g)
[1337]1566 full dEdx = 6.049 MeV/cm (2.24 MeV*cm2/g)
1567
[1342]1568 Leakage : primary = 174.7 GeV +- 246.9 keV secondaries = 0 eV +- 0 eV
[1337]1569 Energy balance : edep + eleak = 174.7 GeV
1570
1571 Total track length (charged) in absorber per event = 8.004 mm +- 0 fm
1572 Total track length (neutral) in absorber per event = 0 fm +- 0 fm
1573
[1342]1574 Number of steps (charged) in absorber per event = 1.266 +- 0
[1337]1575 Number of steps (neutral) in absorber per event = 0 +- 0
1576
[1342]1577 Number of secondaries per event : Gammas = 0; electrons = 0.266; positrons = 0.00074
[1337]1578
1579 Number of events with the primary particle transmitted = 100 %
1580 Number of events with at least 1 particle transmitted (same charge as primary) = 100 %
1581
1582 Number of events with the primary particle reflected = 0 %
1583 Number of events with at least 1 particle reflected (same charge as primary) = 0 %
1584
1585 MultipleScattering:
[1342]1586 rms proj angle of transmit primary particle = 0.02106 mrad (central part only)
[1337]1587 computed theta0 (Highland formula) = 0.02111 mrad
[1342]1588 central part defined as +- 0.06333 mrad; Tail ratio = 1.842 %
[1337]1589
1590--------- Ranecu engine status ---------
1591 Initial seed (index) = 0
[1342]1592 Current couple of seeds = 1430401471, 1601608505
[1337]1593----------------------------------------
1594#
1595/testem/det/setAbsMat Iron
1596/testem/det/setAbsThick 1 m
1597/testem/det/update
1598
[1342]1599 Material: Galactic density: 0.000 kg/m3 RadL: 204727512.315 pc Nucl.Int.Length: 114561548.020 pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm
[1337]1600 ---> Element: Galactic ( ) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
1601
[1342]1602 Material: Iron density: 7.870 g/cm3 RadL: 1.759 cm Nucl.Int.Length: 16.969 cm Imean: 286.000 eV
[1337]1603 ---> Element: Iron ( ) Z = 26.0 N = 55.8 A = 55.85 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
1604
1605 The WORLD is made of 1.5 m of Galactic. The transverse size (YZ) of the world is 6 cm
1606 The ABSORBER is made of 1 m of Iron. The transverse size (YZ) is 5 cm
1607 X position of the middle of the absorber 0 fm
1608Obsolete world logical volume is removed from the default region.
1609World is registered to the default region.
1610/testem/gun/setDefault
1611/gun/particle mu+
1612/gun/energy 100 GeV
1613/run/beamOn 50000
1614
1615phot: for gamma SubType= 12
1616 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1617 PhotoElectric : Emin= 0 eV Emax= 10 TeV
[1337]1618
1619compt: for gamma SubType= 13
1620 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
1621 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1622 Klein-Nishina : Emin= 0 eV Emax= 10 TeV
[1337]1623
1624conv: for gamma SubType= 14
[1342]1625 Lambda tables from 1.022 MeV to 10 TeV in 220 bins, spline: 1
[1337]1626 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1627 BetheHeitler : Emin= 0 eV Emax= 10 TeV
[1337]1628
1629Rayl: for gamma SubType= 11
1630 Lambda tables from 100 eV to 10 TeV in 200 bins, spline: 1
1631 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1632 LivermoreRayleigh : Emin= 0 eV Emax= 100 GeV
[1337]1633
1634msc: for e- SubType= 10
1635 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
1636 RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5
1637 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1638 UrbanMsc93 : Emin= 0 eV Emax= 10 TeV
[1337]1639
1640eIoni: for e- SubType= 2
1641 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1642 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1643 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1644 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1645 MollerBhabha : Emin= 0 eV Emax= 10 TeV
[1337]1646
1647eBrem: for e- SubType= 3
1648 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1649 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1650 LPM flag: 1 for E > 1 GeV
1651 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1652 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban
1653 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban
[1337]1654
1655eIoni: for e+ SubType= 2
1656 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1657 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1658 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1659 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1660 MollerBhabha : Emin= 0 eV Emax= 10 TeV
[1337]1661
1662eBrem: for e+ SubType= 3
1663 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1664 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1665 LPM flag: 1 for E > 1 GeV
1666 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1667 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban
1668 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban
[1337]1669
1670annihil: for e+ SubType= 5
1671 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
1672 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1673 eplus2gg : Emin= 0 eV Emax= 10 TeV
[1337]1674
1675msc: for proton SubType= 10
1676 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
1677 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
1678 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1679 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV
[1337]1680
1681hIoni: for proton SubType= 2
1682 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1683 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1684 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1685 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1686 Bragg : Emin= 0 eV Emax= 2 MeV
1687 BetheBloch : Emin= 2 MeV Emax= 10 TeV
[1337]1688
1689hBrems: for proton SubType= 3
1690 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1691 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1692 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1693 hBrem : Emin= 0 eV Emax= 10 TeV
[1337]1694
1695hPairProd: for proton SubType= 4
1696 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1697 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1698 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1699 hPairProd : Emin= 0 eV Emax= 10 TeV
[1337]1700
1701msc: for GenericIon SubType= 10
1702 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
1703 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1704 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV
[1337]1705
1706ionIoni: for GenericIon SubType= 2
1707 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1708 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1709 finalRange(mm)= 0.01, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
[1342]1710 Stopping Power data for 17 ion/material pairs
[1337]1711 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1712 ParamICRU73 : Emin= 0 eV Emax= 10 TeV
[1337]1713
1714nuclearStopping: for GenericIon SubType= 1
1715 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1716 ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV
[1337]1717
1718nuclearStopping: for alpha SubType= 1
1719 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1720 ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV
[1337]1721
1722hIoni: for anti_proton SubType= 2
1723 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1724 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1725 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1726 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1727 ICRU73QO : Emin= 0 eV Emax= 2 MeV
1728 BetheBloch : Emin= 2 MeV Emax= 10 TeV
[1337]1729
[1342]1730msc: for kaon+ SubType= 10
1731 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
1732 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
1733 ===== EM models for the G4Region DefaultRegionForTheWorld ======
1734 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV
1735
[1337]1736hIoni: for kaon+ SubType= 2
1737 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1738 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1739 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1740 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1741 Bragg : Emin= 0 eV Emax= 1.05231 MeV
1742 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV
[1337]1743
1744hBrems: for kaon+ SubType= 3
1745 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1746 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1747 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1748 hBrem : Emin= 0 eV Emax= 10 TeV
[1337]1749
1750hPairProd: for kaon+ SubType= 4
1751 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1752 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1753 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1754 hPairProd : Emin= 0 eV Emax= 10 TeV
[1337]1755
1756hIoni: for kaon- SubType= 2
1757 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1758 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1759 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1760 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1761 ICRU73QO : Emin= 0 eV Emax= 1.05231 MeV
1762 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV
[1337]1763
1764hBrems: for kaon- SubType= 3
1765 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1766 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1767 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1768 hBrem : Emin= 0 eV Emax= 10 TeV
[1337]1769
1770hPairProd: for kaon- SubType= 4
1771 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1772 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1773 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1774 hPairProd : Emin= 0 eV Emax= 10 TeV
[1337]1775
1776muMsc: for mu+ SubType= 10
1777 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
[1342]1778 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.4592
[1337]1779 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1780 WentzelVIUni : Emin= 0 eV Emax= 10 TeV
[1337]1781
1782muIoni: for mu+ SubType= 2
1783 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1784 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1785 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1786 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1787 Bragg : Emin= 0 eV Emax= 200 keV
1788 BetheBloch : Emin= 200 keV Emax= 1 GeV
1789 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV
[1337]1790
1791muBrems: for mu+ SubType= 3
1792 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1793 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1794 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1795 MuBrem : Emin= 0 eV Emax= 10 TeV
[1337]1796
1797muPairProd: for mu+ SubType= 4
1798 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1799 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1800 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1801 muPairProd : Emin= 0 eV Emax= 10 TeV
[1337]1802
1803CoulombScat: for mu+ SubType= 1
1804 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
[1342]1805 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531
[1337]1806 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1807 eCoulombScattering : Emin= 0 eV Emax= 10 TeV
[1337]1808
1809muIoni: for mu- SubType= 2
1810 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1811 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1812 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1813 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1814 ICRU73QO : Emin= 0 eV Emax= 200 keV
1815 BetheBloch : Emin= 200 keV Emax= 1 GeV
1816 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV
[1337]1817
1818muBrems: for mu- SubType= 3
1819 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1820 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1821 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1822 MuBrem : Emin= 0 eV Emax= 10 TeV
[1337]1823
1824muPairProd: for mu- SubType= 4
1825 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1826 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1827 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1828 muPairProd : Emin= 0 eV Emax= 10 TeV
[1337]1829
1830CoulombScat: for mu- SubType= 1
1831 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
[1342]1832 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531
[1337]1833 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1834 eCoulombScattering : Emin= 0 eV Emax= 10 TeV
[1337]1835
1836hIoni: for pi+ SubType= 2
1837 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1838 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1839 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1840 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1841 Bragg : Emin= 0 eV Emax= 297.505 keV
1842 BetheBloch : Emin= 297.505 keV Emax= 10 TeV
[1337]1843
1844hBrems: for pi+ SubType= 3
1845 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1846 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1847 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1848 hBrem : Emin= 0 eV Emax= 10 TeV
[1337]1849
1850hPairProd: for pi+ SubType= 4
1851 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1852 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1853 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1854 hPairProd : Emin= 0 eV Emax= 10 TeV
[1337]1855
1856msc: for pi- SubType= 10
1857 Lambda tables from 100 eV to 10 TeV in 220 bins, spline: 1
1858 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
1859 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1860 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV
[1337]1861
1862hIoni: for pi- SubType= 2
1863 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1864 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1865 finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
1866 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1867 ICRU73QO : Emin= 0 eV Emax= 297.505 keV
1868 BetheBloch : Emin= 297.505 keV Emax= 10 TeV
[1337]1869
1870hBrems: for pi- SubType= 3
1871 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1872 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1873 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1874 hBrem : Emin= 0 eV Emax= 10 TeV
[1337]1875
1876hPairProd: for pi- SubType= 4
1877 dE/dx and range tables from 100 eV to 10 TeV in 220 bins
1878 Lambda tables from threshold to 10 TeV in 220 bins, spline: 1
1879 ===== EM models for the G4Region DefaultRegionForTheWorld ======
[1342]1880 hPairProd : Emin= 0 eV Emax= 10 TeV
[1337]1881
1882Region <DefaultRegionForTheWorld> -- appears in <World> world volume
1883 Root logical volume(s) : World
1884 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
1885 Materials : Galactic Iron
1886 Production cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm
1887
1888========= Table of registered couples ==============================
1889
1890Index : 0 used in the geometry : Yes recalculation needed : No
1891 Material : Galactic
1892 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm
1893 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 0 eV
1894 Region(s) which use this couple :
1895 DefaultRegionForTheWorld
1896
1897Index : 1 used in the geometry : No recalculation needed : No
1898 Material : Silicon
1899 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm
[1342]1900 Energy thresholds : gamma 6.95018 keV e- 548.291 keV e+ 526.624 keV proton 0 eV
[1337]1901
1902Index : 2 used in the geometry : No recalculation needed : No
1903 Material : Gold
1904 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm
[1342]1905 Energy thresholds : gamma 112.644 keV e- 2.30976 MeV e+ 2.13081 MeV proton 0 eV
[1337]1906
1907Index : 3 used in the geometry : No recalculation needed : No
1908 Material : G4_POLYETHYLENE
1909 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm
[1342]1910 Energy thresholds : gamma 2.26261 keV e- 351.877 keV e+ 342.545 keV proton 0 eV
[1337]1911
1912Index : 4 used in the geometry : No recalculation needed : No
1913 Material : Aluminium
1914 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm
[1342]1915 Energy thresholds : gamma 6.90363 keV e- 598.345 keV e+ 570.85 keV proton 0 eV
[1337]1916
1917Index : 5 used in the geometry : Yes recalculation needed : No
1918 Material : Iron
1919 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm
[1342]1920 Energy thresholds : gamma 20.6438 keV e- 1.29592 MeV e+ 1.21169 MeV proton 0 eV
[1337]1921 Region(s) which use this couple :
1922 DefaultRegionForTheWorld
1923
1924====================================================================
1925
1926Start closing geometry.
1927G4GeometryManager::ReportVoxelStats -- Voxel Statistics
1928
1929 Total memory consumed for geometry optimisation: 0 kByte
1930 Total CPU time elapsed for geometry optimisation: 0 seconds
1931### Run 4 start.
1932
1933--------- Ranecu engine status ---------
1934 Initial seed (index) = 0
[1342]1935 Current couple of seeds = 1430401471, 1601608505
[1337]1936----------------------------------------
1937Start Run processing.
1938
1939---> Begin of Event: 0
1940
1941---> Begin of Event: 10000
1942
1943---> Begin of Event: 20000
1944
1945---> Begin of Event: 30000
1946
1947---> Begin of Event: 40000
1948Run terminated.
1949Run Summary
1950 Number of events processed : 50000
[1342]1951 User=22.69s Real=27.71s Sys=1.28s
[1337]1952
1953 ======================== run summary ======================
1954
1955 The run was 50000 mu+ of 100 GeV through 1 m of Iron (density: 7.87 g/cm3 )
1956
[1342]1957 Total energy deposit in absorber per event = 2.128 GeV +- 14.93 MeV
[1337]1958
[1342]1959 -----> Mean dE/dx = 21.28 MeV/cm (2.704 MeV*cm2/g)
[1337]1960
1961 From formulas :
1962 restricted dEdx = 11.22 MeV/cm (1.425 MeV*cm2/g)
1963 full dEdx = 21.6 MeV/cm (2.745 MeV*cm2/g)
1964
[1342]1965 Leakage : primary = 97.87 GeV +- 14.93 MeV secondaries = 0 eV +- 0 eV
[1337]1966 Energy balance : edep + eleak = 100 GeV
1967
[1342]1968 Total track length (charged) in absorber per event = 99.01 cm +- 52.39 um
[1337]1969 Total track length (neutral) in absorber per event = 0 fm +- 0 fm
1970
[1342]1971 Number of steps (charged) in absorber per event = 47.4 +- 0.001024
[1337]1972 Number of steps (neutral) in absorber per event = 0 +- 0
1973
[1342]1974 Number of secondaries per event : Gammas = 0.04828; electrons = 45.8; positrons = 1.249
[1337]1975
1976 Number of events with the primary particle transmitted = 100 %
1977 Number of events with at least 1 particle transmitted (same charge as primary) = 100 %
1978
1979 Number of events with the primary particle reflected = 0 %
1980 Number of events with at least 1 particle reflected (same charge as primary) = 0 %
1981
1982 MultipleScattering:
[1342]1983 rms proj angle of transmit primary particle = 1.126 mrad (central part only)
[1337]1984 computed theta0 (Highland formula) = 1.182 mrad
[1342]1985 central part defined as +- 3.545 mrad; Tail ratio = 0.308 %
[1337]1986
1987--------- Ranecu engine status ---------
1988 Initial seed (index) = 0
[1342]1989 Current couple of seeds = 1153020195, 500937775
[1337]1990----------------------------------------
1991#
[807]1992UserDetectorConstruction deleted.
1993UserPhysicsList deleted.
1994UserRunAction deleted.
1995UserPrimaryGenerator deleted.
1996G4 kernel has come to Quit state.
1997EventManager deleted.
1998UImanager deleted.
1999Units table cleared.
2000StateManager deleted.
2001RunManagerKernel is deleted.
2002RunManager is deleting.
[1337]2003
Note: See TracBrowser for help on using the repository browser.