lally-10.txt
lally-100.txt
lally-1000.txt
lally-10000.txt
lally-100000.txt
lally-1000000.txt
LESTER_TEST1.cc
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
Link here