(r1) GENIE < Main

Main Web>TWikiUsers>WilliamSeligman>GENIE (revision 1) (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 tthe latest version of !GENIE available at FNAL is v2_12_2b. To set up the environment necessary to run the !GENIE programs for that version:
<pre>
setup genie v2_12_2b -q e10:prof
</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. Typically you'll just copy the C++ compiler version ("eNN") and you'll want the "prof" version for the fastest execution. 

---++ 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 argon target. Note that the numu histogram in this flux file is named "numu". 

<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+ on O, 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 (numu on Ar, 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: r1 - 2017-04-05 - WilliamSeligman
Main
Webs
ATLAS
DOE
Main
TWiki
Veritas