RooEffResModel Class Reference

#include <RooEffResModel.h>

Inheritance diagram for RooEffResModel:

Classes
class	CacheElem

Public Member Functions
	RooEffResModel ()
	RooEffResModel (const char *name, const char *title, RooResolutionModel &model, RooAbsReal &eff)
	RooEffResModel (const RooEffResModel &other, const char *name=0)
virtual RooEffResModel *	clone (const char *newname) const
virtual	~RooEffResModel ()
virtual Int_t	basisCode (const char *name) const
virtual Int_t	getAnalyticalIntegral (RooArgSet &allVars, RooArgSet &analVars, const char *rangeName=0) const
virtual Double_t	analyticalIntegral (Int_t code, const char *rangeName) const
virtual Bool_t	forceAnalyticalInt (const RooAbsArg &dep) const
virtual Int_t	getGenerator (const RooArgSet &directVars, RooArgSet &generateVars, Bool_t staticInitOK=kTRUE) const
virtual void	initGenerator (Int_t code)
virtual void	generateEvent (Int_t code)
virtual RooAbsGenContext *	modelGenContext (const RooAbsAnaConvPdf &convPdf, const RooArgSet &vars, const RooDataSet *prototype=0, const RooArgSet *auxProto=0, Bool_t verbose=kFALSE) const
virtual RooAbsReal *	efficiency () const
	Return pointer to pdf in product. More...
virtual std::vector< RooAbsReal * >	efficiencies () const
	Return pointer to pdf in product. More...
virtual RooResolutionModel &	model () const
virtual RooArgSet *	observables () const
	Return pointer to pdf in product. More...
const RooArgList &	getIntegralRanges (const RooArgSet &iset, const char *rangeName=0) const
Public Member Functions inherited from RooAbsEffResModel
virtual	~RooAbsEffResModel ()

Protected Member Functions
virtual Double_t	evaluate () const
virtual RooEffResModel *	convolution (RooFormulaVar *inBasis, RooAbsArg *owner) const

Private Types
typedef std::map< std::string, RooArgList * >	RangeMap

Private Member Functions
CacheElem *	getCache (const RooArgSet *iset, const TNamed *rangeName=0) const

Private Attributes
RooSetProxy	_observables
RooRealProxy	_model
RooRealProxy	_eff
RangeMap	_ranges
RooObjCacheManager	_cacheMgr

Friends
class	RooMultiEffResModel

Detailed Description

Definition at line 29 of file RooEffResModel.h.

Member Typedef Documentation

RangeMap

typedef std::map<std::string, RooArgList*> RooEffResModel::RangeMap

private

Definition at line 108 of file RooEffResModel.h.

Constructor & Destructor Documentation

RooEffResModel() [1/3]

RooEffResModel::RooEffResModel ( )

inline

Definition at line 32 of file RooEffResModel.h.

{ }

RooEffResModel() [2/3]

RooEffResModel::RooEffResModel	(	const char *	name,
		const char *	title,
		RooResolutionModel &	model,
		RooAbsReal &	eff
	)

Definition at line 188 of file RooEffResModel.cpp.

                                                      :
    RooResolutionModel(name,title,__model.convVar())
    , _observables("observables", "observables", this)
    , _model("!model","Original resolution model",this,__model)
    , _eff("!efficiency","efficiency of convvar", this,eff)
    , _cacheMgr(this)
{
    // FIXME: assert that efficiency is a function of convVar, and there are no overlaps...
    _observables.add(__model.convVar());
}

RooEffResModel() [3/3]

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

Definition at line 201 of file RooEffResModel.cpp.

                                                                            :
    RooResolutionModel(other,name)
    , _observables("observables", this, other._observables)
    , _model("!model", this, other._model)
    , _eff("!efficiency", this, other._eff)
    , _cacheMgr(other._cacheMgr, this)
{
    // copy constructor
}

~RooEffResModel()

RooEffResModel::~RooEffResModel ( )

virtual

Definition at line 215 of file RooEffResModel.cpp.

{
    // Destructor
}

Member Function Documentation

analyticalIntegral()

Double_t RooEffResModel::analyticalIntegral ( Int_t  code, const char *  rangeName  ) const

virtual

Definition at line 327 of file RooEffResModel.cpp.

{
    assert(code > 0);
    CacheElem* cache = static_cast<CacheElem*>(_cacheMgr.getObjByIndex(code - 1));
    if (!cache) {
        myunique_ptr<RooArgSet> vars(getParameters(RooArgSet()));
        myunique_ptr<RooArgSet> iset(
                _cacheMgr.nameSet1ByIndex(code - 1)->select(*vars));
        cache = getCache(iset.get(), RooNameReg::ptr(rangeName));
        assert(cache!=0);
    }
    return cache->getVal();
}

basisCode()

Int_t RooEffResModel::basisCode ( const char *  name ) const

virtual

Definition at line 293 of file RooEffResModel.cpp.

{
    return model().basisCode(name);
}

clone()

RooEffResModel * RooEffResModel::clone ( const char *  newname ) const

virtual

Definition at line 211 of file RooEffResModel.cpp.

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

convolution()

RooEffResModel * RooEffResModel::convolution ( RooFormulaVar *  inBasis, RooAbsArg *  owner  ) const

protectedvirtual

Definition at line 233 of file RooEffResModel.cpp.

{
    // Instantiate a clone of this resolution model representing a convolution with given
    // basis function. The owners object name is incorporated in the clones name
    // to avoid multiple convolution objects with the same name in complex PDF structures.
    //
    // Note: The 'inBasis' formula expression must be a RooFormulaVar that encodes the formula
    // in the title of the object and this expression must be an exact match against the
    // implemented basis function strings (see derived class implementation of method basisCode()
    // for those strings

    // Check that primary variable of basis functions is our convolution variable
    if (inBasis->getParameter(0) != x.absArg()) {
        coutE(InputArguments) << "RooEffResModel::convolution(" << GetName() << "," << this
            << ") convolution parameter of basis function and PDF don't match" << endl
            << "basis->findServer(0) = " << inBasis->findServer(0) << endl
            << "x.absArg()           = " << x.absArg() << endl;
        return 0;
    }

    if (basisCode(inBasis->GetTitle())==0) {
        coutE(InputArguments) << "RooEffResModel::convolution(" << GetName() << "," << this
            << ") basis function '" << inBasis->GetTitle() << "' is not supported." << endl;
        return 0;
    }

    std::string newName(GetName());
    newName += "_conv_";
    newName += inBasis->GetName();
    newName += "_[";
    newName += owner->GetName();
    newName += "]";

    RooResolutionModel *conv = model().convolution(inBasis, owner);

    std::string newTitle(conv->GetTitle());
    newTitle += " convoluted with basis function ";
    newTitle += inBasis->GetName();
    conv->SetTitle(newTitle.c_str());

    RooAbsReal* eff = efficiency();
    RooEffResModel *effConv = new RooEffResModel(newName.c_str(), newTitle.c_str(), *conv, *eff);
    effConv->addOwnedComponents(*conv);
    effConv->changeBasis(inBasis);

    const char* cacheParamsStr = getStringAttribute("CACHEPARAMINT");
    if (cacheParamsStr && strlen(cacheParamsStr)) {
        //cout << "prior has CACHEPARAMINT : " << cacheParamsStr << endl;
        //const char* ecCacheParamsStr = effConv->getStringAttribute("CACHEPARAMINT");
        //if (ecCacheParamsStr && strlen(ecCacheParamsStr)) cout << "bound version has CACHEPARAMINT : " << ecCacheParamsStr << endl;
        effConv->setStringAttribute("CACHEPARAMINT",cacheParamsStr);
        //cout << "2nd time: bound version has CACHEPARAMINT : " << effConv->getStringAttribute("CACHEPARAMINT")  << endl;
        //cout << endl << endl << endl;
    }

    return effConv;
}

efficiencies()

std::vector< RooAbsReal * > RooEffResModel::efficiencies ( ) const

virtual

Return pointer to pdf in product.

Definition at line 226 of file RooEffResModel.cpp.

{ return std::vector<RooAbsReal*>(1, efficiency()); }

efficiency()

RooAbsReal * RooEffResModel::efficiency ( ) const

virtual

Return pointer to pdf in product.

Implements RooAbsEffResModel.

Definition at line 223 of file RooEffResModel.cpp.

{ return static_cast<RooAbsReal*>(_eff.absArg()); }

evaluate()

Double_t RooEffResModel::evaluate ( ) const

protectedvirtual

Definition at line 299 of file RooEffResModel.cpp.

{
    Double_t mod = model().getVal();
    // TODO: replace this by the discretized version, i.e. replace convVar by
    // customized middle of bin...  this in order to ensure evaluate &
    // analyticalIntegral are consistent (in case eff is not discretized!!!)
    Double_t eps = efficiency()->getVal();
    return eps * mod;
}

forceAnalyticalInt()

Bool_t RooEffResModel::forceAnalyticalInt ( const RooAbsArg &  dep ) const

virtual

Definition at line 310 of file RooEffResModel.cpp.

{
    // Force RooRealIntegral to offer all observables for internal integration
    return kTRUE;
}

generateEvent()

void RooEffResModel::generateEvent ( Int_t  code )

virtual

Definition at line 366 of file RooEffResModel.cpp.

{
    // The hit-miss on the efficiency is done at the level of the GenContext.
    model().generateEvent(code);
}

getAnalyticalIntegral()

Int_t RooEffResModel::getAnalyticalIntegral ( RooArgSet &  allVars, RooArgSet &  analVars, const char *  rangeName = 0  ) const

virtual

Definition at line 317 of file RooEffResModel.cpp.

{
    if (_forceNumInt) return 0;
    analVars.add(allVars);
    getCache(&analVars, RooNameReg::ptr(rangeName));
    return 1 + _cacheMgr.lastIndex();
}

getCache()

RooEffResModel::CacheElem * RooEffResModel::getCache ( const RooArgSet *  iset, const TNamed *  rangeName = 0  ) const

private

Definition at line 428 of file RooEffResModel.cpp.

{
    Int_t sterileIndex(-1);
    CacheElem* cache = (CacheElem*) _cacheMgr.getObj(iset, &sterileIndex, rangeName);
    if (cache) return cache;
    _cacheMgr.setObj(iset, new CacheElem( *this,  *iset,  rangeName), rangeName);
    return getCache(iset, rangeName);
}

getGenerator()

Int_t RooEffResModel::getGenerator ( const RooArgSet &  directVars, RooArgSet &  generateVars, Bool_t  staticInitOK = kTRUE  ) const

virtual

Definition at line 352 of file RooEffResModel.cpp.

{
    return model().getGenerator(directVars, generateVars, staticInitOK);
}

getIntegralRanges()

const RooArgList & RooEffResModel::getIntegralRanges	(	const RooArgSet &	iset,
		const char *	rangeName = 0
	)		const

Definition at line 373 of file RooEffResModel.cpp.

{
   rangeName = rangeName ? rangeName : "default";
   RangeMap::const_iterator it = _ranges.find(rangeName);
   if (it != _ranges.end()) return *it->second;

   RooRealVar& x = static_cast<RooRealVar&>(convVar());
   const Double_t xmin = x.getMin(rangeName);
   const Double_t xmax = x.getMax(rangeName);

   RooArgList* ranges = new RooArgList;
   myunique_ptr<std::list<Double_t> > bounds(efficiency()->binBoundaries(x, x.getMin(), x.getMax()));
   if (!bounds.get()) {
       std::string s;
       s += "RooEffResModel("; s += GetName(); s += "): specified efficiency ";
       s += efficiency()->GetName(); s += " does not provide binBoundaries...";
       std::cout << s << std::endl;
       throw s;
   }
   std::list<Double_t>::const_iterator lo, hi = bounds->begin();
   std::stringstream strange;
   for (unsigned int i = 0; i + 1 < bounds->size();++i) {
      lo = hi++;
      if (*hi < xmin) continue; // not there yet...
      if (*lo > xmax) break;    // past the requested interval...
      const Double_t thisxmin = std::max(*lo, xmin);
      const Double_t thisxmax = std::min(*hi, xmax);

      // add eff name, as it specifies the boundaries...
      strange.str("R");
      strange << i << "_" << x.GetName() << "_" << efficiency()->GetName();

      // Add original rangeName if there is one
      if (rangeName) strange << "_" << rangeName;
      strange << "_I_" << RooNameSet(iset).content();
      const std::string trange = strange.str();
      const char* range = trange.c_str();

      // Create a new name for the range
      // check if already exists and matches..
      if (!x.hasRange(range)) {
         x.setRange(range, thisxmin, thisxmax);
      }
      assert(x.getMin(range)==thisxmin);
      assert(x.getMax(range)==thisxmax);
      ranges->addOwned(*new RooStringVar(range, range, range));
   }
   std::pair<RangeMap::iterator, bool> r = _ranges.insert(
       std::make_pair(std::string(rangeName), ranges));
   assert(r.second);
   return *r.first->second;
}

initGenerator()

void RooEffResModel::initGenerator ( Int_t  code )

virtual

Definition at line 360 of file RooEffResModel.cpp.

{
    model().initGenerator(code);
}

model()

RooResolutionModel & RooEffResModel::model ( ) const

virtual

Definition at line 229 of file RooEffResModel.cpp.

{ return dynamic_cast<RooResolutionModel&>(*_model.absArg()); }

modelGenContext()

RooAbsGenContext * RooEffResModel::modelGenContext ( const RooAbsAnaConvPdf &  convPdf, const RooArgSet &  vars, const RooDataSet *  prototype = 0, const RooArgSet *  auxProto = 0, Bool_t  verbose = kFALSE  ) const

virtual

Definition at line 343 of file RooEffResModel.cpp.

{
    return new RooEffConvGenContext(convPdf, vars, prototype, auxProto, verbose);
}

observables()

RooArgSet * RooEffResModel::observables ( ) const

virtual

Return pointer to pdf in product.

Implements RooAbsEffResModel.

Definition at line 220 of file RooEffResModel.cpp.

{ return new RooArgSet(_observables); }

Friends And Related Function Documentation

RooMultiEffResModel

friend class RooMultiEffResModel

friend

Definition at line 73 of file RooEffResModel.h.

Member Data Documentation

_cacheMgr

RooObjCacheManager RooEffResModel::_cacheMgr

mutableprivate

Definition at line 111 of file RooEffResModel.h.

_eff

RooRealProxy RooEffResModel::_eff

private

Definition at line 106 of file RooEffResModel.h.

_model

RooRealProxy RooEffResModel::_model

private

Definition at line 105 of file RooEffResModel.h.

_observables

RooSetProxy RooEffResModel::_observables

private

Definition at line 104 of file RooEffResModel.h.

_ranges

RangeMap RooEffResModel::_ranges

mutableprivate

Definition at line 109 of file RooEffResModel.h.

---

The documentation for this class was generated from the following files:

• Mint/RooEffResModel.h
• src/Mojito/TimeDependent/RooEffResModel.cpp