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

Last change on this file since 812 was 807, checked in by garnier, 16 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.