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

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

update ti head

File size: 20.5 KB
Line 
1
2 ############################################
3 !!! WARNING - FPE detection is activated !!!
4 ############################################
5
6*************************************************************
7 Geant4 version Name: geant4-09-03-ref-09 (25-June-2010)
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/phys/addPhysics emstandard_opt0
16#
17/run/initialize
18userDetector->Construct() start.
19
20-------------------------------------------------------------
21 ---> The calorimeter is 50 layers of:
22 Lead: 2.3 mm
23 liquidArgon: 5.7 mm
24-------------------------------------------------------------
25
26 Material: Galactic density: 0.000 kg/m3 RadL: 204322111.300 pc Nucl.Int.Length: 114334693.470 pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm
27 ---> Element: Galactic ( ) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
28
29 Material: Lead density: 11.350 g/cm3 RadL: 5.613 mm Nucl.Int.Length: 18.261 cm Imean: 823.000 eV
30 ---> Element: Pb (Pb) Z = 82.0 N = 207.2 A = 207.22 g/mole
31 ---> Isotope: Pb204 Z = 82 N = 204 A = 203.97 g/mole abundance: 1.40 %
32 ---> Isotope: Pb206 Z = 82 N = 206 A = 205.97 g/mole abundance: 24.10 %
33 ---> Isotope: Pb207 Z = 82 N = 207 A = 206.98 g/mole abundance: 22.10 %
34 ---> Isotope: Pb208 Z = 82 N = 208 A = 207.98 g/mole abundance: 52.40 % ElmMassFraction: 100.00 % ElmAbundance 100.00 %
35
36 Material: liquidArgon density: 1.390 g/cm3 RadL: 14.064 cm Nucl.Int.Length: 86.001 cm Imean: 188.000 eV
37 ---> Element: Ar (Ar) Z = 18.0 N = 40.0 A = 39.95 g/mole
38 ---> Isotope: Ar36 Z = 18 N = 36 A = 35.97 g/mole abundance: 0.34 %
39 ---> Isotope: Ar38 Z = 18 N = 38 A = 37.96 g/mole abundance: 0.06 %
40 ---> Isotope: Ar40 Z = 18 N = 40 A = 39.96 g/mole abundance: 99.60 % ElmMassFraction: 100.00 % ElmAbundance 100.00 %
41
42-------------------------------------------------------------
43World is registered to the default region.
44physicsList->Construct() start.
45physicsList->Construct() start.
46physicsList->CheckParticleList() start.
47physicsList->setCut() start.
48PhysicsList::SetCuts:CutLength : 1 mm
49/process/eLoss/verbose 1
50/process/msc/Skin 1
51/run/physicsModified
52#
53/gun/particle e-
54/gun/energy 1 GeV
55#
56/testem/histo/setFileName testem3
57/testem/histo/setFileType root
58/testem/histo/setHisto 1 100 700. 1000. MeV
59----> SetHisto 1: Edep in absorber 1 (MeV); 100 bins from 700 MeV to 1000 MeV
60/testem/histo/setHisto 2 100 0. 300. MeV
61----> SetHisto 2: Edep in absorber 2 (MeV); 100 bins from 0 MeV to 300 MeV
62/testem/histo/setHisto 11 52 0. 52. none
63----> SetHisto 11: longit. profile of Edep (MeV/event) in absorber 1; 52 bins from 0 none to 52 none
64/testem/histo/setHisto 12 52 0. 52. none
65----> SetHisto 12: longit. profile of Edep (MeV/event) in absorber 2; 52 bins from 0 none to 52 none
66/testem/histo/setHisto 21 102 0. 102. none
67----> SetHisto 21: energy flow (MeV/event); 102 bins from 0 none to 102 none
68/testem/histo/setHisto 22 102 0. 102. none
69----> SetHisto 22: lateral energy leak (MeV/event); 102 bins from 0 none to 102 none
70#
71/testem/event/printModulo 50
72#
73/run/beamOn 100
74
75phot: for gamma SubType= 12
76 ===== EM models for the G4Region DefaultRegionForTheWorld ======
77 PhotoElectric : Emin= 0 eV Emax= 10 TeV
78
79compt: for gamma SubType= 13
80 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1
81 ===== EM models for the G4Region DefaultRegionForTheWorld ======
82 Klein-Nishina : Emin= 0 eV Emax= 10 TeV
83
84conv: for gamma SubType= 14
85 Lambda tables from 1.022 MeV to 10 TeV in 77 bins, spline: 1
86 ===== EM models for the G4Region DefaultRegionForTheWorld ======
87 BetheHeitler : Emin= 0 eV Emax= 10 TeV
88
89msc: for e- SubType= 10
90 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1
91 RangeFactor= 0.04, stepLimitType: 1, latDisplacement: 1
92 ===== EM models for the G4Region DefaultRegionForTheWorld ======
93 UrbanMsc93 : Emin= 0 eV Emax= 10 TeV
94
95eIoni: for e- SubType= 2
96 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
97 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
98 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
99 ===== EM models for the G4Region DefaultRegionForTheWorld ======
100 MollerBhabha : Emin= 0 eV Emax= 10 TeV
101
102eBrem: for e- SubType= 3
103 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
104 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
105 LPM flag: 1 for E > 1 GeV
106 ===== EM models for the G4Region DefaultRegionForTheWorld ======
107 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban
108 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban
109
110eIoni: for e+ SubType= 2
111 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
112 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
113 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
114 ===== EM models for the G4Region DefaultRegionForTheWorld ======
115 MollerBhabha : Emin= 0 eV Emax= 10 TeV
116
117eBrem: for e+ SubType= 3
118 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
119 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
120 LPM flag: 1 for E > 1 GeV
121 ===== EM models for the G4Region DefaultRegionForTheWorld ======
122 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban
123 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban
124
125annihil: for e+ SubType= 5
126 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1
127 ===== EM models for the G4Region DefaultRegionForTheWorld ======
128 eplus2gg : Emin= 0 eV Emax= 10 TeV
129
130msc: for proton SubType= 10
131 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1
132 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
133 ===== EM models for the G4Region DefaultRegionForTheWorld ======
134 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV
135
136hIoni: for proton SubType= 2
137 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
138 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
139 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
140 ===== EM models for the G4Region DefaultRegionForTheWorld ======
141 Bragg : Emin= 0 eV Emax= 2 MeV
142 BetheBloch : Emin= 2 MeV Emax= 10 TeV
143
144hBrems: for proton SubType= 3
145 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
146 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
147 ===== EM models for the G4Region DefaultRegionForTheWorld ======
148 hBrem : Emin= 0 eV Emax= 10 TeV
149
150hPairProd: for proton SubType= 4
151 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
152 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
153 ===== EM models for the G4Region DefaultRegionForTheWorld ======
154 hPairProd : Emin= 0 eV Emax= 10 TeV
155
156msc: for GenericIon SubType= 10
157 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
158 ===== EM models for the G4Region DefaultRegionForTheWorld ======
159 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV
160
161ionIoni: for GenericIon SubType= 2
162 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
163 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
164 finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
165 Stopping Power data for 17 ion/material pairs
166 ===== EM models for the G4Region DefaultRegionForTheWorld ======
167 BraggIon : Emin= 0 eV Emax= 2 MeV
168 BetheBloch : Emin= 2 MeV Emax= 10 TeV
169
170hIoni: for anti_proton SubType= 2
171 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
172 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
173 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
174 ===== EM models for the G4Region DefaultRegionForTheWorld ======
175 ICRU73QO : Emin= 0 eV Emax= 2 MeV
176 BetheBloch : Emin= 2 MeV Emax= 10 TeV
177
178msc: for kaon+ SubType= 10
179 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1
180 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
181 ===== EM models for the G4Region DefaultRegionForTheWorld ======
182 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV
183
184hIoni: for kaon+ SubType= 2
185 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
186 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
187 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
188 ===== EM models for the G4Region DefaultRegionForTheWorld ======
189 Bragg : Emin= 0 eV Emax= 1.05231 MeV
190 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV
191
192hBrems: for kaon+ SubType= 3
193 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
194 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
195 ===== EM models for the G4Region DefaultRegionForTheWorld ======
196 hBrem : Emin= 0 eV Emax= 10 TeV
197
198hPairProd: for kaon+ SubType= 4
199 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
200 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
201 ===== EM models for the G4Region DefaultRegionForTheWorld ======
202 hPairProd : Emin= 0 eV Emax= 10 TeV
203
204hIoni: for kaon- SubType= 2
205 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
206 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
207 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
208 ===== EM models for the G4Region DefaultRegionForTheWorld ======
209 ICRU73QO : Emin= 0 eV Emax= 1.05231 MeV
210 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV
211
212hBrems: for kaon- SubType= 3
213 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
214 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
215 ===== EM models for the G4Region DefaultRegionForTheWorld ======
216 hBrem : Emin= 0 eV Emax= 10 TeV
217
218hPairProd: for kaon- SubType= 4
219 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
220 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
221 ===== EM models for the G4Region DefaultRegionForTheWorld ======
222 hPairProd : Emin= 0 eV Emax= 10 TeV
223
224muMsc: for mu+ SubType= 10
225 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1
226 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.4592
227 ===== EM models for the G4Region DefaultRegionForTheWorld ======
228 WentzelVIUni : Emin= 0 eV Emax= 10 TeV
229
230muIoni: for mu+ SubType= 2
231 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
232 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
233 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
234 ===== EM models for the G4Region DefaultRegionForTheWorld ======
235 Bragg : Emin= 0 eV Emax= 200 keV
236 BetheBloch : Emin= 200 keV Emax= 1 GeV
237 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV
238
239muBrems: for mu+ SubType= 3
240 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
241 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
242 ===== EM models for the G4Region DefaultRegionForTheWorld ======
243 MuBrem : Emin= 0 eV Emax= 10 TeV
244
245muPairProd: for mu+ SubType= 4
246 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
247 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
248 ===== EM models for the G4Region DefaultRegionForTheWorld ======
249 muPairProd : Emin= 0 eV Emax= 10 TeV
250
251CoulombScat: for mu+ SubType= 1
252 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1
253 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531
254 ===== EM models for the G4Region DefaultRegionForTheWorld ======
255 eCoulombScattering : Emin= 0 eV Emax= 10 TeV
256
257muIoni: for mu- SubType= 2
258 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
259 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
260 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
261 ===== EM models for the G4Region DefaultRegionForTheWorld ======
262 ICRU73QO : Emin= 0 eV Emax= 200 keV
263 BetheBloch : Emin= 200 keV Emax= 1 GeV
264 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV
265
266muBrems: for mu- SubType= 3
267 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
268 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
269 ===== EM models for the G4Region DefaultRegionForTheWorld ======
270 MuBrem : Emin= 0 eV Emax= 10 TeV
271
272muPairProd: for mu- SubType= 4
273 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
274 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
275 ===== EM models for the G4Region DefaultRegionForTheWorld ======
276 muPairProd : Emin= 0 eV Emax= 10 TeV
277
278CoulombScat: for mu- SubType= 1
279 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1
280 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531
281 ===== EM models for the G4Region DefaultRegionForTheWorld ======
282 eCoulombScattering : Emin= 0 eV Emax= 10 TeV
283
284hIoni: for pi+ SubType= 2
285 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
286 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
287 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
288 ===== EM models for the G4Region DefaultRegionForTheWorld ======
289 Bragg : Emin= 0 eV Emax= 297.505 keV
290 BetheBloch : Emin= 297.505 keV Emax= 10 TeV
291
292hBrems: for pi+ SubType= 3
293 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
294 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
295 ===== EM models for the G4Region DefaultRegionForTheWorld ======
296 hBrem : Emin= 0 eV Emax= 10 TeV
297
298hPairProd: for pi+ SubType= 4
299 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
300 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
301 ===== EM models for the G4Region DefaultRegionForTheWorld ======
302 hPairProd : Emin= 0 eV Emax= 10 TeV
303
304msc: for pi- SubType= 10
305 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1
306 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
307 ===== EM models for the G4Region DefaultRegionForTheWorld ======
308 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV
309
310hIoni: for pi- SubType= 2
311 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
312 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
313 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
314 ===== EM models for the G4Region DefaultRegionForTheWorld ======
315 ICRU73QO : Emin= 0 eV Emax= 297.505 keV
316 BetheBloch : Emin= 297.505 keV Emax= 10 TeV
317
318hBrems: for pi- SubType= 3
319 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
320 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
321 ===== EM models for the G4Region DefaultRegionForTheWorld ======
322 hBrem : Emin= 0 eV Emax= 10 TeV
323
324hPairProd: for pi- SubType= 4
325 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
326 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
327 ===== EM models for the G4Region DefaultRegionForTheWorld ======
328 hPairProd : Emin= 0 eV Emax= 10 TeV
329
330Region <DefaultRegionForTheWorld> -- appears in <World> world volume
331 Root logical volume(s) : World
332 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
333 Materials : Galactic Lead liquidArgon
334 Production cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
335
336========= Table of registered couples ==============================
337
338Index : 0 used in the geometry : Yes recalculation needed : No
339 Material : Galactic
340 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
341 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 100 keV
342 Region(s) which use this couple :
343 DefaultRegionForTheWorld
344
345Index : 1 used in the geometry : Yes recalculation needed : No
346 Material : Lead
347 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
348 Energy thresholds : gamma 101.843 keV e- 1.36749 MeV e+ 1.27862 MeV proton 100 keV
349 Region(s) which use this couple :
350 DefaultRegionForTheWorld
351
352Index : 2 used in the geometry : Yes recalculation needed : No
353 Material : liquidArgon
354 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
355 Energy thresholds : gamma 6.19986 keV e- 344.855 keV e+ 333.46 keV proton 100 keV
356 Region(s) which use this couple :
357 DefaultRegionForTheWorld
358
359====================================================================
360
361Start closing geometry.
362G4GeometryManager::ReportVoxelStats -- Voxel Statistics
363
364 Total memory consumed for geometry optimisation: 3 kByte
365 Total CPU time elapsed for geometry optimisation: 0 seconds
366
367 Voxelisation: top CPU users:
368 Percent Total CPU System CPU Memory Volume
369 ------- ---------- ---------- -------- ----------
370 0.00 0.00 0.00 3k Calorimeter
371 0.00 0.00 0.00 0k Layer
372
373 Voxelisation: top memory users:
374 Percent Memory Heads Nodes Pointers Total CPU Volume
375 ------- -------- ------ ------ -------- ---------- ----------
376 92.05 2k 1 50 50 0.00 Calorimeter
377 7.95 0k 1 3 4 0.00 Layer
378### Run 0 start.
379
380--------- Ranecu engine status ---------
381 Initial seed (index) = 0
382 Current couple of seeds = 9876, 54321
383----------------------------------------
384Start Run processing.
385
386---> Begin Of Event: 0
387
388---> Begin Of Event: 50
389Run terminated.
390Run Summary
391 Number of events processed : 100
392 User=5.11s Real=5.11s Sys=0.01s
393
394------------------------------------------------------------
395 material Edep RMS sqrt(E0(GeV))*rmsE/Emean total tracklen
396
397 Lead: 782.46 MeV : 15.18 MeV 1.94 +- 0.194 % 54.6 cm +- 1.11 cm
398 liquidArgon: 206.47 MeV : 13.75 MeV 6.66 +- 0.666 % 1.02 m +- 7.05 cm
399
400------------------------------------------------------------
401 Beam particle e- E = 1 GeV
402 Mean number of gamma 571
403 Mean number of e- 962
404 Mean number of e+ 55.5
405------------------------------------------------------------
406
407 Energy deposition from Energy flow balance :
408 material Total Edep
409
410 Lead: 782.434 MeV
411liquidArgon: 206.478 MeV
412
413------------------------------------------------------------
414
415
416--------- Ranecu engine status ---------
417 Initial seed (index) = 0
418 Current couple of seeds = 225336118, 351963897
419----------------------------------------
420UserDetectorConstruction deleted.
421UserPhysicsList deleted.
422UserRunAction deleted.
423UserPrimaryGenerator deleted.
424G4 kernel has come to Quit state.
425EventManager deleted.
426UImanager deleted.
427Units table cleared.
428StateManager deleted.
429RunManagerKernel is deleted.
430RunManager is deleting.
431
Note: See TracBrowser for help on using the repository browser.