%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)
- 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)
Completed Analyses:
ATLAS Analysis wishlist
Color definition:
Not implemented | | Implementation complete. Unvalidated. | Validation in progress | Validated |
ATLAS analysis wishlist
Here we list ATLAS analyses for which Rivet analyses are wanted:
Analysis |
Physics group |
Contact |
Reference |
Comments |
UE in jet events |
SM |
Andy Buckley |
|
AB will supply |
Z->bb substructure |
SM |
Ben Cooper |
|
Analysis incomplete |
R_jets |
SM |
Gerhardt Brandt |
|
Analysis incomplete |
b-jet masses etc. |
SM |
Javier Merino |
|
Analysis incomplete |
ttbar differential particle-level |
Top |
Will Bell |
|
Analysis incomplete, but communication possible |
ttbar + jets |
Top |
? |
|
Early in analysis |
Two particle angular correlations |
SM |
James Monk |
STDM-2010-11 |
|
Z->tautau |
SM |
Elias Coniavitis |
STDM-2011-18 Plots |
|
Jet shapes and masses |
SM |
Lily Asquith |
STDM-2010- |
|
Event shapes in high momentum transfer events |
SM |
David Miller |
STDM-2010- |
|
ZZ->2l2v |
SM |
|
STDM-2010- |
|
WW->lvlv |
SM |
|
STDM-2010- |
|
W/Z inclusive |
SM |
Uta Klein, Liza Mijovic |
STDM-2011-06 |
|
WZ cross-section |
SM |
Lynn Marx |
http://inspirehep.net/record/1126131/ |
|
ttbar differential parton-level |
Top |
Francesco Spano |
|
Ready for implementation |
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