doc/v634/hf001__example_8C.html

#include <RooStats/HistFactory/MakeModelAndMeasurementsFast.h>

#include <RooStats/HistFactory/Measurement.h>


#include <TFile.h>

#include <TROOT.h>


void hf001_example()

{

   std::string InputFile = "./data/example.root";

   // in case the file is not found

   bool bfile = gSystem->AccessPathName(InputFile.c_str());

   if (bfile) {

      std::cout << "Input file is not found - run prepareHistFactory script " << std::endl;

      gROOT->ProcessLine(".! prepareHistFactory .");

      bfile = gSystem->AccessPathName(InputFile.c_str());

      if (bfile) {

         std::cout << "Still no " << InputFile << ", giving up.\n";

         exit(1);

      }

   }


   // Create the measurement

   RooStats::HistFactory::Measurement meas("meas", "meas");


   meas.SetOutputFilePrefix("./results/example_UsingC");

   meas.SetPOI("SigXsecOverSM");

   meas.AddConstantParam("alpha_syst1");

   meas.AddConstantParam("Lumi");


   meas.SetLumi(1.0);

   meas.SetLumiRelErr(0.10);

   meas.SetBinHigh(2);


   // Create a channel


   RooStats::HistFactory::Channel chan("channel1");

   chan.SetData("data", InputFile);

   chan.SetStatErrorConfig(0.05, "Poisson");


   // Now, create some samples


   // Create the signal sample

   RooStats::HistFactory::Sample signal("signal", "signal", InputFile);

   signal.AddOverallSys("syst1", 0.95, 1.05);

   signal.AddNormFactor("SigXsecOverSM", 1, 0, 3);

   chan.AddSample(signal);


   // Background 1

   RooStats::HistFactory::Sample background1("background1", "background1", InputFile);

   background1.ActivateStatError("background1_statUncert", InputFile);

   background1.AddOverallSys("syst2", 0.95, 1.05);

   chan.AddSample(background1);


   // Background 1

   RooStats::HistFactory::Sample background2("background2", "background2", InputFile);

   background2.ActivateStatError();

   background2.AddOverallSys("syst3", 0.95, 1.05);

   chan.AddSample(background2);


   // Done with this channel

   // Add it to the measurement:

   meas.AddChannel(chan);


   // Collect the histograms from their files,

   // print some output,

   meas.CollectHistograms();

   meas.PrintTree();


   // One can print XML code to an output directory:

   // meas.PrintXML( "xmlFromCCode", meas.GetOutputFilePrefix() );


   // Now, do the measurement

   std::unique_ptr<RooWorkspace> ws{MakeModelAndMeasurementFast(meas)};


   RooStats::ModelConfig *modelConfig = static_cast<RooStats::ModelConfig *>(ws->obj("ModelConfig"));


   // Get probability density function and parameters list from model

   RooAbsPdf *pdf = modelConfig->GetPdf();

   RooArgSet globalObservables{*modelConfig->GetGlobalObservables()};


   // Perform the fit using Minos to get the correct asymmetric uncertainties

   using namespace RooFit;

   std::unique_ptr<RooFitResult> result{

      pdf->fitTo(*ws->data("obsData"), Save(), PrintLevel(-1), GlobalObservables(globalObservables), Minos(true))};


   // Getting list of Parameters of Interest and getting first from them

   RooRealVar *poi = static_cast<RooRealVar *>(modelConfig->GetParametersOfInterest()->first());


   std::unique_ptr<RooAbsReal> nll{pdf->createNLL(*ws->data("obsData"))};

   std::unique_ptr<RooAbsReal> profile{nll->createProfile(*poi)};


   // frame for future plot

   RooPlot *frame = poi->frame();


   frame->SetTitle("");

   frame->GetYaxis()->SetTitle("-log likelihood");

   frame->GetXaxis()->SetTitle(poi->GetTitle());


   TCanvas *profileLikelihoodCanvas = new TCanvas{"combined", "", 800, 600};


   double xmin = poi->getMin();

   double xmax = poi->getMax();

   TLine *line = new TLine(xmin, .5, xmax, .5);

   line->SetLineColor(kGreen);

   TLine *line90 = new TLine(xmin, 2.71 / 2., xmax, 2.71 / 2.);

   line90->SetLineColor(kGreen);

   TLine *line95 = new TLine(xmin, 3.84 / 2., xmax, 3.84 / 2.);

   line95->SetLineColor(kGreen);

   frame->addObject(line);

   frame->addObject(line90);

   frame->addObject(line95);


   nll->plotOn(frame, ShiftToZero(), LineColor(kRed), LineStyle(kDashed));

   profile->plotOn(frame);


   frame->SetMinimum(0);

   frame->SetMaximum(2.);


   frame->Draw();


   // Print fit results to console in verbose mode to see asymmetric uncertainties

   result->Print("v");

}

MakeModelAndMeasurementsFast.h

Measurement.h

kRed
@ kRed
Definition Rtypes.h:66

kGreen
@ kGreen
Definition Rtypes.h:66

kDashed
@ kDashed
Definition TAttLine.h:48

TRangeDynCast
ROOT::Detail::TRangeCast< T, true > TRangeDynCast
TRangeDynCast is an adapter class that allows the typed iteration through a TCollection.
Definition TCollection.h:358

TFile.h

result
Option_t Option_t TPoint TPoint const char GetTextMagnitude GetFillStyle GetLineColor GetLineWidth GetMarkerStyle GetTextAlign GetTextColor GetTextSize void char Point_t Rectangle_t WindowAttributes_t Float_t Float_t Float_t Int_t Int_t UInt_t UInt_t Rectangle_t result
Definition TGWin32VirtualXProxy.cxx:174

xmin
float xmin
Definition THbookFile.cxx:95

xmax
float xmax
Definition THbookFile.cxx:95

TROOT.h

gROOT
#define gROOT
Definition TROOT.h:414

gSystem
R__EXTERN TSystem * gSystem
Definition TSystem.h:572

ROOT::Detail::TRangeCast
Definition TCollection.h:311

RooAbsPdf
Abstract interface for all probability density functions.
Definition RooAbsPdf.h:40

RooAbsPdf::createNLL
RooFit::OwningPtr< RooAbsReal > createNLL(RooAbsData &data, CmdArgs_t const &... cmdArgs)
Construct representation of -log(L) of PDF with given dataset.
Definition RooAbsPdf.h:163

RooAbsPdf::fitTo
RooFit::OwningPtr< RooFitResult > fitTo(RooAbsData &data, CmdArgs_t const &... cmdArgs)
Fit PDF to given dataset.
Definition RooAbsPdf.h:157

RooAbsRealLValue::getMax
virtual double getMax(const char *name=nullptr) const
Get maximum of currently defined range.
Definition RooAbsRealLValue.h:84

RooAbsRealLValue::getMin
virtual double getMin(const char *name=nullptr) const
Get minimum of currently defined range.
Definition RooAbsRealLValue.h:81

RooAbsRealLValue::frame
RooPlot * frame(const RooCmdArg &arg1, const RooCmdArg &arg2={}, const RooCmdArg &arg3={}, const RooCmdArg &arg4={}, const RooCmdArg &arg5={}, const RooCmdArg &arg6={}, const RooCmdArg &arg7={}, const RooCmdArg &arg8={}) const
Create a new RooPlot on the heap with a drawing frame initialized for this object,...
Definition RooAbsRealLValue.cxx:216

RooArgSet
RooArgSet is a container object that can hold multiple RooAbsArg objects.
Definition RooArgSet.h:24

RooPlot
Plot frame and a container for graphics objects within that frame.
Definition RooPlot.h:45

RooPlot::SetTitle
void SetTitle(const char *name) override
Set the title of the RooPlot to 'title'.
Definition RooPlot.cxx:1243

RooPlot::addObject
void addObject(TObject *obj, Option_t *drawOptions="", bool invisible=false)
Add a generic object to this plot.
Definition RooPlot.cxx:366

RooPlot::GetYaxis
TAxis * GetYaxis() const
Definition RooPlot.cxx:1264

RooPlot::Draw
void Draw(Option_t *options=nullptr) override
Draw this plot and all of the elements it contains.
Definition RooPlot.cxx:637

RooPlot::GetXaxis
TAxis * GetXaxis() const
Definition RooPlot.cxx:1262

RooPlot::SetMaximum
SetMaximum(ymax)

RooPlot::SetMinimum
SetMinimum(ymin)

RooRealVar
Variable that can be changed from the outside.
Definition RooRealVar.h:37

RooStats::HistFactory::Channel
This class encapsulates all information for the statistical interpretation of one experiment.
Definition Channel.h:30

RooStats::HistFactory::Measurement
The RooStats::HistFactory::Measurement class can be used to construct a model by combining multiple R...
Definition Measurement.h:33

RooStats::HistFactory::Sample
Definition Sample.h:27

RooStats::ModelConfig
ModelConfig is a simple class that holds configuration information specifying how a model should be u...
Definition ModelConfig.h:35

TAttLine::SetLineColor
virtual void SetLineColor(Color_t lcolor)
Set the line color.
Definition TAttLine.h:40

TCanvas
The Canvas class.
Definition TCanvas.h:23

TLine
Use the TLine constructor to create a simple line.
Definition TLine.h:22

TNamed::SetTitle
virtual void SetTitle(const char *title="")
Set the title of the TNamed.
Definition TNamed.cxx:164

TNamed::GetTitle
const char * GetTitle() const override
Returns title of object.
Definition TNamed.h:48

TSystem::AccessPathName
virtual Bool_t AccessPathName(const char *path, EAccessMode mode=kFileExists)
Returns FALSE if one can access a file using the specified access mode.
Definition TSystem.cxx:1296

line
TLine * line
Definition entrylistblock_figure1.C:235

RooFit::Save
RooCmdArg Save(bool flag=true)
Definition RooGlobalFunc.cxx:649

RooFit::GlobalObservables
RooCmdArg GlobalObservables(Args_t &&... argsOrArgSet)
Definition RooGlobalFunc.h:306

RooFit::Minos
RooCmdArg Minos(bool flag=true)
Definition RooGlobalFunc.cxx:681

RooFit::PrintLevel
RooCmdArg PrintLevel(Int_t code)
Definition RooGlobalFunc.cxx:657

RooFit::ShiftToZero
RooCmdArg ShiftToZero()
Definition RooGlobalFunc.cxx:191

RooFit::LineColor
RooCmdArg LineColor(Color_t color)
Definition RooGlobalFunc.cxx:215

RooFit::LineStyle
RooCmdArg LineStyle(Style_t style)
Definition RooGlobalFunc.cxx:240

RooFit::Detail::MathFuncs::nll
double nll(double pdf, double weight, int binnedL, int doBinOffset)
Definition MathFuncs.h:377

RooFit
The namespace RooFit contains mostly switches that change the behaviour of functions of PDFs (or othe...
Definition JSONIO.h:26

RooStats::HistFactory::MakeModelAndMeasurementFast
RooFit::OwningPtr< RooWorkspace > MakeModelAndMeasurementFast(RooStats::HistFactory::Measurement &measurement, HistoToWorkspaceFactoryFast::Configuration const &cfg={})

hf001_example
Definition hf001_example.py:1