-
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, averaged by
$FLUPRO/flutil/usbsuw
, into a TH3F histogram. -
usxsuw2root: converts USRBDX results, averaged by
$FLUPRO/flutil/usxsuw
, into a TH2F histogram. -
usrsuw2root: converts RESNUCLEI results, averaged by
$FLUPRO/flutil/usrsuw
, into a TH2F histogram and TGraphError. Known bug: isomer data might not be saved in the ROOT file. -
ustsuw2root: converts USRTRACK and USRCOLL results, averaged by
$FLUPRO/flutil/ustsuw
, into a TH1F histogram. - eventdat2root: converts EVENTDAT results into a TTree object.
-
detsuw2root: converts DETECT output, averaged by
$FLUPRO/flutil/detsuw
, into a TH1F histogram.