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

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

update ti head

File size: 19.3 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
14***** Table : Nb of materials = 8 *****
15
16 Material: G4_Al density: 2.699 g/cm3 RadL: 8.896 cm Nucl.Int.Length: 38.877 cm Imean: 166.000 eV
17 ---> Element: Al (Al) Z = 13.0 N = 27.0 A = 26.98 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
18
19 Material: G4_Si density: 2.330 g/cm3 RadL: 9.366 cm Nucl.Int.Length: 45.753 cm Imean: 173.000 eV
20 ---> Element: Si (Si) Z = 14.0 N = 28.1 A = 28.09 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
21
22 Material: G4_Fe density: 7.874 g/cm3 RadL: 1.757 cm Nucl.Int.Length: 16.959 cm Imean: 286.000 eV
23 ---> Element: Fe (Fe) Z = 26.0 N = 55.8 A = 55.85 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
24
25 Material: G4_Ge density: 5.323 g/cm3 RadL: 2.301 cm Nucl.Int.Length: 27.334 cm Imean: 350.000 eV
26 ---> Element: Ge (Ge) Z = 32.0 N = 72.6 A = 72.61 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
27
28 Material: G4_W density: 19.300 g/cm3 RadL: 3.504 mm Nucl.Int.Length: 10.306 cm Imean: 727.000 eV
29 ---> Element: W (W) Z = 74.0 N = 183.8 A = 183.84 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
30
31 Material: G4_Pb density: 11.350 g/cm3 RadL: 5.613 mm Nucl.Int.Length: 18.261 cm Imean: 823.000 eV
32 ---> Element: Pb (Pb) Z = 82.0 N = 207.2 A = 207.22 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
33
34 Material: G4_AIR density: 1.205 mg/cm3 RadL: 303.921 m Nucl.Int.Length: 710.261 m Imean: 85.700 eV temperature: 273.15 K pressure: 1.00 atm
35 ---> Element: C (C) Z = 6.0 N = 12.0 A = 12.01 g/mole ElmMassFraction: 0.01 % ElmAbundance 0.02 %
36 ---> Element: N (N) Z = 7.0 N = 14.0 A = 14.01 g/mole ElmMassFraction: 75.53 % ElmAbundance 78.44 %
37 ---> Element: O (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 23.18 % ElmAbundance 21.07 %
38 ---> Element: Ar (Ar) Z = 18.0 N = 39.9 A = 39.95 g/mole ElmMassFraction: 1.28 % ElmAbundance 0.47 %
39
40 Material: G4_WATER H_2O density: 1.000 g/cm3 RadL: 36.083 cm Nucl.Int.Length: 75.517 cm Imean: 78.000 eV
41 ---> Element: H (H) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 11.19 % ElmAbundance 66.67 %
42 ---> Element: O (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 88.81 % ElmAbundance 33.33 %
43
44
45/run/verbose 2
46#
47/testem/det/setMat G4_WATER
48/testem/det/setRadius 3 cm
49/testem/det/setNbOfLayers 1
50#
51/testem/phys/addPhysics local
52PhysicsList::AddPhysicsList: <local>
53#
54/run/initialize
55userDetector->Construct() start.
56
57---------------------------------------------------------
58---> The Absorber is a sphere of 3 cm radius of G4_WATER divided in 1 slices of 3 cm
59
60---------------------------------------------------------
61Absorber is registered to the default region.
62physicsList->Construct() start.
63physicsList->Construct() start.
64physicsList->CheckParticleList() start.
65physicsList->setCut() start.
66PhysicsList::SetCuts:CutLength : 1 mm
67#
68/testem/phys/setCuts 1 um
69#
70/gun/particle e-
71/gun/energy 4 MeV
72#
73/testem/histo/setFileType root
74/testem/histo/setFileName testem12
75/testem/histo/setHisto 1 120 0. 2.4 cm
76----> SetHisto 1: dE/dr (MeV/mm) along radius (cm); 120 bins from 0 cm to 2.4 cm
77/testem/histo/setHisto 2 180 0. 4.5 MeV
78----> SetHisto 2: total Energy deposited in absorber (MeV); 180 bins from 0 MeV to 4.5 MeV
79/testem/histo/setHisto 3 100 0. 3. cm
80----> SetHisto 3: true track length of the primary particle (cm); 100 bins from 0 cm to 3 cm
81/testem/histo/setHisto 4 100 0. 500. um
82----> SetHisto 4: true step size of the primary particle (um); 100 bins from 0 um to 500 um
83/testem/histo/setHisto 5 200 -3. 3. cm
84----> SetHisto 5: projected range of the primary particle (cm); 200 bins from -3 cm to 3 cm
85/testem/histo/setcsdaRange 2.037 cm
86/testem/histo/setHisto 8 120 0. 1.2 none
87----> SetHisto 8: d(E/E0)/d(r/r0) along r/r0; 120 bins from 0 none to 1.2 none
88#
89/testem/event/printModulo 1000
90#
91/run/beamOn 10000
92
93phot: for gamma SubType= 12
94 ===== EM models for the G4Region DefaultRegionForTheWorld ======
95 PhotoElectric : Emin= 0 eV Emax= 10 TeV
96
97compt: for gamma SubType= 13
98 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1
99 ===== EM models for the G4Region DefaultRegionForTheWorld ======
100 Klein-Nishina : Emin= 0 eV Emax= 10 TeV
101
102conv: for gamma SubType= 14
103 Lambda tables from 1.022 MeV to 10 TeV in 77 bins, spline: 1
104 ===== EM models for the G4Region DefaultRegionForTheWorld ======
105 BetheHeitler : Emin= 0 eV Emax= 10 TeV
106
107msc: for e- SubType= 10
108 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1
109 RangeFactor= 0.04, stepLimitType: 1, latDisplacement: 1
110 ===== EM models for the G4Region DefaultRegionForTheWorld ======
111 UrbanMsc93 : Emin= 0 eV Emax= 10 TeV
112
113eIoni: for e- SubType= 2
114 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
115 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
116 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
117 ===== EM models for the G4Region DefaultRegionForTheWorld ======
118 MollerBhabha : Emin= 0 eV Emax= 10 TeV
119 CSDA range table up to 1 GeV in 160 bins
120
121eBrem: for e- SubType= 3
122 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
123 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
124 LPM flag: 1 for E > 1 GeV
125 ===== EM models for the G4Region DefaultRegionForTheWorld ======
126 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban
127 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban
128
129eIoni: for e+ SubType= 2
130 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
131 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
132 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
133 ===== EM models for the G4Region DefaultRegionForTheWorld ======
134 MollerBhabha : Emin= 0 eV Emax= 10 TeV
135 CSDA range table up to 1 GeV in 160 bins
136
137eBrem: for e+ SubType= 3
138 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
139 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
140 LPM flag: 1 for E > 1 GeV
141 ===== EM models for the G4Region DefaultRegionForTheWorld ======
142 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban
143 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban
144
145annihil: for e+ SubType= 5
146 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1
147 ===== EM models for the G4Region DefaultRegionForTheWorld ======
148 eplus2gg : Emin= 0 eV Emax= 10 TeV
149
150msc: for proton SubType= 10
151 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1
152 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
153 ===== EM models for the G4Region DefaultRegionForTheWorld ======
154 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV
155
156hIoni: for proton SubType= 2
157 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
158 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
159 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
160 ===== EM models for the G4Region DefaultRegionForTheWorld ======
161 Bragg : Emin= 0 eV Emax= 2 MeV
162 BetheBloch : Emin= 2 MeV Emax= 10 TeV
163 CSDA range table up to 1 GeV in 160 bins
164
165hBrems: for proton SubType= 3
166 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
167 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
168 ===== EM models for the G4Region DefaultRegionForTheWorld ======
169 hBrem : Emin= 0 eV Emax= 10 TeV
170
171hPairProd: for proton SubType= 4
172 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
173 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
174 ===== EM models for the G4Region DefaultRegionForTheWorld ======
175 hPairProd : Emin= 0 eV Emax= 10 TeV
176
177msc: for GenericIon SubType= 10
178 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
179 ===== EM models for the G4Region DefaultRegionForTheWorld ======
180 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV
181
182ionIoni: for GenericIon SubType= 2
183 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
184 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
185 finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
186 Stopping Power data for 17 ion/material pairs
187 ===== EM models for the G4Region DefaultRegionForTheWorld ======
188 ParamICRU73 : Emin= 0 eV Emax= 10 TeV
189 CSDA range table up to 1 GeV in 160 bins
190
191nuclearStopping: for GenericIon SubType= 1
192 ===== EM models for the G4Region DefaultRegionForTheWorld ======
193 ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV
194
195nuclearStopping: for alpha SubType= 1
196 ===== EM models for the G4Region DefaultRegionForTheWorld ======
197 ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV
198
199hIoni: for anti_proton SubType= 2
200 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
201 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
202 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
203 ===== EM models for the G4Region DefaultRegionForTheWorld ======
204 ICRU73QO : Emin= 0 eV Emax= 2 MeV
205 BetheBloch : Emin= 2 MeV Emax= 10 TeV
206 CSDA range table up to 1 GeV in 160 bins
207
208msc: for kaon+ SubType= 10
209 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1
210 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
211 ===== EM models for the G4Region DefaultRegionForTheWorld ======
212 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV
213
214hIoni: for kaon+ SubType= 2
215 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
216 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
217 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
218 ===== EM models for the G4Region DefaultRegionForTheWorld ======
219 Bragg : Emin= 0 eV Emax= 1.05231 MeV
220 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV
221 CSDA range table up to 1 GeV in 160 bins
222
223hIoni: for kaon- SubType= 2
224 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
225 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
226 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
227 ===== EM models for the G4Region DefaultRegionForTheWorld ======
228 ICRU73QO : Emin= 0 eV Emax= 1.05231 MeV
229 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV
230 CSDA range table up to 1 GeV in 160 bins
231
232muMsc: for mu+ SubType= 10
233 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1
234 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 0
235 ===== EM models for the G4Region DefaultRegionForTheWorld ======
236 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV
237
238muIoni: for mu+ SubType= 2
239 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
240 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
241 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
242 ===== EM models for the G4Region DefaultRegionForTheWorld ======
243 Bragg : Emin= 0 eV Emax= 200 keV
244 BetheBloch : Emin= 200 keV Emax= 1 GeV
245 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV
246 CSDA range table up to 1 GeV in 160 bins
247
248muBrems: for mu+ SubType= 3
249 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
250 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
251 ===== EM models for the G4Region DefaultRegionForTheWorld ======
252 MuBrem : Emin= 0 eV Emax= 10 TeV
253
254muPairProd: for mu+ SubType= 4
255 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
256 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
257 ===== EM models for the G4Region DefaultRegionForTheWorld ======
258 muPairProd : Emin= 0 eV Emax= 10 TeV
259
260muIoni: for mu- SubType= 2
261 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
262 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
263 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
264 ===== EM models for the G4Region DefaultRegionForTheWorld ======
265 ICRU73QO : Emin= 0 eV Emax= 200 keV
266 BetheBloch : Emin= 200 keV Emax= 1 GeV
267 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV
268 CSDA range table up to 1 GeV in 160 bins
269
270muBrems: for mu- SubType= 3
271 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
272 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
273 ===== EM models for the G4Region DefaultRegionForTheWorld ======
274 MuBrem : Emin= 0 eV Emax= 10 TeV
275
276muPairProd: for mu- SubType= 4
277 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
278 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
279 ===== EM models for the G4Region DefaultRegionForTheWorld ======
280 muPairProd : Emin= 0 eV Emax= 10 TeV
281
282hIoni: for pi+ SubType= 2
283 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
284 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
285 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
286 ===== EM models for the G4Region DefaultRegionForTheWorld ======
287 Bragg : Emin= 0 eV Emax= 297.505 keV
288 BetheBloch : Emin= 297.505 keV Emax= 10 TeV
289 CSDA range table up to 1 GeV in 160 bins
290
291hBrems: for pi+ SubType= 3
292 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
293 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
294 ===== EM models for the G4Region DefaultRegionForTheWorld ======
295 hBrem : Emin= 0 eV Emax= 10 TeV
296
297hPairProd: for pi+ SubType= 4
298 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
299 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
300 ===== EM models for the G4Region DefaultRegionForTheWorld ======
301 hPairProd : Emin= 0 eV Emax= 10 TeV
302
303msc: for pi- SubType= 10
304 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1
305 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
306 ===== EM models for the G4Region DefaultRegionForTheWorld ======
307 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV
308
309hIoni: for pi- SubType= 2
310 dE/dx and range tables from 100 eV to 10 TeV in 77 bins
311 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
312 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
313 ===== EM models for the G4Region DefaultRegionForTheWorld ======
314 ICRU73QO : Emin= 0 eV Emax= 297.505 keV
315 BetheBloch : Emin= 297.505 keV Emax= 10 TeV
316 CSDA range table up to 1 GeV in 160 bins
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 <Absorber> world volume
331 Root logical volume(s) : Absorber
332 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
333 Materials : G4_WATER
334 Production cuts : gamma 1 um e- 1 um e+ 1 um proton 0 fm
335
336========= Table of registered couples ==============================
337
338Index : 0 used in the geometry : Yes recalculation needed : No
339 Material : G4_WATER
340 Range cuts : gamma 1 um e- 1 um e+ 1 um proton 0 fm
341 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 0 eV
342 Region(s) which use this couple :
343 DefaultRegionForTheWorld
344
345====================================================================
346
347Start closing geometry.
348G4GeometryManager::ReportVoxelStats -- Voxel Statistics
349
350 Total memory consumed for geometry optimisation: 0 kByte
351 Total CPU time elapsed for geometry optimisation: 0 seconds
352### Run 0 start.
353
354--------- Ranecu engine status ---------
355 Initial seed (index) = 0
356 Current couple of seeds = 9876, 54321
357----------------------------------------
358
359---> stepMax from HistoManager = 200 um
360Start Run processing.
361
362---> Begin of Event: 0
363
364---> Begin of Event: 1000
365
366---> Begin of Event: 2000
367
368---> Begin of Event: 3000
369
370---> Begin of Event: 4000
371
372---> Begin of Event: 5000
373
374---> Begin of Event: 6000
375
376---> Begin of Event: 7000
377
378---> Begin of Event: 8000
379
380---> Begin of Event: 9000
381Run terminated.
382Run Summary
383 Number of events processed : 10000
384 User=39.49s Real=39.73s Sys=0.22s
385
386 ======================== run summary ======================
387
388 The run consists of 10000 e- of 4.00 MeV through 3.00 cm of G4_WATER (density: 1.00 g/cm3 )
389
390 ============================================================
391
392 Total Energy deposited = 3.937 MeV +- 225.962 keV
393
394 Track length of primary track = 2.038 cm +- 2.841 mm
395 Range from EmCalculator = 2.048 cm (from full dE/dx)
396
397 Projected range = 1.306 cm +- 4.241 mm
398
399 Nb of steps of primary track = 240.82 +- 31.65 Step size= 84.619 um +- 4.544 um
400
401--------- Ranecu engine status ---------
402 Initial seed (index) = 0
403 Current couple of seeds = 1150369497, 2116195687
404----------------------------------------
405UserDetectorConstruction deleted.
406UserPhysicsList deleted.
407UserRunAction deleted.
408UserPrimaryGenerator deleted.
409G4 kernel has come to Quit state.
410EventManager deleted.
411UImanager deleted.
412Units table cleared.
413StateManager deleted.
414RunManagerKernel is deleted.
415RunManager is deleting.
416
Note: See TracBrowser for help on using the repository browser.