Hi Leena,
> I did not find GLOBAL with SDUM=OLDFLAIR option.
>
This is obsolete FLUKA syntax, just removing the !_at_what record should fix
that problem.
> I am uncertain where to find the final result (a reading in Gy/g along
> with its statistical uncertainty).
>
FLUKA reports DOSE in GeV/g. You need to convert into Gy yourself as
explained in note 5 on page 282 of the manual. If you need to normalize Gy
by the region mass, you also do it yourself.
If you merge and plot data with Flair, the data and corresponding
uncertainties are in the last two columns in sim2_21_plot.dat. However,
this file contains only data corresponding to the image plotted with Flair.
To access the whole USRBIN data set, you need to process it:
First, call the $FLUPRO/flutil/usbsuw tool to merge results from different
run cycles and calculate statistical uncertainties. You can run it either
explicitly, e.g.
(ls -1 sim4*.22; echo; echo sim4_22bin) | $FLUPRO/flutil/usbsuw
or with a wrapper script $FLUPRO/flutil/binsum.sh. It needs two arguments:
input file name (without extension) and the estimator unit number:
$FLUPRO/flutil/binsum.sh sim4 22
This generates a binary file sim4_22bin.bnn with two matrices: with results
and corresponding percentage errors.
You can read it directly or first convert into text format. To do the
latter, use the $FUPRO/flutil/usbrea tool, e.g.
(echo sim4_22bin.bnn; echo sim4_22bin.txt) | $FLUPRO/flutil/usbrea
I recommend you to make a test run with a very small number of USRBIN bins
in order to better understand the file format.
> Additionally, I noticed that changing the shape of the beam leads to
> different results. I'm not sure which configuration is best for a
> monoenergetic X-ray beam.
>
What's best depends on your figure of merit, what you would like to
achieve. If your goal is to burn a hole in your target, the pencil beam
looks best, if you want to dissipate the heat load, increasing the beam
footprint area should help.
> One more thing, is there any way to confirm that the mass/volume of the
> region of interest (C2TARGET) is correctly calculated by FLUKA?
>
To my knowledge, FLUKA does not calculate the region masses/volumes itself,
but it is possible to indicate volumes in the input file or in the
post-processing tools.
Konstantin
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Tue Jul 30 2024 - 09:04:40 CEST