%CERTIFY%

Implementation and validation of ATLAS Rivet analyses

Introduction

This page gives a brief overview of the status of Rivet routines for ATLAS analyses. See RivetForAtlas for instructions on using Rivet within Athena jobs.

We want your analysis to be included into Rivet in order to be later used for different MC generators validation/tuning. For that, we need the Rivet routine (.cc code) mimicking your analysis and the data histograms from your plots (.aida file). Once we have that, anyone is able to run different MC generators (with different tunings) to compare to your data results which then allow one to better tune them.

How to write Rivet routine for your analysis

The Rivet routine of your physics analysis is just one .cc file where the analysis flow of your analysis is done. It is then used as plugin to Rivet. You don't need to change anything in Rivet itself. There is tool ('rivet-buildplugin') which makes shared object library from this routine and it can be consequently used within rivet

To start writing the Rivet routine for your analysis you can do this in a coule ways:

• copy one of the existing analysis (.cc file) which is closest to yours and start rewriting it in your favorite editor according to your analysis flow. The list of the analyses can be obtained, e.g. if you are looking for ATLAS dijet analysis by "rivet --list-analyses|grep ATLAS|grep -i dijet". The actual files can then be found in e.g. /afs/cern.ch/sw/lcg/external/MCGenerators_hepmc2.06.05/rivet/1.8.1/share/sources/src/Analyses/

• you can start from scratch by running: 'rivet-mkanalysis MY_ANALYSIS_NAME' which will create a basic template for your code.

Some general information about how to write an analysis can be found here:

http://rivet.hepforge.org/trac/wiki/WritingAnAnalysis

Rules for implementation assistance

We would like to ask those seeking help with the implementation of their analysis in Rivet to consider the following "rules" that will make life much easier. Please

• provide the final corrected histograms (for experimental data only, not MC) in a suitable format, e.g. AIDA-xml (preferable), ROOT or text (each line corresponds to one bin with the following info: xLow xHigh y yErrLow yErrHigh where the errors are total uncertainties)
  • ideally, update your histograms to HepDATA (see e.g. http://hepdata.cedar.ac.uk/view/red4992)

• provide link to latest documentation of your analysis (e.g. paper draft, internal note, etc.)

• provide a table linking those histograms with the figures in your paper or draft

• provide the EVGEN dataset names of those samples that appear on the final plots (see last column of the following table)

Ongoing tasks

Analyses and validation samples Analysis Description Contact Status Validation DS ATLAS_2012_I1083318 http://cdsweb.cern.ch/record/1361687 https://cdsweb.cern.ch/record/1361687/files/ATL-COM-PHYS-2011-769.pdf?version=8 Jets accomp. W Alexander Paramonov Data ready (HS), implementation done by Frank Siegert, in Rivet as unvalidated analysis since r3540 mc11_7TeV.129121.SherpaW5jetstoenu30GeVPDF.evgen.EVNT.e914 mc10_7TeV.107680.AlpgenJimmyWenuNp0_pt20.evgen.EVNT.e600 mc10_7TeV.107681.AlpgenJimmyWenuNp1_pt20.evgen.EVNT.e600 mc10_7TeV.107682.AlpgenJimmyWenuNp2_pt20.evgen.EVNT.e600 mc10_7TeV.107683.AlpgenJimmyWenuNp3_pt20.evgen.EVNT.e600 mc10_7TeV.107684.AlpgenJimmyWenuNp4_pt20.evgen.EVNT.e600 mc10_7TeV.107685.AlpgenJimmyWenuNp5_pt20.evgen.EVNT.e600 mc11_7TeV.129122.SherpaW5jetstomunu30GeVPDF.evgen.EVNT.e914 mc10_7TeV.107690.AlpgenJimmyWmunuNp0_pt20.evgen.EVNT.e600 mc10_7TeV.107691.AlpgenJimmyWmunuNp1_pt20.evgen.EVNT.e600 mc10_7TeV.107692.AlpgenJimmyWmunuNp2_pt20.evgen.EVNT.e600 mc10_7TeV.107693.AlpgenJimmyWmunuNp3_pt20.evgen.EVNT.e600 mc10_7TeV.107694.AlpgenJimmyWmunuNp4_pt20.evgen.EVNT.e600 mc10_7TeV.107695.AlpgenJimmyWmunuNp5_pt20.evgen.EVNT.e600 ATLAS_2011_I945498 http://arxiv.org/abs/ARXIV:1111.2690 plots Z+jets Evelin Meoni In Rivet, branch 2012-06-aidarivet since r3812 (marked unvalidated). mc10_7TeV.119126.SherpaZ5jetstoee30GeVScaleMT.evgen.EVNT.e760 mc10_7TeV.119127.SherpaZ5jetstomumu30GeVScaleMT.evgen.EVNT.e760 ATLAS_2011_I926145 http://arxiv.org/abs/1109.0525 Inclusive lepton analysis Paul Bell Done, in Rivet since r3535 mc09_7TeV.105724.PythiaB_bbe3X.evgen.EVNT.e511 mc09_7TeV.105723.PythiaB_cce3X.evgen.EVNT.e511 mc09_7TeV.105805.filtered_minbias6.evgen.EVNT.e530 mc10_7TeV.105001.pythia_minbias.evgen.EVNT.e574/ ATLAS_2011_I944826 http://arxiv.org/abs/1111.1297 Kshort and lambda Jed Besiada Done. In Rivet since r3577 mc10_7TeV.126209.HerwigppMinBias_K0SLamDec_UE7000.evgen.EVNT.e852/ mc10_7TeV.126214.pythia_minbias_K0SLamDec.evgen.EVNT.e852/ mc10_7TeV.126217.pythia_minbias_K0SLamDec_Perugia2010.evgen.EVNT.e852/ mc10_7TeV.126222.pythia_minbias_Z1_K0SLamDec.evgen.EVNT.e869/ mc10_7TeV.126225.pythia_minbias_AMBT2B_CTEQ6L1_K0SLamDec.evgen.EVNT.e869/ mc10_7TeV.126228.pythia_minbias_Perugia2011_K0SLamDec.evgen.EVNT.e869/ mc10_900GeV.126210.HerwigppMinBias_K0SLamDec_MU900.evgen.EVNT.e854/ mc10_900GeV.126214.pythia_minbias_K0SLamDec.evgen.EVNT.e854/ mc10_900GeV.126217.pythia_minbias_K0SLamDec_Perugia2010.evgen.EVNT.e854/ mc10_900GeV.126222.pythia_minbias_Z1_K0SLamDec.evgen.EVNT.e870/ mc10_900GeV.126225.pythia_minbias_AMBT2B_CTEQ6L1_K0SLamDec.evgen.EVNT.e870/ mc10_900GeV.126228.pythia_minbias_Perugia2011_K0SLamDec.evgen.EVNT.e870/ ATLAS_2012_I1082936 http://arxiv.org/abs/1112.6297 Plots inclusive jets/dijets Eric Feng, Mario Campanelli (Implementation partially by Serguei Yanush), Pavel Starovoitiov Done (HS). In rivet since r3722 mc11_7TeV.126362.SherpaQCDJetsJ0.evgen.EVNT.e969 mc11_7TeV.126363.SherpaQCDJetsJ1.evgen.EVNT.e969_tid584747_00 mc11_7TeV.126364.SherpaQCDJetsJ2.evgen.EVNT.e969 mc11_7TeV.126365.SherpaQCDJetsJ3.evgen.EVNT.e969 mc11_7TeV.126366.SherpaQCDJetsJ4.evgen.EVNT.e969 mc11_7TeV.126367.SherpaQCDJetsJ5.evgen.EVNT.e969 mc11_7TeV.126368.SherpaQCDJetsJ6.evgen.EVNT.e969 mc11_7TeV.126369.SherpaQCDJetsJ7.evgen.EVNT.e969 mc11_7TeV.126370.SherpaQCDJetsJ8.evgen.EVNT.e969 ATLAS_2012_I1084540 http://arxiv.org/abs/1201.2808 plots Rapidity gap cross-section Tim Martin Implementation by Tim Martin. In Rivet since r3548 mc11_7TeV.108316.Pythia8_minbias_ND.evgen.EVNT.e825 mc11_7TeV.108317.Pythia8_minbias_SD.evgen.EVNT.e825 mc11_7TeV.108318.Pythia8_minbias_DD.evgen.EVNT.e825 mc11_7TeV.105001.pythia_minbias.evgen.EVNT.e815 mc11_7TeV.105003.pythia_sdiff.evgen.EVNT.e815 mc11_7TeV.105004.pythia_ddiff.evgen.EVNT.e815 mc10_7TeV.106096.PhojetNdiff.evgen.EVNT.e827 mc10_7TeV.106097.PhojetSdiff.evgen.EVNT.e827 mc10_7TeV.106098.PhojetDdiff.evgen.EVNT.e827 mc10_7TeV.106099.PhojetCdiff.evgen.EVNT.e827 ATLAS_2011_I954993 http://arxiv.org/pdf/1111.5570v2.pdf WZ cross-section Lynn Marx validated, in Rivet since r3940 mc10_7TeV.126053.WZMcAtNloJimmy12870_1_1.evgen.EVNT.e787 mc10_7TeV.126054.WZMcAtNloJimmy12880_1_1.evgen.EVNT.e787 mc10_7TeV.126055.WZMcAtNloJimmy12870_1_3.evgen.EVNT.e787 mc10_7TeV.126056.WZMcAtNloJimmy12880_1_3.evgen.EVNT.e787 mc10_7TeV.126059.WZMcAtNloJimmy12870_2_1.evgen.EVNT.e787 mc10_7TeV.126060.WZMcAtNloJimmy12880_2_1.evgen.EVNT.e787 mc10_7TeV.126061.WZMcAtNloJimmy12870_2_3.evgen.EVNT.e787 mc10_7TeV.126062.WZMcAtNloJimmy12880_2_3.evgen.EVNT.e787 ATLAS_2011_I930220 http://arxiv.org/abs/1109.6833 Plots https://cdsweb.cern.ch/record/1386708 Inclusive and dijets with B-hadrons Stephen Paul Bieniek implemented by Stephen, preparing for validation mc10_7TeV.113893.BBbar_PowHegBOX_ME_MSTW08_Pythia.evgen.EVNT.e735 mc10_7TeV.105009.J0_pythia_jetjet.evgen.EVNT.e574 mc10_7TeV.105010.J1_pythia_jetjet.evgen.EVNT.e913 mc10_7TeV.105011.J2_pythia_jetjet.evgen.EVNT.e913 mc10_7TeV.105012.J3_pythia_jetjet.evgen.EVNT.e913 mc10_7TeV.105013.J4_pythia_jetjet.evgen.EVNT.e913 mc10_7TeV.105014.J5_pythia_jetjet.evgen.EVNT.e913 ATLAS_2012_I1094568 http://arxiv.org/abs/1203.5015 ttbar jet veto Kiran Joshi Done. In Rivet trunk since r3654 mc10_7TeV.105205.AcerMCttbar.evgen.EVNT.e574/ mc10_7TeV.117255.AcerMCttbar_isr_down.evgen.EVNT.e574/ mc10_7TeV.117256.AcerMCttbar_isr_up.evgen.EVNT.e574/ ATLAS_2011_? http://prl.aps.org/abstract/PRL/v108/i4/e041804 ZZ cross-section Shih-Chieh Hsu not implemented mc10_7TeV.109291.PythiaZZ4l.evgen.EVNT.e787 ATLAS_2012_I1093738 http://arxiv.org/pdf/arXiv:1203.3161 photon + jets Giovanni Marchiori Done. In Rivet since r3717. mc10_valid.119901.SherpaY5jets15GeV.evgen.EVNT.e694 mc10_7TeV.113714.SherpaY4JetsPt35.evgen.EVNT.e761 mc10_7TeV.113715.SherpaY4JetsPt70.evgen.EVNT.e761 mc10_7TeV.113716.SherpaY4JetsPt140.evgen.EVNT.e761 mc10_7TeV.113717.SherpaY4JetsPt280.evgen.EVNT.e761 mc10_7TeV.108087.PythiaPhotonJet_Unbinned17.evgen.EVNT.e598 mc10_7TeV.108081.PythiaPhotonJet_Unbinned35.evgen.EVNT.e598 mc10_7TeV.108082.PythiaPhotonJet_Unbinned70.evgen.EVNT.e598 mc10_7TeV.108083.PythiaPhotonJet_Unbinned140.evgen.EVNT.e598 mc10_7TeV.108084.PythiaPhotonJet_Unbinned280.evgen.EVNT.e598 mc10_7TeV.115802.DP17_pythia_photon_filter.evgen.EVNT.e687 mc10_7TeV.115803.DP35_pythia_photon_filter.evgen.EVNT.e615 mc10_7TeV.115804.DP70_pythia_photon_filter.evgen.EVNT.e615 mc10_7TeV.115809.DP140_pythia_photon_filter.evgen.EVNT.e681 mc10_7TeV.115810.DP280_pythia_photon_filter.evgen.EVNT.e681 ATLAS_2012_I1119557 http://inspirehep.net/record/1119557 https://cdsweb.cern.ch/record/1431554/files/ Jet properties Lily Asquith implemented, validating mc11_7TeV.113639.Pythia8JetsJ4.evgen.EVNT.e995 mc11_7TeV.113640.Pythia8JetsJ5.evgen.EVNT.e995 mc11_7TeV.113641.Pythia8JetsJ6.evgen.EVNT.e995 mc11_7TeV.113642.Pythia8JetsJ7.evgen.EVNT.e995 mc11_7TeV.105013.J4_pythia_jetjet.evgen.EVNT.e815 mc11_7TeV.105014.J5_pythia_jetjet.evgen.EVNT.e815 mc11_7TeV.105015.J6_pythia_jetjet.evgen.EVNT.e815 mc11_7TeV.105016.J7_pythia_jetjet.evgen.EVNT.e815 ATLAS_2012_I1094564 http://arxiv.org/abs/1203.4606 plots Jet substructure Adam Davison Implemented, preparing for validation N.N. ATLAS_2012_I1091481 http://arxiv.org/abs/1203.0419 Plots Azimuthal ordering Sharka Todorova Done. In rivet since r3771 (marked unvalidated) ATLAS_2012_I1183818 http://inspirehep.net/record/1183818/ https://cdsweb.cern.ch/record/1435780 ET flow Robindra Prabhu validated, putting to Rivet mc11_7TeV.126921.pythia_AUET2B_CTEQ6L1_minbias.evgen.EVNT.e1250/ mc11_7TeV.113203.HerwigppMinBias.evgen.EVNT.e916/ mc11_7TeV.126922.pythia_AUET2B_CTEQ6L1_minbias_truthJet10.evgen.EVNT.e1250/ mc11_7TeV.126576.HerwigppJetsJ1.evgen.EVNT.e1070/ ATLAS_2012_I1093734 http://inspirehep.net/record/1093734 Plots FB and phi correlations Craig Buttar implemented, validating ATLAS_2012_IXXXXXXX https://cdsweb.cern.ch/record/1455420 ATLAS-COM-CONF-2012-118 flavor composition in dijets Cecile Lapoire implemented, validating mc10_7TeV.105009.J0_pythia_jetjet.evgen.EVNT.e574/ mc10_7TeV.105010.J1_pythia_jetjet.evgen.EVNT.e574/ mc10_7TeV.105011.J2_pythia_jetjet.evgen.EVNT.e574/ mc10_7TeV.105012.J3_pythia_jetjet.evgen.EVNT.e574/ mc10_7TeV.105013.J4_pythia_jetjet.evgen.EVNT.e574/ mc10_7TeV.105014.J5_pythia_jetjet.evgen.EVNT.e574/ mc10_7TeV.105015.J6_pythia_jetjet.evgen.EVNT.e574/ mc10_7TeV.105016.J7_pythia_jetjet.evgen.EVNT.e574/

Color definition:

not implemented

implemented, preparing for validation

validating

validated

ATLAS analysis wishlist

Here we list ATLAS analyses for which Rivet analyses are wanted:

Analysis	Physics group	Contact	Reference	Comments
ttbar + jets	Top	?		Early in analysis
UE in jet events	SM	Andy Buckley		AB will supply
Z->bb substructure	SM	Ben Cooper		Analysis incomplete
Event shapes in high momentum transfer events	SM	David Miller	STDM-2010-
Z->tautau	SM	Elias Coniavitis	STDM-2011-18 Plots
ttbar differential parton-level	Top	Francesco Spano		Ready for implementation
R_jets	SM	Gerhardt Brandt		Analysis incomplete
Two particle angular correlations	SM	James Monk	STDM-2010-11
b-jet masses etc.	SM	Javier Merino		Analysis incomplete
Jet shapes and masses	SM	Lily Asquith	STDM-2010-
WZ cross-section	SM	Lynn Marx	http://inspirehep.net/record/1126131/
ZZ->2l2v	SM		STDM-2010-
WW->lvlv	SM		STDM-2010-
W/Z inclusive	SM	Uta Klein, Liza Mijovic	STDM-2011-06
ttbar differential particle-level	Top	Will Bell		Analysis incomplete, but communication possible

Further projects

Project Description Status ** Tune job GRID submission framework Tesing phase (RL) ** Automated validation Unknown (RL) RivAid Website dedicated for ROOT file conversion to Rivet format (AIDA, YODA) Basically working, no server yet (HS)

Running custom analyses on EVNT datasets

Further explanations can be found at RivetForAtlas.

cd empty
asetup 17.0.4,here
source $SITEROOT/sw/lcg/external/MCGenerators/rivet/1.6.0/i686-slc5-gcc43-opt/rivetenv.sh
rivet-buildplugin RivetATLAS_2011_S9185208.so ATLAS_2011_S9185208.cc -m32

pathena IncLeptonsMB.py --inDS mc09_7TeV.105805.filtered_minbias6.evgen.EVNT.e530/ --outDS user.hschulz.ATLAS_2011_S9185208.filtered_minbias6.Rivet.v1/ --nGBPerJob=MAX  --extOutFile=Analysis.aida --extFile=RivetATLAS_2011_S9185208.so,ATLAS_2011_S9185208.aida

Here's how to run the merging on the Grid:

prun --noBuild --exec "bash ./AIDAmergeWrapper.sh %IN"  --athenaTag=17.0.4,AtlasProduction --outDS=user.hschulz.ATLAS_2011_S9185208.filtered_minbias6.Rivet.merged.v3/  --inDS=user.hschulz.ATLAS_2011_S9185208.filtered_minbias6.Rivet.v2/ --outputs merged.aida --destSE=DESY-HH_LOCALGROUPDISK  --workDir=. --nFilesPerJob 41 --match=aida

Extracting cross-section information from AMI

You can use the attached script after sourcing athena. You need to enter your AMI username and password. The usage is e.g.

./getXSfromAMI.py  mc10_7TeV.126053.WZMcAtNloJimmy12870_1_1.evgen.EVNT.e787  mc10_7TeV.126054.WZMcAtNloJimmy12880_1_1.evgen.EVNT.e787

This will give you the following output:

mc10_7TeV.126053.WZMcAtNloJimmy12870_1_1.evgen.EVNT.e787
    2.359700e-06   [nano barn]
    mc10_7TeV.126054.WZMcAtNloJimmy12880_1_1.evgen.EVNT.e787
    1.266800e-06   [nano barn]

Please do keep in mind that Rivet requires cross-section information to be given in [pb].

---

-- HolgerSchulz - 29-Nov-2011 -- AndyBuckley - 01-Mar-2012 -- RomanLysak - 25-Apr-2012

%RESPONSIBLE% KiranJoshi
%REVIEW% Never reviewed