From: Konstantin Batkov (kbat@dec1.sinp.msu.ru)
Date: Thu Nov 16 2006 - 14:22:40 CET
I am trying to realize how to get the energy deposition in lattice cells
using the EVENTDAT estimator.
I am using an example of input file from the FLUKA Course.
This file is cited below (example.inp).
The output of EVENTDAT shows the energy deposition only for 17
non-lattice regions, but there is no any information about energy
deposition in lattice cells.
Is it possible to get this information using the EVENTDAT estimator?
Best regards,
Konstantin Batkov
****************** example001_energy *************************
**** Event-Data ****
FLUKA Course Exercise
DATE: 11/16/ 6, TIME: 13:47:32
No. of regions. 17 No. of distr. scored 1
Event # 1
Primary Weight 1. Primary Energy 3.5 GeV
Contributions to the energy deposition
(GeV not normalized to the weight):
1.67194021 0.675703585 0.139613718 0. 0.297916889
0. 0. 0.10776981 0.
0. 0. 0.607055843
Generalized scoring distribution # 208
from first to last region:
0.297916889 0. 0.00736697856 0.00672934763 0.795559287 0.
0.000258351269 0. 0. 0.000264887116 0.130817473 1.32444572
0.195700422
0.109640747 0.0242440086 0. 0.
Seeds after event # 1
*** 19CB2 0 0 0 0 0 181CD 3039
0 0***
************************ end of example001_energy ******************
*********************** example.inp ********************************
TITLE
FLUKA Course Exercise
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7..
DEFAULTS NEW-DEFA
BEAM -3.5 -0.082425 -1.7 0.0 0.0 1.0PROTON
BEAMPOS 0.0 0.0 0.1 0.0 0.0 0.0
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7..
GEOBEGIN COMBNAME
0 0 Cylindrical Target
**AAA*IIII_________+_________+_________+_________+_________+_________+
SPH BLK 0.0 0.0 0.0 10000.
* vacuum box
RPP VOI -1000. 1000. -1000. 1000. -1000. 1000.
* Lead target
ZCC TARG 0.0 0.0 5.0
* planes limiting the target
XYP ZTlow 0.0
XYP ZThigh 10.
**
* planes segmenting the target
XYP T1seg 1.
XYP T2seg 2.
* shield inner cylinder
ZCC ShInCy 0.0 0.0 10.
* Shield 1st seg
ZCC Sh1seg 0.0 0.0 50.
* Shield 2nd seg
ZCC Sh2seg 0.0 0.0 90.
* Shield 3nd seg
ZCC Sh3seg 0.0 0.0 140.
* Shield 4th seg
ZCC Sh4seg 0.0 0.0 180.
* Shield 5th seg
ZCC Sh5seg 0.0 0.0 220.
* Shield outer surface
ZCC ShOuC 0.0 0.0 260.
* shield bakwrd outer surface
XYP ShZminO -55.
* shield bakwrd inner surface
XYP ShZminI -50.
* shield forw. inner surface
XYP ShZmaxI 50.
* shield forw. outer surface
XYP ShZmaxO 60.
**
*Plane for lattice
XYP Zlat 20.
END
* Regions
*_AAA.....AA.....AA.....AA.....AA.....AA.....AA.....AA.....AA.....AA.....
* Black Hole
BLKHOLE 5 +BLK -VOI
* Void around
VAC 5 | +VOI -ShOuC | +VOI +ShOuC +ShZminO | +VOI +ShOuC -ShZmaxO
* Target seg. 1
TARGS1 5 +TARG -ZTlow +T1seg
* Target seg. 2
TARGS2 5 +TARG -T1seg +T2seg
* Target seg. 3
TARGS3 5 +TARG -T2seg +ZThigh
* Target replica
TARGRP 5 +TARG -ZThigh +Zlat
* Air inside shield
INAIR 5 +ShInCy -ShZminI +ShZmaxI
- ( +TARG -ZTlow +ZThigh ) - ( +TARG -ZThigh +Zlat )
*Shield back hole
SHBHOL 5 +TARG -ShZminO +ShZminI
*Shield front hole
SHFHOL 5 +TARG -ShZmaxI +ShZmaxO
*shield back
SHBACK 5 +ShOuC -TARG -ShZminO +ShZminI
*shield front
SHFRONT 5 +ShOuC -TARG -ShZmaxI +ShZmaxO
* shield 1st seg
SH1 5 -ShZminI +ShZmaxI +Sh1seg -ShInCy
* shield 2nd seg
SH2 5 -ShZminI +ShZmaxI +Sh2seg -Sh1seg
* shield 3rd seg
SH3 5 -ShZminI +ShZmaxI +Sh3seg -Sh2seg
* shield 4th seg
SH4 5 -ShZminI +ShZmaxI +Sh4seg -Sh3seg
* shield 5th seg
SH5 5 -ShZminI +ShZmaxI +Sh5seg -Sh4seg
* shield 6th seg
SH6 5 -ShZminI +ShZmaxI +ShOuC -Sh5seg
END
* define region TARGRP as a replica, with name Target2, associated
* to roto-traslation number 1
LATTICE TARGRP Target2 Rot# 1
*AAAAAAAxx_________+_________+_________+_________+_________+_________+
GEOEND
*** switch on to debug this geometry
*GEOEND 250.0 0.0 70.0 0.0 0.0 -60.0DEBUG
*GEOEND 2501. 1. 1301. &
*** switch on to plot this geometry
PLOTGEOM 1.0 -2000.0 1.
paw
0.0 -300.0 -100.0 0.0 300.0 100.0
+0.0 +0.0 1.0 +0.0 +1.0 +0.0
+1.0 +1.0 +0.0
*
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
MATERIAL 18.0 39.948 1.40 ARGON
MATERIAL 19.0 39.0983 0.862 POTASSIU
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
MATERIAL 1.00 WATER
COMPOUND 2.0 HYDROGEN 1.0 OXYGEN WATER
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
MATERIAL 0.001225 AIR
COMPOUND -.9256E-03 NITROGEN-.2837E-03 OXYGEN-.01572E-3 ARGON AIR
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
MATERIAL 2.42 CONCRETE
COMPOUND -0.01 HYDROGEN -0.001 CARBON -0.5290 OXYGEN CONCRETE
COMPOUND -0.016 SODIUM -0.002 MAGNESIU -0.034 ALUMINUM CONCRETE
COMPOUND -0.3370 SILICON -0.013 POTASSIU -0.044 CALCIUM CONCRETE
COMPOUND -0.014 IRON CONCRETE
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7..
ASSIGNMA BLCKHOLE BLKHOLE
ASSIGNMA CARBON TARGS1
ASSIGNMA ALUMINUM TARGS2
ASSIGNMA LEAD TARGS3
ASSIGNMA VACUUM VAC
ASSIGNMA AIR INAIR
ASSIGNMAT AIR SHBHOL
ASSIGNMAT AIR SHFHOL
ASSIGNMAT CONCRETE SHBACK SH6 1.0
*ASSIGNMAT AIR TARGRP
* bias: set 1 everywhere
BIASING 1.0 1.0 1.0 1.0 50.0 PRINT
* now bias in steps of 2
BIASING 1.0 1.0 2.0 SH2 SH2
BIASING 1.0 1.0 4.0 SH3 SH3
BIASING 1.0 1.0 8.0 SH4 SH4
BIASING 1.0 1.0 16.0 SH5 SH5
BIASING 1.0 1.0 32.0 SH6 SH6
*AAAAAAAxx_________+_________+_________+_________+_________+_________+
*define the transformation for lattice cell Target2
* rot.numb. theta phi dx dy dz
ROT-DEFI 1.0 0.0 0.0 0.0 0.0 -10.0
SCORE ENERGY
EVENTDAT 75.0 0.0 0.0 0.0 0.0 0.0 energy
RANDOMIZ 1.0
START 10.0 0.0
STOP
**************************** end of example.inp **********************
This archive was generated by hypermail 2.1.6 : Thu Nov 16 2006 - 15:48:41 CET