From: paola sala (paola.sala@cern.ch)
Date: Thu Nov 16 2006 - 16:40:25 CET
Hi
You cannot do it with eventdat,
but you can use the EVENTBIN card instead,
with the "special" binning type 8
As described in the manual ( actually in the USRBIN description..)
what(1)= 8.0 : Special user-defined 3D binning. Two variables are
discontinuous (e.g. region number), the third one is continuous,
but not necessarily a space coordinate. See 12} for instructions
on how to write, compile and link the user routines.
Variable # type default Override routine
1st integer region number MUSRBR
2nd integer lattice cell number LUSRBL
3rd continuous pseudorapidity FUSRBV
Paola
On Thu, 2006-11-16 at 16:22 +0300, Konstantin Batkov wrote:
> 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 - 17:53:56 CET