Herwig info on rivet analyses
See this
Other Rivet info
See also ./LESTER

Installed like this:


wget http://svn.hepforge.org/rivet/bootstrap/rivet-bootstrap
chmod +x rivet-bootstrap
export PATH=/afs/cern.ch/sw/lcg/external/Python/2.6.5/x86_64-slc5-gcc43-opt/bin:$PATH
export LD_LIBRARY_PATH=/afs/cern.ch/sw/lcg/external/Python/2.6.5/x86_64-slc5-gcc43-opt/lib:$LD_LIBRARY_PATH
export PATH=/afs/cern.ch/sw/lcg/external/swig/1.3.40/x86_64-slc5-gcc43-opt/bin:$PATH
source /afs/cern.ch/sw/lcg/contrib/gcc/4.3/x86_64-slc5-gcc43-opt/setup.sh
./rivet-bootstrap --prefix=/atlas/lester/usr/local
./rivet-bootstrap --prefix=/atlas/lester/usr/local --install-agile
mv {rivet,agile}*sh /atlas/lester/usr/local/rivet

following instructions on http://projects.hepforge.org/rivet/trac/wiki/GettingStarted

Then I copied the build directoreies to

/usera/lester/proj/utils/MCGen/Rivet-1.6.0



Use with:

First 4 lines make us use the compiler and python versions that LCG is expecting, not slc5 defaults:

export PATH=/afs/cern.ch/sw/lcg/external/Python/2.6.5/x86_64-slc5-gcc43-opt/bin:$PATH
export LD_LIBRARY_PATH=/afs/cern.ch/sw/lcg/external/Python/2.6.5/x86_64-slc5-gcc43-opt/lib:$LD_LIBRARY_PATH
export PATH=/afs/cern.ch/sw/lcg/external/swig/1.3.40/x86_64-slc5-gcc43-opt/bin:$PATH
source /afs/cern.ch/sw/lcg/contrib/gcc/4.3/x86_64-slc5-gcc43-opt/setup.sh
source /usera/lester/proj/utils/MCGen/Rivet-1.6.0/rivetenv.sh
source /usera/lester/proj/utils/MCGen/AGILe-1.2.2/agileenv.sh
#source /atlas/lester/usr/local/rivet/rivetenv.sh
#source /atlas/lester/usr/local/rivet/agileenv.sh

and then follow instructions on http://projects.hepforge.org/rivet/trac/wiki/GettingStarted and http://projects.hepforge.org/rivet/trac/wiki/FirstRivetRun

agile-runmc --list-gens
agile-runmc Pythia6:423 --beams pp:14TeV

etc

and

My test games are all in

cd /usera/lester/proj/utils/MCGen/Rivet-1.6.0/LESTER

mkfifo /tmp/fifo.hepmc
agile-runmc Pythia6:423 --beams=TVT:1800 -n 500 -o /tmp/fifo.hepmc &
rivet --analysis=CDF_2001_S4751469 /tmp/fifo.hepmc
rivet-mkhtml Rivet.aida
# makes 21 plots with names like
#  ./plots/CDF_2001_S4751469/d01-x01-y01.dat
#  ./plots/CDF_2001_S4751469/d01-x01-y01.pdf
# viewable with:
firefox file:///usera/lester/proj/utils/MCGen/Rivet-1.6.0/LESTER/plots/CDF_2001_S4751469/index.html &

agile-runmc Pythia6:423 --beams=pp:900 -n 500 -o /tmp/fifo.hepmc &
rivet --analysis=ATLAS_2010_S8591806 /tmp/fifo.hepmc
rivet-mkhtml Rivet.aida

###### START ASIDE
#Note that there are MT2 rivet analyses:
grep -i mt2 ../src/Analyses/*.cc
../src/Analyses/ATLAS_2011_S8983313.cc:#include "Rivet/Tools/RivetMT2.hh"
../src/Analyses/ATLAS_2011_S8983313.cc:      _hist_mT2_B   = bookHistogram1D("m_T2", 25, 0., 1000.);
../src/Analyses/ATLAS_2011_S8983313.cc:      double m_T2 = mT2::mT2( a,
../src/Analyses/ATLAS_2011_S8983313.cc:      _hist_mT2_B->fill(m_T2, weight);
####### END ASIDE

cp ../src/Analyses/ATLAS_2010_S8591806.cc ./RENNER.cc
sed -i -e 's/ATLAS_2010_S8591806/RENNER/g' RENNER.cc
# Change class name to RENNER, Analsis name to "RENNER", 
# MUST THEN MAKE THE HISTO BOOKING CONCRETE in the init() method of RENNER.cc ...
      _h_dNch_deta = bookHistogram1D("eta",40, -3, 3);
      _h_dNch_dpT = bookHistogram1D("pT",30, 0, 100*GeV);
      _h_dNevt_dNch = bookHistogram1D("Nch",40, 0, 200);
      _p_meanpT_Nch = bookProfile1D("meanpT_Nch",50, 0, 100*GeV);
# THEN
rivet-buildplugin RivetRENNER.so RENNER.cc
rivet --analysis-path-append=. --list-analyses 'R.*'
agile-runmc Pythia6:423 --beams=pp:900 -n 500 -o /tmp/fifo.hepmc &
rivet --analysis=RENNER --analysis-path-append=. -H=RENNER.aida /tmp/fifo.hepmc
rivet-mkhtml RENNER.aida



# To make web pages under ./plots/ANALYSIS_NAME/ from Rivet.aida
rivet-mkhtml --mc-errs Rivet.aida
# In addition, this makes lots of dat files in ./plots/ANALYSIS_NAME/*.dat (one for each plot) that contain both styles and data.  These can be edited (eg to remove plot styles).  Then you can "re-run" the plot creation with

make-plots --pdfpng plots/LESTER_TEST1/*.dat

# Alternatively, you could put somthing like this in a file LESTER_TEAS1.plot:
cat << EOF > LESTER_TEST1.plot
# BEGIN PLOT ./plots/LESTER_TEST1/abc.dat
XTwosidedTicks=1
Title=Moo
LogY=0
RatioPlotYLabel=Ratio
XLabel=
YLabel=
YTwosidedTicks=1
Legend=1
# END PLOT

# ... add more histograms as you need them ...
EOF

and then do 

rivet-mkhtml --mc-errs -c LESTER_TEST1.plot Rivet.aida

and it will get the web-page right first time.

You can control herwig pt cuts (at generation) by adding lines like


set /Herwig/Cuts/JetKtCut:MinKT 30.0*GeV

to the WHATEVER.in file that makes the WHATEVER.run file ...

See B.5 of  http://arxiv.org/abs/0803.0883 and examples in

/usera/lester/proj/utils/MCGen/herwig-fun-2/../Herwig++-2.5.1/src/defaults/Cuts.in

and LESTERs own LHC.in file: /usera/lester/proj/utils/MCGen/herwig-fun-2/LHC.in ... 

export NUM=500000
( cd /usera/lester/proj/utils/MCGen/herwig-fun-2 ; Herwig++ run LHC.run -N$NUM -d1 >& Herwig++.log ) &


Turning on hepMC writing in herwig++:

insert LHCGenerator:AnalysisHandlers 0 /Herwig/Analysis/HepMCFile
set /Herwig/Analysis/HepMCFile:PrintEvent 10000000
set /Herwig/Analysis/HepMCFile:Format GenEvent
set /Herwig/Analysis/HepMCFile:Units GeV_mm
set /Herwig/Analysis/HepMCFile:Filename /tmp/fifo.hepmc


and then running with -N50000 (or whatever)

Lally archive