Difference: GavinHeskethSandbox (2 vs. 3)

Revision 32012-05-06 - GavinHesketh

  META TOPICPARENT 
 name="Main.GavinHesketh"
- META TOPICPARENT
+ name="Main.GavinHesketh"
-<
<
+ Running Professor on Alpgen
->
>
+ 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
->
>
+export PATH=`pwd`/local/bin/:$PATH
 export LD_LIBRARY_PATH=`pwd`/local/lib
-<
<
+Then apply some patches:
->
>
+./rivet-bootstrap --install-boost
-<
<
+patch -p0 -i XXX.patch
->
>
+./rivet-bootstrap --install-agile
-<
<
+) HepMC.patch
->
>
+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.
-<
<
+) rivet.patch
->
>
+This can be done eg with this patch rivet.patch
-<
<
+These fix a crash when the last event is empty, and save the analysis normalisations to allow combination of algen samples. Agile has now been fixed to allow PYTUNE calls, so no need to patch that any more!
  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
-<
<
+Now some more patches:
->
>
+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
->
>
+ alpgen.patch
-<
<
+) pythia.patch

to automate the alpgen matching stuff, and to allow the latest tunes in pytune. Remake the packages after patching...
->
>
+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









 META FILEATTACHMENT 
 attachment="HepMC.patch" attr="" comment="HepMC patch" date="1293810097" name="HepMC.patch" path="HepMC.patch" size="729" user="GavinHesketh" version="1" 


 META FILEATTACHMENT 
 attachment="agile.patch" attr="" comment="Agile patch" date="1293810004" name="agile.patch" path="agile.patch" size="2417" user="GavinHesketh" version="2" 


 META FILEATTACHMENT 
 attachment="pythia.patch" attr="" comment="" date="1293814942" name="pythia.patch" path="pythia.patch" size="23418" user="GavinHesketh" version="1" 


 META FILEATTACHMENT 
 attachment="alpgen.patch" attr="" comment="" date="1293814934" name="alpgen.patch" path="alpgen.patch" size="7113" user="GavinHesketh" version="1" 


 META FILEATTACHMENT 
 attachment="rivet.patch" attr="" comment="Rivet patch" date="1293811164" name="rivet.patch" path="rivet.patch" size="12115" user="GavinHesketh" version="1" 


 META FILEATTACHMENT 
 attachment="agile-genser-bootstrap" attr="" comment="agile genser script" date="1293814689" name="agile-genser-bootstrap" path="agile-genser-bootstrap" size="15411" user="GavinHesketh" version="1"
- META FILEATTACHMENT
+ attachment="HepMC.patch" attr="" comment="HepMC patch" date="1293810097" name="HepMC.patch" path="HepMC.patch" size="729" user="GavinHesketh" version="1"
- META FILEATTACHMENT
+ attachment="agile.patch" attr="" comment="Agile patch" date="1293810004" name="agile.patch" path="agile.patch" size="2417" user="GavinHesketh" version="2"
- META FILEATTACHMENT
+ attachment="pythia.patch" attr="" comment="" date="1293814942" name="pythia.patch" path="pythia.patch" size="23418" user="GavinHesketh" version="1"
- META FILEATTACHMENT
+ attachment="alpgen.patch" attr="" comment="" date="1293814934" name="alpgen.patch" path="alpgen.patch" size="7113" user="GavinHesketh" version="1"
- META FILEATTACHMENT
+ attachment="rivet.patch" attr="" comment="Rivet patch" date="1293811164" name="rivet.patch" path="rivet.patch" size="12115" user="GavinHesketh" version="1"
- META FILEATTACHMENT
+ attachment="agile-genser-bootstrap" attr="" comment="agile genser script" date="1293814689" name="agile-genser-bootstrap" path="agile-genser-bootstrap" size="15411" user="GavinHesketh" version="1"

View topic | History: r3 < r2 < r1 | More topic actions...