source: trunk/examples/advanced/microbeam/microbeam.out@ 1197

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

update

File size: 18.2 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
9Visualization Manager instantiating...
10Visualization Manager initialising...
11Registering graphics systems...
12
13You have successfully registered the following graphics systems.
14Current available graphics systems are:
15 ASCIITree (ATree)
16 DAWNFILE (DAWNFILE)
17 G4HepRep (HepRepXML)
18 G4HepRepFile (HepRepFile)
19 RayTracer (RayTracer)
20 VRML1FILE (VRML1FILE)
21 VRML2FILE (VRML2FILE)
22 OpenGLImmediateX (OGLIX)
23 OpenGLStoredX (OGLSX)
24
25Registering model factories...
26
27You have successfully registered the following model factories.
28Registered model factories:
29 generic
30 drawByCharge
31 drawByParticleID
32 drawByOriginVolume
33 drawByAttribute
34
35Registered filter factories:
36 chargeFilter
37 particleFilter
38 originVolumeFilter
39 attributeFilter
40
41
42
43***** Table : Nb of materials = 16 *****
44
45 Material: Vacuum density: 0.000 mg/cm3 RadL: 204727576.737 pc Imean: 21.800 eV temperature: 273.15 K pressure: 1.00 atm
46 ---> Element: Vacuum ( ) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
47
48 Material: H2O density: 1.000 g/cm3 RadL: 36.092 cm Imean: 70.893 eV
49 ---> Element: Hydrogen (H) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 11.21 % ElmAbundance 66.67 %
50 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 88.79 % ElmAbundance 33.33 %
51
52 Material: Air density: 1.290 mg/cm3 RadL: 285.161 m Imean: 85.684 eV temperature: 293.16 K pressure: 1.00 atm
53 ---> Element: Nitrogen (N) Z = 7.0 N = 14.0 A = 14.01 g/mole ElmMassFraction: 70.00 % ElmAbundance 72.71 %
54 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 30.00 % ElmAbundance 27.29 %
55
56 Material: LPAir density: 0.000 mg/cm3 RadL: 56273252.573 km Imean: 87.308 eV temperature: 293.16 K pressure: 1.00 atm
57 ---> Element: Nitrogen (N) Z = 7.0 N = 14.0 A = 14.01 g/mole ElmMassFraction: 71.50 % ElmAbundance 75.57 %
58 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 25.00 % ElmAbundance 23.14 %
59 ---> Element: Argon (Ar) Z = 18.0 N = 39.9 A = 39.95 g/mole ElmMassFraction: 3.50 % ElmAbundance 1.30 %
60
61 Material: Pl density: 21.400 g/cm3 RadL: 3.058 mm Imean: 787.800 eV
62 ---> Element: Pl ( ) Z = 78.0 N = 195.1 A = 195.09 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 %
63
64 Material: Butane density: 0.026 mg/cm3 RadL: 17.729 km Imean: 53.612 eV temperature: 293.16 K pressure: 0.01 atm
65 ---> Element: Carbon (C) Z = 6.0 N = 12.0 A = 12.01 g/mole ElmMassFraction: 82.63 % ElmAbundance 28.57 %
66 ---> Element: Hydrogen (H) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 17.37 % ElmAbundance 71.43 %
67
68 Material: Polyprop density: 900.000 mg/cm3 RadL: 49.764 cm Imean: 56.713 eV
69 ---> Element: Carbon (C) Z = 6.0 N = 12.0 A = 12.01 g/mole ElmMassFraction: 85.60 % ElmAbundance 33.33 %
70 ---> Element: Hydrogen (H) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 14.40 % ElmAbundance 66.67 %
71
72 Material: Si3N4 density: 3.440 g/cm3 RadL: 7.644 cm Imean: 128.542 eV
73 ---> Element: Silicon (Si) Z = 14.0 N = 28.1 A = 28.09 g/mole ElmMassFraction: 60.06 % ElmAbundance 42.86 %
74 ---> Element: Nitrogen (N) Z = 7.0 N = 14.0 A = 14.01 g/mole ElmMassFraction: 39.94 % ElmAbundance 57.14 %
75
76 Material: SiO2 density: 2.500 g/cm3 RadL: 10.819 cm Imean: 126.007 eV
77 ---> Element: Silicon (Si) Z = 14.0 N = 28.1 A = 28.09 g/mole ElmMassFraction: 46.74 % ElmAbundance 33.33 %
78 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 53.26 % ElmAbundance 66.67 %
79
80 Material: Laiton density: 8.500 g/cm3 RadL: 1.487 cm Imean: 325.993 eV
81 ---> Element: Cuivre (Cu) Z = 29.0 N = 63.5 A = 63.55 g/mole ElmMassFraction: 49.28 % ElmAbundance 50.00 %
82 ---> Element: Zinc (Zn) Z = 30.0 N = 65.4 A = 65.41 g/mole ElmMassFraction: 50.72 % ElmAbundance 50.00 %
83
84 Material: Cytoplasm1 density: 1.000 g/cm3 RadL: 48.413 cm Imean: 31.293 eV
85 ---> Element: Hydrogen (H) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 59.60 % ElmAbundance 95.53 %
86 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 24.24 % ElmAbundance 2.45 %
87 ---> Element: Carbon (C) Z = 6.0 N = 12.0 A = 12.01 g/mole ElmMassFraction: 11.11 % ElmAbundance 1.50 %
88 ---> Element: Nitrogen (N) Z = 7.0 N = 14.0 A = 14.01 g/mole ElmMassFraction: 4.04 % ElmAbundance 0.47 %
89 ---> Element: Phosphorus (P) Z = 15.0 N = 31.0 A = 30.97 g/mole ElmMassFraction: 1.01 % ElmAbundance 0.05 %
90
91 Material: Cytoplasm2 density: 10.000 g/cm3 RadL: 3.619 cm Imean: 71.338 eV
92 ---> Element: Hydrogen (H) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 10.64 % ElmAbundance 64.80 %
93 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 74.50 % ElmAbundance 28.64 %
94 ---> Element: Carbon (C) Z = 6.0 N = 12.0 A = 12.01 g/mole ElmMassFraction: 9.04 % ElmAbundance 4.63 %
95 ---> Element: Nitrogen (N) Z = 7.0 N = 14.0 A = 14.01 g/mole ElmMassFraction: 3.21 % ElmAbundance 1.41 %
96 ---> Element: Phosphorus (P) Z = 15.0 N = 31.0 A = 30.97 g/mole ElmMassFraction: 2.61 % ElmAbundance 0.52 %
97
98 Material: Cytoplasm3 density: 1.000 g/cm3 RadL: 48.413 cm Imean: 31.293 eV
99 ---> Element: Hydrogen (H) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 59.60 % ElmAbundance 95.53 %
100 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 24.24 % ElmAbundance 2.45 %
101 ---> Element: Carbon (C) Z = 6.0 N = 12.0 A = 12.01 g/mole ElmMassFraction: 11.11 % ElmAbundance 1.50 %
102 ---> Element: Nitrogen (N) Z = 7.0 N = 14.0 A = 14.01 g/mole ElmMassFraction: 4.04 % ElmAbundance 0.47 %
103 ---> Element: Phosphorus (P) Z = 15.0 N = 31.0 A = 30.97 g/mole ElmMassFraction: 1.01 % ElmAbundance 0.05 %
104
105 Material: Nucleus1 density: 1.000 g/cm3 RadL: 36.185 cm Imean: 71.338 eV
106 ---> Element: Hydrogen (H) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 10.64 % ElmAbundance 64.80 %
107 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 74.50 % ElmAbundance 28.64 %
108 ---> Element: Carbon (C) Z = 6.0 N = 12.0 A = 12.01 g/mole ElmMassFraction: 9.04 % ElmAbundance 4.63 %
109 ---> Element: Nitrogen (N) Z = 7.0 N = 14.0 A = 14.01 g/mole ElmMassFraction: 3.21 % ElmAbundance 1.41 %
110 ---> Element: Phosphorus (P) Z = 15.0 N = 31.0 A = 30.97 g/mole ElmMassFraction: 2.61 % ElmAbundance 0.52 %
111
112 Material: Nucleus2 density: 1.100 g/cm3 RadL: 32.896 cm Imean: 71.338 eV
113 ---> Element: Hydrogen (H) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 10.64 % ElmAbundance 64.80 %
114 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 74.50 % ElmAbundance 28.64 %
115 ---> Element: Carbon (C) Z = 6.0 N = 12.0 A = 12.01 g/mole ElmMassFraction: 9.04 % ElmAbundance 4.63 %
116 ---> Element: Nitrogen (N) Z = 7.0 N = 14.0 A = 14.01 g/mole ElmMassFraction: 3.21 % ElmAbundance 1.41 %
117 ---> Element: Phosphorus (P) Z = 15.0 N = 31.0 A = 30.97 g/mole ElmMassFraction: 2.61 % ElmAbundance 0.52 %
118
119 Material: Nucleus3 density: 1.000 g/cm3 RadL: 36.185 cm Imean: 71.338 eV
120 ---> Element: Hydrogen (H) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 10.64 % ElmAbundance 64.80 %
121 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 74.50 % ElmAbundance 28.64 %
122 ---> Element: Carbon (C) Z = 6.0 N = 12.0 A = 12.01 g/mole ElmMassFraction: 9.04 % ElmAbundance 4.63 %
123 ---> Element: Nitrogen (N) Z = 7.0 N = 14.0 A = 14.01 g/mole ElmMassFraction: 3.21 % ElmAbundance 1.41 %
124 ---> Element: Phosphorus (P) Z = 15.0 N = 31.0 A = 30.97 g/mole ElmMassFraction: 2.61 % ElmAbundance 0.52 %
125
126
127 ==========> The phantom contains 53480 voxels
128
129Thank you for using G4BinaryCascade.
130MicrobeamPhysicsList::SetCuts:CutLength : 10 nm
131
132msc: Model variant of multiple scattering for e-
133 Lambda tables from 100 eV to 100 TeV in 120 bins.
134 LateralDisplacementFlag= 1 Skin= 0
135 Boundary/stepping algorithm is active with RangeFactor= 0.02 Step limit type 1
136G4AugerData for Element no. 6 are loaded
137G4AugerData for Element no. 7 are loaded
138G4AugerData for Element no. 8 are loaded
139G4AugerData for Element no. 14 are loaded
140G4AugerData for Element no. 15 are loaded
141G4AugerData for Element no. 18 are loaded
142G4AugerData for Element no. 29 are loaded
143G4AugerData for Element no. 30 are loaded
144G4AugerData for Element no. 78 are loaded
145AugerTransitionTable complete
146
147IONI: Total cross sections from EEDL database.
148 Gamma energy sampled from a parametrised formula.
149 Implementation of the continuous dE/dx part.
150 At present it can be used for electrons in the energy range [250eV,100GeV].
151 The process must work with G4LowEnergyBremsstrahlung.
152
153LowEnBrem: Total cross sections from EEDL database.
154 Gamma energy sampled from a parameterised formula.
155 Implementation of the continuous dE/dx part.
156 At present it can be used for electrons in the energy range [250eV,100GeV].
157 The process must work with G4LowEnergyIonisation.
158
159eIoni: tables are built for e+
160 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
161 Lambda tables from threshold to 100 TeV in 120 bins.
162 Delta cross sections and sampling from MollerBhabha model
163 Good description from 1 KeV to 100 GeV.
164 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1
165
166eBrem: tables are built for e+
167 dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
168 Lambda tables from threshold to 100 TeV in 120 bins.
169 Total cross sections and sampling from StandBrem model (based on the EEDL data library)
170 Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. LPM flag 1
171
172annihil: Sampling according eplus2gg model
173 tables are built for e+
174 Lambda tables from 100 eV to 100 TeV in 120 bins.
175
176msc: Model variant of multiple scattering for proton
177 Lambda tables from 100 eV to 100 TeV in 120 bins.
178 LateralDisplacementFlag= 1 Skin= 0
179 Boundary/stepping algorithm is active with RangeFactor= 0.2 Step limit type 0
180
181hLowEIoni: Knock-on electron cross sections .
182 Good description above the mean excitation energy.
183 Delta ray energy sampled from differential Xsection.
184 PhysicsTables from 10 eV to 0.1 TeV in 360 bins.
185 Electronic stopping power model is ICRU_R49He
186 from 1 keV to 2.5 MeV .
187
188 Parametrisation model for antiprotons is ICRU_R49p
189 from 25 keV to 2 MeV .
190 Parametrization of the Barkas effect is switched on.
191 Nuclear stopping power model is ICRU_R49
192
193msc: Model variant of multiple scattering for GenericIon
194 LateralDisplacementFlag= 0 Skin= 0
195 Boundary/stepping algorithm is active with RangeFactor= 0.2 Step limit type 1
196
197hLowEIoni: Knock-on electron cross sections .
198 Good description above the mean excitation energy.
199 Delta ray energy sampled from differential Xsection.
200 PhysicsTables from 10 eV to 0.1 TeV in 360 bins.
201 Electronic stopping power model is ICRU_R49p
202 from 1 keV to 2 MeV .
203
204 Parametrisation model for antiprotons is ICRU_R49He
205 from 25 keV to 2 MeV .
206 Parametrization of the Barkas effect is switched on.
207 Nuclear stopping power model is ICRU_R49
208
209msc: Model variant of multiple scattering for pi-
210 Lambda tables from 100 eV to 100 TeV in 120 bins.
211 LateralDisplacementFlag= 1 Skin= 0
212 Boundary/stepping algorithm is active with RangeFactor= 0.2 Step limit type 0
213
214========= Table of registered couples ==============================
215
216Index : 0 used in the geometry : Yes recalculation needed : No
217 Material : Vacuum
218 Range cuts : gamma 10 nm e- 10 nm e+ 10 nm
219 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV
220 Region(s) which use this couple :
221 DefaultRegionForTheWorld
222
223Index : 1 used in the geometry : Yes recalculation needed : No
224 Material : Pl
225 Range cuts : gamma 10 nm e- 10 nm e+ 10 nm
226 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV
227 Region(s) which use this couple :
228 DefaultRegionForTheWorld
229
230Index : 2 used in the geometry : Yes recalculation needed : No
231 Material : Butane
232 Range cuts : gamma 10 nm e- 10 nm e+ 10 nm
233 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV
234 Region(s) which use this couple :
235 DefaultRegionForTheWorld
236
237Index : 3 used in the geometry : Yes recalculation needed : No
238 Material : Laiton
239 Range cuts : gamma 10 nm e- 10 nm e+ 10 nm
240 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV
241 Region(s) which use this couple :
242 DefaultRegionForTheWorld
243
244Index : 4 used in the geometry : Yes recalculation needed : No
245 Material : Si3N4
246 Range cuts : gamma 10 nm e- 10 nm e+ 10 nm
247 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV
248 Region(s) which use this couple :
249 DefaultRegionForTheWorld
250
251Index : 5 used in the geometry : Yes recalculation needed : No
252 Material : Air
253 Range cuts : gamma 10 nm e- 10 nm e+ 10 nm
254 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV
255 Region(s) which use this couple :
256 DefaultRegionForTheWorld
257
258Index : 6 used in the geometry : Yes recalculation needed : No
259 Material : Polyprop
260 Range cuts : gamma 10 nm e- 10 nm e+ 10 nm
261 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV
262 Region(s) which use this couple :
263 DefaultRegionForTheWorld
264
265Index : 7 used in the geometry : Yes recalculation needed : No
266 Material : H2O
267 Range cuts : gamma 10 nm e- 10 nm e+ 10 nm
268 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV
269 Region(s) which use this couple :
270 DefaultRegionForTheWorld
271
272Index : 8 used in the geometry : Yes recalculation needed : No
273 Material : Cytoplasm1
274 Range cuts : gamma 10 nm e- 10 nm e+ 10 nm
275 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV
276 Region(s) which use this couple :
277 DefaultRegionForTheWorld
278
279Index : 9 used in the geometry : Yes recalculation needed : No
280 Material : Nucleus1
281 Range cuts : gamma 10 nm e- 10 nm e+ 10 nm
282 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV
283 Region(s) which use this couple :
284 DefaultRegionForTheWorld
285
286Index : 10 used in the geometry : Yes recalculation needed : No
287 Material : Nucleus2
288 Range cuts : gamma 10 nm e- 10 nm e+ 10 nm
289 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV
290 Region(s) which use this couple :
291 DefaultRegionForTheWorld
292
293Index : 11 used in the geometry : Yes recalculation needed : No
294 Material : Cytoplasm2
295 Range cuts : gamma 10 nm e- 10 nm e+ 10 nm
296 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV
297 Region(s) which use this couple :
298 DefaultRegionForTheWorld
299
300Index : 12 used in the geometry : Yes recalculation needed : No
301 Material : SiO2
302 Range cuts : gamma 10 nm e- 10 nm e+ 10 nm
303 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV
304 Region(s) which use this couple :
305 DefaultRegionForTheWorld
306
307====================================================================
308
309-> Event # 1 generated
310 ===> Sorry, the incident alpha particle has missed the targeted cell !
311
312G4VisManager: Using G4TrajectoryDrawByCharge as default trajectory model.
313See commands in /vis/modeling/trajectories/ for other options.
314Trajectory drawing configuration will be based on imode value of 50
315WARNING: G4VisManager::IsValidView(): Attempt to draw when no graphics system
316 has been instantiated. Use "/vis/open" or "/vis/sceneHandler/create".
317 Alternatively, to avoid this message, suppress instantiation of vis
318 manager (G4VisExecutive), possibly by setting G4VIS_NONE, and ensure
319 drawing code is executed only if G4VVisManager::GetConcreteInstance()
320 is non-zero.
321-> Event # 2 generated
322 ===> The incident alpha particle has reached the targeted cell :
323 -----> total absorbed dose within Nucleus is (Gy) = 0.40389574
324 -----> total absorbed dose within Cytoplasm is (Gy) = 0.04089877
325
326-> Event # 3 generated
327 ===> The incident alpha particle has reached the targeted cell :
328 -----> total absorbed dose within Nucleus is (Gy) = 0.30514622
329 -----> total absorbed dose within Cytoplasm is (Gy) = 0.090157568
330
331-> Event # 4 generated
332 ===> Sorry, the incident alpha particle has missed the targeted cell !
333
334-> Event # 5 generated
335 ===> Sorry, the incident alpha particle has missed the targeted cell !
336
337-> Event # 6 generated
338 ===> Sorry, the incident alpha particle has missed the targeted cell !
339
340-> Event # 7 generated
341 ===> The incident alpha particle has reached the targeted cell :
342 -----> total absorbed dose within Nucleus is (Gy) = 0.37321898
343 -----> total absorbed dose within Cytoplasm is (Gy) = 0.070893854
344
345-> Event # 8 generated
346 ===> Sorry, the incident alpha particle has missed the targeted cell !
347
348-> Event # 9 generated
349 ===> Sorry, the incident alpha particle has missed the targeted cell !
350
351-> Event # 10 generated
352 ===> Sorry, the incident alpha particle has missed the targeted cell !
353
354ERROR: G4VisCommandsViewerUpdate::SetNewValue: no current viewer.
355-> Total number of particles detected by the gas detector : 3
356
357Graphics systems deleted.
358Visualization Manager deleting...
Note: See TracBrowser for help on using the repository browser.