Logo ROOT  

Reference Guide
 
Loading...
Searching...
No Matches
df102_NanoAODDimuonAnalysis.C File Reference
Tutorials » Data analysis tutorials » RDataFrame analysis tutorials

Detailed Description

View in nbviewer Open in SWAN
Show how NanoAOD files can be processed with RDataFrame.

This tutorial illustrates how NanoAOD files can be processed with ROOT dataframes. The NanoAOD-like input files are filled with 66 mio. events from CMS OpenData containing muon candidates part of 2012 dataset (DOI: 10.7483/OPENDATA.CMS.YLIC.86ZZ and DOI: 10.7483/OPENDATA.CMS.M5AD.Y3V3). The macro matches muon pairs and produces an histogram of the dimuon mass spectrum showing resonances up to the Z mass. Note that the bump at 30 GeV is not a resonance but a trigger effect.

More details about the dataset can be found on the CERN Open Data portal.

#include "ROOT/RDataFrame.hxx"
#include "ROOT/RDFHelpers.hxx"
#include "ROOT/RVec.hxx"
#include "TCanvas.h"
#include "TH1D.h"
#include "TLatex.h"
#include "TStyle.h"
using namespace ROOT::VecOps;
void df102_NanoAODDimuonAnalysis()
{
// Enable multi-threading
ROOT::EnableImplicitMT();
// Create dataframe from NanoAOD files
ROOT::RDataFrame df("Events", "root://eospublic.cern.ch//eos/opendata/cms/derived-data/AOD2NanoAODOutreachTool/"
"Run2012BC_DoubleMuParked_Muons.root");
// Add ProgressBar
ROOT::RDF::Experimental::AddProgressBar(df);
// For simplicity, select only events with exactly two muons and require opposite charge
auto df_2mu = df.Filter("nMuon == 2", "Events with exactly two muons");
auto df_os = df_2mu.Filter("Muon_charge[0] != Muon_charge[1]", "Muons with opposite charge");
// Compute invariant mass of the dimuon system
auto df_mass = df_os.Define("Dimuon_mass", InvariantMass<float>, {"Muon_pt", "Muon_eta", "Muon_phi", "Muon_mass"});
// Make histogram of dimuon mass spectrum. Note how we can set title and axis labels in one go
auto h = df_mass.Histo1D({"Dimuon_mass", "Dimuon mass;m_{#mu#mu} (GeV);N_{Events}", 30000, 0.25, 300}, "Dimuon_mass");
// Request cut-flow report
auto report = df.Report();
// Produce plot
gStyle->SetOptStat(0); gStyle->SetTextFont(42);
auto c = new TCanvas("c", "", 800, 700);
c->SetLogx(); c->SetLogy();
h->GetXaxis()->SetTitleSize(0.04);
h->GetYaxis()->SetTitleSize(0.04);
h->DrawClone();
TLatex label; label.SetNDC(true);
label.DrawLatex(0.175, 0.740, "#eta");
label.DrawLatex(0.205, 0.775, "#rho,#omega");
label.DrawLatex(0.270, 0.740, "#phi");
label.DrawLatex(0.400, 0.800, "J/#psi");
label.DrawLatex(0.415, 0.670, "#psi'");
label.DrawLatex(0.485, 0.700, "Y(1,2,3S)");
label.DrawLatex(0.755, 0.680, "Z");
label.SetTextSize(0.040); label.DrawLatex(0.100, 0.920, "#bf{CMS Open Data}");
label.SetTextSize(0.030); label.DrawLatex(0.630, 0.920, "#sqrt{s} = 8 TeV, L_{int} = 11.6 fb^{-1}");
c->SaveAs("dimuon_spectrum.pdf");
// Print cut-flow report
report->Print();
}
int main()
{
df102_NanoAODDimuonAnalysis();
}
main
int main()
Definition Prototype.cxx:12
c
#define c(i)
Definition RSha256.hxx:101
h
#define h(i)
Definition RSha256.hxx:106
gStyle
R__EXTERN TStyle * gStyle
Definition TStyle.h:442
ROOT::RDataFrame
ROOT's RDataFrame offers a modern, high-level interface for analysis of data stored in TTree ,...
Definition RDataFrame.hxx:41
TAttText::SetTextFont
virtual void SetTextFont(Font_t tfont=62)
Set the text font.
Definition TAttText.h:48
TAttText::SetTextSize
virtual void SetTextSize(Float_t tsize=1)
Set the text size.
Definition TAttText.h:49
TCanvas
The Canvas class.
Definition TCanvas.h:23
TLatex
To draw Mathematical Formula.
Definition TLatex.h:18
TLatex::DrawLatex
TLatex * DrawLatex(Double_t x, Double_t y, const char *text)
Make a copy of this object with the new parameters And copy object attributes.
Definition TLatex.cxx:1943
TObject::SaveAs
virtual void SaveAs(const char *filename="", Option_t *option="") const
Save this object in the file specified by filename.
Definition TObject.cxx:705
TStyle::SetOptStat
void SetOptStat(Int_t stat=1)
The type of information printed in the histogram statistics box can be selected via the parameter mod...
Definition TStyle.cxx:1642
TText::SetNDC
virtual void SetNDC(Bool_t isNDC=kTRUE)
Set NDC mode on if isNDC = kTRUE, off otherwise.
Definition TText.cxx:817
ROOT::RDF::Experimental::AddProgressBar
void AddProgressBar(ROOT::RDF::RNode df)
Add ProgressBar to a ROOT::RDF::RNode.
Definition RDFHelpers.cxx:375
ROOT::EnableImplicitMT
void EnableImplicitMT(UInt_t numthreads=0)
Enable ROOT's implicit multi-threading for all objects and methods that provide an internal paralleli...
Definition TROOT.cxx:539
df102_NanoAODDimuonAnalysis
Definition df102_NanoAODDimuonAnalysis.py:1
|> | [Elapsed time: 0:01m processing file: 1 / 1 processed evts: 1000 / 50123395 5.33e+02 evt/s 26:08:17h remaining time (per file being processed)]
|================================================================> | [Elapsed time: 0:02m processing file: 1 / 1 processed evts: 5332000 / 54231870 2.67e+06 evt/s 0:18m remaining time (per file being processed)]
|===================================================================================================================================================> | [Elapsed time: 0:03m processing file: 1 / 1 processed evts: 12374000 / 55053565 4.12e+06 evt/s 0:10m remaining time (per file being processed)]
|=======================================================================================================================================================================================================================================> | [Elapsed time: 0:04m processing file: 1 / 1 processed evts: 19603000 / 55875260 4.90e+06 evt/s 0:07m remaining time (per file being processed)]
|=============================================================================================================================================================================================================================================================================================================================> | [Elapsed time: 0:05m processing file: 1 / 1 processed evts: 27317000 / 56696955 5.46e+06 evt/s 0:05m remaining time (per file being processed)]
|==========================================================================================================================================================================================================================================================================================================================================================================================================> | [Elapsed time: 0:06m processing file: 1 / 1 processed evts: 34362000 / 57518650 5.73e+06 evt/s 0:04m remaining time (per file being processed)]
|===============================================================================================================================================================================================================================================================================================================================================================================================================================================================================> | [Elapsed time: 0:07m processing file: 1 / 1 processed evts: 42124000 / 59983735 6.02e+06 evt/s 0:02m remaining time (per file being processed)]
|============================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================> | [Elapsed time: 0:08m processing file: 1 / 1 processed evts: 49732000 / 60805430 6.22e+06 evt/s 0:01m remaining time (per file being processed)]
|==================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================> | [Elapsed time: 0:09m processing file: 1 / 1 processed evts: 56877000 / 61540413 6.32e+06 evt/s 0:00m remaining time (per file being processed)]
[Total elapsed time: 0:10m processed files: 1 / 1 processed evts: 61540413 / 61540413]
Events with exactly two muons: pass=31104343 all=61540413 -- eff=50.54 % cumulative eff=50.54 %
Muons with opposite charge: pass=24067843 all=31104343 -- eff=77.38 % cumulative eff=39.11 %
Date
August 2018
Author
Stefan Wunsch (KIT, CERN)

Definition in file df102_NanoAODDimuonAnalysis.C.