105 _pdfList(
"!pdfs",
"List of PDFs",this),
140 _pdfList(
"!pdfs",
"List of PDFs",this),
161 <<
") multiple components with extended terms detected,"
162 <<
" product will not be extendible." << endl ;
195 _pdfList(
"!pdfs",
"List of PDFs",this),
203 Int_t numExtended(0) ;
208 << arg->
GetName() <<
" is not a PDF, ignored" << endl ;
225 <<
") WARNING: multiple components with extended terms detected,"
226 <<
" 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 ;
492 cxcoutD(
Eval) <<
"RooProdPdf::calculate(" <<
GetName() <<
") rearranged product calculation"
506 const auto normSet = cache.
_normList[i].get();
508 const Double_t piVal = partInt.getVal(normSet->getSize() > 0 ? normSet :
nullptr);
527 if (cache->_isRearranged) {
529 cxcoutD(
Eval) <<
"RooProdPdf::calculate(" <<
GetName() <<
") rearranged product calculation"
530 <<
" calculate: num = " << cache->_rearrangedNum->GetName() <<
" = " << cache->_rearrangedNum->getVal() << endl ;
531 cxcoutD(
Eval) <<
"calculate: den = " << cache->_rearrangedDen->GetName() <<
" = " << cache->_rearrangedDen->getVal() << endl ;
535 auto numerator = cache->_rearrangedNum->getValBatch(begin, size);
536 auto denominator = cache->_rearrangedDen->getValBatch(begin, size);
538 for (std::size_t i=0; i < outputs.size(); ++i) {
539 outputs[i] = numerator[i] / denominator[i];
546 assert(cache->_normList.size() == cache->_partList.size());
547 for (std::size_t i = 0; i < cache->_partList.size(); ++i) {
548 const auto& partInt =
static_cast<const RooAbsReal&
>(cache->_partList[i]);
549 const auto normSet = cache->_normList[i].get();
551 const auto partialInts = partInt.
getValBatch(begin, size, normSet->getSize() > 0 ? normSet :
nullptr);
552 for (std::size_t j=0; j < outputs.size(); ++j) {
553 outputs[j] *= partialInts[j];
587 (pdfNSetOrig = (
RooArgSet*) nIter.next(),
592 if (0 == strcmp(
"nset", pdfNSetOrig->
GetName())) {
595 }
else if (0 == strcmp(
"cset", pdfNSetOrig->
GetName())) {
598 pdfCSet = pdfNSetOrig;
612 pdfAllDeps.
add(*tmp);
622 pdfNormDeps.
add(*tmp2);
626 pdfNormDeps.
add(pdfAllDeps);
649 (termNormDeps = (
RooArgSet*) ldIter.next(),
650 termAllDeps = (
RooArgSet*) laIter.next(),
668 termIntDeps =
new RooArgSet(
"termIntDeps");
671 if (!termIntNoNormDeps) {
672 termIntNoNormDeps =
new RooArgSet(
"termIntNoNormDeps");
674 termIntNoNormDeps->
add(pdfIntNoNormDeps,
kFALSE);
685 termNormDeps =
new RooArgSet(
"termNormDeps");
686 termAllDeps =
new RooArgSet(
"termAllDeps");
687 termIntDeps =
new RooArgSet(
"termIntDeps");
688 termIntNoNormDeps =
new RooArgSet(
"termIntNoNormDeps");
694 termIntNoNormDeps->
add(pdfIntNoNormDeps,
kFALSE);
697 normList.
Add(termNormDeps);
698 depAllList.
Add(termAllDeps);
699 intList.
Add(termIntDeps);
700 depIntNoNormList.
Add(termIntNoNormDeps);
707 if (pdfCSet != pdfNSetOrig) {
713 RooArgSet *normDeps, *allDeps, *intNoNormDeps;
720 intNoNormDeps = (
RooArgSet*) innIter.next(),
737 depIntNoNormList.
Delete();
759 Int_t sterileIdx(-1);
793 map<string, RooArgSet> ratioTerms;
795 if (1 ==
group->GetSize()) {
803 RooArgSet termNSet(*norm), termImpSet(*imps);
814 RooFIter niter = termNSet.fwdIterator();
830 if (!rangeIdentical || 1) {
833 ostringstream str; termImpSet.
printValue(str);
835 ratioTerms[str.str()].add(*ratio);
850 RooArgSet termNSet(*norm), termImpSet(*imps);
856 RooFIter niter = termNSet.fwdIterator();
871 if (!rangeIdentical || 1) {
874 ostringstream str; termImpSet.
printValue(str);
875 ratioTerms[str.str()].add(*ratio);
887 if (1 ==
group->GetSize()) {
893 RooArgSet termNSet(*norm), termImpSet(*imps);
896 ostringstream str; termNSet.printValue(str);
897 if (ratioTerms[str.str()].getSize() > 0) {
899 term->
add(ratioTerms[str.str()]);
909 RooArgSet termNSet(*norm), termImpSet(*imps);
912 ostringstream str; termNSet.printValue(str);
913 if (ratioTerms[str.str()].getSize() > 0) {
915 term->
add(ratioTerms[str.str()]);
926 if (1 ==
group->GetSize()) {
936 RooArgSet termNSet, termISet, termXSet, termImpSet;
939 termISet.
add(*integ);
943 termXSet.
add(*xdeps);
944 termImpSet.
add(*imps);
953 vector<RooAbsReal*> func =
processProductTerm(nset, iset, isetRangeName, term, termNSet, termISet, isOwned);
968 RooArgSet compTermSet, compTermNorm, compTermNum, compTermDen;
979 RooArgSet termNSet, termISet, termXSet, termImpSet;
980 termISet.
add(*integ);
982 termXSet.
add(*xdeps);
983 termImpSet.
add(*imps);
995 vector<RooAbsReal*> func =
processProductTerm(nset, iset, isetRangeName, term, termNSet, termISet, isOwned,
kTRUE);
998 compTermSet.
add(*func[0]);
1002 compTermNum.
add(*func[1]);
1003 compTermDen.
add(*func[2]);
1022 const std::string prodname =
makeRGPPName(
"SPECPROD", compTermSet, outerIntDeps,
RooArgSet(), isetRangeName);
1026 const std::string intname =
makeRGPPName(
"SPECINT", compTermSet, outerIntDeps,
RooArgSet(), isetRangeName);
1035 RooProduct* prodtmp_num =
new RooProduct(prodname_num.c_str(), prodname_num.c_str(), compTermNum);
1041 RooProduct* prodtmp_den =
new RooProduct(prodname_den.c_str(), prodname_den.c_str(), compTermDen);
1046 string name =
Form(
"SPEC_RATIO(%s,%s)", prodname_num.c_str(), prodname_den.c_str());
1127 list<string> rangeComps ;
1134 rangeComps.push_back(token) ;
1141 map<string,RooArgSet> denListList ;
1143 string specIntRange ;
1186 nomList.
add(*funcCust) ;
1188 nomList.
add(*func) ;
1209 origNumTerm.
add(*arg) ;
1217 nomList.
add(origNumTerm) ;
1224 for (list<string>::iterator iter = rangeComps.begin() ; iter != rangeComps.end() ; ++iter) {
1270 partCust = cust.
build() ;
1275 partCust = cust.
build() ;
1285 string name =
Form(
"%s_divided_by_ratio",specializedPartCust->
GetName()) ;
1288 denListList[*iter].add(*specIntFinal) ;
1309 tmp.
add(*specRatio) ;
1321 throw string(
"this should not happen") ;
1327 denListList[*iter].add(*specIntFinal) ;
1349 for (map<string,RooArgSet>::iterator iter = denListList.begin() ; iter != denListList.end() ; ++iter) {
1354 products.
add(*prod_comp) ;
1360 if (specIntDeps.
getSize()>0) {
1369 numerator = numtmp ;
1469 outerIntDeps.
add(*tmp) ;
1496 if (termNormDeps->
contains(*outerIntDep) ||
1497 termIntDeps->
contains(*outerIntDep) ||
1498 termImpDeps->
contains(*outerIntDep)) {
1511 tIter2 =
group->fwdIterator() ;
1513 newGroup->
Add(term2) ;
1523 groupedTerms.
Add(newGroup) ;
1549 vector<RooAbsReal*> ret(3) ; ret[0] = 0 ; ret[1] = 0 ; ret[2] = 0 ;
1565 if (nset && termNSet.
getSize()==0) {
1573 if (iset && termISet.
getSize()>0) {
1604 const std::string
name =
makeRGPPName(
"GENPROJ_",*term,termISet,termNSet,isetRangeName) ;
1631 const std::string
name =
makeRGPPName(
"GENPROJ_",*term,termISet,termNSet,isetRangeName) ;
1664 name.Append(
"_NORM[") ;
1711 coutE(
Eval) <<
"RooProdPdf::processProductTerm(" <<
GetName() <<
") unidentified term!!!" << endl ;
1722 const RooArgSet& nset,
const char* isetRangeName)
const
1726 std::ostringstream os(pfx);
1734 if (!
first) os <<
"_X_";
1771 const RooArgSet* normSet,
const char* rangeName)
const
1776 analVars.
add(allVars) ;
1857 coutF(
Generation) <<
"Requesting expected number of events from a RooProdPdf that does not contain an extended p.d.f" << endl ;
1858 throw std::logic_error(std::string(
"RooProdPdf ") +
GetName() +
" could not be extended.");
1897 std::vector<Int_t> code;
1903 code.push_back(pdfCode);
1905 generateVars.
add(pdfDirect) ;
1910 if (generateVars.
getSize()>0) {
1912 return masterCode+1 ;
1933 if (codeList[i]!=0) {
1956 if (codeList[i]!=0) {
1988 os <<
indent <<
"RooProdPdf begin partial integral cache" << endl ;
1991 RooFIter iter = _partList.fwdIterator() ;
1994 indent2 +=
Form(
"[%d] ",curElem) ;
1999 if (curElem==maxElem) {
2000 os <<
indent <<
"RooProdPdf end partial integral cache" << endl ;
2024 if (thePdf)
return kFALSE ;
2041 if (idx<0)
return 0 ;
2065 constraints.
add(*pdf) ;
2081 finalConstraints->
add(*pdf) ;
2084 <<
" as constraint term as it does not share any parameters with the other pdfs in product. "
2085 <<
"To force inclusion in likelihood, add an explicit Constrain() argument for the target parameter" << endl ;
2096 return finalConstraints ;
2113 if (arg->dependsOn(observables)) {
2114 RooArgSet* tmp = arg->getParameters(observables) ;
2115 connectedPars->
add(*tmp) ;
2119 return connectedPars ;
2129 if (!stripDisconnected) return ;
2130 if (!nset || nset->
getSize()==0) return ;
2138 for (
auto param : *params) {
2140 for (
auto term : plist) {
2141 if (term->dependsOnValue(*param)) {
2146 tostrip.
add(*param) ;
2250 if (parg->canNodeBeCached()==
Always) {
2251 trackNodes.
add(*parg) ;
2258 if (
string(
"nset")==pdf_nset->
GetName() && pdf_nset->
getSize()>0) {
2260 parg->setStringAttribute(
"CATNormSet",
n.content()) ;
2262 if (
string(
"cset")==pdf_nset->
GetName()) {
2264 parg->setStringAttribute(
"CATCondSet",
c.content()) ;
2267 coutW(
Optimization) <<
"RooProdPdf::setCacheAndTrackHints(" <<
GetName() <<
") WARNING product pdf does not specify a normalization set for component " << parg->GetName() << endl ;
2283 if (i>0) os <<
" * " ;
2287 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 > makeWritableBatchInit(std::size_t begin, std::size_t batchSize, double value)
Make a batch and return a span pointing to the pdf-local memory.
RooSpan< double > makeWritableBatchUnInit(std::size_t begin, std::size_t batchSize)
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
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.
RooFIter fwdIterator() const R__SUGGEST_ALTERNATIVE("begin()
One-time forward iterator.
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
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.
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
virtual Double_t getValV(const RooArgSet *set=0) const
Return current value, normalized by integrating over the observables in nset.
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.
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.
RooArgSet * _curNormSet
Registry of composite direct generator codes.
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.
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
virtual Double_t getValV(const RooArgSet *set=0) const
Overload getVal() to intercept normalization set for use in evaluate()
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
R__EXTERN PyObject * gIter
RooConstVar & RooConst(Double_t val)