Logo ROOT  
Reference Guide
 
Loading...
Searching...
No Matches
PiecewiseInterpolation.h
Go to the documentation of this file.
1/*****************************************************************************
2 * Project: RooFit *
3 * Package: RooFitCore *
4 * File: $Id: PiecewiseInterpolation.h,v 1.3 2007/05/11 09:11:30 verkerke 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_PIECEWISEINTERPOLATION
17#define ROO_PIECEWISEINTERPOLATION
18
19#include "RooAbsReal.h"
20#include "RooRealProxy.h"
21#include "RooListProxy.h"
22
23#include "RooObjCacheManager.h"
24#include <vector>
25#include <list>
26
27class RooRealVar;
28class RooArgList;
29
31public:
32
34 PiecewiseInterpolation(const char *name, const char *title, const RooAbsReal& nominal, const RooArgList& lowSet, const RooArgList& highSet, const RooArgList& paramSet, Bool_t takeOwnerShip=kFALSE) ;
35 virtual ~PiecewiseInterpolation() ;
36
37 PiecewiseInterpolation(const PiecewiseInterpolation& other, const char* name = 0);
38 virtual TObject* clone(const char* newname) const { return new PiecewiseInterpolation(*this, newname); }
39
40 // virtual Double_t defaultErrorLevel() const ;
41
42 // void printMetaArgs(std::ostream& os) const ;
43
44 const RooArgList& lowList() const { return _lowSet ; }
45 const RooArgList& highList() const { return _highSet ; }
46 const RooArgList& paramList() const { return _paramSet ; }
47
48 //virtual Bool_t forceAnalyticalInt(const RooAbsArg&) const { return kTRUE ; }
50
51 Int_t getAnalyticalIntegralWN(RooArgSet& allVars, RooArgSet& analVars, const RooArgSet* normSet,const char* rangeName=0) const ;
52 Double_t analyticalIntegralWN(Int_t code, const RooArgSet* normSet, const char* rangeName=0) const ;
53
54 void setPositiveDefinite(bool flag=true){_positiveDefinite=flag;}
55
56 void setInterpCode(RooAbsReal& param, int code);
57 void setAllInterpCodes(int code);
59
60 virtual std::list<Double_t>* binBoundaries(RooAbsRealLValue& /*obs*/, Double_t /*xlo*/, Double_t /*xhi*/) const ;
61 virtual std::list<Double_t>* plotSamplingHint(RooAbsRealLValue& obs, Double_t xlo, Double_t xhi) const ;
62 virtual Bool_t isBinnedDistribution(const RooArgSet& obs) const ;
63
64protected:
65
67 public:
68 CacheElem() {} ;
69 virtual ~CacheElem() {} ;
72 ret.add(_lowIntList);
73 ret.add(_highIntList);
74 return ret ;
75 }
79 // will want std::vector<RooRealVar*> for low and high also
80 } ;
81 mutable RooObjCacheManager _normIntMgr ; // The integration cache manager
82
83 RooRealProxy _nominal; // The nominal value
84 RooArgList _ownedList ; // List of owned components
85 RooListProxy _lowSet ; // Low-side variation
86 RooListProxy _highSet ; // High-side varaition
87 RooListProxy _paramSet ; // interpolation parameters
88 RooListProxy _normSet ; // interpolation parameters
89 Bool_t _positiveDefinite; // protect against negative and 0 bins.
90
91 std::vector<int> _interpCode;
92
93 Double_t evaluate() const;
94
95 ClassDef(PiecewiseInterpolation,3) // Sum of RooAbsReal objects
96};
97
98#endif
const Bool_t kFALSE
Definition RtypesCore.h:92
double Double_t
Definition RtypesCore.h:59
#define ClassDef(name, id)
Definition Rtypes.h:325
char name[80]
Definition TGX11.cxx:110
virtual RooArgList containedArgs(Action)
The PiecewiseInterpolation is a class that can morph distributions into each other,...
const RooArgList & highList() const
virtual std::list< Double_t > * binBoundaries(RooAbsRealLValue &, Double_t, Double_t) const
WVE note: assumes nominal and alternates have identical structure, must add explicit check.
Double_t analyticalIntegralWN(Int_t code, const RooArgSet *normSet, const char *rangeName=0) const
Implement analytical integrations by doing appropriate weighting from component integrals functions t...
const RooArgList & lowList() const
Bool_t setBinIntegrator(RooArgSet &allVars)
void setPositiveDefinite(bool flag=true)
Int_t getAnalyticalIntegralWN(RooArgSet &allVars, RooArgSet &analVars, const RooArgSet *normSet, const char *rangeName=0) const
Advertise that all integrals can be handled internally.
const RooArgList & paramList() const
Double_t evaluate() const
Calculate and return current value of self.
virtual std::list< Double_t > * plotSamplingHint(RooAbsRealLValue &obs, Double_t xlo, Double_t xhi) const
Interface for returning an optional hint for initial sampling points when constructing a curve projec...
virtual TObject * clone(const char *newname) const
void setInterpCode(RooAbsReal &param, int code)
virtual ~PiecewiseInterpolation()
Destructor.
virtual Bool_t isBinnedDistribution(const RooArgSet &obs) const
WVE note: assumes nominal and alternates have identical structure, must add explicit check.
RooAbsCacheElement is the abstract base class for objects to be stored in RooAbsCache cache manager o...
virtual Bool_t add(const RooAbsArg &var, Bool_t silent=kFALSE)
Add the specified argument to list.
RooAbsRealLValue is the common abstract base class for objects that represent a real value that may a...
RooAbsReal is the common abstract base class for objects that represent a real value and implements f...
Definition RooAbsReal.h:61
RooArgList is a container object that can hold multiple RooAbsArg objects.
Definition RooArgList.h:21
RooArgSet is a container object that can hold multiple RooAbsArg objects.
Definition RooArgSet.h:29
RooListProxy is the concrete proxy for RooArgList objects.
Class RooObjCacheManager is an implementation of class RooCacheManager<RooAbsCacheElement> and specia...
RooRealVar represents a variable that can be changed from the outside.
Definition RooRealVar.h:39
Mother of all ROOT objects.
Definition TObject.h:37