15 auto ret = dataFrame.Filter(
"TMath::Abs(md0_d - 1.8646) < 0.04")
16 .Filter(
"ptds_d > 2.5")
17 .Filter(
"TMath::Abs(etads_d) < 1.5")
18 .Filter([](
int ik,
int ipi,
RVecI& nhitrp) {
return nhitrp[ik - 1] * nhitrp[ipi - 1] > 1; },
19 {
"ik",
"ipi",
"nhitrp"})
20 .
Filter([](
int ik,
RVecF& rstart,
RVecF& rend) {
return rend[ik - 1] - rstart[ik - 1] > 22; },
21 {
"ik",
"rstart",
"rend"})
22 .
Filter([](
int ipi,
RVecF& rstart,
RVecF& rend) {
return rend[ipi - 1] - rstart[ipi - 1] > 22; },
23 {
"ipi",
"rstart",
"rend"})
24 .
Filter([](
int ik,
RVecF& nlhk) {
return nlhk[ik - 1] > 0.1; }, {
"ik",
"nlhk"})
25 .
Filter([](
int ipi,
RVecF& nlhpi) {
return nlhpi[ipi - 1] > 0.1; }, {
"ipi",
"nlhpi"})
26 .
Filter([](
int ipis,
RVecF& nlhpi) {
return nlhpi[ipis - 1] > 0.1; }, {
"ipis",
"nlhpi"})
41 dxbin * (par[0] *
pow(
x - 0.13957, par[1]) + par[2] / 2.5066 / par[4] *
exp(-xp3 / 2 / par[4] / par[4]));
56void FitAndPlotHdmd(
TH1 &hdmd)
60 auto c1 =
new TCanvas(
"c1",
"h1analysis analysis", 10, 10, 800, 600);
64 auto f5 =
new TF1(
"f5",
fdm5, 0.139, 0.17, 5);
65 f5->SetParameters(1000000, .25, 2000, .1454, .001);
71void FitAndPlotH2(
TH2 &h2)
74 auto c2 =
new TCanvas(
"c2",
"tauD0", 100, 100, 800, 600);
77 c2->SetBottomMargin(0.15);
83 auto f2 =
new TF1(
"f2",
fdm2, 0.139, 0.17, 2);
84 f2->SetParameters(10000, 10);
89 h2_1->GetXaxis()->SetTitle(
"#tau [ps]");
90 h2_1->SetMarkerStyle(21);
98void df101_h1Analysis()
101 chain.Add(
"root://eospublic.cern.ch//eos/root-eos/h1/dstarmb.root");
102 chain.Add(
"root://eospublic.cern.ch//eos/root-eos/h1/dstarp1a.root");
103 chain.Add(
"root://eospublic.cern.ch//eos/root-eos/h1/dstarp1b.root");
104 chain.Add(
"root://eospublic.cern.ch//eos/root-eos/h1/dstarp2.root");
109 auto selected = Select(dataFrame);
111 auto hdmdARP = selected.Histo1D({
"hdmd",
"Dm_d;m_{K#pi#pi} - m_{K#pi}[GeV/c^{2}]", 40, 0.13, 0.17},
"dm_d");
112 auto selectedAddedBranch = selected.Define(
"h2_y",
"rpd0_t / 0.029979f * 1.8646f / ptd0_d");
113 auto h2ARP = selectedAddedBranch.Histo2D({
"h2",
"ptD0 vs Dm_d", 30, 0.135, 0.165, 30, -3, 6},
"dm_d",
"h2_y");
115 FitAndPlotHdmd(*hdmdARP);
116 FitAndPlotH2(*h2ARP);
R__EXTERN TStyle * gStyle
ROOT's RDataFrame offers a modern, high-level interface for analysis of data stored in TTree ,...
virtual void SetTitleOffset(Float_t offset=1)
Set distance between the axis and the axis title.
A chain is a collection of files containing TTree objects.
1-D histogram with a double per channel (see TH1 documentation)}
TH1 is the base class of all histogram classes in ROOT.
virtual TFitResultPtr Fit(const char *formula, Option_t *option="", Option_t *goption="", Double_t xmin=0, Double_t xmax=0)
Fit histogram with function fname.
Service class for 2-D histogram classes.
virtual void FitSlicesX(TF1 *f1=nullptr, Int_t firstybin=0, Int_t lastybin=-1, Int_t cut=0, Option_t *option="QNR", TObjArray *arr=nullptr)
Project slices along X in case of a 2-D histogram, then fit each slice with function f1 and make a hi...
Use the TLine constructor to create a simple line.
virtual TObject * DrawClone(Option_t *option="") const
Draw a clone of this object in the current selected pad with: gROOT->SetSelectedPad(c1).
virtual void Draw(Option_t *option="")
Default Draw method for all objects.
void SetOptFit(Int_t fit=1)
The type of information about fit parameters printed in the histogram statistics box can be selected ...
RVec< PromoteTypes< T0, T1 > > pow(const T0 &x, const RVec< T1 > &v)
RVec< PromoteType< T > > exp(const RVec< T > &v)
RVec< T > Filter(const RVec< T > &v, F &&f)
Create a new collection with the elements passing the filter expressed by the predicate.
Double_t fdm5(Double_t *xx, Double_t *par)
Double_t fdm2(Double_t *xx, Double_t *par)
This file contains a specialised ROOT message handler to test for diagnostic in unit tests.
void EnableImplicitMT(UInt_t numthreads=0)
Enable ROOT's implicit multi-threading for all objects and methods that provide an internal paralleli...