FastVerticalInterpHistPdf Class Reference

#include <VerticalInterpHistPdf.h>

Inheritance diagram for FastVerticalInterpHistPdf:

Public Member Functions
virtual TObject *	clone (const char *newname) const
Double_t	evaluate () const
	FastVerticalInterpHistPdf ()
	FastVerticalInterpHistPdf (const char *name, const char *title, const RooRealVar &x, const RooArgList &funcList, const RooArgList &coefList, Double_t smoothRegion=1., Int_t smoothAlgo=1)
	FastVerticalInterpHistPdf (const FastVerticalInterpHistPdf &other, const char *name=0)
Bool_t	hasCache () const
Bool_t	isCacheReady () const
virtual	~FastVerticalInterpHistPdf ()
Public Member Functions inherited from FastVerticalInterpHistPdfBase
const RooArgList &	coefList () const
	FastVerticalInterpHistPdfBase ()
	FastVerticalInterpHistPdfBase (const char *name, const char *title, const RooArgSet &obs, const RooArgList &funcList, const RooArgList &coefList, Double_t smoothRegion=1., Int_t smoothAlgo=1)
	FastVerticalInterpHistPdfBase (const FastVerticalInterpHistPdfBase &other, const char *name=0)
const RooArgList &	funcList () const
Bool_t	selfNormalized () const
virtual	~FastVerticalInterpHistPdfBase ()

Protected Member Functions
void	setupCaches () const
	not to be serialized More...
void	syncComponent (int which) const
void	syncComponents (int dimension) const
void	syncNominal () const
void	syncTotal () const
Protected Member Functions inherited from FastVerticalInterpHistPdfBase
double	smoothStepFunc (double x) const
void	syncMorph (Morph &out, const FastTemplate &nominal, FastTemplate &lo, FastTemplate &hi) const
	not to be serialized More...
void	syncTotal (FastTemplate &cache, const FastTemplate &cacheNominal, const FastTemplate &cacheNominalLog) const

Protected Attributes
FastHisto	_cache
	Cache of the result. More...
FastHisto	_cacheNominal
	Cache of nominal pdf (additive morphing) and its bin-by-bin logarithm (multiplicative) More...
FastHisto	_cacheNominalLog
	not to be serialized More...
RooRealProxy	_x
Protected Attributes inherited from FastVerticalInterpHistPdfBase
TIterator *	_coefIter
	Iterator over FUNC list. More...
RooListProxy	_coefList
TIterator *	_funcIter
RooListProxy	_funcList
bool	_init
	Iterator over coefficient list. More...
std::vector< RooAbsReal * >	_morphParams
	not to be serialized More...
std::vector< Morph >	_morphs
	not to be serialized More...
SimpleCacheSentry	_sentry
	not to be serialized More...
Int_t	_smoothAlgo
Double_t	_smoothRegion
RooRealProxy	_x

Detailed Description

Definition at line 125 of file VerticalInterpHistPdf.h.

Constructor & Destructor Documentation

FastVerticalInterpHistPdf::FastVerticalInterpHistPdf ( )

inline

Definition at line 128 of file VerticalInterpHistPdf.h.

Referenced by clone().

: FastVerticalInterpHistPdfBase() {}

FastVerticalInterpHistPdf::FastVerticalInterpHistPdf ( const char *  name, const char *  title, const RooRealVar &  x, const RooArgList &  funcList, const RooArgList &  coefList, Double_t  smoothRegion = 1., Int_t  smoothAlgo = 1  )

inline

Definition at line 129 of file VerticalInterpHistPdf.h.

                                                                                                                                                                                            :
    FastVerticalInterpHistPdfBase(name,title,RooArgSet(x),funcList,coefList,smoothRegion,smoothAlgo),
    _x("x","Independent variable",this,const_cast<RooRealVar&>(x)),
    _cache(), _cacheNominal(), _cacheNominalLog()  {}

FastVerticalInterpHistPdf::FastVerticalInterpHistPdf ( const FastVerticalInterpHistPdf &  other, const char *  name = 0  )

inline

Definition at line 133 of file VerticalInterpHistPdf.h.

                                                                                        :
    FastVerticalInterpHistPdfBase(other, name),
    _x("x",this,other._x),
#ifdef FastVerticalInterpHistPdf_CopyConstructor
    _cache(other._cache), _cacheNominal(other._cacheNominal), _cacheNominalLog(other._cacheNominalLog)  {}

virtual FastVerticalInterpHistPdf::~FastVerticalInterpHistPdf ( )

inlinevirtual

Definition at line 142 of file VerticalInterpHistPdf.h.

{}

Member Function Documentation

virtual TObject* FastVerticalInterpHistPdf::clone ( const char *  newname ) const

inlinevirtual

Implements FastVerticalInterpHistPdfBase.

Definition at line 141 of file VerticalInterpHistPdf.h.

References FastVerticalInterpHistPdf().

{ return new FastVerticalInterpHistPdf(*this,newname) ; }

Double_t FastVerticalInterpHistPdf::evaluate ( ) const

virtual

Implements FastVerticalInterpHistPdfBase.

Definition at line 324 of file VerticalInterpHistPdf.cc.

References _cache, FastVerticalInterpHistPdfBase::_init, FastVerticalInterpHistPdfBase::_sentry, _x, FastHisto::GetAt(), SimpleCacheSentry::good(), setupCaches(), FastTemplate::size(), syncTotal(), and TRACEME.

{
  TRACEME()
  if (_cache.size() == 0) setupCaches();
  if (!_sentry.good() || !_init) syncTotal();
  //return std::max<double>(1e-9, _cache.GetAt(_x));
  return _cache.GetAt(_x);
}

Bool_t FastVerticalInterpHistPdf::hasCache ( ) const

inline

Definition at line 146 of file VerticalInterpHistPdf.h.

References _cache, and FastTemplate::size().

{ return _cache.size() > 0; }

Bool_t FastVerticalInterpHistPdf::isCacheReady ( ) const

inline

Definition at line 147 of file VerticalInterpHistPdf.h.

References _cache, FastVerticalInterpHistPdfBase::_init, and FastTemplate::size().

{ return _cache.size() > 0 && _init; }

void FastVerticalInterpHistPdf::setupCaches ( ) const

protected

not to be serialized

Definition at line 489 of file VerticalInterpHistPdf.cc.

References _cache, _cacheNominal, FastVerticalInterpHistPdfBase::_coefIter, FastVerticalInterpHistPdfBase::_coefList, FastVerticalInterpHistPdfBase::_morphParams, FastVerticalInterpHistPdfBase::_morphs, FastVerticalInterpHistPdfBase::_sentry, SimpleCacheSentry::addVars(), SimpleCacheSentry::empty(), i, FastTemplate::size(), syncComponents(), syncNominal(), syncTotal(), and TRACEME.

Referenced by evaluate().

                                                  {
    TRACEME()
    int ndim = _coefList.getSize();

    _morphs.resize(ndim);
    _morphParams.resize(ndim);
    syncNominal();
    //printf("Nominal template has been set up: \n");  _cacheNominal.Dump();
    _coefIter->Reset();
    for (int i = 0; i < ndim; ++i) {
        _morphParams[i] = dynamic_cast<RooAbsReal *>(_coefIter->Next());
        _morphs[i].sum.Resize(_cacheNominal.size());
        _morphs[i].diff.Resize(_cacheNominal.size());
        syncComponents(i);
    } 
    _cache = FastHisto(_cacheNominal);

    if (_sentry.empty()) _sentry.addVars(_coefList); 
    syncTotal();
}

void FastVerticalInterpHistPdf::syncComponent ( int  which ) const

protected

void FastVerticalInterpHistPdf::syncComponents ( int  dimension ) const

protected

Definition at line 392 of file VerticalInterpHistPdf.cc.

References _cacheNominal, FastVerticalInterpHistPdfBase::_funcList, FastVerticalInterpHistPdfBase::_morphs, _x, FastHisto::Normalize(), FastVerticalInterpHistPdfBase::syncMorph(), TRACEME, and vdt::x.

Referenced by setupCaches().

                                                            {
    TRACEME()
    RooAbsPdf *pdfHi = dynamic_cast<RooAbsPdf *>(_funcList.at(2*dim+1));
    RooAbsPdf *pdfLo = dynamic_cast<RooAbsPdf *>(_funcList.at(2*dim+2));
    const RooRealVar &x = dynamic_cast<const RooRealVar &>(_x.arg());
    std::auto_ptr<TH1> histHi(pdfHi->createHistogram("",x)); histHi->SetDirectory(0); 
    std::auto_ptr<TH1> histLo(pdfLo->createHistogram("",x)); histLo->SetDirectory(0);

    FastHisto hi(*histHi), lo(*histLo); 
    //printf("Un-normalized templates for dimension %d: \n", dim);  hi.Dump(); lo.Dump();
    hi.Normalize(); lo.Normalize();
    //printf("Normalized templates for dimension %d: \n", dim);  hi.Dump(); lo.Dump();
    syncMorph(_morphs[dim], _cacheNominal, lo, hi);
}

void FastVerticalInterpHistPdf::syncNominal ( ) const

protected

Definition at line 346 of file VerticalInterpHistPdf.cc.

References _cacheNominal, _cacheNominalLog, FastVerticalInterpHistPdfBase::_funcList, FastVerticalInterpHistPdfBase::_smoothAlgo, _x, estimatePileup::hist, FastTemplate::Log(), FastHisto::Normalize(), TRACEME, and vdt::x.

Referenced by setupCaches().

                                                  {
    TRACEME()
    RooAbsPdf *pdf = dynamic_cast<RooAbsPdf *>(_funcList.at(0));
    const RooRealVar &x = dynamic_cast<const RooRealVar &>(_x.arg());
    std::auto_ptr<TH1> hist(pdf->createHistogram("",x));
    hist->SetDirectory(0); 
    _cacheNominal = FastHisto(*hist);
    _cacheNominal.Normalize();
    if (_smoothAlgo < 0) {
        _cacheNominalLog = _cacheNominal;
        _cacheNominalLog.Log();
    }
}

void FastVerticalInterpHistPdf::syncTotal ( ) const

protected

Definition at line 471 of file VerticalInterpHistPdf.cc.

References _cache, _cacheNominal, _cacheNominalLog, FastHisto::Normalize(), and FastVerticalInterpHistPdfBase::syncTotal().

Referenced by evaluate(), and setupCaches().

                                                {
    FastVerticalInterpHistPdfBase::syncTotal(_cache, _cacheNominal, _cacheNominalLog);

    // normalize the result
    _cache.Normalize(); 
    //printf("Normalized result\n");  _cache.Dump();
}

Member Data Documentation

FastHisto FastVerticalInterpHistPdf::_cache

mutableprotected

Cache of the result.

Definition at line 152 of file VerticalInterpHistPdf.h.

Referenced by evaluate(), hasCache(), isCacheReady(), setupCaches(), and syncTotal().

FastHisto FastVerticalInterpHistPdf::_cacheNominal

mutableprotected

Cache of nominal pdf (additive morphing) and its bin-by-bin logarithm (multiplicative)

not to be serialized

Definition at line 154 of file VerticalInterpHistPdf.h.

Referenced by setupCaches(), syncComponents(), syncNominal(), and syncTotal().

FastHisto FastVerticalInterpHistPdf::_cacheNominalLog

mutableprotected

not to be serialized

Definition at line 155 of file VerticalInterpHistPdf.h.

Referenced by syncNominal(), and syncTotal().

RooRealProxy FastVerticalInterpHistPdf::_x

protected

Definition at line 149 of file VerticalInterpHistPdf.h.

Referenced by evaluate(), syncComponents(), and syncNominal().