Logo ROOT  
Reference Guide
RooNLLVar.h
Go to the documentation of this file.
1/*****************************************************************************
2 * Project: RooFit *
3 * Package: RooFitCore *
4 * File: $Id: RooNLLVar.h,v 1.10 2007/07/21 21:32:52 wouter Exp $
5 * Authors: *
6 * WV, Wouter Verkerke, UC Santa Barbara, verkerke@slac.stanford.edu *
7 * DK, David Kirkby, UC Irvine, dkirkby@uci.edu *
8 * *
9 * Copyright (c) 2000-2005, Regents of the University of California *
10 * and Stanford University. All rights reserved. *
11 * *
12 * Redistribution and use in source and binary forms, *
13 * with or without modification, are permitted according to the terms *
14 * listed in LICENSE (http://roofit.sourceforge.net/license.txt) *
15 *****************************************************************************/
16#ifndef ROO_NLL_VAR
17#define ROO_NLL_VAR
18
20#include "RooCmdArg.h"
21#include "RooAbsPdf.h"
22#include <vector>
23#include <utility>
24
25class RooRealSumPdf ;
26namespace RooBatchCompute {
27struct RunContext;
28}
29
31public:
32
33 // Constructors, assignment etc
34 RooNLLVar();
35 RooNLLVar(const char *name, const char* title, RooAbsPdf& pdf, RooAbsData& data,
36 const RooCmdArg& arg1=RooCmdArg::none(), const RooCmdArg& arg2=RooCmdArg::none(),const RooCmdArg& arg3=RooCmdArg::none(),
37 const RooCmdArg& arg4=RooCmdArg::none(), const RooCmdArg& arg5=RooCmdArg::none(),const RooCmdArg& arg6=RooCmdArg::none(),
38 const RooCmdArg& arg7=RooCmdArg::none(), const RooCmdArg& arg8=RooCmdArg::none(),const RooCmdArg& arg9=RooCmdArg::none()) ;
39
40 RooNLLVar(const char *name, const char *title, RooAbsPdf& pdf, RooAbsData& data,
41 bool extended,
43
44 RooNLLVar(const char *name, const char *title, RooAbsPdf& pdf, RooAbsData& data,
45 const RooArgSet& projDeps, bool extended = false,
47
48 RooNLLVar(const RooNLLVar& other, const char* name=0);
49 virtual TObject* clone(const char* newname) const { return new RooNLLVar(*this,newname); }
50
51 virtual RooAbsTestStatistic* create(const char *name, const char *title, RooAbsReal& pdf, RooAbsData& adata,
52 const RooArgSet& projDeps, RooAbsTestStatistic::Configuration const& cfg);
53
54 virtual ~RooNLLVar();
55
56 void applyWeightSquared(Bool_t flag) ;
57
58 virtual Double_t defaultErrorLevel() const { return 0.5 ; }
59
60 void batchMode(bool on = true) {
62 }
63
64 using ComputeResult = std::pair<ROOT::Math::KahanSum<double>, double>;
65
66 static RooNLLVar::ComputeResult computeBatchedFunc(const RooAbsPdf *pdfClone, RooAbsData *dataClone,
67 std::unique_ptr<RooBatchCompute::RunContext> &evalData,
68 RooArgSet *normSet, bool weightSq, std::size_t stepSize,
69 std::size_t firstEvent, std::size_t lastEvent);
70 static RooNLLVar::ComputeResult computeScalarFunc(const RooAbsPdf *pdfClone, RooAbsData *dataClone, RooArgSet *normSet,
71 bool weightSq, std::size_t stepSize, std::size_t firstEvent,
72 std::size_t lastEvent);
73
74protected:
75
76 virtual Bool_t processEmptyDataSets() const { return _extended ; }
77 virtual Double_t evaluatePartition(std::size_t firstEvent, std::size_t lastEvent, std::size_t stepSize) const;
78
79 static RooArgSet _emptySet ; // Supports named argument constructor
80
81private:
82 ComputeResult computeBatched(std::size_t stepSize, std::size_t firstEvent, std::size_t lastEvent) const;
83 ComputeResult computeScalar(std::size_t stepSize, std::size_t firstEvent, std::size_t lastEvent) const;
84
86 bool _batchEvaluations{false};
87 Bool_t _weightSq{false}; // Apply weights squared?
88 mutable Bool_t _first{true}; //!
90
91 mutable std::vector<Double_t> _binw ; //!
92 mutable RooRealSumPdf* _binnedPdf{nullptr}; //!
93 mutable std::unique_ptr<RooBatchCompute::RunContext> _evalData; //! Struct to store function evaluation workspaces.
94
95 ClassDef(RooNLLVar,3) // Function representing (extended) -log(L) of p.d.f and dataset
96};
97
98#endif
99
bool Bool_t
Definition: RtypesCore.h:63
double Double_t
Definition: RtypesCore.h:59
#define ClassDef(name, id)
Definition: Rtypes.h:325
char name[80]
Definition: TGX11.cxx:110
The Kahan summation is a compensated summation algorithm, which significantly reduces numerical error...
Definition: Util.h:122
RooAbsData is the common abstract base class for binned and unbinned datasets.
Definition: RooAbsData.h:79
RooAbsOptTestStatistic is the abstract base class for test statistics objects that evaluate a functio...
RooAbsReal is the common abstract base class for objects that represent a real value and implements f...
Definition: RooAbsReal.h:61
RooAbsTestStatistic is the abstract base class for all test statistics.
RooArgSet is a container object that can hold multiple RooAbsArg objects.
Definition: RooArgSet.h:35
RooCmdArg is a named container for two doubles, two integers two object points and three string point...
Definition: RooCmdArg.h:27
static const RooCmdArg & none()
Return reference to null argument.
Definition: RooCmdArg.cxx:52
Class RooNLLVar implements a -log(likelihood) calculation from a dataset and a PDF.
Definition: RooNLLVar.h:30
ComputeResult computeScalar(std::size_t stepSize, std::size_t firstEvent, std::size_t lastEvent) const
Definition: RooNLLVar.cxx:520
void applyWeightSquared(Bool_t flag)
Definition: RooNLLVar.cxx:250
virtual TObject * clone(const char *newname) const
Definition: RooNLLVar.h:49
RooRealSumPdf * _binnedPdf
Definition: RooNLLVar.h:92
static RooNLLVar::ComputeResult computeScalarFunc(const RooAbsPdf *pdfClone, RooAbsData *dataClone, RooArgSet *normSet, bool weightSq, std::size_t stepSize, std::size_t firstEvent, std::size_t lastEvent)
Definition: RooNLLVar.cxx:526
ROOT::Math::KahanSum< double > _offsetSaveW2
Definition: RooNLLVar.h:89
static RooNLLVar::ComputeResult computeBatchedFunc(const RooAbsPdf *pdfClone, RooAbsData *dataClone, std::unique_ptr< RooBatchCompute::RunContext > &evalData, RooArgSet *normSet, bool weightSq, std::size_t stepSize, std::size_t firstEvent, std::size_t lastEvent)
Definition: RooNLLVar.cxx:417
Bool_t _extended
Definition: RooNLLVar.h:85
bool _batchEvaluations
Definition: RooNLLVar.h:86
static RooArgSet _emptySet
Definition: RooNLLVar.h:79
virtual ~RooNLLVar()
Destructor.
Definition: RooNLLVar.cxx:241
Bool_t _first
Definition: RooNLLVar.h:88
void batchMode(bool on=true)
Definition: RooNLLVar.h:60
virtual Double_t defaultErrorLevel() const
Definition: RooNLLVar.h:58
virtual Bool_t processEmptyDataSets() const
Definition: RooNLLVar.h:76
std::pair< ROOT::Math::KahanSum< double >, double > ComputeResult
Definition: RooNLLVar.h:64
std::vector< Double_t > _binw
Definition: RooNLLVar.h:91
virtual RooAbsTestStatistic * create(const char *name, const char *title, RooAbsReal &pdf, RooAbsData &adata, const RooArgSet &projDeps, RooAbsTestStatistic::Configuration const &cfg)
Create a test statistic using several properties of the current instance.
Definition: RooNLLVar.cxx:226
virtual Double_t evaluatePartition(std::size_t firstEvent, std::size_t lastEvent, std::size_t stepSize) const
Calculate and return likelihood on subset of data.
Definition: RooNLLVar.cxx:278
std::unique_ptr< RooBatchCompute::RunContext > _evalData
Definition: RooNLLVar.h:93
Bool_t _weightSq
Definition: RooNLLVar.h:87
ComputeResult computeBatched(std::size_t stepSize, std::size_t firstEvent, std::size_t lastEvent) const
Compute probabilites of all data events.
Definition: RooNLLVar.cxx:410
The class RooRealSumPdf implements a PDF constructed from a sum of functions:
Definition: RooRealSumPdf.h:24
Mother of all ROOT objects.
Definition: TObject.h:37
Namespace for dispatching RooFit computations to various backends.