Logo ROOT  
Reference Guide
rf503_wspaceread.C
Go to the documentation of this file.
1 /// \file
2 /// \ingroup tutorial_roofit
3 /// \notebook -js
4 ///
5 /// Organisation and simultaneous fits: reading and using a workspace
6 ///
7 /// The input file for this macro is generated by rf502_wspacewrite.C
8 ///
9 /// \macro_image
10 /// \macro_output
11 /// \macro_code
12 ///
13 /// \date 07/2008
14 /// \author Wouter Verkerke
15 
16 #include "RooRealVar.h"
17 #include "RooDataSet.h"
18 #include "RooGaussian.h"
19 #include "RooConstVar.h"
20 #include "RooChebychev.h"
21 #include "RooAddPdf.h"
22 #include "RooWorkspace.h"
23 #include "RooPlot.h"
24 #include "TCanvas.h"
25 #include "TAxis.h"
26 #include "TFile.h"
27 #include "TH1.h"
28 using namespace RooFit;
29 
30 void rf503_wspaceread()
31 {
32  // R e a d w o r k s p a c e f r o m f i l e
33  // -----------------------------------------------
34 
35  // Open input file with workspace (generated by rf14_wspacewrite)
36  TFile *f = new TFile("rf502_workspace.root");
37 
38  // Retrieve workspace from file
39  RooWorkspace *w = (RooWorkspace *)f->Get("w");
40 
41  // R e t r i e v e p d f , d a t a f r o m w o r k s p a c e
42  // -----------------------------------------------------------------
43 
44  // Retrieve x,model and data from workspace
45  RooRealVar *x = w->var("x");
46  RooAbsPdf *model = w->pdf("model");
47  RooAbsData *data = w->data("modelData");
48 
49  // Print structure of composite p.d.f.
50  model->Print("t");
51 
52  // F i t m o d e l t o d a t a , p l o t m o d e l
53  // ---------------------------------------------------------
54 
55  // Fit model to data
56  model->fitTo(*data);
57 
58  // Plot data and PDF overlaid
59  RooPlot *xframe = x->frame(Title("Model and data read from workspace"));
60  data->plotOn(xframe);
61  model->plotOn(xframe);
62 
63  // Overlay the background component of model with a dashed line
64  model->plotOn(xframe, Components("bkg"), LineStyle(kDashed));
65 
66  // Overlay the background+sig2 components of model with a dotted line
67  model->plotOn(xframe, Components("bkg,sig2"), LineStyle(kDotted));
68 
69  // Draw the frame on the canvas
70  new TCanvas("rf503_wspaceread", "rf503_wspaceread", 600, 600);
71  gPad->SetLeftMargin(0.15);
72  xframe->GetYaxis()->SetTitleOffset(1.4);
73  xframe->Draw();
74 }
RooWorkspace::data
RooAbsData * data(const char *name) const
Retrieve dataset (binned or unbinned) with given name. A null pointer is returned if not found.
Definition: RooWorkspace.cxx:1368
RooWorkspace.h
RooPlot::Draw
virtual void Draw(Option_t *options=0)
Draw this plot and all of the elements it contains.
Definition: RooPlot.cxx:691
RooChebychev.h
f
#define f(i)
Definition: RSha256.hxx:122
RooAbsData
Definition: RooAbsData.h:46
RooGaussian.h
x
Double_t x[n]
Definition: legend1.C:17
RooAddPdf.h
TCanvas.h
RooDataSet.h
TFile.h
RooFit
Definition: RooCFunction1Binding.h:29
RooAbsData::plotOn
virtual RooPlot * plotOn(RooPlot *frame, const RooCmdArg &arg1=RooCmdArg::none(), const RooCmdArg &arg2=RooCmdArg::none(), const RooCmdArg &arg3=RooCmdArg::none(), const RooCmdArg &arg4=RooCmdArg::none(), const RooCmdArg &arg5=RooCmdArg::none(), const RooCmdArg &arg6=RooCmdArg::none(), const RooCmdArg &arg7=RooCmdArg::none(), const RooCmdArg &arg8=RooCmdArg::none()) const
Definition: RooAbsData.cxx:547
RooPlot.h
RooPlot::GetYaxis
TAxis * GetYaxis() const
Definition: RooPlot.cxx:1256
RooPlot
Definition: RooPlot.h:44
RooWorkspace::pdf
RooAbsPdf * pdf(const char *name) const
Retrieve p.d.f (RooAbsPdf) with given name. A null pointer is returned if not found.
Definition: RooWorkspace.cxx:1277
RooRealVar.h
TFile
Definition: TFile.h:54
RooConstVar.h
RooWorkspace
Definition: RooWorkspace.h:43
TCanvas
Definition: TCanvas.h:23
TAxis.h
kDashed
@ kDashed
Definition: TAttLine.h:48
gPad
#define gPad
Definition: TVirtualPad.h:287
RooWorkspace::var
RooRealVar * var(const char *name) const
Retrieve real-valued variable (RooRealVar) with given name. A null pointer is returned if not found.
Definition: RooWorkspace.cxx:1295
make_cnn_model.model
model
Definition: make_cnn_model.py:6
RooAbsPdf
Definition: RooAbsPdf.h:40
RooFit::LineStyle
RooCmdArg LineStyle(Style_t style)
Definition: RooGlobalFunc.cxx:57
TAttAxis::SetTitleOffset
virtual void SetTitleOffset(Float_t offset=1)
Set distance between the axis and the axis title.
Definition: TAttAxis.cxx:293
RooRealVar
Definition: RooRealVar.h:35
RooFit::Components
RooCmdArg Components(const RooArgSet &compSet)
Definition: RooGlobalFunc.cxx:74
TH1.h
RooFit::Title
RooCmdArg Title(const char *name)
Definition: RooGlobalFunc.cxx:173
kDotted
@ kDotted
Definition: TAttLine.h:48