(r4) GENIE < Main

Main Web>TWikiUsers>WilliamSeligman>GENIE (revision 4) (raw view)~~EditAttach~~
---+ Running !GENIE at Nevis

To run the [[http://www.genie-mc.org/][GENIE Neutrino Monte Carlo]] at Nevis, the simplest thing to do is to access the version maintained at [[https://web.fnal.gov/project/GENIE/SitePages/Home.aspx][Fermilab]]. I recommend against trying to install it on your laptop or in your own work area on the cluster; !GENIE requires quite a few external packages, is finicky about those packages' versions, and the complete build/installation process has a number of pitfalls. Even the supposedly pre-configured [[https://github.com/GENIEMC/lamp][lamp]] scripts, which are supposed to install all of !GENIE's dependencies for you, does not work reliably. 

---++ Setting up !GENIE

First, set up the Fermilab [[http://www.fnal.gov/docs/products/ups/ReferenceManual/][UPS]] environment:
<pre>
source /a/share/amsterdam/seligman/microboone/setup_uboone.sh
</pre>

You can look over the available versions of GENIE:
<pre>
ups list -aK+ genie
</pre>

As of 5-Apr-2017, some of the output of that command includes:
<pre>
"genie" "v2_12_4" "Linux64bit+2.6-2.12" "debug:e10:r6" "" 
"genie" "v2_12_4" "Linux64bit+2.6-2.12" "e10:prof:r6" "" 
"genie" "v2_12_2b" "Linux64bit+2.6-2.12" "e10:prof" "" 
"genie" "v2_12_2b" "Linux64bit+2.6-2.12" "debug:e10" "" 
</pre>


It appears that the latest version of !GENIE available at FNAL is v2_12_4. To set up the environment necessary to run the !GENIE programs for that version:
<pre>
setup genie v2_12_4 -q e10:prof:r6
</pre>
Note that you don't include the OS version in the <tt>setup</tt> command, but you do have to include the "qualifier" to select the options that were used to compile !GENIE. Just copy whatever is in the quotes and put it after the =-q= option. 

---++ Running !GENIE

The !GENIE manual is [[https://arxiv.org/pdf/1510.05494][here]] (PDF file). The practical examples start on page 59. Here's the first example, adapted for running at Nevis, and with an Ar<sup>40</sup> target. Note that the &nu;<sub>&mu;</sub>  histogram in the flux file is named "numu"; that name follows the comma after the .root file in the =-f= option. . 

<pre>
# Create a work area for transient output files
mkdir -p /a/data/westside/seligman/work
cd /a/data/westside/seligman/work

# Define variables to point at directories with flux and cross-section files.
export GENIE_FLUXES=/a/data/westside/uboone/flux
export GENIE_XSECS=/a/data/westside/uboone/genie/genie_xsec

# Generate 10,000 numu+Ar events with energies between 0 and 10 GeV
gevgen -n 10000 -p 14 -t 1000180400 -e 0,10 --run 100 \
   -f $GENIE_FLUXES/bnb/current/histograms/ub470_r200_bnb_fhc/numode_bnb_470m_r200.root,numu 
   --seed 2989819 --cross-sections $GENIE_XSECS/v2_12_0/NULL/DefaultPlusMECWithNC/data/gxspl-FNALsmall.xml
</pre>

Archeology: on 5-Apr-2017, I copied the flux files from =uboonegpvm01.fnal.gov:/uboone/data/flux=, and the !GENIE spline files from [[https://www.hepforge.org/archive/genie/data/]]. 

This example is straight from the manual: &pi;<sup>+</sup> on O<sup>16</sup>, with kinetic energy 200 !MeV:
<pre>
gevgen_hadron -n 10000 -p 211 -t 1000080160 -k 0.2 --seed 9839389
</pre>

---++ What's next?

As described in the !GENIE manual, to turn various models on/off you need to generate new spline files. The program for doing that is =gmkspl=, described on page 63 of the manual. You'll want to copy the directory =$GENIE/config/= to your work area, then copy and edit file =EventGeneratorListAssembler.xml=; that revised file will become an input list to =gmkspl=. An possible example (&nu;<sub>&mu;</sub> on Ar<sup>40</sup>, up to 10 !GeV):
<pre>
gmkspl -p 14 -t 1000180400 -e 10 \
   --event-generator-list my-edited-process-list.xml \
   --output-cross-sections my-splines.xml
</pre>

Note that =gmkspl= typically takes several hours to run.
Topic revision: r4 - 2017-06-21 - WilliamSeligman
Main
Webs
ATLAS
DOE
Main
TWiki
Veritas