104 _pdfList(
"!pdfs",
"List of PDFs",this),
139 _pdfList(
"!pdfs",
"List of PDFs",this),
160 <<
") multiple components with extended terms detected,"
161 <<
" product will not be extendible." << endl ;
194 _pdfList(
"!pdfs",
"List of PDFs",this),
202 Int_t numExtended(0) ;
207 << arg->
GetName() <<
" is not a PDF, ignored" << endl ;
224 <<
") WARNING: multiple components with extended terms detected,"
225 <<
" product will not be extendible." << endl ;
276 _pdfList(
"!pdfs",
"List of PDFs",this),
306 _pdfList(
"!pdfList",
"List of PDFs",this),
332 _pdfList(
"!pdfs",
"List of PDFs",this),
349 _cacheMgr(other._cacheMgr,this),
350 _genCode(other._genCode),
351 _cutOff(other._cutOff),
352 _pdfList(
"!pdfs",this,other._pdfList),
353 _extendedIndex(other._extendedIndex),
354 _useDefaultGen(other._useDefaultGen),
355 _refRangeName(other._refRangeName),
356 _selfNorm(other._selfNorm),
357 _defNormSet(other._defNormSet)
378 Int_t numExtended(0) ;
400 if (0 == strcmp(carg->
GetName(),
"Conditional")) {
412 tmp->
setName(0 == argType ?
"nset" :
"cset") ;
422 }
else if (0 != strlen(carg->
GetName())) {
430 <<
") WARNING: multiple components with extended terms detected,"
431 <<
" product will not be extendible." << endl ;
481 cxcoutD(
Eval) <<
"RooProdPdf::calculate(" <<
GetName() <<
") rearranged product calculation"
495 const auto normSet = cache.
_normList[i].get();
497 const Double_t piVal = partInt.getVal(normSet->getSize() > 0 ? normSet :
nullptr);
516 if (cache->_isRearranged) {
518 cxcoutD(
Eval) <<
"RooProdPdf::calculate(" <<
GetName() <<
") rearranged product calculation"
519 <<
" calculate: num = " << cache->_rearrangedNum->GetName() <<
" = " << cache->_rearrangedNum->getVal() << endl ;
520 cxcoutD(
Eval) <<
"calculate: den = " << cache->_rearrangedDen->GetName() <<
" = " << cache->_rearrangedDen->getVal() << endl ;
524 auto numerator = cache->_rearrangedNum->getValBatch(begin, size);
525 auto denominator = cache->_rearrangedDen->getValBatch(begin, size);
527 for (std::size_t i=0; i < outputs.size(); ++i) {
528 outputs[i] = numerator[i] / denominator[i];
535 assert(cache->_normList.size() == cache->_partList.size());
536 for (std::size_t i = 0; i < cache->_partList.size(); ++i) {
537 const auto& partInt =
static_cast<const RooAbsReal&
>(cache->_partList[i]);
538 const auto normSet = cache->_normList[i].get();
540 const auto partialInts = partInt.
getValBatch(begin, size, normSet->getSize() > 0 ? normSet :
nullptr);
541 for (std::size_t j=0; j < outputs.size(); ++j) {
542 outputs[j] *= partialInts[j];
576 (pdfNSetOrig = (
RooArgSet*) nIter.next(),
581 if (0 == strcmp(
"nset", pdfNSetOrig->
GetName())) {
584 }
else if (0 == strcmp(
"cset", pdfNSetOrig->
GetName())) {
587 pdfCSet = pdfNSetOrig;
601 pdfAllDeps.
add(*tmp);
611 pdfNormDeps.
add(*tmp2);
615 pdfNormDeps.
add(pdfAllDeps);
638 (termNormDeps = (
RooArgSet*) ldIter.next(),
639 termAllDeps = (
RooArgSet*) laIter.next(),
657 termIntDeps =
new RooArgSet(
"termIntDeps");
660 if (!termIntNoNormDeps) {
661 termIntNoNormDeps =
new RooArgSet(
"termIntNoNormDeps");
663 termIntNoNormDeps->
add(pdfIntNoNormDeps,
kFALSE);
674 termNormDeps =
new RooArgSet(
"termNormDeps");
675 termAllDeps =
new RooArgSet(
"termAllDeps");
676 termIntDeps =
new RooArgSet(
"termIntDeps");
677 termIntNoNormDeps =
new RooArgSet(
"termIntNoNormDeps");
683 termIntNoNormDeps->
add(pdfIntNoNormDeps,
kFALSE);
686 normList.
Add(termNormDeps);
687 depAllList.
Add(termAllDeps);
688 intList.
Add(termIntDeps);
689 depIntNoNormList.
Add(termIntNoNormDeps);
696 if (pdfCSet != pdfNSetOrig) {
702 RooArgSet *normDeps, *allDeps, *intNoNormDeps;
709 intNoNormDeps = (
RooArgSet*) innIter.next(),
726 depIntNoNormList.
Delete();
748 Int_t sterileIdx(-1);
782 map<string, RooArgSet> ratioTerms;
784 if (1 ==
group->GetSize()) {
792 RooArgSet termNSet(*norm), termImpSet(*imps);
803 RooFIter niter = termNSet.fwdIterator();
819 if (!rangeIdentical || 1) {
822 ostringstream str; termImpSet.
printValue(str);
824 ratioTerms[str.str()].add(*ratio);
839 RooArgSet termNSet(*norm), termImpSet(*imps);
845 RooFIter niter = termNSet.fwdIterator();
860 if (!rangeIdentical || 1) {
863 ostringstream str; termImpSet.
printValue(str);
864 ratioTerms[str.str()].add(*ratio);
876 if (1 ==
group->GetSize()) {
882 RooArgSet termNSet(*norm), termImpSet(*imps);
885 ostringstream str; termNSet.printValue(str);
886 if (ratioTerms[str.str()].getSize() > 0) {
888 term->
add(ratioTerms[str.str()]);
898 RooArgSet termNSet(*norm), termImpSet(*imps);
901 ostringstream str; termNSet.printValue(str);
902 if (ratioTerms[str.str()].getSize() > 0) {
904 term->
add(ratioTerms[str.str()]);
915 if (1 ==
group->GetSize()) {
925 RooArgSet termNSet, termISet, termXSet, termImpSet;
928 termISet.
add(*integ);
932 termXSet.
add(*xdeps);
933 termImpSet.
add(*imps);
942 vector<RooAbsReal*> func =
processProductTerm(nset, iset, isetRangeName, term, termNSet, termISet, isOwned);
957 RooArgSet compTermSet, compTermNorm, compTermNum, compTermDen;
968 RooArgSet termNSet, termISet, termXSet, termImpSet;
969 termISet.
add(*integ);
971 termXSet.
add(*xdeps);
972 termImpSet.
add(*imps);
984 vector<RooAbsReal*> func =
processProductTerm(nset, iset, isetRangeName, term, termNSet, termISet, isOwned,
kTRUE);
987 compTermSet.
add(*func[0]);
991 compTermNum.
add(*func[1]);
992 compTermDen.
add(*func[2]);
1011 const std::string prodname =
makeRGPPName(
"SPECPROD", compTermSet, outerIntDeps,
RooArgSet(), isetRangeName);
1015 const std::string intname =
makeRGPPName(
"SPECINT", compTermSet, outerIntDeps,
RooArgSet(), isetRangeName);
1024 RooProduct* prodtmp_num =
new RooProduct(prodname_num.c_str(), prodname_num.c_str(), compTermNum);
1030 RooProduct* prodtmp_den =
new RooProduct(prodname_den.c_str(), prodname_den.c_str(), compTermDen);
1035 string name =
Form(
"SPEC_RATIO(%s,%s)", prodname_num.c_str(), prodname_den.c_str());
1116 list<string> rangeComps ;
1123 rangeComps.push_back(token) ;
1130 map<string,RooArgSet> denListList ;
1132 string specIntRange ;
1175 nomList.
add(*funcCust) ;
1177 nomList.
add(*func) ;
1198 origNumTerm.
add(*arg) ;
1206 nomList.
add(origNumTerm) ;
1213 for (list<string>::iterator iter = rangeComps.begin() ; iter != rangeComps.end() ; ++iter) {
1259 partCust = cust.
build() ;
1264 partCust = cust.
build() ;
1274 string name =
Form(
"%s_divided_by_ratio",specializedPartCust->
GetName()) ;
1277 denListList[*iter].add(*specIntFinal) ;
1298 tmp.
add(*specRatio) ;
1310 throw string(
"this should not happen") ;
1316 denListList[*iter].add(*specIntFinal) ;
1338 for (map<string,RooArgSet>::iterator iter = denListList.begin() ; iter != denListList.end() ; ++iter) {
1343 products.
add(*prod_comp) ;
1349 if (specIntDeps.
getSize()>0) {
1358 numerator = numtmp ;
1458 outerIntDeps.
add(*tmp) ;
1485 if (termNormDeps->
contains(*outerIntDep) ||
1486 termIntDeps->
contains(*outerIntDep) ||
1487 termImpDeps->
contains(*outerIntDep)) {
1500 tIter2 =
group->fwdIterator() ;
1502 newGroup->
Add(term2) ;
1512 groupedTerms.
Add(newGroup) ;
1538 vector<RooAbsReal*> ret(3) ; ret[0] = 0 ; ret[1] = 0 ; ret[2] = 0 ;
1554 if (nset && termNSet.
getSize()==0) {
1562 if (iset && termISet.
getSize()>0) {
1593 const std::string
name =
makeRGPPName(
"GENPROJ_",*term,termISet,termNSet,isetRangeName) ;
1620 const std::string
name =
makeRGPPName(
"GENPROJ_",*term,termISet,termNSet,isetRangeName) ;
1653 name.Append(
"_NORM[") ;
1700 coutE(
Eval) <<
"RooProdPdf::processProductTerm(" <<
GetName() <<
") unidentified term!!!" << endl ;
1711 const RooArgSet& nset,
const char* isetRangeName)
const
1715 std::ostringstream os(pfx);
1723 if (!
first) os <<
"_X_";
1760 const RooArgSet* normSet,
const char* rangeName)
const
1765 analVars.
add(allVars) ;
1846 coutF(
Generation) <<
"Requesting expected number of events from a RooProdPdf that does not contain an extended p.d.f" << endl ;
1847 throw std::logic_error(std::string(
"RooProdPdf ") +
GetName() +
" could not be extended.");
1886 std::vector<Int_t> code;
1892 code.push_back(pdfCode);
1894 generateVars.
add(pdfDirect) ;
1899 if (generateVars.
getSize()>0) {
1901 return masterCode+1 ;
1922 if (codeList[i]!=0) {
1945 if (codeList[i]!=0) {
1977 os <<
indent <<
"RooProdPdf begin partial integral cache" << endl ;
1980 RooFIter iter = _partList.fwdIterator() ;
1983 indent2 +=
Form(
"[%d] ",curElem) ;
1988 if (curElem==maxElem) {
1989 os <<
indent <<
"RooProdPdf end partial integral cache" << endl ;
2013 if (thePdf)
return kFALSE ;
2030 if (idx<0)
return 0 ;
2054 constraints.
add(*pdf) ;
2070 finalConstraints->
add(*pdf) ;
2073 <<
" as constraint term as it does not share any parameters with the other pdfs in product. "
2074 <<
"To force inclusion in likelihood, add an explicit Constrain() argument for the target parameter" << endl ;
2085 return finalConstraints ;
2102 if (arg->dependsOn(observables)) {
2103 RooArgSet* tmp = arg->getParameters(observables) ;
2104 connectedPars->
add(*tmp) ;
2108 return connectedPars ;
2118 if (!stripDisconnected) return ;
2119 if (!nset || nset->
getSize()==0) return ;
2127 for (
auto param : *params) {
2129 for (
auto term : plist) {
2130 if (term->dependsOnValue(*param)) {
2135 tostrip.
add(*param) ;
2239 if (parg->canNodeBeCached()==
Always) {
2240 trackNodes.
add(*parg) ;
2247 if (
string(
"nset")==pdf_nset->
GetName() && pdf_nset->
getSize()>0) {
2249 parg->setStringAttribute(
"CATNormSet",
n.content()) ;
2251 if (
string(
"cset")==pdf_nset->
GetName()) {
2253 parg->setStringAttribute(
"CATCondSet",
c.content()) ;
2256 coutW(
Optimization) <<
"RooProdPdf::setCacheAndTrackHints(" <<
GetName() <<
") WARNING product pdf does not specify a normalization set for component " << parg->GetName() << endl ;
2272 if (i>0) os <<
" * " ;
2276 if (
string(
"nset")==ncset->
GetName()) {
char * R__STRTOK_R(char *str, const char *delim, char **saveptr)
static void indent(ostringstream &buf, int indent_level)
char * Form(const char *fmt,...)
RooSpan< double > makeWritableBatchUnInit(std::size_t begin, std::size_t batchSize, const RooArgSet *const normSet=nullptr, Tag_t ownerTag=kUnspecified)
Make a batch and return a span pointing to the pdf-local memory.
RooSpan< double > makeWritableBatchInit(std::size_t begin, std::size_t batchSize, double value, const RooArgSet *const normSet=nullptr, Tag_t ownerTag=kUnspecified)
Make a batch and return a span pointing to the pdf-local memory.
const std::vector< Int_t > & retrieve(Int_t masterCode) const
Retrieve the array of integer codes associated with the given master code.
Int_t store(const std::vector< Int_t > &codeList, RooArgSet *set1=0, RooArgSet *set2=0, RooArgSet *set3=0, RooArgSet *set4=0)
Store given arrays of integer codes, and up to four RooArgSets in the registry (each setX pointer may...
RooAbsArg is the common abstract base class for objects that represent a value (of arbitrary type) an...
RooExpensiveObjectCache & expensiveObjectCache() const
RooArgSet * getObservables(const RooArgSet &set, Bool_t valueOnly=kTRUE) const
Return the observables of this pdf given a set of observables.
void printCompactTree(const char *indent="", const char *fileName=0, const char *namePat=0, RooAbsArg *client=0)
Print tree structure of expression tree on stdout, or to file if filename is specified.
Bool_t dependsOn(const RooAbsCollection &serverList, const RooAbsArg *ignoreArg=0, Bool_t valueOnly=kFALSE) const
Test whether we depend on (ie, are served by) any object in the specified collection.
void setStringAttribute(const Text_t *key, const Text_t *value)
Associate string 'value' to this object under key 'key'.
virtual TObject * Clone(const char *newname=0) const
Make a clone of an object using the Streamer facility.
virtual void setExpensiveObjectCache(RooExpensiveObjectCache &cache)
const Text_t * getStringAttribute(const Text_t *key) const
Get string attribute mapped under key 'key'.
RooArgSet * getParameters(const RooAbsData *data, Bool_t stripDisconnected=kTRUE) const
Create a list of leaf nodes in the arg tree starting with ourself as top node that don't match any of...
void setAttribute(const Text_t *name, Bool_t value=kTRUE)
Set (default) or clear a named boolean attribute of this object.
Bool_t addOwnedComponents(const RooArgSet &comps)
Take ownership of the contents of 'comps'.
Bool_t getAttribute(const Text_t *name) const
Check if a named attribute is set. By default, all attributes are unset.
Bool_t dependsOnValue(const RooAbsCollection &serverList, const RooAbsArg *ignoreArg=0) const
Check whether this object depends on values from an element in the serverList.
RooArgSet * getComponents() const
RooAbsCacheElement is the abstract base class for objects to be stored in RooAbsCache cache manager o...
RooAbsCollection is an abstract container object that can hold multiple RooAbsArg objects.
virtual void removeAll()
Remove all arguments from our set, deleting them if we own them.
RooAbsCollection * selectCommon(const RooAbsCollection &refColl) const
Create a subset of the current collection, consisting only of those elements that are contained as we...
Bool_t contains(const RooAbsArg &var) const
Check if collection contains an argument with the same name as var.
virtual Bool_t addOwned(RooAbsArg &var, Bool_t silent=kFALSE)
Add the specified argument to list.
Int_t index(const RooAbsArg *arg) const
Returns index of given arg, or -1 if arg is not in the collection.
RooFIter fwdIterator() const
One-time forward iterator.
virtual Bool_t add(const RooAbsArg &var, Bool_t silent=kFALSE)
Add the specified argument to list.
Storage_t::size_type size() const
RooAbsArg * first() const
Bool_t overlaps(const RooAbsCollection &otherColl) const
Check if this and other collection have common entries.
const char * GetName() const
Returns name of object.
void setName(const char *name)
virtual Bool_t remove(const RooAbsArg &var, Bool_t silent=kFALSE, Bool_t matchByNameOnly=kFALSE)
Remove the specified argument from our list.
RooAbsArg * find(const char *name) const
Find object with given name in list.
virtual void printValue(std::ostream &os) const
Print value of collection, i.e.
RooAbsGenContext is the abstract base class for generator contexts of RooAbsPdf objects.
virtual void generateEvent(Int_t code)
Interface for generation of an event using the algorithm corresponding to the specified code.
virtual RooAbsGenContext * genContext(const RooArgSet &vars, const RooDataSet *prototype=0, const RooArgSet *auxProto=0, Bool_t verbose=kFALSE) const
Interface function to create a generator context from a p.d.f.
TString _normRange
MC generator configuration specific for this object.
virtual Int_t getGenerator(const RooArgSet &directVars, RooArgSet &generateVars, Bool_t staticInitOK=kTRUE) const
Load generatedVars with the subset of directVars that we can generate events for, and return a code t...
Bool_t canBeExtended() const
friend class RooRealIntegral
const char * normRange() const
virtual Bool_t isDirectGenSafe(const RooAbsArg &arg) const
Check if given observable can be safely generated using the pdfs internal generator mechanism (if tha...
virtual void initGenerator(Int_t code)
Interface for one-time initialization to setup the generator for the specified code.
RooArgSet * _normSet
Normalization integral (owned by _normMgr)
RooAbsRealLValue is the common abstract base class for objects that represent a real value that may a...
virtual Double_t getMax(const char *name=0) const
Get maximum of currently defined range.
virtual Double_t getMin(const char *name=0) const
Get miniminum of currently defined range.
RooAbsReal is the common abstract base class for objects that represent a real value and implements f...
virtual std::list< Double_t > * plotSamplingHint(RooAbsRealLValue &, Double_t, Double_t) const
TString integralNameSuffix(const RooArgSet &iset, const RooArgSet *nset=0, const char *rangeName=0, Bool_t omitEmpty=kFALSE) const
Construct string with unique suffix name to give to integral object that encodes integrated observabl...
virtual RooSpan< const double > getValBatch(std::size_t begin, std::size_t maxSize, const RooArgSet *normSet=nullptr) const
Return value of object for all data events in the batch.
virtual std::list< Double_t > * binBoundaries(RooAbsRealLValue &, Double_t, Double_t) const
Double_t getVal(const RooArgSet *normalisationSet=nullptr) const
Evaluate object.
RooAbsReal * createIntegral(const RooArgSet &iset, const RooCmdArg &arg1, const RooCmdArg &arg2=RooCmdArg::none(), const RooCmdArg &arg3=RooCmdArg::none(), const RooCmdArg &arg4=RooCmdArg::none(), const RooCmdArg &arg5=RooCmdArg::none(), const RooCmdArg &arg6=RooCmdArg::none(), const RooCmdArg &arg7=RooCmdArg::none(), const RooCmdArg &arg8=RooCmdArg::none()) const
Create an object that represents the integral of the function over one or more observables listed in ...
BatchHelpers::BatchData _batchData
virtual Bool_t isBinnedDistribution(const RooArgSet &) const
RooAddition calculates the sum of a set of RooAbsReal terms, or when constructed with two sets,...
const RooArgList & list1() const
RooArgList is a container object that can hold multiple RooAbsArg objects.
RooAbsArg * at(Int_t idx) const
Return object at given index, or nullptr if index is out of range.
RooArgSet is a container object that can hold multiple RooAbsArg objects.
RooArgSet * snapshot(bool deepCopy=true) const
Use RooAbsCollection::snapshot(), but return as RooArgSet.
virtual Bool_t add(const RooAbsCollection &col, Bool_t silent=kFALSE)
Add a collection of arguments to this collection by calling add() for each element in the source coll...
T * getObj(const RooArgSet *nset, Int_t *sterileIndex=0, const TNamed *isetRangeName=0)
const RooNameSet * nameSet2ByIndex(Int_t index) const
Retrieve RooNameSet associated with slot at given index.
T * getObjByIndex(Int_t index) const
Retrieve payload object by slot index.
void reset()
Clear the cache.
const RooNameSet * nameSet1ByIndex(Int_t index) const
Retrieve RooNameSet associated with slot at given index.
Int_t setObj(const RooArgSet *nset, T *obj, const TNamed *isetRangeName=0)
RooCmdArg is a named container for two doubles, two integers two object points and three string point...
const RooArgSet * getSet(Int_t idx) const
Return RooArgSet stored in slot idx.
Int_t getInt(Int_t idx) const
virtual TObject * clone(const char *newname) const
RooCustomizer is a factory class to produce clones of a prototype composite PDF object with the same ...
void replaceArg(const RooAbsArg &orig, const RooAbsArg &subst)
Replace any occurence of arg 'orig' with arg 'subst'.
RooAbsArg * build(const char *masterCatState, Bool_t verbose=kFALSE)
Build a clone of the prototype executing all registered 'replace' rules and 'split' rules for the mas...
RooDataSet is a container class to hold unbinned data.
A one-time forward iterator working on RooLinkedList or RooAbsCollection.
RooAbsArg * next()
Return next element or nullptr if at end.
RooGenProdProj is an auxiliary class for RooProdPdf that calculates a general normalised projection o...
RooLinkedList is an collection class for internal use, storing a collection of RooAbsArg pointers in ...
RooFIter fwdIterator() const
Create a one-time-use forward iterator for this list.
TObject * At(Int_t index) const
Return object stored in sequential position given by index.
void Delete(Option_t *o=0)
Remove all elements in collection and delete all elements NB: Collection does not own elements,...
virtual void Add(TObject *arg)
virtual Bool_t Remove(TObject *arg)
Remove object from collection.
Int_t IndexOf(const char *name) const
Return position of given object in list.
virtual Bool_t add(const RooAbsArg &var, Bool_t silent=kFALSE)
Reimplementation of standard RooArgList::add()
virtual Bool_t remove(const RooAbsArg &var, Bool_t silent=kFALSE, Bool_t matchByNameOnly=kFALSE)
Reimplementation of standard RooArgList::remove()
static const char * str(const TNamed *ptr)
Return C++ string corresponding to given TNamed pointer.
static const TNamed * ptr(const char *stringPtr)
Return a unique TNamed pointer for given C++ string.
RooNameSet is a utility class that stores the names the objects in a RooArget.
RooArgSet * select(const RooArgSet &list) const
Construct a RooArgSet of objects in input 'list' whose names match to those in the internal name list...
virtual void printCompactTreeHook(std::ostream &, const char *, Int_t, Int_t)
Hook function to print cache contents in tree printing of RooProdPdf.
std::unique_ptr< RooAbsReal > _rearrangedNum
std::vector< std::unique_ptr< RooArgSet > > _normList
std::unique_ptr< RooAbsReal > _rearrangedDen
virtual RooArgList containedArgs(Action)
Return RooAbsArg components contained in the cache.
RooProdPdf is an efficient implementation of a product of PDFs of the form.
void rearrangeProduct(CacheElem &) const
virtual Double_t expectedEvents(const RooArgSet *nset) const
Return the expected number of events associated with the extendable input PDF in the product.
void factorizeProduct(const RooArgSet &normSet, const RooArgSet &intSet, RooLinkedList &termList, RooLinkedList &normList, RooLinkedList &impDepList, RooLinkedList &crossDepList, RooLinkedList &intList) const
Factorize product in irreducible terms for given choice of integration/normalization.
void groupProductTerms(RooLinkedList &groupedTerms, RooArgSet &outerIntDeps, const RooLinkedList &terms, const RooLinkedList &norms, const RooLinkedList &imps, const RooLinkedList &ints, const RooLinkedList &cross) const
Group product into terms that can be calculated independently.
Int_t getPartIntList(const RooArgSet *nset, const RooArgSet *iset, const char *isetRangeName=0) const
Return list of (partial) integrals of product terms for integration of p.d.f over observables iset wh...
RooProdPdf()
Default constructor.
virtual Bool_t checkObservables(const RooArgSet *nset) const
Obsolete.
Double_t _cutOff
Registry of composite direct generator codes.
void printMetaArgs(std::ostream &os) const
Customized printing of arguments of a RooProdPdf to more intuitively reflect the contents of the prod...
Double_t evaluate() const
Calculate current value of object.
Bool_t isBinnedDistribution(const RooArgSet &obs) const
If all components that depend on obs are binned that so is the product.
virtual Bool_t isDirectGenSafe(const RooAbsArg &arg) const
Forward determination of safety of internal generator code to component p.d.f that would generate the...
Int_t getAnalyticalIntegralWN(RooArgSet &allVars, RooArgSet &numVars, const RooArgSet *normSet, const char *rangeName=0) const
Determine which part (if any) of given integral can be performed analytically.
RooAbsReal * specializeIntegral(RooAbsReal &orig, const char *targetRangeName) const
virtual ExtendMode extendMode() const
If this product contains exactly one extendable p.d.f return the extension abilities of that p....
RooObjCacheManager _cacheMgr
virtual void initGenerator(Int_t code)
Forward one-time initialization call to component generation initialization methods.
std::string makeRGPPName(const char *pfx, const RooArgSet &term, const RooArgSet &iset, const RooArgSet &nset, const char *isetRangeName) const
Make an appropriate automatic name for a RooGenProdProj object in getPartIntList()
Double_t analyticalIntegralWN(Int_t code, const RooArgSet *normSet, const char *rangeName=0) const
Return analytical integral defined by given scenario code.
virtual ~RooProdPdf()
Destructor.
virtual void getParametersHook(const RooArgSet *, RooArgSet *, Bool_t stripDisconnected) const
RooAbsReal * makeCondPdfRatioCorr(RooAbsReal &term, const RooArgSet &termNset, const RooArgSet &termImpSet, const char *normRange, const char *refRange) const
For single normalization ranges.
RooArgSet * findPdfNSet(RooAbsPdf &pdf) const
Look up user specified normalization set for given input PDF component.
virtual RooSpan< double > evaluateBatch(std::size_t begin, std::size_t size) const
Evaluate function for a batch of input data points.
virtual void setCacheAndTrackHints(RooArgSet &)
Label OK'ed components of a RooProdPdf with cache-and-track, and label all RooProdPdf descendants wit...
virtual RooArgSet * getConstraints(const RooArgSet &observables, RooArgSet &constrainedParams, Bool_t stripDisconnected) const
Return all parameter constraint p.d.f.s on parameters listed in constrainedParams.
RooLinkedList _pdfNSetList
friend class RooProdGenContext
virtual Int_t getGenerator(const RooArgSet &directVars, RooArgSet &generateVars, Bool_t staticInitOK=kTRUE) const
Query internal generation capabilities of component p.d.f.s and aggregate capabilities into master co...
virtual void selectNormalizationRange(const char *rangeName=0, Bool_t force=kFALSE)
Interface function used by test statistics to freeze choice of range for interpretation of conditiona...
RooArgSet * getConnectedParameters(const RooArgSet &observables) const
Return all parameter constraint p.d.f.s on parameters listed in constrainedParams.
Double_t calculate(const RooProdPdf::CacheElem &cache, Bool_t verbose=kFALSE) const
Calculate running product of pdfs terms, using the supplied normalization set in 'normSetList' for ea...
virtual std::list< Double_t > * plotSamplingHint(RooAbsRealLValue &obs, Double_t xlo, Double_t xhi) const
Forward the plot sampling hint from the p.d.f. that defines the observable obs.
RooAbsReal * specializeRatio(RooFormulaVar &input, const char *targetRangeName) const
Bool_t redirectServersHook(const RooAbsCollection &, Bool_t, Bool_t, Bool_t)
Implement support for node removal.
void fixRefRange(const char *rangeName)
virtual Bool_t forceAnalyticalInt(const RooAbsArg &dep) const
Force RooRealIntegral to offer all observables for internal integration.
virtual TObject * clone(const char *newname) const
virtual RooAbsGenContext * genContext(const RooArgSet &vars, const RooDataSet *prototype=0, const RooArgSet *auxProto=0, Bool_t verbose=kFALSE) const
Return generator context optimized for generating events from product p.d.f.s.
void initializeFromCmdArgList(const RooArgSet &fullPdfSet, const RooLinkedList &l)
Initialize RooProdPdf configuration from given list of RooCmdArg configuration arguments and set of '...
virtual void generateEvent(Int_t code)
Generate a single event with configuration specified by 'code' Defer internal generation to component...
std::vector< RooAbsReal * > processProductTerm(const RooArgSet *nset, const RooArgSet *iset, const char *isetRangeName, const RooArgSet *term, const RooArgSet &termNSet, const RooArgSet &termISet, Bool_t &isOwned, Bool_t forceWrap=kFALSE) const
Calculate integrals of factorized product terms over observables iset while normalized to observables...
virtual std::list< Double_t > * binBoundaries(RooAbsRealLValue &, Double_t, Double_t) const
Forward the plot sampling hint from the p.d.f. that defines the observable obs.
A RooProduct represents the product of a given set of RooAbsReal objects.
RooRealIntegral performs hybrid numerical/analytical integrals of RooAbsReal objects.
RooArgSet intVars() const
const RooAbsReal & integrand() const
RooRealVar represents a variable that can be changed from the outside.
A simple container to hold a batch of data values.
The TNamed class is the base class for all named ROOT classes.
virtual const char * GetName() const
Returns name of object.
Mother of all ROOT objects.
virtual Bool_t InheritsFrom(const char *classname) const
Returns kTRUE if object inherits from class "classname".
const char * Data() const
Bool_t Contains(const char *pat, ECaseCompare cmp=kExact) const
RooConstVar & RooConst(Double_t val)