Running Professor on Alpgen (December 2010)

Part 1: running alpgen with agile & rivet

First, some setups. Grab the rivet bootstrap script:

http://projects.hepforge.org/rivet/trac/wiki/GettingStarted

cd /scratch/rivet/
wget http://svn.hepforge.org/rivet/bootstrap/rivet-bootstrap
chmod +x rivet-bootstrap

Then:

./rivet-bootstrap

It probably won't work. Now, install gsl and a recent python:

ftp://ftp.gnu.org/gnu/gsl/ http://www.python.org/ftp/python/2.6.6/Python-2.6.6.tgz

Download to the rivet "downoads" directory, unpack to the build directory. Configure with:

./configure --prefix=`pwd`/../../local --exec-prefix=`pwd`/../../local make make install

Hack

export PATH=`pwd`/local/bin/:$PATH
export LD_LIBRARY_PATH=`pwd`/local/lib

./rivet-bootstrap --install-boost

./rivet-bootstrap --install-agile

As of 2012, there should only be one fix needed, and this is to handle the normalisation needed to combine the individual npX alpgen samples. If your rivet analysis does any kind of normalisation to cross section or number of events, you have to change this. I generally just turn it off and sort out the normalisation later.

This can be done eg with this patch rivet.patch

Install the generators

Get the rivet genser script:

http://svn.hepforge.org/agile/genser/agile-genser-bootstrap

This is my hacked version to build alpgen:

agile-genser-bootstrap

and set AGILE_GEN_PATH=`pwd`/genser

The main thing you need to check is that NMAXHEP is the same in agile and algpen. It can happen that agile defaults to 10000, but alpgen has 4000 hard-coded. You can just grep for 4000 to catch all instances, or this patch should work. It also fixes some other things, like automation of file names etc rather than having to type them in whenever you run.

alpgen.patch

Remake the packages after patching...

Running

That should be it. Assuming you have some unweighted alpgen events ready for the matching step, it should then be easy:

export PATH=`pwd`/local/bin/:$PATH

export LD_LIBRARY_PATH=`pwd`/local/lib

export AGILE_GEN_PATH=`pwd`/genser

agile-runmc AlpGenPythia6:423:2.1.3d --beams=TVT:1960 -n 100 -P params.txt -o hepmc.fifo > agile.log &
rivet --analysis=D0_2010_S8671338 --analysis=D0_2009_S8349509 --analysis=D0_2008_S7863608 --analysis=CDF_2008_NOTE_9351 hepmc.fif
o

TODO

Correct scale alteration in pythia - ktfac**2?

-- GavinHesketh - 2010-12-31