FLUKA

This section describes how to use the FLUKA-related tools.

fluka2root

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 averaged by $FLUPRO/flutil/usbsuw into a TH3F histogram.
• usxsuw2root: converts USRBDX results into a TH2F histogram. (See the notes for usbsuw2root above.)
• usrsuw2root: converts RESNUCLEI results into a TH2F histogram and TGraphError. (See the notes for usbsuw2root above.) Known bug: isomer data might not be saved in the ROOT file.
• ustsuw2root: converts USRTRACK and USRCOLL results into a TH1F histogram. (Again, see the notes for usbsuw2root.)
• eventdat2root: converts EVENTDAT results into a TTree object.
• detsuw2root: converts DETECT output into a TH1F histogram.