Logo ROOT   6.30.04
Reference Guide
 All Namespaces Files Pages
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
38  RooAbsAnaConvPdf() ;
39  RooAbsAnaConvPdf(const char *name, const char *title,
40  const RooResolutionModel& model,
41  RooRealVar& convVar) ;
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 
83  Bool_t _isCopy ;
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 
100 
101  class CacheElem : public RooAbsCacheElement {
102  public:
103  virtual ~CacheElem() {} ;
104 
105  RooArgList containedArgs(Action) {
106  RooArgList l(_coefVarList) ;
107  l.add(_normList) ;
108  return l ;
109  }
110 
111  RooArgList _coefVarList ;
112  RooArgList _normList ;
113  } ;
114  mutable RooObjCacheManager _coefNormMgr ; // Coefficient normalization manager
115 
116  mutable RooAICRegistry _codeReg ; //! Registry of analytical integration codes
117 
118  ClassDef(RooAbsAnaConvPdf,2) // Abstract Composite Convoluted PDF
119 };
120 
121 #endif