Re: Can't score the energy deposition blackhole region by USRBIN?

From: Paola Sala <>
Date: Fri, 13 May 2011 15:27:54 +0200 (CEST)

Dear Hantao Jing
Why do you put blackhole? You should NOT do it to estimate energy deposition.
To score energy deposition, just put the material corresponding to your
real geometry, you will get the energy deposited there.
By the way, the energy "seen" in the blackhole is NOT a correct estimate
of energy deposition, since all the effects of nuclear interactions, non
ionizing particles, escaping particles and so on and so fort are not
accounted for.
Be also careful : you are asking DOSE ( generalized particle 228) in a
USRBIN by region. The results will NOT be GeV/g because the usrbin by
region does NOT normalize to the voulume. (..besides the difficulty in
assigning a density to blackhole..)

> Dear FLUKA experts,
> I want score the energy deposition along the particle outgoing
> direction from a target, so I setup several
> blackhole region to score the energy disposition. From the USRBIN output
> results, the reulsts is zero:
> 1
> Region binning n. 2 "magnetcD " , generalized particle n. 228
> 4 bins corresponding to the region sets:
> from region 37 to region 40 in step of 1 regions, or
> from region 0 to region 0 in step of 1 regions, or
> from region 0 to region 0 in step of 1 regions
> Data follow in an array A(ir), format (1(5x,1p,10(1x,e11.4)))
> accurate deposition along the tracks requested
> this is a track-length binning
> 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
> Percentage errors follow in an array A(ir), format
> (1(5x,1p,10(1x,e11.4)))
> 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
> But, from the xxx.out file, we can find the black regions have the
> output results:
> 1Region # name volume ALL-PART Star Density BEAMPART Star
> Density ENERGY Density EM-ENRGY Density
> in cubic cm Stars/cm**3
> Stars/cm**3 GeV/cm**3 GeV/cm**3
> /one beam particle /one beam
> particle /one beam particle /one beam particle
> 37 Q30DC 1.000000000E+00 0.000000000E+00
> 0.000000000E+00 8.835180331E-06 0.000000000E+00
> 38 Q29DC 1.000000000E+00 0.000000000E+00
> 0.000000000E+00 3.756331285E-06 0.000000000E+00
> 39 Q28AC 1.000000000E+00 0.000000000E+00
> 0.000000000E+00 1.983465558E-06 0.000000000E+00
> 40 RTBHC 1.000000000E+00 0.000000000E+00
> 0.000000000E+00 5.267937506E-08 0.000000000E+00
> In fact, I run many jobs at the same time on servers. If I can't use
> USRBIN to score the blackhole region energy deposition, how to merge all
> output results and obtain the statistical errors quickly?
> Thank you very much!
> Hantao Jing
> 2011-5-13

Paola Sala
INFN Milano
tel. Milano +39-0250317374
tel. CERN +41-227679148
Received on Fri May 13 2011 - 15:59:31 CEST

This archive was generated by hypermail 2.2.0 : Fri May 13 2011 - 15:59:32 CEST