Top MC validation (How To)

Rivet

Rivet is a software toolkit for data analysis of MC event generator samples, providing tools (jet finders, event shapes, etc.) for easy writing of analyses which are not generator-specific, and a large (250+) collection of MC versions of published experimental analyses. http://rivet.hepforge.org.

Producing EVGEN files

Running Rivet 1.X

First, log on lxplus5 (I haven't had luck with SLC6 yet) and create a temporal directory to work there:

tmpdir=$(mktemp -d)
mkdir ${tmpdir}/evgen

Rivet will run over EVGEN files, so one can produce these files or they can be retrieved from the GRID if they already exists. Let's do the latter:

cd ${tmpdir}/evgen
export ATLAS_LOCAL_ROOT_BASE=/cvmfs/atlas.cern.ch/repo/ATLASLocalRootBase
alias setupATLAS='source ${ATLAS_LOCAL_ROOT_BASE}/user/atlasLocalSetup.sh'
export ALRB_localConfigDir=$HOME/localConfig
setupATLAS --quiet
localSetupDQ2Client --quiet --skipConfirm
voms-proxy-init -voms atlas

Take for example:

dq2-get -f "*.root" -n 2 mc12_8TeV.110101.AcerMCPythia_P2011CCTEQ6L1_singletop_tchan_l.evgen.EVNT.e1731/

Now, one can easily produce the AIDA files just running the runRivetValidation.sh script which is available here. It can be retrieved simply doing:

cd ${tmpdir}
svn export svn+ssh://svn.cern.ch/reps/atlas-cescobar/cescobar/TopMC/validation/scripts/runRivetValidation.sh

Just run this script in the temporal directory that you just produce:

./runRivetValidation.sh

This script will run over all the directories you have in ${tmpdir}/evgen and it will run locally one athena job per directory in parallel. Once all the jobs have finished you will find a AIDA file per DSID in ${tmpdir} directory and also a copy in your EOS area. Now let's produce some plots.

Producing plots using Rivet

Setup rivet:

asetup 17.2.10,here
source /afs/cern.ch/sw/lcg/external/MCGenerators_hepmc2.06.05/rivet/1.8.1/${CMTCONFIG}/rivetenv.sh

Now, one can use compare-histos utility to produce all the plots that are defined in the plug-in you selected.

mkdir plots
compare-histos mc12_8TeV.110101.aida:'Title=AcerMC+Pythia,P2011C,110101' mc12_8TeV.110105.aida:'Title=AcerMC+Pythia,P2011C,MorePS,110105' mc12_8TeV.110106.aida:'Title=AcerMC+Pythia,P2011C,LessPS,110106'
make-plots --png *.dat
mv MC_* plots

Or one can use the rivet-mkhtml utility to produce the plots together with a (ugly) webpage.

rivet-mkhtml mc12_8TeV.110101.aida:'Title=AcerMC+Pythia,P2011C,110101' mc12_8TeV.110105.aida:'Title=AcerMC+Pythia,P2011C,MorePS,110105' mc12_8TeV.110106.aida:'Title=AcerMC+Pythia,P2011C,LessPS,110106' --mc-errs -s

-- CarlosEscobar - 14-May-2013

Edit | Attach | Watch | Print version | History: r3 < r2 < r1 | Backlinks | Raw View | Raw edit | More topic actions...
Topic revision: r2 - 2014-03-07 - CarlosEscobar
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Main All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2019 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback