Mini-Analysis on CMS OpenData with RDataFrame.
This tutorial illustrates that analyzing data with RDataFrame works the same for both TTree data and RNTuple data. The RNTuple data are converted from the Events tree in http://root.cern.ch/files/NanoAOD_DoubleMuon_CMS2011OpenData.root Based on RDataFrame's df102_NanoAODDimuonAnalysis.C
#include <ROOT/RDataFrame.hxx>
#include <cassert>
#include <cmath>
#include <iostream>
#include <memory>
#include <string>
#include <vector>
#include <utility>
constexpr char const* kNTupleFileName = "http://root.cern.ch/files/tutorials/ntpl004_dimuon_v1rc1.root";
void ntpl004_dimuon() {
auto df_2mu = df.Filter("#Muon == 2", "Events with exactly two muons");
auto df_os = df_2mu.Filter("Muon.charge[0] != Muon.charge[1]", "Muons with opposite charge");
auto df_mass = df_os.Define("Dimuon_mass", InvariantMass<float>, {"Muon.pt", "Muon.eta", "Muon.phi", "Muon.mass"});
auto h = df_mass.Histo1D({
"Dimuon_mass",
"Dimuon_mass", 30000, 0.25, 300},
"Dimuon_mass");
auto report = df_mass.Report();
auto c =
new TCanvas(
"c",
"", 800, 700);
c->SetLogx();
c->SetLogy();
h->GetXaxis()->SetTitle(
"m_{#mu#mu} (GeV)");
h->GetXaxis()->SetTitleSize(0.04);
h->GetYaxis()->SetTitle(
"N_{Events}");
h->GetYaxis()->SetTitleSize(0.04);
label.
DrawLatex(0.205, 0.775,
"#rho,#omega");
label.
SetTextSize(0.030); label.
DrawLatex(0.630, 0.920,
"#sqrt{s} = 8 TeV, L_{int} = 11.6 fb^{-1}");
}
#define R__LOAD_LIBRARY(LIBRARY)
R__EXTERN TStyle * gStyle
The RDataSource implementation for RNTuple.
An RNTuple that is used to read data from storage.
virtual void SetTextFont(Font_t tfont=62)
Set the text font.
virtual void SetTextSize(Float_t tsize=1)
Set the text size.
TLatex * DrawLatex(Double_t x, Double_t y, const char *text)
void SetOptStat(Int_t stat=1)
The type of information printed in the histogram statistics box can be selected via the parameter mod...
virtual void SetNDC(Bool_t isNDC=kTRUE)
void Print(Option_t *option="") const override
Print TNamed name and title.
RDataFrame FromRNTuple(std::string_view ntupleName, std::string_view fileName)
void EnableImplicitMT(UInt_t numthreads=0)
Enable ROOT's implicit multi-threading for all objects and methods that provide an internal paralleli...
- Date
- April 2019
- Author
- The ROOT Team
Definition in file ntpl004_dimuon.C.