-
Notifications
You must be signed in to change notification settings - Fork 17
FLUKA
This section describes how to use the FLUKA-related tools.
Copy $FLUPRO/exfixed.inp
into an empty folder and run:
$FLUPRO/flutil/rfluka exfixed.inp
As expected, this generates output files from 5 run cycles. We can use the fluka2root
script to call the standard FLUKA tools to average the results, and then convert them into a single ROOT file:
fluka2root exfixed.inp
This produces the exfixed.root
file with TH1, TH2 and TH3 histograms for the USRBIN, USRBDX and USRTRACK estimators defined in exfixed.inp
:
lsroot exfixed.root
KEY: TH3F piFluBin;1 PIONS#pm
KEY: TH3F Edeposit;1 ENERGY
KEY: TH2F piFluenUD;1 PIONS#pm fluence #diamond reg 3 #rightarrow 4 #diamond 400 cm^{2}
KEY: TH2F piCurrUD;1 PIONS#pm current #diamond reg 3 #rightarrow 4 #diamond 400 cm^{2}
KEY: TH1F piFluenU;1 PIONS#pm #diamond reg 3 #diamond 1000 cm^{3} #diamond 0.001 < E < 50 GeV
KEY: TH1F piFluenD;1 PIONS#pm #diamond reg 4 #diamond 1000 cm^{3} #diamond 0.001 < E < 50 GeV
The histograms can be plotted and analysed with standard ROOT tools. For example, this script demonstrates how to plot these histograms. Note that it creates a PDF file with multiple pages of plots.
If you run several input files with different random seeds, fluka2root
can average over all data by specifying the corresponding input files in the argument list:
fluka2root files*.inp
Note that the free format FLUKA input is not supported.
Internally, fluka2root
runs the following convertes, which can technically be called individually, though it's not recommended:
-
usbsuw2root: converts the USRBIN results into a TH3F histogram. Note that this tool does not directly convert the files produced by the USRBIN card. These files must first be averaged using $FLUPRO/flutil/usbsuw. Once avaraged, the resulting file can be converted to ROOT format using
usbsuw2root
. When using the general the fluka2root converter, this $FLUPRO/flutil/usbsuw call is performed automatically. - usxsuw2root script to convert the USRBDX results into a TH2F histogram. + see the comments for the previous item.
-
usrsuw2root script to convertert the
RESNUCLEI results into a TH2F histogram and TGraphError + see the comments for
usbsuw2root
above. Known bug: isomer data might be not saved into the ROOT file. -
ustsuw2root script to convert the USRTRACK and USRCOLL results into a TH1F histogram. + see the comments for
usbsuw2root
above. - eventdat2root script to convert the EVENTDAT results into a TTree object.
- detsuw2root script to convert DETECT output into TH1F histogram.