Project

General

Profile

How to Run LOAF » History » Version 5

Version 4 (Christopher Marshall, 02/17/2017 02:53 PM) → Version 5/14 (Christopher Marshall, 02/17/2017 03:02 PM)

h1. How to Run LOAF

* Check out the Fast MC as described on https://cdcvs.fnal.gov/redmine/projects/lbne-beamsim/wiki
* Go to the LOAF subdirectory directory
<pre>
$ cd trunk/LOAF
</pre>
* Build LOAF:
<pre>
$ gmake
</pre>
* Make inputs, or copy them from Dan's direcotry:
<pre>
cp /dune/data/users/lblpwg_tools/FastMC_Data/globes_files/cherdack/Globes_Inputs/* /dune/data/users/lblpwg_tools/FastMC_Data/globes_files/<username>/Globes_Inputs/
</pre>
* Make your output directory (if not running in batch mode, in which case this is done automatically):
<pre>
mkdir -p /dune/data/users/lblpwg_tools/FastMC_Data/globes_files/ljf26/Globes_Output/<configuration_name>/
</pre>
* Run LOAF.
<pre>
./Analysis_tools <configuration_name> <delta_CP> <IH or NH> <outputfile>
</pre>
The first argument is a configuration name, for example CDR_BASE_FullFD. Look in FastMC_GlobesParams.xml to see the available configurations, or add new ones. The configuration names appear between quotation marks on lines such as: <param_set name="Default">. The second argument is the true value of delta_CP in radians. The third argunment is the true mass hierarchy. The final argument is the name of a text output file that will be created.

* The output text file is a series of numbers. The identity of the numbers is given in:
/dune/data/users/lblpwg_tools/FastMC_Data/globes_files/<username>/Globes_Output/<configuration_name>/TTreeReadHeader.dat
That directory also has an image file comparing far detector energy spectra Markov data and Monte Carlo assuming delta_CP = 0 (this means the data and Monte Carlo will not match unless you have specified the true value of delta_CP to be zero.

* Run on the grid using

*
Create the inputs using split.C, in FullMCProcessing. The input files are MVASelect files, which should be merged with AnaTree eventually. The path is hard-coded into split.C. It needs to be run on the nue and numu near samples, and the commands are
<pre>
$ root -l -q -b split.C\(\"nue\"\)
$ root -l -q -b split.C\(\"numu\"\)
</pre>
The stupid forward slashes are if you want to do it from the command line, otherwise you can open root and use .L. Then run the scripts run_dis.sh and run_app.sh:
<pre>
$ . run_dis.sh
$ . run_app.sh
</pre>
This will take a long time. It will put stuff in LOAF_Inputs and LOAF_Inputs2. LOAF_Inputs2 needs to be copied to the Globes inputs area to run, which is /dune/data/users/lblpwg_tools/FastMC_Data/globes_files/$USER/Globes_Inputs