source: trunk/examples/extended/electromagnetic/TestEm9/TestEm9.out@ 1230

Last change on this file since 1230 was 807, checked in by garnier, 17 years ago

update

File size: 15.6 KB
Line 
1
2*************************************************************
3 Geant4 version Name: global-V09-00-03 (9-May-2008)
4 Copyright : Geant4 Collaboration
5 Reference : NIM A 506 (2003), 250-303
6 WWW : http://cern.ch/geant4
7*************************************************************
8
9### NIST DataBase for Materials is used
10G4Material: Mean excitation energy is changed for G4_AIR Iold= 85.627eV; Inew= 85.7 eV;
11G4Material: Mean excitation energy is changed for G4_Al Iold= 166.4eV; Inew= 166 eV;
12G4Material: Mean excitation energy is changed for G4_Si Iold= 173.6eV; Inew= 173 eV;
13G4Material: Mean excitation energy is changed for G4_Fe Iold= 286eV; Inew= 286 eV;
14G4Material: Mean excitation energy is changed for G4_CESIUM_IODIDE Iold= 491.166eV; Inew= 553.1 eV;
15/run/verbose 1
16#/tracking/verbose 1
17#
18# select EM PhysicsList
19#
20/testem/phys/addPhysics emstandard_opt1
21PhysicsList::Set emstandard_opt1 EM physics
22#
23/testem/phys/setCuts 0.3 mm
24/testem/phys/VertexCuts 0.01 mm
25/testem/phys/MuonCuts 10.0 mm
26/testem/stepMax 100 mm
27#
28/testem/det/acceptance1 0.8271 0.02974 2.5
29/testem/det/acceptance9 0.9570 0.01286 2.5
30/testem/det/acceptance25 0.9793 0.00884 2.5
31#
32#/process/eLoss/preciseRange true
33#/process/eLoss/integral true
34/run/initialize
35
36***** Table : Nb of materials = 5 *****
37
38 Material: G4_AIR density: 1.205 mg/cm3 RadL: 303.921 m Imean: 85.700 eV temperature: 273.15 K pressure: 1.00 atm
39 ---> Element: C (C) Z = 6.0 N = 12.0 A = 12.01 g/mole
40 ---> Isotope: C Z = 6 N = 12 A = 12.00 g/mole abundance: 98.93 %
41 ---> Isotope: C Z = 6 N = 13 A = 13.00 g/mole abundance: 1.07 % ElmMassFraction: 0.01 % ElmAbundance 0.02 %
42 ---> Element: N (N) Z = 7.0 N = 14.0 A = 14.01 g/mole
43 ---> Isotope: N Z = 7 N = 14 A = 14.00 g/mole abundance: 99.63 %
44 ---> Isotope: N Z = 7 N = 15 A = 15.00 g/mole abundance: 0.37 % ElmMassFraction: 75.53 % ElmAbundance 78.44 %
45 ---> Element: O (O) Z = 8.0 N = 16.0 A = 16.00 g/mole
46 ---> Isotope: O Z = 8 N = 16 A = 15.99 g/mole abundance: 99.76 %
47 ---> Isotope: O Z = 8 N = 17 A = 17.00 g/mole abundance: 0.04 %
48 ---> Isotope: O Z = 8 N = 18 A = 18.00 g/mole abundance: 0.20 % ElmMassFraction: 23.18 % ElmAbundance 21.07 %
49 ---> Element: Ar (Ar) Z = 18.0 N = 40.0 A = 39.95 g/mole
50 ---> Isotope: Ar Z = 18 N = 36 A = 35.97 g/mole abundance: 0.34 %
51 ---> Isotope: Ar Z = 18 N = 38 A = 37.96 g/mole abundance: 0.06 %
52 ---> Isotope: Ar Z = 18 N = 40 A = 39.96 g/mole abundance: 99.60 % ElmMassFraction: 1.28 % ElmAbundance 0.47 %
53
54 Material: G4_Al density: 2.699 g/cm3 RadL: 8.896 cm Imean: 166.000 eV
55 ---> Element: Al (Al) Z = 13.0 N = 27.0 A = 26.98 g/mole
56 ---> Isotope: Al Z = 13 N = 27 A = 26.98 g/mole abundance: 100.00 % ElmMassFraction: 100.00 % ElmAbundance 100.00 %
57
58 Material: G4_Si density: 2.330 g/cm3 RadL: 9.366 cm Imean: 173.000 eV
59 ---> Element: Si (Si) Z = 14.0 N = 28.1 A = 28.09 g/mole
60 ---> Isotope: Si Z = 14 N = 28 A = 27.98 g/mole abundance: 92.23 %
61 ---> Isotope: Si Z = 14 N = 29 A = 28.98 g/mole abundance: 4.68 %
62 ---> Isotope: Si Z = 14 N = 30 A = 29.97 g/mole abundance: 3.09 % ElmMassFraction: 100.00 % ElmAbundance 100.00 %
63
64 Material: G4_Fe density: 7.874 g/cm3 RadL: 1.757 cm Imean: 286.000 eV
65 ---> Element: Fe (Fe) Z = 26.0 N = 55.9 A = 55.85 g/mole
66 ---> Isotope: Fe Z = 26 N = 54 A = 53.94 g/mole abundance: 5.84 %
67 ---> Isotope: Fe Z = 26 N = 56 A = 55.93 g/mole abundance: 91.75 %
68 ---> Isotope: Fe Z = 26 N = 57 A = 56.94 g/mole abundance: 2.12 %
69 ---> Isotope: Fe Z = 26 N = 58 A = 57.93 g/mole abundance: 0.28 % ElmMassFraction: 100.00 % ElmAbundance 100.00 %
70
71 Material: G4_CESIUM_IODIDE density: 4.510 g/cm3 RadL: 1.860 cm Imean: 553.100 eV
72 ---> Element: I (I) Z = 53.0 N = 127.0 A = 126.90 g/mole
73 ---> Isotope: I Z = 53 N = 127 A = 126.90 g/mole abundance: 100.00 % ElmMassFraction: 48.85 % ElmAbundance 50.00 %
74 ---> Element: Cs (Cs) Z = 55.0 N = 133.0 A = 132.91 g/mole
75 ---> Isotope: Cs Z = 55 N = 133 A = 132.91 g/mole abundance: 100.00 % ElmMassFraction: 51.15 % ElmAbundance 50.00 %
76
77
78Ecal is 36 cm of G4_CESIUM_IODIDE
79Absorber is 2 mm of G4_Al
80Vertex is 3 cm of 3 layers of Si of 100 um npads= 601
81### New geometry is constructed
82PhysicsList: world cuts are set cutG= 0.3 mm cutE= 0.3 mm
83Vertex cuts are set
84Muon cuts are set
85#
86#/testem/histo/fileName e_pr.root
87#/testem/histo/fileType root
88/gun/particle e-
89/gun/energy 1 GeV
90/run/beamOn 100
91
92phot: Total cross sections from Sandia parametrisation.
93 Sampling according PhotoElectric model
94
95compt: Total cross sections has a good parametrisation from 10 KeV to (100/Z) GeV
96 Sampling according Klein-Nishina model
97 tables are built for gamma
98 Lambda tables from 100 eV to 100 GeV in 90 bins.
99
100conv: Total cross sections has a good parametrisation from 1.5 MeV to 100 GeV for all Z;
101 sampling secondary e+e- according Bethe-Heitler model
102 tables are built for gamma
103 Lambda tables from 1.022 MeV to 100 GeV in 100 bins.
104
105msc: Model variant of multiple scattering for e-
106 Lambda tables from 100 eV to 100 TeV in 120 bins.
107 LateralDisplacementFlag= 1 Skin= 0
108 Boundary/stepping algorithm is active with RangeFactor= 0.2 Step limit type 0
109
110eIoni: tables are built for e-
111 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
112 Lambda tables from threshold to 100 TeV in 120 bins.
113 Delta cross sections and sampling from MollerBhabha model
114 Good description from 1 KeV to 100 GeV.
115 Step function: finalRange(mm)= 1, dRoverRange= 0.8, integral: 1, fluct: 1
116
117eBrem: tables are built for e-
118 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
119 Lambda tables from threshold to 100 TeV in 120 bins.
120 Total cross sections and sampling from StandBrem model (based on the EEDL data library)
121 Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. LPM flag 1
122
123eIoni: tables are built for e+
124 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
125 Lambda tables from threshold to 100 TeV in 120 bins.
126 Delta cross sections and sampling from MollerBhabha model
127 Good description from 1 KeV to 100 GeV.
128 Step function: finalRange(mm)= 1, dRoverRange= 0.8, integral: 1, fluct: 1
129
130eBrem: tables are built for e+
131 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
132 Lambda tables from threshold to 100 TeV in 120 bins.
133 Total cross sections and sampling from StandBrem model (based on the EEDL data library)
134 Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. LPM flag 1
135
136annihil: Sampling according eplus2gg model
137 tables are built for e+
138 Lambda tables from 100 eV to 100 TeV in 120 bins.
139
140msc: Model variant of multiple scattering for proton
141 Lambda tables from 100 eV to 100 TeV in 120 bins.
142 LateralDisplacementFlag= 1 Skin= 0
143 Boundary/stepping algorithm is active with RangeFactor= 0.2 Step limit type 0
144
145hIoni: tables are built for proton
146 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
147 Lambda tables from threshold to 100 TeV in 120 bins.
148 Scaling relation is used from proton dE/dx and range.
149 Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV
150 Parametrisation from Bragg for protons below. NuclearStopping= 1
151 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
152
153msc: Model variant of multiple scattering for GenericIon
154 LateralDisplacementFlag= 0 Skin= 0
155 Boundary/stepping algorithm is active with RangeFactor= 0.2 Step limit type 1
156
157ionIoni: tables are built for GenericIon
158 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
159 Lambda tables from threshold to 100 TeV in 120 bins.
160 Scaling relation is used from proton dE/dx and range.
161 Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV
162 Parametrisation from BraggIon for protons below. NuclearStopping= 1
163
164 Stopping Power data for 8 ion/material pairs are used.
165 Step function: finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1, fluct: 1
166
167hIoni: tables are built for anti_proton
168 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
169 Lambda tables from threshold to 100 TeV in 120 bins.
170 Scaling relation is used from proton dE/dx and range.
171 Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV
172 Parametrisation from Bragg for protons below. NuclearStopping= 1
173 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
174
175msc: Model variant of multiple scattering for mu+
176 Lambda tables from 100 eV to 100 TeV in 120 bins.
177 LateralDisplacementFlag= 1 Skin= 0
178 Boundary/stepping algorithm is active with RangeFactor= 0.2 Step limit type 0
179
180muIoni: tables are built for mu+
181 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
182 Lambda tables from threshold to 100 TeV in 120 bins.
183 Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below,
184 radiative corrections for E > 1 GeV
185 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
186
187muBrems: tables are built for mu+
188 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
189 Lambda tables from threshold to 100 TeV in 120 bins.
190 Parametrised model
191
192muPairProd: tables are built for mu+
193 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
194 Lambda tables from threshold to 100 TeV in 120 bins.
195 Parametrised model
196
197muIoni: tables are built for mu-
198 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
199 Lambda tables from threshold to 100 TeV in 120 bins.
200 Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below,
201 radiative corrections for E > 1 GeV
202 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
203
204muBrems: tables are built for mu-
205 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
206 Lambda tables from threshold to 100 TeV in 120 bins.
207 Parametrised model
208
209muPairProd: tables are built for mu-
210 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
211 Lambda tables from threshold to 100 TeV in 120 bins.
212 Parametrised model
213
214hIoni: tables are built for pi+
215 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
216 Lambda tables from threshold to 100 TeV in 120 bins.
217 Scaling relation is used from proton dE/dx and range.
218 Delta cross sections and sampling from BetheBloch model for scaled energy > 0.297504 MeV
219 Parametrisation from Bragg for protons below. NuclearStopping= 1
220 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
221
222msc: Model variant of multiple scattering for pi-
223 Lambda tables from 100 eV to 100 TeV in 120 bins.
224 LateralDisplacementFlag= 1 Skin= 0
225 Boundary/stepping algorithm is active with RangeFactor= 0.2 Step limit type 0
226
227hIoni: tables are built for pi-
228 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
229 Lambda tables from threshold to 100 TeV in 120 bins.
230 Scaling relation is used from proton dE/dx and range.
231 Delta cross sections and sampling from BetheBloch model for scaled energy > 0.297504 MeV
232 Parametrisation from Bragg for protons below. NuclearStopping= 1
233 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
234
235========= Table of registered couples ==============================
236
237Index : 0 used in the geometry : Yes recalculation needed : No
238 Material : G4_AIR
239 Range cuts : gamma 300 um e- 300 um e+ 300 um
240 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV
241 Region(s) which use this couple :
242 DefaultRegionForTheWorld
243
244Index : 1 used in the geometry : Yes recalculation needed : No
245 Material : G4_CESIUM_IODIDE
246 Range cuts : gamma 300 um e- 300 um e+ 300 um
247 Energy thresholds : gamma 18.1187 keV e- 288.607 keV e+ 278.143 keV
248 Region(s) which use this couple :
249 DefaultRegionForTheWorld
250
251Index : 2 used in the geometry : Yes recalculation needed : No
252 Material : G4_Al
253 Range cuts : gamma 300 um e- 300 um e+ 300 um
254 Energy thresholds : gamma 3.89739 keV e- 258.339 keV e+ 252.056 keV
255 Region(s) which use this couple :
256 DefaultRegionForTheWorld
257
258Index : 3 used in the geometry : Yes recalculation needed : No
259 Material : G4_AIR
260 Range cuts : gamma 10 um e- 10 um e+ 10 um
261 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV
262 Region(s) which use this couple :
263 VertexDetector
264
265Index : 4 used in the geometry : Yes recalculation needed : No
266 Material : G4_Al
267 Range cuts : gamma 10 um e- 10 um e+ 10 um
268 Energy thresholds : gamma 990 eV e- 33.8876 keV e+ 33.4729 keV
269 Region(s) which use this couple :
270 VertexDetector
271
272Index : 5 used in the geometry : Yes recalculation needed : No
273 Material : G4_Si
274 Range cuts : gamma 10 um e- 10 um e+ 10 um
275 Energy thresholds : gamma 990 eV e- 31.4747 keV e+ 31.0896 keV
276 Region(s) which use this couple :
277 VertexDetector
278
279Index : 6 used in the geometry : Yes recalculation needed : No
280 Material : G4_Fe
281 Range cuts : gamma 1 cm e- 1 cm e+ 1 cm
282 Energy thresholds : gamma 60.4019 keV e- 12.9518 MeV e+ 12.0297 MeV
283 Region(s) which use this couple :
284 MuonDetector
285
286Index : 7 used in the geometry : Yes recalculation needed : No
287 Material : G4_Al
288 Range cuts : gamma 1 cm e- 1 cm e+ 1 cm
289 Energy thresholds : gamma 20.4949 keV e- 4.54874 MeV e+ 4.2772 MeV
290 Region(s) which use this couple :
291 MuonDetector
292
293====================================================================
294
295### Run 0 start
296HistoManager: Histograms are booked and run has been started
297Start Run processing.
298Run terminated.
299Run Summary
300 Number of events processed : 100
301 User=4.19s Real=4.21s Sys=0s
302RunAction: End of run actions are started
303HistoManager: End of run actions are started
304=================================================================
305Number of events 100
306Average number of e- 1320
307Average number of gamma 701.3
308Average number of e+ 48.12
309Average number of steps 3980
310Edep 1x1 = 0.8336 +- 0.002518 res= 3.021 %
311Edep 3x3 = 0.9601 +- 0.001163 res= 1.212 %
312Edep 5x5 = 0.9807 +- 0.0007853 res= 0.8007 %
313=========== Mean values without trancating =====================
314Edep 1x1 = 0.8327 +- 0.002681 res= 3.219 %
315Edep 3x3 = 0.9596 +- 0.001241 res= 1.293 %
316Edep 5x5 = 0.9807 +- 0.0007853 res= 0.8007 %
317==================================================================
318
319
320<<<<<ACCEPTANCE>>>>> 100 events for Crystal Calorimeter
321Edep1x1: 0.8336 delEdep1x1= 0.006515 nrms= 2.191
322Erms1x1: 0.02518 delErms1x1= -0.00456 nrms= -1.533
323Edep3x3: 0.9601 delEdep3x3= 0.003087 nrms= 2.401
324Erms3x3: 0.01163 delErms3x3= -0.001226 nrms= -0.9531
325Edep5x5: 0.9807 delEdep5x5= 0.001428 nrms= 1.615
326Erms5x5: 0.007853 delErms5x5= -0.0009872 nrms= -1.117
327<<<<<END>>>>> IS ACCEPTED
328
329#
330G4 kernel has come to Quit state.
Note: See TracBrowser for help on using the repository browser.