Logo ROOT  

Reference Guide
RooAbsPdf.h
Go to the documentation of this file.
1 /*****************************************************************************
2  * Project: RooFit *
3  * Package: RooFitCore *
4  * File: $Id: RooAbsPdf.h,v 1.90 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_ABS_PDF
17 #define ROO_ABS_PDF
18 
19 #include "RooAbsReal.h"
20 //#include "RooRealIntegral.h"
21 #include "RooNameSet.h"
22 #include "RooObjCacheManager.h"
23 #include "RooCmdArg.h"
24 
25 class RooDataSet;
26 class RooDataHist ;
27 class RooArgSet ;
28 class RooAbsGenContext ;
29 class RooFitResult ;
30 class RooExtendPdf ;
31 class RooCategory ;
32 class TPaveText;
33 class TH1F;
34 class TH2F;
35 class TList ;
36 class RooLinkedList ;
37 class RooNumGenConfig ;
38 class RooRealIntegral ;
39 
40 class RooAbsPdf : public RooAbsReal {
41 public:
42 
43  // Constructors, assignment etc
44  RooAbsPdf() ;
45  RooAbsPdf(const char *name, const char *title=0) ;
46  RooAbsPdf(const char *name, const char *title, Double_t minVal, Double_t maxVal) ;
47  // RooAbsPdf(const RooAbsPdf& other, const char* name=0);
48  virtual ~RooAbsPdf();
49 
50  // Toy MC generation
51 
52  ////////////////////////////////////////////////////////////////////////////////
53  /// See RooAbsPdf::generate(const RooArgSet&,const RooCmdArg&,const RooCmdArg&,const RooCmdArg&,const RooCmdArg&,const RooCmdArg&,const RooCmdArg&)
54  /// \param[in] nEvents How many events to generate
55  RooDataSet *generate(const RooArgSet &whatVars, Int_t nEvents, const RooCmdArg& arg1,
56  const RooCmdArg& arg2=RooCmdArg::none(), const RooCmdArg& arg3=RooCmdArg::none(),
57  const RooCmdArg& arg4=RooCmdArg::none(), const RooCmdArg& arg5=RooCmdArg::none()) {
58  return generate(whatVars,RooFit::NumEvents(nEvents),arg1,arg2,arg3,arg4,arg5) ;
59  }
60  RooDataSet *generate(const RooArgSet &whatVars,
61  const RooCmdArg& arg1=RooCmdArg::none(),const RooCmdArg& arg2=RooCmdArg::none(),
62  const RooCmdArg& arg3=RooCmdArg::none(),const RooCmdArg& arg4=RooCmdArg::none(),
63  const RooCmdArg& arg5=RooCmdArg::none(),const RooCmdArg& arg6=RooCmdArg::none()) ;
64  RooDataSet *generate(const RooArgSet &whatVars, Double_t nEvents = 0, Bool_t verbose=kFALSE, Bool_t autoBinned=kTRUE,
65  const char* binnedTag="", Bool_t expectedData=kFALSE, Bool_t extended = kFALSE) const;
66  RooDataSet *generate(const RooArgSet &whatVars, const RooDataSet &prototype, Int_t nEvents= 0,
67  Bool_t verbose=kFALSE, Bool_t randProtoOrder=kFALSE, Bool_t resampleProto=kFALSE) const;
68 
69 
70  class GenSpec {
71  public:
72  virtual ~GenSpec() ;
73  GenSpec() { _genContext = 0 ; _protoData = 0 ; _init = kFALSE ; _extended=kFALSE, _nGen=0 ; _randProto = kFALSE ; _resampleProto=kFALSE ; }
74  private:
75  GenSpec(RooAbsGenContext* context, const RooArgSet& whatVars, RooDataSet* protoData, Int_t nGen, Bool_t extended,
76  Bool_t randProto, Bool_t resampleProto, TString dsetName, Bool_t init=kFALSE) ;
77  GenSpec(const GenSpec& other) ;
78 
79  friend class RooAbsPdf ;
80  RooAbsGenContext* _genContext ;
81  RooArgSet _whatVars ;
82  RooDataSet* _protoData ;
83  Int_t _nGen ;
84  Bool_t _extended ;
85  Bool_t _randProto ;
86  Bool_t _resampleProto ;
87  TString _dsetName ;
88  Bool_t _init ;
89  ClassDef(GenSpec,0) // Generation specification
90  } ;
91 
92  ///Prepare GenSpec configuration object for efficient generation of multiple datasets from identical specification.
93  GenSpec* prepareMultiGen(const RooArgSet &whatVars,
94  const RooCmdArg& arg1=RooCmdArg::none(),const RooCmdArg& arg2=RooCmdArg::none(),
95  const RooCmdArg& arg3=RooCmdArg::none(),const RooCmdArg& arg4=RooCmdArg::none(),
96  const RooCmdArg& arg5=RooCmdArg::none(),const RooCmdArg& arg6=RooCmdArg::none()) ;
97  ///Generate according to GenSpec obtained from prepareMultiGen().
98  RooDataSet* generate(GenSpec&) const ;
99 
100 
101  ////////////////////////////////////////////////////////////////////////////////
102  /// As RooAbsPdf::generateBinned(const RooArgSet&, const RooCmdArg&,const RooCmdArg&, const RooCmdArg&,const RooCmdArg&, const RooCmdArg&,const RooCmdArg&)
103  /// \param[in] nEvents How many events to generate
104  virtual RooDataHist *generateBinned(const RooArgSet &whatVars, Double_t nEvents, const RooCmdArg& arg1,
105  const RooCmdArg& arg2=RooCmdArg::none(), const RooCmdArg& arg3=RooCmdArg::none(),
106  const RooCmdArg& arg4=RooCmdArg::none(), const RooCmdArg& arg5=RooCmdArg::none()) const {
107  return generateBinned(whatVars,RooFit::NumEvents(nEvents),arg1,arg2,arg3,arg4,arg5);
108  }
109  virtual RooDataHist *generateBinned(const RooArgSet &whatVars,
110  const RooCmdArg& arg1=RooCmdArg::none(),const RooCmdArg& arg2=RooCmdArg::none(),
111  const RooCmdArg& arg3=RooCmdArg::none(),const RooCmdArg& arg4=RooCmdArg::none(),
112  const RooCmdArg& arg5=RooCmdArg::none(),const RooCmdArg& arg6=RooCmdArg::none()) const;
113  virtual RooDataHist *generateBinned(const RooArgSet &whatVars, Double_t nEvents, Bool_t expectedData=kFALSE, Bool_t extended=kFALSE) const;
114 
115  virtual RooDataSet* generateSimGlobal(const RooArgSet& whatVars, Int_t nEvents) ;
116 
117  ///Helper calling plotOn(RooPlot*, RooLinkedList&) const
118  virtual RooPlot* plotOn(RooPlot* frame,
119  const RooCmdArg& arg1=RooCmdArg::none(), const RooCmdArg& arg2=RooCmdArg::none(),
120  const RooCmdArg& arg3=RooCmdArg::none(), const RooCmdArg& arg4=RooCmdArg::none(),
121  const RooCmdArg& arg5=RooCmdArg::none(), const RooCmdArg& arg6=RooCmdArg::none(),
122  const RooCmdArg& arg7=RooCmdArg::none(), const RooCmdArg& arg8=RooCmdArg::none(),
123  const RooCmdArg& arg9=RooCmdArg::none(), const RooCmdArg& arg10=RooCmdArg::none()
124  ) const {
125  return RooAbsReal::plotOn(frame,arg1,arg2,arg3,arg4,arg5,arg6,arg7,arg8,arg9,arg10) ;
126  }
127  virtual RooPlot* plotOn(RooPlot* frame, RooLinkedList& cmdList) const ;
128 
129  /// Add a box with parameter values (and errors) to the specified frame
130  virtual RooPlot* paramOn(RooPlot* frame,
131  const RooCmdArg& arg1=RooCmdArg::none(), const RooCmdArg& arg2=RooCmdArg::none(),
132  const RooCmdArg& arg3=RooCmdArg::none(), const RooCmdArg& arg4=RooCmdArg::none(),
133  const RooCmdArg& arg5=RooCmdArg::none(), const RooCmdArg& arg6=RooCmdArg::none(),
134  const RooCmdArg& arg7=RooCmdArg::none(), const RooCmdArg& arg8=RooCmdArg::none()) ;
135 
136  virtual RooPlot* paramOn(RooPlot* frame, const RooAbsData* data, const char *label= "", Int_t sigDigits = 2,
137  Option_t *options = "NELU", Double_t xmin=0.50,
138  Double_t xmax= 0.99,Double_t ymax=0.95) ;
139 
140  // Built-in generator support
141  virtual Int_t getGenerator(const RooArgSet& directVars, RooArgSet &generateVars, Bool_t staticInitOK=kTRUE) const;
142  virtual void initGenerator(Int_t code) ;
143  virtual void generateEvent(Int_t code);
144  virtual Bool_t isDirectGenSafe(const RooAbsArg& arg) const ;
145 
146  // Configuration of MC generators used for this pdf
147  const RooNumGenConfig* getGeneratorConfig() const ;
148  static RooNumGenConfig* defaultGeneratorConfig() ;
149  RooNumGenConfig* specialGeneratorConfig() const ;
150  RooNumGenConfig* specialGeneratorConfig(Bool_t createOnTheFly) ;
151  void setGeneratorConfig() ;
152  void setGeneratorConfig(const RooNumGenConfig& config) ;
153 
154  // -log(L) fits to binned and unbinned data
155  virtual RooFitResult* fitTo(RooAbsData& data, const RooCmdArg& arg1=RooCmdArg::none(), const RooCmdArg& arg2=RooCmdArg::none(),
156  const RooCmdArg& arg3=RooCmdArg::none(), const RooCmdArg& arg4=RooCmdArg::none(), const RooCmdArg& arg5=RooCmdArg::none(),
157  const RooCmdArg& arg6=RooCmdArg::none(), const RooCmdArg& arg7=RooCmdArg::none(), const RooCmdArg& arg8=RooCmdArg::none()) ;
158  virtual RooFitResult* fitTo(RooAbsData& data, const RooLinkedList& cmdList) ;
159 
160  virtual RooAbsReal* createNLL(RooAbsData& data, const RooLinkedList& cmdList) ;
161  virtual RooAbsReal* createNLL(RooAbsData& data, const RooCmdArg& arg1=RooCmdArg::none(), const RooCmdArg& arg2=RooCmdArg::none(),
162  const RooCmdArg& arg3=RooCmdArg::none(), const RooCmdArg& arg4=RooCmdArg::none(), const RooCmdArg& arg5=RooCmdArg::none(),
163  const RooCmdArg& arg6=RooCmdArg::none(), const RooCmdArg& arg7=RooCmdArg::none(), const RooCmdArg& arg8=RooCmdArg::none()) ;
164 
165  // Chi^2 fits to histograms
166  using RooAbsReal::chi2FitTo ;
167  using RooAbsReal::createChi2 ;
168  virtual RooFitResult* chi2FitTo(RooDataHist& data, const RooLinkedList& cmdList) ;
169  virtual RooAbsReal* createChi2(RooDataHist& data, const RooCmdArg& arg1=RooCmdArg::none(), const RooCmdArg& arg2=RooCmdArg::none(),
170  const RooCmdArg& arg3=RooCmdArg::none(), const RooCmdArg& arg4=RooCmdArg::none(), const RooCmdArg& arg5=RooCmdArg::none(),
171  const RooCmdArg& arg6=RooCmdArg::none(), const RooCmdArg& arg7=RooCmdArg::none(), const RooCmdArg& arg8=RooCmdArg::none()) ;
172 
173  // Chi^2 fits to X-Y datasets
174  virtual RooAbsReal* createChi2(RooDataSet& data, const RooLinkedList& cmdList) ;
175 
176 
177 
178 
179 
180  // Constraint management
181  virtual RooArgSet* getConstraints(const RooArgSet& /*observables*/, RooArgSet& /*constrainedParams*/, Bool_t /*stripDisconnected*/) const {
182  // Interface to retrieve constraint terms on this pdf. Default implementation returns null
183  return 0 ;
184  }
185  virtual RooArgSet* getAllConstraints(const RooArgSet& observables, RooArgSet& constrainedParams, Bool_t stripDisconnected=kTRUE) const ;
186 
187  // Project p.d.f into lower dimensional p.d.f
188  virtual RooAbsPdf* createProjection(const RooArgSet& iset) ;
189 
190  // Create cumulative density function from p.d.f
191  RooAbsReal* createCdf(const RooArgSet& iset, const RooArgSet& nset=RooArgSet()) ;
192  RooAbsReal* createCdf(const RooArgSet& iset, const RooCmdArg& arg1, const RooCmdArg& arg2=RooCmdArg::none(),
193  const RooCmdArg& arg3=RooCmdArg::none(), const RooCmdArg& arg4=RooCmdArg::none(),
194  const RooCmdArg& arg5=RooCmdArg::none(), const RooCmdArg& arg6=RooCmdArg::none(),
195  const RooCmdArg& arg7=RooCmdArg::none(), const RooCmdArg& arg8=RooCmdArg::none()) ;
196  RooAbsReal* createScanCdf(const RooArgSet& iset, const RooArgSet& nset, Int_t numScanBins, Int_t intOrder) ;
197 
198  // Function evaluation support
199  virtual Double_t getValV(const RooArgSet* set=0) const ;
200  virtual Double_t getLogVal(const RooArgSet* set=0) const ;
201 
202  RooSpan<const double> getValBatch(std::size_t begin, std::size_t batchSize,
203  const RooArgSet* normSet = nullptr) const final;
204  RooSpan<const double> getLogValBatch(std::size_t begin, std::size_t batchSize,
205  const RooArgSet* normSet = nullptr) const;
206 
207  /// \copydoc getNorm(const RooArgSet*) const
208  Double_t getNorm(const RooArgSet& nset) const {
209  return getNorm(&nset) ;
210  }
211  virtual Double_t getNorm(const RooArgSet* set=0) const ;
212 
213  virtual void resetErrorCounters(Int_t resetValue=10) ;
214  void setTraceCounter(Int_t value, Bool_t allNodes=kFALSE) ;
215 private:
216  Bool_t traceEvalPdf(Double_t value) const;
217 
218 public:
219 
220  Double_t analyticalIntegralWN(Int_t code, const RooArgSet* normSet, const char* rangeName=0) const ;
221 
222  /// Shows if a PDF is self-normalized, which means that no attempt is made to add a normalization term.
223  /// Always returns false, unless a PDF overrides this function.
224  virtual Bool_t selfNormalized() const {
225  return kFALSE ;
226  }
227 
228  // Support for extended maximum likelihood, switched off by default
229  enum ExtendMode { CanNotBeExtended, CanBeExtended, MustBeExtended } ;
230  virtual ExtendMode extendMode() const {
231  // Returns ability of p.d.f to provided extended likelihood terms. Possible
232  // answers are CanNotBeExtended, CanBeExtended or MustBeExtended. This
233  // default implementation always return CanNotBeExtended
234  return CanNotBeExtended ;
235  }
236  inline Bool_t canBeExtended() const {
237  // If true p.d.f can provide extended likelihood term
238  return (extendMode() != CanNotBeExtended) ;
239  }
240  inline Bool_t mustBeExtended() const {
241  // If true p.d.f must extended likelihood term
242  return (extendMode() == MustBeExtended) ;
243  }
244  virtual Double_t expectedEvents(const RooArgSet* nset) const ;
245  virtual Double_t expectedEvents(const RooArgSet& nset) const {
246  // Return expecteded number of p.d.fs to be used in calculated of extended likelihood
247  return expectedEvents(&nset) ;
248  }
249 
250  // Printing interface (human readable)
251  virtual void printValue(std::ostream& os) const ;
252  virtual void printMultiline(std::ostream& os, Int_t contents, Bool_t verbose=kFALSE, TString indent="") const ;
253 
254  static void verboseEval(Int_t stat) ;
255  static int verboseEval() ;
256 
257  virtual Double_t extendedTerm(Double_t observedEvents, const RooArgSet* nset=0) const ;
258 
259  void setNormRange(const char* rangeName) ;
260  const char* normRange() const {
261  return _normRange.Length()>0 ? _normRange.Data() : 0 ;
262  }
263  void setNormRangeOverride(const char* rangeName) ;
264 
265  const RooAbsReal* getNormIntegral(const RooArgSet& nset) const { return getNormObj(0,&nset,0) ; }
266 
267 protected:
268 
269 public:
270  virtual const RooAbsReal* getNormObj(const RooArgSet* set, const RooArgSet* iset, const TNamed* rangeName=0) const ;
271 
272  virtual RooAbsGenContext* binnedGenContext(const RooArgSet &vars, Bool_t verbose= kFALSE) const ;
273 
274  virtual RooAbsGenContext* genContext(const RooArgSet &vars, const RooDataSet *prototype=0,
275  const RooArgSet* auxProto=0, Bool_t verbose= kFALSE) const ;
276 
277  virtual RooAbsGenContext* autoGenContext(const RooArgSet &vars, const RooDataSet* prototype=0, const RooArgSet* auxProto=0,
278  Bool_t verbose=kFALSE, Bool_t autoBinned=kTRUE, const char* binnedTag="") const ;
279 
280 private:
281 
282  RooDataSet *generate(RooAbsGenContext& context, const RooArgSet& whatVars, const RooDataSet* prototype,
283  Double_t nEvents, Bool_t verbose, Bool_t randProtoOrder, Bool_t resampleProto, Bool_t skipInit=kFALSE,
284  Bool_t extended=kFALSE) const ;
285 
286  // Implementation version
287  virtual RooPlot* paramOn(RooPlot* frame, const RooArgSet& params, Bool_t showConstants=kFALSE,
288  const char *label= "", Int_t sigDigits = 2, Option_t *options = "NELU", Double_t xmin=0.65,
289  Double_t xmax= 0.99,Double_t ymax=0.95, const RooCmdArg* formatCmd=0) ;
290 
291  void logBatchComputationErrors(RooSpan<const double>& outputs, std::size_t begin) const;
292 
293 protected:
294  virtual RooPlot *plotOn(RooPlot *frame, PlotOpt o) const;
295 
296  friend class RooEffGenContext ;
297  friend class RooAddGenContext ;
298  friend class RooProdGenContext ;
299  friend class RooSimGenContext ;
300  friend class RooSimSplitGenContext ;
301  friend class RooConvGenContext ;
302  friend class RooSimultaneous ;
303  friend class RooAddGenContextOrig ;
304  friend class RooProdPdf ;
305  friend class RooMCStudy ;
306 
307  Int_t* randomizeProtoOrder(Int_t nProto,Int_t nGen,Bool_t resample=kFALSE) const ;
308 
309  friend class RooExtendPdf ;
310  // This also forces the definition of a copy ctor in derived classes
311  RooAbsPdf(const RooAbsPdf& other, const char* name = 0);
312 
313  friend class RooRealIntegral ;
314  static Int_t _verboseEval ;
315 
316  virtual Bool_t syncNormalization(const RooArgSet* dset, Bool_t adjustProxies=kTRUE) const ;
317 
318  friend class RooAbsAnaConvPdf ;
319  mutable Double_t _rawValue ;
320  mutable RooAbsReal* _norm ; //! Normalization integral (owned by _normMgr)
321  mutable RooArgSet* _normSet ; //! Normalization set with for above integral
322 
323  class CacheElem : public RooAbsCacheElement {
324  public:
325  CacheElem(RooAbsReal& norm) : _norm(&norm) {} ;
326  void operModeHook(RooAbsArg::OperMode) ;
327  virtual ~CacheElem() ;
328  virtual RooArgList containedArgs(Action) { return RooArgList(*_norm) ; }
329  RooAbsReal* _norm ;
330  } ;
331  mutable RooObjCacheManager _normMgr ; // The cache manager
332 
333  friend class CacheElem ; // Cache needs to be able to clear _norm pointer
334 
335  virtual Bool_t redirectServersHook(const RooAbsCollection&, Bool_t, Bool_t, Bool_t) {
336  // Hook function intercepting redirectServer calls. Discard current normalization
337  // object if any server is redirected
338 
339  // Object is own by _normCacheManager that will delete object as soon as cache
340  // is sterilized by server redirect
341  _norm = 0 ;
342  return kFALSE ;
343  } ;
344 
345 
346  mutable Int_t _errorCount ; // Number of errors remaining to print
347  mutable Int_t _traceCount ; // Number of traces remaining to print
348  mutable Int_t _negCount ; // Number of negative probablities remaining to print
349 
350  Bool_t _selectComp ; // Component selection flag for RooAbsPdf::plotCompOn
351 
352  RooNumGenConfig* _specGeneratorConfig ; //! MC generator configuration specific for this object
353 
354  TString _normRange ; // Normalization range
355  static TString _normRangeOverride ;
356 
357  ClassDef(RooAbsPdf,4) // Abstract PDF with normalization support
358 };
359 
360 
361 #endif
RooAbsPdf::createScanCdf
RooAbsReal * createScanCdf(const RooArgSet &iset, const RooArgSet &nset, Int_t numScanBins, Int_t intOrder)
Definition: RooAbsPdf.cxx:3484
RooAbsPdf::_normSet
RooArgSet * _normSet
Normalization integral (owned by _normMgr)
Definition: RooAbsPdf.h:321
RooAbsPdf::_normRange
TString _normRange
MC generator configuration specific for this object.
Definition: RooAbsPdf.h:354
RooAbsPdf::getNormIntegral
const RooAbsReal * getNormIntegral(const RooArgSet &nset) const
Definition: RooAbsPdf.h:265
RooCmdArg
Definition: RooCmdArg.h:27
RooAbsPdf::genContext
virtual RooAbsGenContext * genContext(const RooArgSet &vars, const RooDataSet *prototype=0, const RooArgSet *auxProto=0, Bool_t verbose=kFALSE) const
Interface function to create a generator context from a p.d.f.
Definition: RooAbsPdf.cxx:2005
ymax
float ymax
Definition: THbookFile.cxx:95
RooAbsPdf::generateSimGlobal
virtual RooDataSet * generateSimGlobal(const RooArgSet &whatVars, Int_t nEvents)
Special generator interface for generation of 'global observables' – for RooStats tools.
Definition: RooAbsPdf.cxx:2671
kTRUE
const Bool_t kTRUE
Definition: RtypesCore.h:91
RooAbsReal.h
RooAbsPdf::traceEvalPdf
Bool_t traceEvalPdf(Double_t value) const
Check that passed value is positive and not 'not-a-number'.
Definition: RooAbsPdf.cxx:428
RooSimSplitGenContext
Definition: RooSimSplitGenContext.h:27
RooFit::NumEvents
RooCmdArg NumEvents(Int_t numEvents)
Definition: RooGlobalFunc.cxx:230
RooAbsPdf::autoGenContext
virtual RooAbsGenContext * autoGenContext(const RooArgSet &vars, const RooDataSet *prototype=0, const RooArgSet *auxProto=0, Bool_t verbose=kFALSE, Bool_t autoBinned=kTRUE, const char *binnedTag="") const
Definition: RooAbsPdf.cxx:2014
RooAbsData
Definition: RooAbsData.h:46
RooAbsPdf::setNormRangeOverride
void setNormRangeOverride(const char *rangeName)
Definition: RooAbsPdf.cxx:3636
RooAbsPdf::generateEvent
virtual void generateEvent(Int_t code)
Interface for generation of an event using the algorithm corresponding to the specified code.
Definition: RooAbsPdf.cxx:2435
RooAbsPdf::getLogValBatch
RooSpan< const double > getLogValBatch(std::size_t begin, std::size_t batchSize, const RooArgSet *normSet=nullptr) const
Compute the log-likelihoods for all events in the requested batch.
Definition: RooAbsPdf.cxx:740
RooAbsPdf::extendMode
virtual ExtendMode extendMode() const
Definition: RooAbsPdf.h:230
TH2F
2-D histogram with a float per channel (see TH1 documentation)}
Definition: TH2.h:251
RooProdGenContext
Definition: RooProdGenContext.h:31
TString::Data
const char * Data() const
Definition: TString.h:369
RooAbsPdf::CanBeExtended
@ CanBeExtended
Definition: RooAbsPdf.h:229
xmax
float xmax
Definition: THbookFile.cxx:95
RooMCStudy
Definition: RooMCStudy.h:32
RooAbsPdf::binnedGenContext
virtual RooAbsGenContext * binnedGenContext(const RooArgSet &vars, Bool_t verbose=kFALSE) const
Return a binned generator context.
Definition: RooAbsPdf.cxx:1995
RooAbsPdf::CacheElem::containedArgs
virtual RooArgList containedArgs(Action)
Definition: RooAbsPdf.h:328
RooAbsPdf::createCdf
RooAbsReal * createCdf(const RooArgSet &iset, const RooArgSet &nset=RooArgSet())
Create a cumulative distribution function of this p.d.f in terms of the observables listed in iset.
Definition: RooAbsPdf.cxx:3406
RooArgList
Definition: RooArgList.h:21
BatchHelpers::init
EvaluateInfo init(std::vector< RooRealProxy > parameters, std::vector< ArrayWrapper * > wrappers, std::vector< double * > arrays, size_t begin, size_t batchSize)
RooAbsPdf::getValV
virtual Double_t getValV(const RooArgSet *set=0) const
Return current value, normalized by integrating over the observables in nset.
Definition: RooAbsPdf.cxx:278
RooAbsPdf::printMultiline
virtual void printMultiline(std::ostream &os, Int_t contents, Bool_t verbose=kFALSE, TString indent="") const
Print multi line detailed information of this RooAbsPdf.
Definition: RooAbsPdf.cxx:1978
RooAbsPdf::RooAbsPdf
RooAbsPdf()
Default constructor.
Definition: RooAbsPdf.cxx:202
RooAbsPdf::GenSpec::_protoData
RooDataSet * _protoData
Definition: RooAbsPdf.h:82
RooAbsPdf::RooAddGenContextOrig
friend class RooAddGenContextOrig
Definition: RooAbsPdf.h:303
RooExtendPdf
Definition: RooExtendPdf.h:22
RooAbsPdf::specialGeneratorConfig
RooNumGenConfig * specialGeneratorConfig() const
Returns the specialized integrator configuration for this RooAbsReal.
Definition: RooAbsPdf.cxx:3533
RooAbsPdf::CacheElem::_norm
RooAbsReal * _norm
Definition: RooAbsPdf.h:329
RooAbsPdf::GenSpec
Definition: RooAbsPdf.h:70
RooAbsPdf::getGenerator
virtual Int_t getGenerator(const RooArgSet &directVars, RooArgSet &generateVars, Bool_t staticInitOK=kTRUE) const
Load generatedVars with the subset of directVars that we can generate events for, and return a code t...
Definition: RooAbsPdf.cxx:2413
TString::Length
Ssiz_t Length() const
Definition: TString.h:410
RooAbsPdf::CacheElem
Normalization set with for above integral.
Definition: RooAbsPdf.h:323
indent
static void indent(ostringstream &buf, int indent_level)
Definition: TClingCallFunc.cxx:87
RooNameSet.h
RooAbsPdf::defaultGeneratorConfig
static RooNumGenConfig * defaultGeneratorConfig()
Returns the default numeric MC generator configuration for all RooAbsReals.
Definition: RooAbsPdf.cxx:3523
RooAbsReal
Definition: RooAbsReal.h:61
RooAbsPdf::initGenerator
virtual void initGenerator(Int_t code)
Interface for one-time initialization to setup the generator for the specified code.
Definition: RooAbsPdf.cxx:2423
RooAbsPdf::_normMgr
RooObjCacheManager _normMgr
Definition: RooAbsPdf.h:331
RooAbsPdf::chi2FitTo
virtual RooFitResult * chi2FitTo(RooDataHist &data, const RooLinkedList &cmdList)
Calls RooAbsPdf::createChi2(RooDataSet& data, const RooLinkedList& cmdList) and returns fit result.
Definition: RooAbsPdf.cxx:1812
RooAbsPdf::GenSpec::_randProto
Bool_t _randProto
Definition: RooAbsPdf.h:85
RooAbsCacheElement
Definition: RooAbsCacheElement.h:26
TString
Definition: TString.h:136
RooSimGenContext
Definition: RooSimGenContext.h:27
RooCmdArg::none
static const RooCmdArg & none()
Return reference to null argument.
Definition: RooCmdArg.cxx:52
RooAbsPdf::_traceCount
Int_t _traceCount
Definition: RooAbsPdf.h:347
RooAbsPdf::createChi2
virtual RooAbsReal * createChi2(RooDataHist &data, 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())
Create a from a histogram and this function.
Definition: RooAbsPdf.cxx:1857
RooAbsPdf::isDirectGenSafe
virtual Bool_t isDirectGenSafe(const RooAbsArg &arg) const
Check if given observable can be safely generated using the pdfs internal generator mechanism (if tha...
Definition: RooAbsPdf.cxx:2448
RooAbsPdf::GenSpec::GenSpec
GenSpec()
Definition: RooAbsPdf.h:73
RooAbsPdf::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 RooCmdArg &arg9=RooCmdArg::none(), const RooCmdArg &arg10=RooCmdArg::none()) const
Helper calling plotOn(RooPlot*, RooLinkedList&) const.
Definition: RooAbsPdf.h:118
bool
RooFitResult
Definition: RooFitResult.h:40
RooAbsPdf::verboseEval
static int verboseEval()
Return global level of verbosity for p.d.f. evaluations.
Definition: RooAbsPdf.cxx:3329
RooAbsPdf::setTraceCounter
void setTraceCounter(Int_t value, Bool_t allNodes=kFALSE)
Reset trace counter to given value, limiting the number of future trace messages for this pdf to 'val...
Definition: RooAbsPdf.cxx:633
RooRealIntegral
Definition: RooRealIntegral.h:34
RooAbsPdf::getNormObj
virtual const RooAbsReal * getNormObj(const RooArgSet *set, const RooArgSet *iset, const TNamed *rangeName=0) const
Return pointer to RooAbsReal object that implements calculation of integral over observables iset in ...
Definition: RooAbsPdf.cxx:487
RooAbsPdf::setGeneratorConfig
void setGeneratorConfig()
Remove the specialized numeric MC generator configuration associated with this object.
Definition: RooAbsPdf.cxx:3588
RooAbsPdf::_errorCount
Int_t _errorCount
Definition: RooAbsPdf.h:343
RooAbsPdf::extendedTerm
virtual Double_t extendedTerm(Double_t observedEvents, const RooArgSet *nset=0) const
Return the extended likelihood term ( ) of this PDF for the given number of observed events.
Definition: RooAbsPdf.cxx:777
RooAbsPdf::GenSpec::_whatVars
RooArgSet _whatVars
Definition: RooAbsPdf.h:81
RooAbsPdf::ExtendMode
ExtendMode
Definition: RooAbsPdf.h:229
RooAbsPdf::CacheElem::CacheElem
CacheElem(RooAbsReal &norm)
Definition: RooAbsPdf.h:325
RooAbsPdf::MustBeExtended
@ MustBeExtended
Definition: RooAbsPdf.h:229
RooAbsPdf::~RooAbsPdf
virtual ~RooAbsPdf()
Destructor.
Definition: RooAbsPdf.cxx:260
RooAbsPdf::setNormRange
void setNormRange(const char *rangeName)
Definition: RooAbsPdf.cxx:3619
RooDataHist
Definition: RooDataHist.h:39
RooAbsPdf::_rawValue
Double_t _rawValue
Definition: RooAbsPdf.h:319
Option_t
const typedef char Option_t
Definition: RtypesCore.h:66
RooAbsPdf::CanNotBeExtended
@ CanNotBeExtended
Definition: RooAbsPdf.h:229
RooAbsPdf::getNorm
Double_t getNorm(const RooArgSet &nset) const
Get normalisation term needed to normalise the raw values returned by getVal().
Definition: RooAbsPdf.h:208
RooAbsPdf::GenSpec::~GenSpec
virtual ~GenSpec()
Definition: RooAbsPdf.cxx:3600
RooAbsPdf::GenSpec::_genContext
RooAbsGenContext * _genContext
Definition: RooAbsPdf.h:80
RooAbsPdf::getValBatch
RooSpan< const double > getValBatch(std::size_t begin, std::size_t batchSize, const RooArgSet *normSet=nullptr) const final
Compute batch of values for given range, and normalise by integrating over the observables in nset.
Definition: RooAbsPdf.cxx:345
xmin
float xmin
Definition: THbookFile.cxx:95
RooAbsPdf::paramOn
virtual RooPlot * paramOn(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())
Add a box with parameter values (and errors) to the specified frame.
Definition: RooAbsPdf.cxx:3152
TNamed
Definition: TNamed.h:29
RooAbsReal::createChi2
virtual RooAbsReal * createChi2(RooDataHist &data, const RooLinkedList &cmdList)
Create a variable from a histogram and this function.
Definition: RooAbsReal.cxx:4449
RooAbsPdf::prepareMultiGen
GenSpec * prepareMultiGen(const RooArgSet &whatVars, 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())
Prepare GenSpec configuration object for efficient generation of multiple datasets from identical spe...
Definition: RooAbsPdf.cxx:2177
RooConvGenContext
Definition: RooConvGenContext.h:31
kFALSE
const Bool_t kFALSE
Definition: RtypesCore.h:92
RooAbsPdf::analyticalIntegralWN
Double_t analyticalIntegralWN(Int_t code, const RooArgSet *normSet, const char *rangeName=0) const
Analytical integral with normalization (see RooAbsReal::analyticalIntegralWN() for further informatio...
Definition: RooAbsPdf.cxx:408
RooAbsPdf::createNLL
virtual RooAbsReal * createNLL(RooAbsData &data, const RooLinkedList &cmdList)
Construct representation of -log(L) of PDFwith given dataset.
Definition: RooAbsPdf.cxx:917
RooAbsPdf::getGeneratorConfig
const RooNumGenConfig * getGeneratorConfig() const
Return the numeric MC generator configuration used for this object.
Definition: RooAbsPdf.cxx:3561
RooLinkedList
Definition: RooLinkedList.h:35
RooAbsPdf::_normRangeOverride
static TString _normRangeOverride
Definition: RooAbsPdf.h:355
RooAbsPdf::CacheElem::~CacheElem
virtual ~CacheElem()
Destructor of normalization cache element.
Definition: RooAbsPdf.cxx:3350
RooAbsPdf::normRange
const char * normRange() const
Definition: RooAbsPdf.h:260
RooAbsCollection
Definition: RooAbsCollection.h:30
RooPlot
Definition: RooPlot.h:44
RooAbsArg::RooArgSet
friend class RooArgSet
Definition: RooAbsArg.h:579
RooAbsGenContext
Definition: RooAbsGenContext.h:26
RooAbsPdf::getAllConstraints
virtual RooArgSet * getAllConstraints(const RooArgSet &observables, RooArgSet &constrainedParams, Bool_t stripDisconnected=kTRUE) const
This helper function finds and collects all constraints terms of all component p.d....
Definition: RooAbsPdf.cxx:3501
RooAbsReal::chi2FitTo
virtual RooFitResult * chi2FitTo(RooDataHist &data, 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())
Perform a fit to given histogram.
Definition: RooAbsReal.cxx:4381
RooAddGenContext
Definition: RooAddGenContext.h:32
RooAbsPdf::logBatchComputationErrors
void logBatchComputationErrors(RooSpan< const double > &outputs, std::size_t begin) const
Scan through outputs and fix+log all nans and negative values.
Definition: RooAbsPdf.cxx:717
RooAbsPdf::CacheElem::operModeHook
void operModeHook(RooAbsArg::OperMode)
Dummy implementation.
Definition: RooAbsPdf.cxx:3339
RooCmdArg.h
RooObjCacheManager.h
RooAbsReal::plotOn
virtual RooPlot * plotOn(RooPlot *frame, const RooCmdArg &arg1=RooCmdArg(), 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(), const RooCmdArg &arg9=RooCmdArg(), const RooCmdArg &arg10=RooCmdArg()) const
Plot (project) PDF on specified frame.
Definition: RooAbsReal.cxx:1714
RooAbsReal::maxVal
virtual Double_t maxVal(Int_t code) const
Return maximum value for set of observables identified by code assigned in getMaxVal.
Definition: RooAbsReal.cxx:3662
RooAbsPdf::getConstraints
virtual RooArgSet * getConstraints(const RooArgSet &, RooArgSet &, Bool_t) const
Definition: RooAbsPdf.h:181
RooAbsPdf::printValue
virtual void printValue(std::ostream &os) const
Print value of p.d.f, also print normalization integral that was last used, if any.
Definition: RooAbsPdf.cxx:1962
RooAbsPdf::_negCount
Int_t _negCount
Definition: RooAbsPdf.h:348
RooAbsArg::OperMode
OperMode
Definition: RooAbsArg.h:491
RooAbsPdf::_selectComp
Bool_t _selectComp
Definition: RooAbsPdf.h:350
Double_t
double Double_t
Definition: RtypesCore.h:59
RooAbsPdf::mustBeExtended
Bool_t mustBeExtended() const
Definition: RooAbsPdf.h:240
RooNumGenConfig
Definition: RooNumGenConfig.h:25
RooAbsPdf::GenSpec::_dsetName
TString _dsetName
Definition: RooAbsPdf.h:87
RooAbsPdf::createProjection
virtual RooAbsPdf * createProjection(const RooArgSet &iset)
Return a p.d.f that represent a projection of this p.d.f integrated over given observables.
Definition: RooAbsPdf.cxx:3368
RooAbsPdf::GenSpec::_nGen
Int_t _nGen
Definition: RooAbsPdf.h:83
RooAbsPdf::GenSpec::_init
Bool_t _init
Definition: RooAbsPdf.h:88
RooCategory
Definition: RooCategory.h:27
RooAbsPdf::_norm
RooAbsReal * _norm
Definition: RooAbsPdf.h:320
RooAbsPdf::_specGeneratorConfig
RooNumGenConfig * _specGeneratorConfig
Definition: RooAbsPdf.h:352
RooAbsPdf::expectedEvents
virtual Double_t expectedEvents(const RooArgSet &nset) const
Definition: RooAbsPdf.h:245
RooAbsPdf::generate
RooDataSet * generate(const RooArgSet &whatVars, Int_t nEvents, const RooCmdArg &arg1, const RooCmdArg &arg2=RooCmdArg::none(), const RooCmdArg &arg3=RooCmdArg::none(), const RooCmdArg &arg4=RooCmdArg::none(), const RooCmdArg &arg5=RooCmdArg::none())
See RooAbsPdf::generate(const RooArgSet&,const RooCmdArg&,const RooCmdArg&,const RooCmdArg&,...
Definition: RooAbsPdf.h:55
TH1F
1-D histogram with a float per channel (see TH1 documentation)}
Definition: TH1.h:572
RooObjCacheManager
Definition: RooObjCacheManager.h:32
RooAbsCacheElement::Action
Action
Definition: RooAbsCacheElement.h:39
RooAbsPdf::GenSpec::_extended
Bool_t _extended
Definition: RooAbsPdf.h:84
ClassDef
#define ClassDef(name, id)
Definition: Rtypes.h:325
RooAbsPdf::_verboseEval
static Int_t _verboseEval
Definition: RooAbsPdf.h:314
name
char name[80]
Definition: TGX11.cxx:110
RooAbsPdf::redirectServersHook
virtual Bool_t redirectServersHook(const RooAbsCollection &, Bool_t, Bool_t, Bool_t)
Definition: RooAbsPdf.h:335
genreflex::verbose
bool verbose
Definition: rootcling_impl.cxx:133
RooDataSet
Definition: RooDataSet.h:33
RooAbsPdf::GenSpec::_resampleProto
Bool_t _resampleProto
Definition: RooAbsPdf.h:86
RooAbsPdf::generateBinned
virtual RooDataHist * generateBinned(const RooArgSet &whatVars, Double_t nEvents, const RooCmdArg &arg1, const RooCmdArg &arg2=RooCmdArg::none(), const RooCmdArg &arg3=RooCmdArg::none(), const RooCmdArg &arg4=RooCmdArg::none(), const RooCmdArg &arg5=RooCmdArg::none()) const
As RooAbsPdf::generateBinned(const RooArgSet&, const RooCmdArg&,const RooCmdArg&, const RooCmdArg&,...
Definition: RooAbsPdf.h:104
TPaveText
Definition: TPaveText.h:21
RooAbsArg
Definition: RooAbsArg.h:73
RooAbsPdf
Definition: RooAbsPdf.h:40
RooAbsPdf::selfNormalized
virtual Bool_t selfNormalized() const
Shows if a PDF is self-normalized, which means that no attempt is made to add a normalization term.
Definition: RooAbsPdf.h:224
RooSimultaneous
Definition: RooSimultaneous.h:37
RooAbsPdf::syncNormalization
virtual Bool_t syncNormalization(const RooArgSet *dset, Bool_t adjustProxies=kTRUE) const
Verify that the normalization integral cached with this PDF is valid for given set of normalization o...
Definition: RooAbsPdf.cxx:522
RooAbsPdf::getLogVal
virtual Double_t getLogVal(const RooArgSet *set=0) const
Return the log of the current value with given normalization An error message is printed if the argum...
Definition: RooAbsPdf.cxx:659
RooProdPdf
Definition: RooProdPdf.h:33
RooAbsPdf::resetErrorCounters
virtual void resetErrorCounters(Int_t resetValue=10)
Reset error counter to given value, limiting the number of future error messages for this pdf to 'res...
Definition: RooAbsPdf.cxx:621
RooAbsPdf::randomizeProtoOrder
Int_t * randomizeProtoOrder(Int_t nProto, Int_t nGen, Bool_t resample=kFALSE) const
Return lookup table with randomized access order for prototype events, given nProto prototype data ev...
Definition: RooAbsPdf.cxx:2369
RooEffGenContext
Definition: RooEffGenContext.h:23
RooAbsPdf::canBeExtended
Bool_t canBeExtended() const
Definition: RooAbsPdf.h:236
RooAbsAnaConvPdf
Definition: RooAbsAnaConvPdf.h:34
TList
Definition: TList.h:44
RooSpan
A simple container to hold a batch of data values.
Definition: RooSpan.h:33
RooArgSet
Definition: RooArgSet.h:28
int
RooAbsPdf::expectedEvents
virtual Double_t expectedEvents(const RooArgSet *nset) const
Return expected number of events from this p.d.f for use in extended likelihood calculations.
Definition: RooAbsPdf.cxx:3309
RooAbsPdf::fitTo
virtual RooFitResult * fitTo(RooAbsData &data, 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())
Fit PDF to given dataset.
Definition: RooAbsPdf.cxx:1261