Logo ROOT   6.18/05
Reference Guide
rf402_datahandling.C File Reference

Detailed Description

View in nbviewer Open in SWAN Data and categories: tools for manipulation of (un)binned datasets

␛[1mRooFit v3.60 -- Developed by Wouter Verkerke and David Kirkby␛[0m
Copyright (C) 2000-2013 NIKHEF, University of California & Stanford University
All rights reserved, please read http://roofit.sourceforge.net/license.txt
DataStore d (d)
Contains 1000 entries
Observables:
1) x = 9 L(-10 - 10) "x"
2) y = 31.607 L(0 - 40) "y"
3) c = Plus(idx = 1)
"c"
1) 0x55e550426020 RooRealVar:: x = 9 L(-10 - 10) "x"
2) 0x55e550796a50 RooRealVar:: y = 31.607 L(0 - 40) "y"
3) 0x55e550772440 RooCategory:: c = Plus(idx = 1)
"c"
1) 0x55e550426020 RooRealVar:: x = 8 L(-10 - 10) "x"
2) 0x55e550796a50 RooRealVar:: y = 30 L(0 - 40) "y"
3) 0x55e550772440 RooCategory:: c = Minus(idx = -1)
"c"
>> d1 has only columns x,c
DataStore d (d)
Contains 1000 entries
Observables:
1) x = 9 L(-10 - 10) "x"
2) c = Plus(idx = 1)
"c"
>> d2 has only column y
DataStore d (d)
Contains 1000 entries
Observables:
1) y = 31.607 L(0 - 40) "y"
>> d3 has only the points with y>5.17
DataStore d (d)
Contains 973 entries
Observables:
1) x = 9 L(-10 - 10) "x"
2) y = 31.607 L(0 - 40) "y"
3) c = Plus(idx = 1)
"c"
>> d4 has only columns x,c for data points with y>5.17
DataStore d (d)
Contains 973 entries
Observables:
1) x = 9 L(-10 - 10) "x"
2) c = Plus(idx = 1)
"c"
>> merge d2(y) with d1(x,c) to form d1(x,c,y)
DataStore d (d)
Contains 1000 entries
Observables:
1) x = 9 L(-10 - 10) "x"
2) c = Plus(idx = 1)
"c"
3) y = 31.607 L(0 - 40) "y"
>> append data points of d3 to d1
DataStore d (d)
Contains 1973 entries
Observables:
1) x = 9 L(-10 - 10) "x"
2) c = Plus(idx = 1)
"c"
3) y = 31.607 L(0 - 40) "y"
>> construct dh (binned) from d(unbinned) but only take the x and y dimensions,
>> the category 'c' will be projected in the filling process
DataStore dh (binned version of d)
Contains 100 entries
Observables:
1) x = 9 L(-10 - 10) B(10) "x"
2) y = 31.607 L(0 - 40) B(10) "y"
Binned Dataset dh (binned version of d)
Contains 100 bins with a total weight of 1000
Observables: 1) x = 9 L(-10 - 10) B(10) "x"
2) y = 31.607 L(0 - 40) B(10) "y"
>> number of bins in dh : 100
>> sum of weights in dh : 1000
>> integral over histogram: 8000
>> retrieving the properties of the bin enclosing coordinate (x,y) = (0.3,20.5)
bin center:
1) 0x55e55059ef60 RooRealVar:: x = 1 L(-10 - 10) B(10) "x"
2) 0x55e55087c170 RooRealVar:: y = 22 L(0 - 40) B(10) "y"
weight = 76
>> Creating 1-dimensional projection on y of dh for bins with x>0
DataStore dh (binned version of d)
Contains 10 entries
Observables:
1) y = 38 L(0 - 40) B(10) "y"
Binned Dataset dh (binned version of d)
Contains 10 bins with a total weight of 500
Observables: 1) y = 38 L(0 - 40) B(10) "y"
[#1] INFO:Plotting -- RooPlot::updateFitRangeNorm: New event count of 500 will supercede previous event count of 1000 for normalization of PDF projections
>> Persisting d via ROOT I/O
TFile** rf402_datahandling.root
TFile* rf402_datahandling.root
KEY: RooDataSet d;1 d
KEY: TProcessID ProcessID0;1 c9599a9c-bac5-11ec-b0d2-942c8a89beef
#include "RooRealVar.h"
#include "RooDataSet.h"
#include "RooDataHist.h"
#include "RooGaussian.h"
#include "RooConstVar.h"
#include "RooCategory.h"
#include "TCanvas.h"
#include "TAxis.h"
#include "RooPlot.h"
#include "TFile.h"
using namespace RooFit;
// WVE Add reduction by range
{
// Binned (RooDataHist) and unbinned datasets (RooDataSet) share
// many properties and inherit from a common abstract base class
// (RooAbsData), that provides an interface for all operations
// that can be performed regardless of the data format
RooRealVar x("x", "x", -10, 10);
RooRealVar y("y", "y", 0, 40);
RooCategory c("c", "c");
c.defineType("Plus", +1);
c.defineType("Minus", -1);
// B a s i c O p e r a t i o n s o n u n b i n n e d d a t a s e t s
// --------------------------------------------------------------
// RooDataSet is an unbinned dataset (a collection of points in N-dimensional space)
RooDataSet d("d", "d", RooArgSet(x, y, c));
// Unlike RooAbsArgs (RooAbsPdf,RooFormulaVar,....) datasets are not attached to
// the variables they are constructed from. Instead they are attached to an internal
// clone of the supplied set of arguments
// Fill d with dummy values
Int_t i;
for (i = 0; i < 1000; i++) {
x = i / 50 - 10;
y = sqrt(1.0 * i);
c.setLabel((i % 2) ? "Plus" : "Minus");
// We must explicitly refer to x,y,c here to pass the values because
// d is not linked to them (as explained above)
d.add(RooArgSet(x, y, c));
}
d.Print("v");
cout << endl;
// The get() function returns a pointer to the internal copy of the RooArgSet(x,y,c)
// supplied in the constructor
const RooArgSet *row = d.get();
row->Print("v");
cout << endl;
// Get with an argument loads a specific data point in row and returns
// a pointer to row argset. get() always returns the same pointer, unless
// an invalid row number is specified. In that case a null ptr is returned
d.get(900)->Print("v");
cout << endl;
// R e d u c i n g , A p p e n d i n g a n d M e r g i n g
// -------------------------------------------------------------
// The reduce() function returns a new dataset which is a subset of the original
cout << endl << ">> d1 has only columns x,c" << endl;
RooDataSet *d1 = (RooDataSet *)d.reduce(RooArgSet(x, c));
d1->Print("v");
cout << endl << ">> d2 has only column y" << endl;
RooDataSet *d2 = (RooDataSet *)d.reduce(RooArgSet(y));
d2->Print("v");
cout << endl << ">> d3 has only the points with y>5.17" << endl;
RooDataSet *d3 = (RooDataSet *)d.reduce("y>5.17");
d3->Print("v");
cout << endl << ">> d4 has only columns x,c for data points with y>5.17" << endl;
RooDataSet *d4 = (RooDataSet *)d.reduce(RooArgSet(x, c), "y>5.17");
d4->Print("v");
// The merge() function adds two data set column-wise
cout << endl << ">> merge d2(y) with d1(x,c) to form d1(x,c,y)" << endl;
d1->merge(d2);
d1->Print("v");
// The append() function addes two datasets row-wise
cout << endl << ">> append data points of d3 to d1" << endl;
d1->append(*d3);
d1->Print("v");
// O p e r a t i o n s o n b i n n e d d a t a s e t s
// ---------------------------------------------------------
// A binned dataset can be constructed empty, from an unbinned dataset, or
// from a ROOT native histogram (TH1,2,3)
cout << ">> construct dh (binned) from d(unbinned) but only take the x and y dimensions," << endl
<< ">> the category 'c' will be projected in the filling process" << endl;
// The binning of real variables (like x,y) is done using their fit range
// 'get/setRange()' and number of specified fit bins 'get/setBins()'.
// Category dimensions of binned datasets get one bin per defined category state
x.setBins(10);
y.setBins(10);
RooDataHist dh("dh", "binned version of d", RooArgSet(x, y), d);
dh.Print("v");
RooPlot *yframe = y.frame(Bins(10), Title("Operations on binned datasets"));
dh.plotOn(yframe); // plot projection of 2D binned data on y
// Examine the statistics of a binned dataset
cout << ">> number of bins in dh : " << dh.numEntries() << endl;
cout << ">> sum of weights in dh : " << dh.sum(kFALSE) << endl;
cout << ">> integral over histogram: " << dh.sum(kTRUE) << endl; // accounts for bin volume
// Locate a bin from a set of coordinates and retrieve its properties
x = 0.3;
y = 20.5;
cout << ">> retrieving the properties of the bin enclosing coordinate (x,y) = (0.3,20.5) " << endl;
cout << " bin center:" << endl;
dh.get(RooArgSet(x, y))->Print("v"); // load bin center coordinates in internal buffer
cout << " weight = " << dh.weight() << endl; // return weight of last loaded coordinates
// Reduce the 2-dimensional binned dataset to a 1-dimensional binned dataset
//
// All reduce() methods are interfaced in RooAbsData. All reduction techniques
// demonstrated on unbinned datasets can be applied to binned datasets as well.
cout << ">> Creating 1-dimensional projection on y of dh for bins with x>0" << endl;
RooDataHist *dh2 = (RooDataHist *)dh.reduce(y, "x>0");
dh2->Print("v");
// Add dh2 to yframe and redraw
dh2->plotOn(yframe, LineColor(kRed), MarkerColor(kRed));
// S a v i n g a n d l o a d i n g f r o m f i l e
// -------------------------------------------------------
// Datasets can be persisted with ROOT I/O
cout << endl << ">> Persisting d via ROOT I/O" << endl;
TFile f("rf402_datahandling.root", "RECREATE");
d.Write();
f.ls();
// To read back in future session:
// > TFile f("rf402_datahandling.root") ;
// > RooDataSet* d = (RooDataSet*) f.FindObject("d") ;
new TCanvas("rf402_datahandling", "rf402_datahandling", 600, 600);
gPad->SetLeftMargin(0.15);
yframe->GetYaxis()->SetTitleOffset(1.4);
yframe->Draw();
}
#define d(i)
Definition: RSha256.hxx:102
#define f(i)
Definition: RSha256.hxx:104
#define c(i)
Definition: RSha256.hxx:101
int Int_t
Definition: RtypesCore.h:41
const Bool_t kFALSE
Definition: RtypesCore.h:88
const Bool_t kTRUE
Definition: RtypesCore.h:87
@ kRed
Definition: Rtypes.h:64
double sqrt(double)
#define gPad
Definition: TVirtualPad.h:286
virtual void Print(Option_t *options=0) const
This method must be overridden when a class wants to print itself.
RooAbsData * reduce(const RooCmdArg &arg1, const RooCmdArg &arg2=RooCmdArg(), const RooCmdArg &arg3=RooCmdArg(), const RooCmdArg &arg4=RooCmdArg(), const RooCmdArg &arg5=RooCmdArg(), const RooCmdArg &arg6=RooCmdArg(), const RooCmdArg &arg7=RooCmdArg(), const RooCmdArg &arg8=RooCmdArg())
Create a reduced copy of this dataset.
Definition: RooAbsData.cxx:381
virtual void Print(Option_t *options=0) const
Print TNamed name and title.
Definition: RooAbsData.h:162
RooArgSet is a container object that can hold multiple RooAbsArg objects.
Definition: RooArgSet.h:28
RooCategory represents a fundamental (non-derived) discrete value object.
Definition: RooCategory.h:24
The RooDataHist is a container class to hold N-dimensional binned data.
Definition: RooDataHist.h:40
virtual RooPlot * plotOn(RooPlot *frame, PlotOpt o) const
Back end function to plotting functionality.
RooDataSet is a container class to hold unbinned data.
Definition: RooDataSet.h:31
void append(RooDataSet &data)
Add all data points of given data set to this data set.
Bool_t merge(RooDataSet *data1, RooDataSet *data2=0, RooDataSet *data3=0, RooDataSet *data4=0, RooDataSet *data5=0, RooDataSet *data6=0)
A RooPlot is a plot frame and a container for graphics objects within that frame.
Definition: RooPlot.h:41
TAxis * GetYaxis() const
Definition: RooPlot.cxx:1123
virtual void Draw(Option_t *options=0)
Draw this plot and all of the elements it contains.
Definition: RooPlot.cxx:558
RooRealVar represents a fundamental (non-derived) real valued object.
Definition: RooRealVar.h:36
virtual void SetTitleOffset(Float_t offset=1)
Set distance between the axis and the axis title Offset is a correction factor with respect to the "s...
Definition: TAttAxis.cxx:294
The Canvas class.
Definition: TCanvas.h:31
A ROOT file is a suite of consecutive data records (TKey instances) with a well defined format.
Definition: TFile.h:48
Double_t y[n]
Definition: legend1.C:17
Double_t x[n]
Definition: legend1.C:17
Template specialisation used in RooAbsArg:
RooCmdArg MarkerColor(Color_t color)
RooCmdArg LineColor(Color_t color)
RooCmdArg Bins(Int_t nbin)
const char * Title
Definition: TXMLSetup.cxx:67
Author
07/2008 - Wouter Verkerke

Definition in file rf402_datahandling.C.