GavinHeskethSandbox < Sandbox

Tags: view all tags

---+ *Running Professor on Alpgen (December 2010)*

---++ Part 1: running alpgen with agile & rivet

First, some setups. Grab the rivet bootstrap script:

<a target="_blank" href="http://projects.hepforge.org/rivet/trac/wiki/GettingStarted" title="Getting Started with Rivet">http://projects.hepforge.org/rivet/trac/wiki/GettingStarted</a>

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

Then:

./rivet-bootstrap

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

<a target="_blank" href="ftp://ftp.gnu.org/gnu/gsl/" title="gsl">ftp://ftp.gnu.org/gnu/gsl/</a> <a target="_blank" href="http://www.python.org/ftp/python/2.6.6/Python-2.6.6.tgz" title="Python 2.6.6">http://www.python.org/ftp/python/2.6.6/Python-2.6.6.tgz</a>

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<br /> 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 [[%ATTACHURL%/rivet.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:

[[%ATTACHURL%/agile-genser-bootstrap][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.

 [[%ATTACHURL%/alpgen.patch][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 &gt; agile.log &<br />rivet --analysis=D0_2010_S8671338 --analysis=D0_2009_S8349509 --analysis=D0_2008_S7863608 --analysis=CDF_2008_NOTE_9351 hepmc.fif<br />o
---++ TODO

Correct scale alteration in pythia - ktfac**2?

-- Main.GavinHesketh - 2010-12-31

Attachments

Topic attachments
I	Attachment	History	Action	Size	Date	Who	Comment
patch	HepMC.patch	r1	manage	0.7 K	2010-12-31 - 15:41	GavinHesketh	HepMC patch
ext	agile-genser-bootstrap	r1	manage	15.0 K	2010-12-31 - 16:58	GavinHesketh	agile genser script
patch	agile.patch	r2 r1	manage	2.4 K	2010-12-31 - 15:40	GavinHesketh	Agile patch
patch	alpgen.patch	r1	manage	6.9 K	2010-12-31 - 17:02	GavinHesketh
patch	pythia.patch	r1	manage	22.9 K	2010-12-31 - 17:02	GavinHesketh
patch	rivet.patch	r1	manage	11.8 K	2010-12-31 - 15:59	GavinHesketh	Rivet patch

Topic revision: r3 - 2012-05-06 - GavinHesketh