ROOT  6.06/09
Reference Guide
RooAbsAnaConvPdf.h
Go to the documentation of this file.
1 /*****************************************************************************
2  * Project: RooFit *
3  * Package: RooFitCore *
4  * File: $Id: RooAbsAnaConvPdf.h,v 1.8 2007/07/16 21:04:28 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_ANA_CONV_PDF
17 #define ROO_ABS_ANA_CONV_PDF
18 
19 
20 class TIterator ;
21 #include "RooAbsPdf.h"
22 #include "RooRealProxy.h"
23 #include "RooListProxy.h"
24 #include "RooDataSet.h"
25 #include "RooAICRegistry.h"
26 #include "RooObjCacheManager.h"
27 #include "RooAbsCacheElement.h"
28 
29 class RooResolutionModel ;
30 class RooRealVar ;
31 class RooAbsGenContext ;
32 class RooConvGenContext ;
33 
34 class RooAbsAnaConvPdf : public RooAbsPdf {
35 public:
36 
37  // Constructors, assignment etc
39  RooAbsAnaConvPdf(const char *name, const char *title,
40  const RooResolutionModel& model,
42 
43  RooAbsAnaConvPdf(const RooAbsAnaConvPdf& other, const char* name=0);
44  virtual ~RooAbsAnaConvPdf();
45 
46  Int_t declareBasis(const char* expression, const RooArgList& params) ;
47  virtual void printMultiline(std::ostream& stream, Int_t contents, Bool_t verbose=kFALSE, TString indent= "") const ;
48 
49  // Coefficient normalization access
50  inline Double_t getCoefNorm(Int_t coefIdx, const RooArgSet& nset, const char* rangeName) const {
51  // Returns normalization integral for coefficient coefIdx for observables nset in range rangeNae
52  return getCoefNorm(coefIdx,&nset,rangeName) ;
53  }
54  Double_t getCoefNorm(Int_t coefIdx, const RooArgSet* nset=0, const char* rangeName=0) const {
55  return getCoefNorm(coefIdx,nset,RooNameReg::ptr(rangeName));
56  }
57 
58  // Analytical integration support
59  virtual Int_t getAnalyticalIntegralWN(RooArgSet& allVars, RooArgSet& analVars, const RooArgSet* normSet, const char* rangeName=0) const ;
60  virtual Double_t analyticalIntegralWN(Int_t code, const RooArgSet* normSet, const char* rangeName=0) const ;
61 
62  // Coefficient Analytical integration support
63  virtual Int_t getCoefAnalyticalIntegral(Int_t coef, RooArgSet& allVars, RooArgSet& analVars, const char* rangeName=0) const ;
64  virtual Double_t coefAnalyticalIntegral(Int_t coef, Int_t code, const char* rangeName=0) const ;
65  virtual Bool_t forceAnalyticalInt(const RooAbsArg& dep) const ;
66 
67  virtual Double_t coefficient(Int_t basisIndex) const = 0 ;
68  virtual RooArgSet* coefVars(Int_t coefIdx) const ;
69 
70  virtual Bool_t isDirectGenSafe(const RooAbsArg& arg) const ;
71 
72  virtual void setCacheAndTrackHints(RooArgSet&) ;
73 
74  virtual RooAbsGenContext* genContext(const RooArgSet &vars, const RooDataSet *prototype=0,
75  const RooArgSet* auxProto=0, Bool_t verbose= kFALSE) const ;
76  virtual Bool_t changeModel(const RooResolutionModel& newModel) ;
77 
78  const RooRealVar* convVar() const ; // Convolution variable
79 
80 protected:
81  Double_t getCoefNorm(Int_t coefIdx, const RooArgSet* nset, const TNamed* rangeName) const ;
82 
84 
85  virtual Double_t evaluate() const ;
86 
87  void makeCoefVarList(RooArgList&) const ;
88 
89  friend class RooConvGenContext ;
90 
91  RooRealProxy _model ; // Original model
92  RooRealProxy _convVar ; // Convolution variable
93 
94  RooArgSet* parseIntegrationRequest(const RooArgSet& intSet, Int_t& coefCode, RooArgSet* analVars=0) const ;
95 
96  RooListProxy _convSet ; // Set of (resModel (x) basisFunc) convolution objects
97  RooArgList _basisList ; //! List of created basis functions
98  mutable RooArgSet* _convNormSet ; //! Subset of last normalization that applies to convolutions
99  mutable TIterator* _convSetIter ; //! Iterator over _convNormSet
100 
101 
102  class CacheElem : public RooAbsCacheElement {
103  public:
104  virtual ~CacheElem() {} ;
105 
108  l.add(_normList) ;
109  return l ;
110  }
111 
114  } ;
115  mutable RooObjCacheManager _coefNormMgr ; // Coefficient normalization manager
116 
117  mutable RooAICRegistry _codeReg ; //! Registry of analytical integration codes
118 
119  ClassDef(RooAbsAnaConvPdf,2) // Abstract Composite Convoluted PDF
120 };
121 
122 #endif
RooArgList _basisList
virtual Bool_t add(const RooAbsArg &var, Bool_t silent=kFALSE)
Add the specified argument to list.
void makeCoefVarList(RooArgList &) const
Build complete list of coefficient variables.
virtual Double_t coefficient(Int_t basisIndex) const =0
RooAbsAnaConvPdf()
Default constructor, required for persistence.
virtual ~RooAbsAnaConvPdf()
Destructor.
const RooRealVar * convVar() const
Return a pointer to the convolution variable instance used in the resolution model.
Basic string class.
Definition: TString.h:137
int Int_t
Definition: RtypesCore.h:41
bool Bool_t
Definition: RtypesCore.h:59
virtual Bool_t forceAnalyticalInt(const RooAbsArg &dep) const
This function forces RooRealIntegral to offer all integration dependents to RooAbsAnaConvPdf::getAnal...
const Bool_t kFALSE
Definition: Rtypes.h:92
virtual Int_t getCoefAnalyticalIntegral(Int_t coef, RooArgSet &allVars, RooArgSet &analVars, const char *rangeName=0) const
Default implementation of function advertising integration capabilities.
Iterator abstract base class.
Definition: TIterator.h:32
virtual Int_t getAnalyticalIntegralWN(RooArgSet &allVars, RooArgSet &analVars, const RooArgSet *normSet, const char *rangeName=0) const
Advertise capability to perform (analytical) integrals internally.
virtual void setCacheAndTrackHints(RooArgSet &)
Label OK'ed components with cache-and-track.
#define ClassDef(name, id)
Definition: Rtypes.h:254
The TNamed class is the base class for all named ROOT classes.
Definition: TNamed.h:33
RooListProxy _convSet
RooArgSet * _convNormSet
List of created basis functions.
virtual Bool_t changeModel(const RooResolutionModel &newModel)
Change the current resolution model to newModel.
RooObjCacheManager _coefNormMgr
virtual Double_t coefAnalyticalIntegral(Int_t coef, Int_t code, const char *rangeName=0) const
Default implementation of function implementing advertised integrals.
static const TNamed * ptr(const char *stringPtr)
Return a unique TNamed pointer for given C++ string.
Definition: RooNameReg.cxx:124
RooRealProxy _model
virtual RooAbsGenContext * genContext(const RooArgSet &vars, const RooDataSet *prototype=0, const RooArgSet *auxProto=0, Bool_t verbose=kFALSE) const
Create a generator context for this p.d.f.
virtual void printMultiline(std::ostream &stream, Int_t contents, Bool_t verbose=kFALSE, TString indent="") const
Print info about this object to the specified stream.
bool verbose
TLine * l
Definition: textangle.C:4
static void indent(ostringstream &buf, int indent_level)
RooAICRegistry _codeReg
virtual Double_t analyticalIntegralWN(Int_t code, const RooArgSet *normSet, const char *rangeName=0) const
Return analytical integral defined by given code, which is returned by getAnalyticalIntegralWN() ...
virtual RooArgSet * coefVars(Int_t coefIdx) const
Return set of parameters with are used exclusively by the coefficient functions.
Iterator over _convNormSet.
double Double_t
Definition: RtypesCore.h:55
Int_t declareBasis(const char *expression, const RooArgList &params)
Declare a basis function for use in this physics model.
TIterator * _convSetIter
Subset of last normalization that applies to convolutions.
RooArgList containedArgs(Action)
#define name(a, b)
Definition: linkTestLib0.cpp:5
virtual Double_t evaluate() const
Calculate the current unnormalized value of the PDF.
RooAbsPdf is the abstract interface for all probability density functions The class provides hybrid a...
Definition: RooAbsPdf.h:41
virtual Bool_t isDirectGenSafe(const RooAbsArg &arg) const
Return true if it is safe to generate the convolution observable from the internal generator (this is...
RooAbsArg is the common abstract base class for objects that represent a value (of arbitrary type) an...
Definition: RooAbsArg.h:66
Double_t getCoefNorm(Int_t coefIdx, const RooArgSet *nset=0, const char *rangeName=0) const
Double_t getCoefNorm(Int_t coefIdx, const RooArgSet &nset, const char *rangeName) const
RooRealProxy _convVar
RooArgSet * parseIntegrationRequest(const RooArgSet &intSet, Int_t &coefCode, RooArgSet *analVars=0) const