138 RooAbsArg(
name,title), _plotMin(0), _plotMax(0), _plotBins(100),
139 _value(0), _unit(unit), _forceNumInt(
kFALSE), _specIntegratorConfig(0), _selectComp(
kTRUE), _lastNSet(0)
153 RooAbsArg(
name,title), _plotMin(inMinVal), _plotMax(inMaxVal), _plotBins(100),
154 _value(0), _unit(unit), _forceNumInt(
kFALSE), _specIntegratorConfig(0), _selectComp(
kTRUE), _lastNSet(0)
167 _plotBins(other._plotBins), _value(other._value), _unit(other._unit), _label(other._label),
168 _forceNumInt(other._forceNumInt), _selectComp(other._selectComp), _lastNSet(0)
173 _specIntegratorConfig = 0 ;
219 return (
getVal()==value) ;
239 if (!assumeSameType) {
255 if(appendUnit && 0 != strlen(
getUnit())) {
312 maxSize = ret.size();
352 <<
"): validation failed: " << value << endl ;
372 const RooArgSet* ,
const char* rangeName)
const
403 if (code==0)
return getVal(normSet) ;
417 coutF(
Eval) <<
"RooAbsReal::analyticalIntegral(" <<
GetName() <<
") code " << code <<
" not implemented" << endl ;
478 os <<
indent <<
"--- RooAbsReal ---" << endl;
568 pc.defineString(
"rangeName",
"RangeWithName",0,
"",
kTRUE) ;
569 pc.defineObject(
"normSet",
"NormSet",0,0) ;
570 pc.defineObject(
"numIntConfig",
"NumIntConfig",0,0) ;
573 pc.process(arg1,arg2,arg3,arg4,arg5,arg6,arg7,arg8) ;
579 const char* rangeName =
pc.getString(
"rangeName",0,
kTRUE) ;
604 if (!rangeName || strchr(rangeName,
',')==0) {
614 for (
const std::string& token : tokens) {
616 components.
add(*compIntegral);
620 title.
Prepend(
"Integral of ") ;
648 title.
Prepend(
"Integral of ") ;
673 title.
Prepend(
"Integral of ") ;
682 if (integrand !=
this) {
690 if (integrand ==
this && iset.
getSize()>0) {
691 coutI(
Integration) <<
GetName() <<
" : multidimensional integration over observables with parameterized ranges in terms of other integrated observables detected, using recursive integration strategy to construct final integral" << endl ;
695 integrand = integral ;
703 coutE(
Integration) <<
GetName() <<
" : ERROR while defining recursive integral over observables with parameterized integration ranges, please check that integration rangs specify uniquely defined integral " << endl;
712 if (cacheParamsStr && strlen(cacheParamsStr)) {
720 if (cacheParams->
getSize()>0) {
722 <<
"-dim value cache for integral over " << iset2 <<
" as a function of " << *cacheParams <<
" in range " << (rangeName?rangeName:
"<none>") << endl ;
732 integral = cachedIntegral ;
763 for (
const auto aarg : allObs) {
776 obsWithParamRange.
add(*aarg) ;
777 obsWithFixedRange.
remove(*aarg) ;
778 obsServingAsRangeParams.
add(*loBoundObs,
kFALSE) ;
779 obsServingAsRangeParams.
add(*hiBoundObs,
kFALSE) ;
788 RooArgSet obsWithFixedRangeNP(obsWithFixedRange) ;
789 obsWithFixedRangeNP.
remove(obsServingAsRangeParams) ;
792 RooArgSet obsWithParamRangeNP(obsWithParamRange) ;
793 obsWithParamRangeNP.
remove(obsServingAsRangeParams) ;
797 innerObs.
add(obsWithFixedRangeNP) ;
798 innerObs.
add(obsWithParamRangeNP) ;
815 name.Append(
"_Int[") ;
830 name.Append(rangeName) ;
833 }
else if (!omitEmpty) {
834 name.Append(
"_Int[]") ;
837 if (nset && nset->
getSize()>0 ) {
842 name.Append(
"_Norm[") ;
929 for (
const auto arg : dependentVars) {
930 if(!arg->isFundamental() && !
dynamic_cast<const RooAbsLValue*
>(arg)) {
932 <<
"\" of wrong type: " << arg->ClassName() << endl;
939 <<
"\" is not a dependent and will be ignored." << endl;
944 leafNodes.
replace(*found,*arg);
946 leafNodes.
add(*arg) ;
949 RooArgSet* lvDep = arg->getObservables(&leafNodes) ;
950 for (
const auto lvs : *lvDep) {
954 leafNodes.
add(*lvs) ;
961 if(0 != projectedVars && projectedVars->
find(arg->GetName())) {
963 <<
"\" cannot be both a dependent and a projected variable." << endl;
969 if(0 != projectedVars) leafNodes.
remove(*projectedVars,
kTRUE);
974 coutE(
Plotting) <<
"RooAbsPdf::createPlotProjection(" <<
GetName() <<
") Couldn't deep-clone PDF, abort," << endl ;
987 delete plotLeafNodes ;
991 if(0 != projectedVars) normSet.
add(*projectedVars);
1000 if(0 == projectedVars) projectedVars= ∅
1006 title.
Prepend(
"Projection of ");
1011 if(0 == projected || !projected->
isValid()) {
1015 if(0 != projected)
delete projected;
1020 static_cast<RooRealIntegral*
>(projected)->setAllowComponentSelection(
true);
1063 if(hdim != plotVars.
getSize()) {
1072 for(
Int_t index= 0; index < plotVars.
getSize(); index++) {
1077 <<
"\" of type " << var->
ClassName() << endl;
1082 <<
":fillHistogram: WARNING: variable is not an explicit dependent: " << realVar->
GetName() << endl;
1098 if (projectedVars) {
1099 allDeps.
add(*projectedVars) ;
1113 Int_t xbins(0),ybins(1),zbins(1);
1125 assert(0 != zvar && 0 != zaxis);
1126 if (scaleForDensity) {
1134 assert(0 != yvar && 0 != yaxis);
1135 if (scaleForDensity) {
1143 assert(0 != xvar && 0 != xaxis);
1144 if (scaleForDensity) {
1150 << hdim <<
" dimensions" << endl;
1157 Int_t xbin(0),ybin(0),zbin(0);
1158 Int_t bins= xbins*ybins*zbins;
1159 for(
Int_t bin= 0; bin < bins; bin++) {
1162 if(bin % (xbins*ybins) == 0) {
1168 if(bin % xbins == 0) {
1169 ybin= (ybin%ybins) + 1;
1174 xbin= (xbin%xbins) + 1;
1184 coutW(
Plotting) <<
"WARNING: Function evaluation error(s) at coordinates [x]=" << xvar->
getVal() ;
1257 if (showProgress && (i%onePct==0)) {
1261 Double_t binVal = theClone->
getVal(normSet?normSet:obs)*scaleFactor ;
1262 if (correctForBinSize) {
1294 strlcpy(buf,varNameList,1024) ;
1295 char* varName = strtok(buf,
",:") ;
1300 varName = strtok(0,
",") ;
1302 varName = strtok(0,
",") ;
1394 pc.defineInt(
"scaling",
"Scaling",0,1) ;
1395 pc.defineInt(
"intBinning",
"IntrinsicBinning",0,2) ;
1396 pc.defineInt(
"extended",
"Extended",0,2) ;
1398 pc.defineObject(
"compSet",
"SelectCompSet",0) ;
1399 pc.defineString(
"compSpec",
"SelectCompSpec",0) ;
1400 pc.defineSet(
"projObs",
"ProjectedObservables",0,0) ;
1401 pc.defineObject(
"yvar",
"YVar",0,0) ;
1402 pc.defineObject(
"zvar",
"ZVar",0,0) ;
1403 pc.defineMutex(
"SelectCompSet",
"SelectCompSpec") ;
1404 pc.defineMutex(
"IntrinsicBinning",
"Binning") ;
1405 pc.defineMutex(
"IntrinsicBinning",
"BinningName") ;
1406 pc.defineMutex(
"IntrinsicBinning",
"BinningSpec") ;
1407 pc.allowUndefined() ;
1410 pc.process(argList) ;
1428 Bool_t doScaling =
pc.getInt(
"scaling") ;
1429 Int_t doIntBinning =
pc.getInt(
"intBinning") ;
1430 Int_t doExtended =
pc.getInt(
"extended") ;
1434 if (!pdfSelf && doExtended>0) {
1435 coutW(
InputArguments) <<
"RooAbsReal::createHistogram(" <<
GetName() <<
") WARNING extended mode requested for a non-pdf object, ignored" << endl ;
1439 coutW(
InputArguments) <<
"RooAbsReal::createHistogram(" <<
GetName() <<
") WARNING extended mode requested for a non-extendable pdf, ignored" << endl ;
1442 if (pdfSelf && doExtended==2) {
1446 const char* compSpec =
pc.getString(
"compSpec") ;
1448 Bool_t haveCompSel = ( (compSpec && strlen(compSpec)>0) || compSet) ;
1451 if (doIntBinning>0) {
1456 if (doIntBinning==1) {
1458 <<
") WARNING, intrinsic model binning requested for histogram, but model does not define bin boundaries, reverting to default binning"<< endl ;
1461 if (doIntBinning==2) {
1463 <<
") INFO: Model has intrinsic binning definition, selecting that binning for the histogram"<< endl ;
1466 for (list<double>::iterator it=bl->begin() ; it!=bl->end() ; ++it) { ba[i++] = *it ; }
1467 intBinning =
new RooBinning(bl->size()-1,ba) ;
1473 pc.stripCmdList(argListCreate,
"Scaling,ProjectedObservables,IntrinsicBinning,SelectCompSet,SelectCompSpec,Extended") ;
1478 argListCreate.
Add(&tmp) ;
1497 branchNodeSet.
remove(*arg) ;
1508 if (dirSelNodes->
getSize()>0) {
1509 coutI(
Plotting) <<
"RooAbsPdf::createHistogram(" <<
GetName() <<
") directly selected PDF components: " << *dirSelNodes << endl ;
1515 coutE(
Plotting) <<
"RooAbsPdf::createHistogram(" <<
GetName() <<
") ERROR: component selection set " << *compSet <<
" does not match any components of p.d.f." << endl ;
1517 coutE(
Plotting) <<
"RooAbsPdf::createHistogram(" <<
GetName() <<
") ERROR: component selection expression '" << compSpec <<
"' does not select any components of p.d.f." << endl ;
1521 delete dirSelNodes ;
1557 for (
unsigned int i = 0; i < branchNodeSet.
size(); ++i) {
1558 const auto arg = branchNodeSet[i];
1560 branchNodeSet.
remove(*arg) ;
1567 for (
const auto arg : branchNodeSet) {
1576 for (
const auto arg : branchNodeSet) {
1577 for (
const auto selNode : *selNodes) {
1578 if (selNode->dependsOn(*arg)) {
1585 for (
const auto arg : branchNodeSet) {
1586 if (arg->dependsOn(*selNodes)) {
1591 tmp.
remove(*selNodes,
true);
1594 coutI(
Plotting) <<
"RooAbsPdf::plotOn(" <<
GetName() <<
") indirectly selected PDF components: " << tmp << endl ;
1597 for (
const auto arg : branchNodeSet) {
1598 Bool_t select = selNodes->
find(arg->GetName()) !=
nullptr;
1743 argList.
Add(&rnorm) ;
1745 std::vector<string> rlist;
1749 rlist.emplace_back(rangeNameToken);
1752 for (
const auto& rangeString : rlist) {
1754 rcmd->
setString(0, rangeString.c_str());
1763 pc.defineString(
"drawOption",
"DrawOption",0,
"L") ;
1764 pc.defineString(
"projectionRangeName",
"ProjectionRange",0,
"",
kTRUE) ;
1765 pc.defineString(
"curveNameSuffix",
"CurveNameSuffix",0,
"") ;
1766 pc.defineString(
"sliceCatState",
"SliceCat",0,
"",
kTRUE) ;
1767 pc.defineDouble(
"scaleFactor",
"Normalization",0,1.0) ;
1768 pc.defineInt(
"scaleType",
"Normalization",0,
Relative) ;
1769 pc.defineObject(
"sliceSet",
"SliceVars",0) ;
1770 pc.defineObject(
"sliceCatList",
"SliceCat",0,0,
kTRUE) ;
1771 pc.defineObject(
"projSet",
"Project",0) ;
1772 pc.defineObject(
"asymCat",
"Asymmetry",0) ;
1773 pc.defineDouble(
"precision",
"Precision",0,1
e-3) ;
1774 pc.defineDouble(
"evalErrorVal",
"EvalErrorValue",0,0) ;
1775 pc.defineInt(
"doEvalError",
"EvalErrorValue",0,0) ;
1776 pc.defineInt(
"shiftToZero",
"ShiftToZero",0,0) ;
1777 pc.defineObject(
"projDataSet",
"ProjData",0) ;
1778 pc.defineObject(
"projData",
"ProjData",1) ;
1779 pc.defineObject(
"errorFR",
"VisualizeError",0) ;
1780 pc.defineDouble(
"errorZ",
"VisualizeError",0,1.) ;
1781 pc.defineSet(
"errorPars",
"VisualizeError",0) ;
1782 pc.defineInt(
"linearMethod",
"VisualizeError",0,0) ;
1783 pc.defineInt(
"binProjData",
"ProjData",0,0) ;
1784 pc.defineDouble(
"rangeLo",
"Range",0,-999.) ;
1785 pc.defineDouble(
"rangeHi",
"Range",1,-999.) ;
1786 pc.defineInt(
"numee",
"PrintEvalErrors",0,10) ;
1787 pc.defineInt(
"rangeAdjustNorm",
"Range",0,0) ;
1788 pc.defineInt(
"rangeWNAdjustNorm",
"RangeWithName",0,0) ;
1789 pc.defineInt(
"VLines",
"VLines",0,2) ;
1790 pc.defineString(
"rangeName",
"RangeWithName",0,
"") ;
1791 pc.defineString(
"normRangeName",
"NormRange",0,
"") ;
1792 pc.defineInt(
"markerColor",
"MarkerColor",0,-999) ;
1793 pc.defineInt(
"markerStyle",
"MarkerStyle",0,-999) ;
1794 pc.defineDouble(
"markerSize",
"MarkerSize",0,-999) ;
1795 pc.defineInt(
"lineColor",
"LineColor",0,-999) ;
1796 pc.defineInt(
"lineStyle",
"LineStyle",0,-999) ;
1797 pc.defineInt(
"lineWidth",
"LineWidth",0,-999) ;
1798 pc.defineInt(
"fillColor",
"FillColor",0,-999) ;
1799 pc.defineInt(
"fillStyle",
"FillStyle",0,-999) ;
1800 pc.defineString(
"curveName",
"Name",0,
"") ;
1801 pc.defineInt(
"curveInvisible",
"Invisible",0,0) ;
1802 pc.defineInt(
"showProg",
"ShowProgress",0,0) ;
1803 pc.defineInt(
"numCPU",
"NumCPU",0,1) ;
1804 pc.defineInt(
"interleave",
"NumCPU",1,0) ;
1805 pc.defineString(
"addToCurveName",
"AddTo",0,
"") ;
1806 pc.defineDouble(
"addToWgtSelf",
"AddTo",0,1.) ;
1807 pc.defineDouble(
"addToWgtOther",
"AddTo",1,1.) ;
1808 pc.defineInt(
"moveToBack",
"MoveToBack",0,0) ;
1809 pc.defineMutex(
"SliceVars",
"Project") ;
1810 pc.defineMutex(
"AddTo",
"Asymmetry") ;
1811 pc.defineMutex(
"Range",
"RangeWithName") ;
1812 pc.defineMutex(
"VisualizeError",
"VisualizeErrorData") ;
1815 pc.process(argList) ;
1821 TString drawOpt(
pc.getString(
"drawOption"));
1826 Bool_t linMethod =
pc.getInt(
"linearMethod") ;
1827 if (!drawOpt.
Contains(
"P") && errFR) {
1834 o.
numee =
pc.getInt(
"numee") ;
1844 o.
eeval =
pc.getDouble(
"evalErrorVal") ;
1854 const char* sliceCatState =
pc.getString(
"sliceCatState",0,
kTRUE) ;
1855 const RooLinkedList& sliceCatList =
pc.getObjectList(
"sliceCatList") ;
1856 if (sliceCatState) {
1865 strlcpy(buf,sliceCatState,1024) ;
1866 const char* slabel = strtok(buf,
",") ;
1878 slabel = strtok(0,
",") ;
1885 Int_t vlines =
pc.getInt(
"VLines");
1886 if (
pc.hasProcessed(
"Range")) {
1890 if (vlines==2) vlines=0 ;
1891 }
else if (
pc.hasProcessed(
"RangeWithName")) {
1896 if (vlines==2) vlines=0 ;
1901 if (
pc.hasProcessed(
"NormRange")) {
1922 cxcoutD(
Plotting) <<
"RooAbsReal::plotOn(" <<
GetName() <<
") Preprocessing: have slice " << *sliceSet << endl ;
1932 projectedVars.
remove(*arg) ;
1935 << sliceArg->
GetName() <<
" was not projected anyway" << endl ;
1939 }
else if (projSet) {
1940 cxcoutD(
Plotting) <<
"RooAbsReal::plotOn(" <<
GetName() <<
") Preprocessing: have projSet " << *projSet << endl ;
1943 cxcoutD(
Plotting) <<
"RooAbsReal::plotOn(" <<
GetName() <<
") Preprocessing: have neither sliceSet nor projSet " << endl ;
1948 cxcoutD(
Plotting) <<
"RooAbsReal::plotOn(" <<
GetName() <<
") Preprocessing: projectedVars = " << projectedVars << endl ;
1963 Int_t lineColor =
pc.getInt(
"lineColor") ;
1964 Int_t lineStyle =
pc.getInt(
"lineStyle") ;
1965 Int_t lineWidth =
pc.getInt(
"lineWidth") ;
1966 Int_t markerColor =
pc.getInt(
"markerColor") ;
1967 Int_t markerStyle =
pc.getInt(
"markerStyle") ;
1968 Size_t markerSize =
pc.getDouble(
"markerSize") ;
1969 Int_t fillColor =
pc.getInt(
"fillColor") ;
1970 Int_t fillStyle =
pc.getInt(
"fillStyle") ;
1980 if ((fillColor != -999 || fillStyle != -999) && !drawOpt.
Contains(
"F")) {
1982 <<
"\", but these only have an effect when 'DrawOption(\"F\")' for fill is used at the same time." << std::endl;
1986 if (
pc.getInt(
"moveToBack") && frame->
numItems()>1) {
2025 projDataVars.
add(*tmp) ;
2048 cxcoutD(
Plotting) <<
"RooAbsReal::plotOn(" <<
GetName() <<
") frame->getNormVars() that are also observables = " << *sliceSetTmp << endl ;
2061 << frame->
getPlotVar()->
GetName() <<
" represents a slice in " << *sliceSetTmp << endl ;
2063 sliceSet.
add(*sliceSetTmp) ;
2064 delete sliceSetTmp ;
2070 cxcoutD(
Plotting) <<
"RooAbsReal::plotOn(" <<
GetName() <<
") projectedVars = " << projectedVars <<
" sliceSet = " << sliceSet << endl ;
2083 if (!plotCloneSet) {
2084 coutE(
Plotting) <<
"RooAbsReal::plotOn(" <<
GetName() <<
") Couldn't deep-clone self, abort," << endl ;
2090 if (projectedVars.
getSize()) {
2092 <<
" integrates over variables " << projectedVars
2095 if (projDataNeededVars && projDataNeededVars->
getSize()>0) {
2097 <<
" averages using data variables " << *projDataNeededVars << endl ;
2105 if (projDataNeededVars) {
2109 deps->
add(*plotVar) ;
2115 coutE(
Plotting) <<
"RooAbsReal::plotOn(" <<
GetName() <<
") error in checkObservables, abort" << endl ;
2117 delete plotCloneSet ;
2118 if (projDataNeededVars)
delete projDataNeededVars ;
2133 fullNormSet.
add(projectedVars) ;
2134 if (projDataNeededVars && projDataNeededVars->
getSize()>0) {
2135 fullNormSet.
add(*projDataNeededVars) ;
2151 if (o.
projData && projDataNeededVars && projDataNeededVars->
getSize()>0) {
2161 if (sliceDataSet->
getSize()>0) {
2182 delete sliceDataSet ;
2184 if (!cutString.
IsNull()) {
2186 coutI(
Plotting) <<
"RooAbsReal::plotOn(" <<
GetName() <<
") reducing given projection dataset to entries with " << cutString << endl ;
2191 <<
") only the following components of the projection data will be used: " << *projDataNeededVars << endl ;
2207 coutI(
Plotting) <<
"RooAbsReal::plotOn(" <<
GetName() <<
") unbinned projection dataset consist only of discrete variables,"
2208 <<
" performing projection with binned copy for optimization." << endl ;
2216 if (projDataSel!=o.
projData)
delete projDataSel ;
2223 projection->
getVal(projDataSel->
get()) ;
2281 if (projDataSel!=o.
projData)
delete projDataSel ;
2364 if (projDataNeededVars)
delete projDataNeededVars ;
2366 delete projectionCompList ;
2367 delete plotCloneSet ;
2389 projectedVars.
remove(*arg) ;
2392 << sliceArg->
GetName() <<
" was not projected anyway" << endl ;
2435 projDataVars.
add(*tmp) ;
2445 <<
") function doesn't depend on asymmetry category " << asymCat.
GetName() << endl ;
2452 <<
") asymmetry category must have 2 or 3 states with index values -1,0,1" << endl ;
2476 << frame->
getPlotVar()->
GetName() <<
" represents a slice in " << *sliceSetTmp << endl ;
2478 sliceSet.
add(*sliceSetTmp) ;
2479 delete sliceSetTmp ;
2503 if (projectedVars.
getSize()) {
2505 <<
" projects variables " << projectedVars << endl ;
2507 if (projDataNeededVars && projDataNeededVars->
getSize()>0) {
2509 <<
" averages using data variables "<< *projDataNeededVars << endl ;
2528 RooArgSet *posProjCompList, *negProjCompList ;
2534 depPos.
add(projDataVars) ;
2535 depNeg.
add(projDataVars) ;
2539 if (!posProj || !negProj) {
2540 coutE(
Plotting) <<
"RooAbsReal::plotAsymOn(" <<
GetName() <<
") Unable to create projections, abort" << endl ;
2546 asymName.
Append(
"_Asym[") ;
2550 asymTitle.
Append(
" Asymmetry of ") ;
2563 if (sliceDataSet->
getSize()>0) {
2584 delete sliceDataSet ;
2586 if (!cutString.
IsNull()) {
2589 <<
") reducing given projection dataset to entries with " << cutString << endl ;
2594 <<
") only the following components of the projection data will be used: " << *projDataNeededVars << endl ;
2604 ((
RooAbsReal*)posProj)->attachDataSet(*projDataSel) ;
2605 ((
RooAbsReal*)negProj)->attachDataSet(*projDataSel) ;
2639 if (projDataSel!=o.
projData)
delete projDataSel ;
2678 delete posProjCompList ;
2679 delete negProjCompList ;
2711 fpf_stripped.
add(*frv);
2725 vector<int> fpf_idx;
2729 paramList.
add(*par);
2730 fpf_idx.push_back(i);
2734 vector<Double_t> plusVar, minusVar ;
2741 for (
Int_t ivar=0 ; ivar<paramList.
getSize() ; ivar++) {
2749 ((
RooRealVar*)paramList.
at(ivar))->setVal(cenVal+errVal) ;
2750 plusVar.push_back(cloneFunc->
getVal(nset)) ;
2753 ((
RooRealVar*)paramList.
at(ivar))->setVal(cenVal-errVal) ;
2754 minusVar.push_back(cloneFunc->
getVal(nset)) ;
2760 vector<double> errVec(paramList.
getSize()) ;
2761 for (
int i=0 ; i<paramList.
getSize() ; i++) {
2762 errVec[i] =
sqrt(V(i,i)) ;
2763 for (
int j=i ; j<paramList.
getSize() ; j++) {
2764 C(i,j) = V(i,j)/
sqrt(V(i,i)*V(j,j)) ;
2771 for (
unsigned int j=0 ; j<plusVar.size() ; j++) {
2772 F[j] = (plusVar[j]-minusVar[j])/2 ;
2779 delete errorParams ;
2831 pc.stripCmdList(plotArgListTmp,
"VisualizeError,MoveToBack") ;
2838 if (std::string(
"Normalization")==cmd->
GetName()) {
2841 plotArgList.
Add(cmd) ;
2844 plotArgList.
Add(cmd) ;
2875 coutI(
Plotting) <<
"RooAbsReal::plotOn(" <<
GetName() <<
") INFO: visualizing " << Z <<
"-sigma uncertainties in parameters "
2876 << *errorParams <<
" from fit result " << fr.
GetName() <<
" using " <<
n <<
" samplings." << endl ;
2882 vector<RooCurve*> cvec ;
2883 for (
int i=0 ; i<
d->numEntries() ; i++) {
2884 *cloneParams = (*
d->get(i)) ;
2886 cloneFunc->
plotOn(frame,tmp2) ;
2887 cvec.push_back(frame->
getCurve()) ;
2900 for (vector<RooCurve*>::iterator i=cvec.begin() ; i!=cvec.end() ; ++i) {
2922 fpf_stripped.
add(*frv);
2935 vector<int> fpf_idx ;
2939 paramList.add(*par) ;
2940 fpf_idx.push_back(i) ;
2944 vector<RooCurve*> plusVar, minusVar ;
2953 for (
Int_t ivar=0 ; ivar<paramList.getSize() ; ivar++) {
2961 ((
RooRealVar*)paramList.at(ivar))->setVal(cenVal+Z*errVal) ;
2965 cloneFunc->
plotOn(frame,tmp2) ;
2966 plusVar.push_back(frame->
getCurve()) ;
2971 ((
RooRealVar*)paramList.at(ivar))->setVal(cenVal-Z*errVal) ;
2973 cloneFunc->
plotOn(frame,tmp3) ;
2974 minusVar.push_back(frame->
getCurve()) ;
2977 ((
RooRealVar*)paramList.at(ivar))->setVal(cenVal) ;
2981 vector<double> errVec(paramList.getSize()) ;
2982 for (
int i=0 ; i<paramList.getSize() ; i++) {
2983 errVec[i] =
sqrt(V(i,i)) ;
2984 for (
int j=i ; j<paramList.getSize() ; j++) {
2985 C(i,j) = V(i,j)/
sqrt(V(i,i)*V(j,j)) ;
2995 for (vector<RooCurve*>::iterator i=plusVar.begin() ; i!=plusVar.end() ; ++i) {
2998 for (vector<RooCurve*>::iterator i=minusVar.begin() ; i!=minusVar.end() ; ++i) {
3005 if (!band)
return frame ;
3008 pc.defineString(
"drawOption",
"DrawOption",0,
"F") ;
3009 pc.defineString(
"curveNameSuffix",
"CurveNameSuffix",0,
"") ;
3010 pc.defineInt(
"lineColor",
"LineColor",0,-999) ;
3011 pc.defineInt(
"lineStyle",
"LineStyle",0,-999) ;
3012 pc.defineInt(
"lineWidth",
"LineWidth",0,-999) ;
3013 pc.defineInt(
"markerColor",
"MarkerColor",0,-999) ;
3014 pc.defineInt(
"markerStyle",
"MarkerStyle",0,-999) ;
3015 pc.defineDouble(
"markerSize",
"MarkerSize",0,-999) ;
3016 pc.defineInt(
"fillColor",
"FillColor",0,-999) ;
3017 pc.defineInt(
"fillStyle",
"FillStyle",0,-999) ;
3018 pc.defineString(
"curveName",
"Name",0,
"") ;
3019 pc.defineInt(
"curveInvisible",
"Invisible",0,0) ;
3020 pc.defineInt(
"moveToBack",
"MoveToBack",0,0) ;
3021 pc.allowUndefined() ;
3024 pc.process(argList) ;
3030 frame->
addPlotable(band,
pc.getString(
"drawOption"),
pc.getInt(
"curveInvisible")) ;
3033 Int_t lineColor =
pc.getInt(
"lineColor") ;
3034 Int_t lineStyle =
pc.getInt(
"lineStyle") ;
3035 Int_t lineWidth =
pc.getInt(
"lineWidth") ;
3036 Int_t markerColor =
pc.getInt(
"markerColor") ;
3037 Int_t markerStyle =
pc.getInt(
"markerStyle") ;
3038 Size_t markerSize =
pc.getDouble(
"markerSize") ;
3039 Int_t fillColor =
pc.getInt(
"fillColor") ;
3040 Int_t fillStyle =
pc.getInt(
"fillStyle") ;
3051 if (
pc.getString(
"curveName",0,
kTRUE)) {
3053 }
else if (
pc.getString(
"curveNameSuffix",0,
kTRUE)) {
3055 name.Append(
pc.getString(
"curveNameSuffix",0,
kTRUE)) ;
3060 if (
pc.getInt(
"moveToBack") && frame->
numItems()>1) {
3086 <<
":plotOn: frame does not specify a plot variable" << endl;
3121 <<
" allVars = " << (allVars?(*allVars):
RooArgSet()) << endl ;
3124 if (!allVars) return ;
3127 projectedVars.
add(*allVars) ;
3132 projectedVars.
remove(*found);
3142 <<
" from projection set because it a server of " << plotVar->
GetName() << endl ;
3143 projectedVars.
remove(*tmp) ;
3147 delete plotServers ;
3151 <<
") WARNING: cannot project out frame variable ("
3152 << found->
GetName() <<
"), ignoring" << endl ;
3164 <<
") function doesn't depend on projection variable "
3165 << arg->
GetName() <<
", ignoring" << endl ;
3204 if(binding && !binding->
isValid()) {
3214struct TreeReadBuffer {
3215 virtual ~TreeReadBuffer() =
default;
3216 virtual operator double() = 0;
3227 auto other =
static_cast<const RooAbsReal*
>(source);
3228 assert(
dynamic_cast<const RooAbsReal*
>(source));
3230 _value = other->_treeReadBuffer ? other->_treeReadBuffer->operator
double() : other->_value;
3253struct TypedTreeReadBuffer final :
public TreeReadBuffer {
3254 operator double()
override {
3265std::unique_ptr<TreeReadBuffer> createTreeReadBuffer(
const TString& branchName,
TTree&
tree) {
3266 auto buf =
new TypedTreeReadBuffer<T>();
3267 tree.SetBranchAddress(branchName.
Data(), &buf->_value);
3268 return std::unique_ptr<TreeReadBuffer>(buf);
3299 <<
" is an array and cannot be attached to a RooAbsReal" << endl ;
3309 std::map<std::string, std::pair<std::string, std::function<std::unique_ptr<TreeReadBuffer>()>>> typeMap {
3310 {
"Float_t", {
"FLOAT_TREE_BRANCH", [&](){
return createTreeReadBuffer<Float_t >(cleanName, t); }}},
3311 {
"Int_t", {
"INTEGER_TREE_BRANCH", [&](){
return createTreeReadBuffer<Int_t >(cleanName, t); }}},
3312 {
"UChar_t", {
"BYTE_TREE_BRANCH", [&](){
return createTreeReadBuffer<UChar_t >(cleanName, t); }}},
3313 {
"Bool_t", {
"BOOL_TREE_BRANCH", [&](){
return createTreeReadBuffer<Bool_t >(cleanName, t); }}},
3314 {
"Char_t", {
"SIGNEDBYTE_TREE_BRANCH", [&](){
return createTreeReadBuffer<Char_t >(cleanName, t); }}},
3315 {
"UInt_t", {
"UNSIGNED_INTEGER_TREE_BRANCH", [&](){
return createTreeReadBuffer<UInt_t >(cleanName, t); }}},
3316 {
"Long64_t", {
"LONG_TREE_BRANCH", [&](){
return createTreeReadBuffer<Long64_t >(cleanName, t); }}},
3317 {
"ULong64_t", {
"UNSIGNED_LONG_TREE_BRANCH", [&](){
return createTreeReadBuffer<ULong64_t>(cleanName, t); }}},
3318 {
"Short_t", {
"SHORT_TREE_BRANCH", [&](){
return createTreeReadBuffer<Short_t >(cleanName, t); }}},
3319 {
"UShort_t", {
"UNSIGNED_SHORT_TREE_BRANCH", [&](){
return createTreeReadBuffer<UShort_t >(cleanName, t); }}},
3322 auto typeDetails = typeMap.find(typeName.
Data());
3323 if (typeDetails != typeMap.end()) {
3325 <<
" will be converted to double precision." << endl ;
3335 coutE(
InputArguments) <<
"RooAbsReal::attachToTree(" <<
GetName() <<
") data type " << typeName <<
" is not supported." << endl ;
3498 const TList &nameList)
const
3507 matched.
add(*found);
3521 if(isMatched) matchedArgs.
add(matched);
3570 if (config)
return config ;
3583 if (config)
return config ;
3685 if (inLogEvalError) {
3688 inLogEvalError =
kTRUE ;
3693 if (serverValueString) {
3698 oocoutE((
TObject*)0,
Eval) <<
"RooAbsReal::logEvalError(" <<
"<STATIC>" <<
") evaluation error, " << endl
3699 <<
" origin : " << origName << endl
3700 <<
" message : " << ee.
_msg << endl
3701 <<
" server values: " << ee.
_srvval << endl ;
3708 inLogEvalError =
kFALSE ;
3740 if (inLogEvalError) {
3743 inLogEvalError =
kTRUE ;
3748 if (serverValueString) {
3768 ostringstream oss2 ;
3772 coutE(
Eval) <<
"RooAbsReal::logEvalError(" <<
GetName() <<
") evaluation error, " << endl
3773 <<
" origin : " << oss2.str() << endl
3774 <<
" message : " << ee.
_msg << endl
3775 <<
" server values: " << ee.
_srvval << endl ;
3784 <<
") delayed evaluation error, " << endl
3785 <<
" origin : " << oss2.str() << endl
3786 <<
" message : " << oee.
_msg << endl
3787 <<
" server values: " << oee.
_srvval << endl ;
3794 inLogEvalError =
kFALSE ;
3830 if (maxPerNode<0) return ;
3832 map<const RooAbsArg*,pair<string,list<EvalError> > >::iterator iter =
_evalErrorList.begin() ;
3835 if (maxPerNode==0) {
3838 os << iter->second.first ;
3840 os <<
" has " << iter->second.second.size() <<
" errors" << endl ;
3845 os << iter->second.first << endl ;
3849 std::list<EvalError>::iterator iter2 = iter->second.second.begin() ;
3850 for(;iter2!=iter->second.second.end() ; ++iter2, i++) {
3851 os <<
" " << iter2->_msg <<
" @ " << iter2->_srvval << endl ;
3853 os <<
" ... (remaining " << iter->second.second.size() - maxPerNode <<
" messages suppressed)" << endl ;
3873 map<const RooAbsArg*,pair<string,list<EvalError> > >::iterator iter =
_evalErrorList.begin() ;
3875 ntot += iter->second.second.size() ;
3950 orderedObs.
add(obs) ;
4008 pc.defineObject(
"supNormSet",
"SupNormSet",0,0) ;
4009 pc.defineInt(
"numScanBins",
"ScanParameters",0,1000) ;
4010 pc.defineInt(
"intOrder",
"ScanParameters",1,2) ;
4011 pc.defineInt(
"doScanNum",
"ScanNum",0,1) ;
4012 pc.defineInt(
"doScanAll",
"ScanAll",0,0) ;
4013 pc.defineInt(
"doScanNon",
"ScanNone",0,0) ;
4014 pc.defineMutex(
"ScanNum",
"ScanAll",
"ScanNone") ;
4017 pc.process(arg1,arg2,arg3,arg4,arg5,arg6,arg7,arg8) ;
4028 Int_t numScanBins =
pc.getInt(
"numScanBins") ;
4029 Int_t intOrder =
pc.getInt(
"intOrder") ;
4030 Int_t doScanNum =
pc.getInt(
"doScanNum") ;
4031 Int_t doScanAll =
pc.getInt(
"doScanAll") ;
4032 Int_t doScanNon =
pc.getInt(
"doScanNon") ;
4047 coutI(
NumIntegration) <<
"RooAbsPdf::createRunningIntegral(" <<
GetName() <<
") integration over observable(s) " << iset <<
" involves numeric integration," << endl
4048 <<
" constructing cdf though numeric integration of sampled pdf in " << numScanBins <<
" bins and applying order "
4049 << intOrder <<
" interpolation on integrated histogram." << endl
4050 <<
" To override this choice of technique use argument ScanNone(), to change scan parameters use ScanParameters(nbins,order) argument" << endl ;
4068 ivar->
setBins(numScanBins,
"numcdf") ;
4112 cloneList.
add(*cloneArg) ;
4118 loList.
add(*cloneLo) ;
4151 coutE(
InputArguments) <<
"RooAbsReal::functor(" <<
GetName() <<
") ERROR: one or more specified observables are not variables of this p.d.f" << endl ;
4157 coutE(
InputArguments) <<
"RooAbsReal::functor(" <<
GetName() <<
") ERROR: one or more specified parameters are not variables of this p.d.f" << endl ;
4178 coutE(
InputArguments) <<
"RooAbsReal::functor(" <<
GetName() <<
") ERROR: one or more specified observables are not variables of this p.d.f" << endl ;
4184 coutE(
InputArguments) <<
"RooAbsReal::functor(" <<
GetName() <<
") ERROR: one or more specified parameters are not variables of this p.d.f" << endl ;
4192 for (
int i=0 ; i<obs.
getSize() ; i++) {
4198 for (
int i=0 ; i<pars.
getSize() ; i++) {
4232 <<
" observables specified, but a ROOT TFx can only have 1,2 or 3 observables" << endl ;
4237 for (
int i=0 ; i<pars.
getSize() ; i++) {
4267 return new RooDerivative(
name.c_str(),title.c_str(),*
this,obs,normSet,order,eps) ;
4282 string title=
Form(
"%sMoment of order %d of %s w.r.t %s ",(central?
"Central ":
""),order,
GetName(),obs.
GetName()) ;
4284 if (order==2)
return new RooSecondMoment(
name.c_str(),title.c_str(),*
this,obs,central,takeRoot) ;
4285 return new RooMoment(
name.c_str(),title.c_str(),*
this,obs,order,central,takeRoot) ;
4301 string title=
Form(
"%sMoment of order %d of %s w.r.t %s ",(central?
"Central ":
""),order,
GetName(),obs.
GetName()) ;
4303 if (order==1)
return new RooFirstMoment(
name.c_str(),title.c_str(),*
this,obs,normObs,intNormObs) ;
4304 if (order==2)
return new RooSecondMoment(
name.c_str(),title.c_str(),*
this,obs,normObs,central,takeRoot,intNormObs) ;
4305 return new RooMoment(
name.c_str(),title.c_str(),*
this,obs,normObs,order,central,takeRoot,intNormObs) ;
4405 RooLinkedList chi2CmdList =
pc.filterCmdList(fitCmdList,
"Range,RangeWithName,NumCPU,Optimize") ;
4438 return new RooChi2Var(
name.c_str(),
name.c_str(),*
this,data,arg1,arg2,arg3,arg4,arg5,arg6,arg7,arg8) ;
4463 return createChi2(data,*cmds[0],*cmds[1],*cmds[2],*cmds[3],*cmds[4],*cmds[5],*cmds[6],*cmds[7]) ;
4529 RooLinkedList chi2CmdList =
pc.filterCmdList(fitCmdList,
"YVar,Integrate") ;
4578 pc.defineInt(
"integrate",
"Integrate",0,0) ;
4579 pc.defineObject(
"yvar",
"YVar",0,0) ;
4582 pc.process(cmdList) ;
4588 Bool_t integrate =
pc.getInt(
"integrate") ;
4613 pc.defineString(
"fitOpt",
"FitOptions",0,
"") ;
4615 pc.defineInt(
"optConst",
"Optimize",0,1) ;
4616 pc.defineInt(
"verbose",
"Verbose",0,0) ;
4617 pc.defineInt(
"doSave",
"Save",0,0) ;
4618 pc.defineInt(
"doTimer",
"Timer",0,0) ;
4619 pc.defineInt(
"plevel",
"PrintLevel",0,1) ;
4620 pc.defineInt(
"strat",
"Strategy",0,1) ;
4621 pc.defineInt(
"initHesse",
"InitialHesse",0,0) ;
4622 pc.defineInt(
"hesse",
"Hesse",0,1) ;
4623 pc.defineInt(
"minos",
"Minos",0,0) ;
4624 pc.defineInt(
"ext",
"Extended",0,2) ;
4625 pc.defineInt(
"numee",
"PrintEvalErrors",0,10) ;
4626 pc.defineInt(
"doWarn",
"Warnings",0,1) ;
4627 pc.defineString(
"mintype",
"Minimizer",0,
"Minuit") ;
4628 pc.defineString(
"minalg",
"Minimizer",1,
"minuit") ;
4629 pc.defineObject(
"minosSet",
"Minos",0,0) ;
4631 pc.defineMutex(
"FitOptions",
"Verbose") ;
4632 pc.defineMutex(
"FitOptions",
"Save") ;
4633 pc.defineMutex(
"FitOptions",
"Timer") ;
4634 pc.defineMutex(
"FitOptions",
"Strategy") ;
4635 pc.defineMutex(
"FitOptions",
"InitialHesse") ;
4636 pc.defineMutex(
"FitOptions",
"Hesse") ;
4637 pc.defineMutex(
"FitOptions",
"Minos") ;
4640 pc.process(cmdList) ;
4646 const char* fitOpt =
pc.getString(
"fitOpt",0,
kTRUE) ;
4647#ifdef __ROOFIT_NOROOMINIMIZER
4648 const char* minType =0 ;
4650 const char* minType =
pc.getString(
"mintype",
"Minuit") ;
4651 const char* minAlg =
pc.getString(
"minalg",
"minuit") ;
4653 Int_t optConst =
pc.getInt(
"optConst") ;
4655 Int_t doSave =
pc.getInt(
"doSave") ;
4656 Int_t doTimer =
pc.getInt(
"doTimer") ;
4657 Int_t plevel =
pc.getInt(
"plevel") ;
4658 Int_t strat =
pc.getInt(
"strat") ;
4659 Int_t initHesse=
pc.getInt(
"initHesse") ;
4660 Int_t hesse =
pc.getInt(
"hesse") ;
4661 Int_t minos =
pc.getInt(
"minos") ;
4662 Int_t numee =
pc.getInt(
"numee") ;
4663 Int_t doWarn =
pc.getInt(
"doWarn") ;
4668#ifdef __ROOFIT_NOROOMINIMIZER
4671 if (
"OldMinuit" ==
string(minType)) {
4680 m.setPrintEvalErrors(numee) ;
4682 m.setPrintLevel(plevel) ;
4687 m.optimizeConst(optConst);
4693 ret =
m.fit(fitOpt) ;
4708 m.setStrategy(strat) ;
4727 m.minos(*minosSet) ;
4736 string title =
Form(
"Result of fit of %s ",
GetName()) ;
4737 ret =
m.save(
name.c_str(),title.c_str()) ;
4742#ifndef __ROOFIT_NOROOMINIMIZER
4745 m.setMinimizerType(minType);
4751 m.setPrintEvalErrors(numee) ;
4753 m.setPrintLevel(plevel) ;
4758 m.optimizeConst(optConst);
4764 ret =
m.fit(fitOpt) ;
4779 m.setStrategy(strat) ;
4788 m.minimize(minType, minAlg) ;
4798 m.minos(*minosSet) ;
4807 string title =
Form(
"Result of fit of %s ",
GetName()) ;
4808 ret =
m.save(
name.c_str(),title.c_str()) ;
4850 while((arg=iter.
next())) {
4853 <<
") function does not depend on listed parameter " << arg->
GetName() <<
", ignoring" << endl ;
4856 if (plist.size()>0) plist +=
":" ;
4878 coutI(
Optimization) <<
"The class " << IsA()->GetName() <<
" could benefit from implementing the faster batch evaluation interface."
4879 <<
" If it is part of ROOT, consider requesting this on https://root.cern." << std::endl;
4883 std::vector<std::tuple<RooRealVar*, RooSpan<const double>,
double>> batchLeafs;
4884 for (
auto leaf : allLeafs) {
4885 auto leafRRV =
dynamic_cast<RooRealVar*
>(leaf);
4889 auto leafBatch = leafRRV->
getValBatch(begin, maxSize);
4890 if (leafBatch.empty())
4893 maxSize = std::min(maxSize, leafBatch.size());
4894 batchLeafs.emplace_back(leafRRV, leafBatch, leafRRV->_value);
4897 if (batchLeafs.empty() || maxSize == 0)
4910 for (std::size_t i = 0; i <
output.size(); ++i) {
4911 for (
auto& tup : batchLeafs) {
4913 auto batch = std::get<1>(tup);
4924 for (
auto& tup : batchLeafs) {
4925 std::get<0>(tup)->setVal(std::get<2>(tup));
4942 const bool tmpFast =
_fast;
4943 const double tmp =
_value;
4945 double fullEval = 0.;
4947 fullEval =
getValV(normalisationSet);
4956 if (std::isfinite(ret) && ( ret != 0. ? (ret - fullEval)/ret : ret - fullEval) > 1.E-9) {
4959 formatter <<
"--> (Scalar computation wrong here:)\n"
4960 <<
GetName() <<
" " <<
this <<
" _fast=" << tmpFast
4961 <<
"\n\tcached _value=" << std::setprecision(16) << tmp
4962 <<
"\n\treturning =" << ret
4963 <<
"\n\trecomputed =" << fullEval
4964 <<
"\n\tnew _value =" <<
_value <<
"] ";
4965 formatter <<
"\nServers:";
4981 auto realServer =
dynamic_cast<RooAbsReal*
>(server);
4993 const double batchVal = batch[0];
4998 formatter <<
"--> (Batch computation wrong here:)\n";
5000 formatter << std::setprecision(17)
5001 <<
"\n _batch[" << std::setw(7) << evtNo-1 <<
"]= " << (enclosingBatch.
empty() ? 0 : enclosingBatch[0])
5002 <<
"\n _batch[" << std::setw(7) << evtNo <<
"]= " << batchVal <<
" !!!"
5003 <<
"\n expected ('_value'): " <<
_value
5004 <<
"\n delta " <<
" = " <<
_value - batchVal
5005 <<
"\n rel delta " <<
" = " << relDiff
5006 <<
"\n _batch[" << std::setw(7) << evtNo+1 <<
"]= " << (enclosingBatch.
empty() ? 0 : enclosingBatch[2]);
5008 formatter <<
"\n" << std::left << std::setw(24) <<
"evaluate(unnorm.)" <<
'=' <<
evaluate();
5010 formatter <<
"\nServers: ";
5012 formatter <<
"\n - ";
5014 formatter << std::setprecision(17);
5016 auto serverAsReal =
dynamic_cast<RooAbsReal*
>(server);
5020 if (!theBatch.
empty()) {
5021 formatter <<
"\n _batch[" << evtNo-1 <<
"]=" << theBatch[0]
5022 <<
"\n _batch[" << evtNo <<
"]=" << theBatch[1]
5023 <<
"\n _batch[" << evtNo+1 <<
"]=" << theBatch[2];
5026 formatter << std::setprecision(17)
5027 <<
"\n getVal()=" << serverAsReal->getVal(normSet);
static RooMathCoreReg dummy
static void indent(ostringstream &buf, int indent_level)
char * Form(const char *fmt,...)
R__EXTERN TSystem * gSystem
A class to store batches of data points that can be accessed via RooSpan.
bool setStatus(std::size_t begin, std::size_t size, Status_t stat, const RooArgSet *const normSet=nullptr, Tag_t ownerTag=kUnspecified)
Set the status of a batch with the given properties.
RooSpan< const double > getBatch(std::size_t begin, std::size_t maxSize, const RooArgSet *const normSet=nullptr, Tag_t ownerTag=kUnspecified) const
Retrieve an existing batch.
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.
void print(std::ostream &os, const std::string &indent) const
Print to given output stream.
void markDirty()
Mark all batches dirty. This will trigger recomputations.
Status_t status(std::size_t begin, const RooArgSet *const normSet=nullptr, Tag_t ownerTag=kUnspecified) const
Return the status of the batch starting at begin.
void attachForeignStorage(const std::vector< double > &vec)
Attach a foreign storage. Batches coming from this storage will be read only.
RooAbsArg is the common abstract base class for objects that represent a value (of arbitrary type) an...
virtual Bool_t checkObservables(const RooArgSet *nset) const
Overloadable function in which derived classes can implement consistency checks of the variables.
virtual RooAbsArg * cloneTree(const char *newname=0) const
Clone tree expression of objects.
RooArgSet * getObservables(const RooArgSet &set, Bool_t valueOnly=kTRUE) const
Return the observables of this pdf given a set of observables.
void setShapeDirty()
Notify that a shape-like property (e.g. binning) has changed.
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'.
void leafNodeServerList(RooAbsCollection *list, const RooAbsArg *arg=0, Bool_t recurseNonDerived=kFALSE) const
Fill supplied list with all leaf nodes of the arg tree, starting with ourself as top node.
virtual void printMultiline(std::ostream &os, Int_t contents, Bool_t verbose=kFALSE, TString indent="") const
Implement multi-line detailed printing.
virtual TObject * Clone(const char *newname=0) const
Make a clone of an object using the Streamer facility.
static void setDirtyInhibit(Bool_t flag)
Control global dirty inhibit mode.
Bool_t inhibitDirty() const
Delete watch flag.
const Text_t * getStringAttribute(const Text_t *key) const
Get string attribute mapped under key 'key'.
void treeNodeServerList(RooAbsCollection *list, const RooAbsArg *arg=0, Bool_t doBranch=kTRUE, Bool_t doLeaf=kTRUE, Bool_t valueOnly=kFALSE, Bool_t recurseNonDerived=kFALSE) const
Fill supplied list with nodes of the arg tree, following all server links, starting with ourself as t...
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.
RooArgSet * getVariables(Bool_t stripDisconnected=kTRUE) const
Return RooArgSet with all variables (tree leaf nodes of expresssion tree)
static Bool_t _inhibitDirty
void setOperMode(OperMode mode, Bool_t recurseADirty=kTRUE)
Change cache operation mode to given mode.
OperMode _operMode
Mark batches as dirty (only meaningful for RooAbsReal).
void setValueDirty()
Mark the element dirty. This forces a re-evaluation when a value is requested.
Bool_t addOwnedComponents(const RooArgSet &comps)
Take ownership of the contents of 'comps'.
TString cleanBranchName() const
Construct a mangled name from the actual name that is free of any math symbols that might be interpre...
Int_t numProxies() const
Return the number of registered proxies.
friend class RooCustomizer
Bool_t dependsOnValue(const RooAbsCollection &serverList, const RooAbsArg *ignoreArg=0) const
Check whether this object depends on values from an element in the serverList.
void setProxyNormSet(const RooArgSet *nset)
Forward a change in the cached normalization argset to all the registered proxies.
RooAbsProxy * getProxy(Int_t index) const
Return the nth proxy from the proxy list.
RooAbsArg & operator=(const RooAbsArg &other)
Assign all boolean and string properties of the original object.
Bool_t isValueDirtyAndClear() const
RefCountList_t _serverList
RooArgSet * getComponents() const
void branchNodeServerList(RooAbsCollection *list, const RooAbsArg *arg=0, Bool_t recurseNonDerived=kFALSE) const
Fill supplied list with all branch nodes of the arg tree starting with ourself as top node.
Bool_t recursiveRedirectServers(const RooAbsCollection &newServerList, Bool_t mustReplaceAll=kFALSE, Bool_t nameChange=kFALSE, Bool_t recurseInNewSet=kTRUE)
void attachDataSet(const RooAbsData &set)
Replace server nodes with names matching the dataset variable names with those data set variables,...
void SetName(const char *name)
Set the name of the TNamed.
OperMode operMode() const
RooAbsBinning is the abstract base class for RooRealVar binning definitions This class defines the in...
virtual Bool_t isParameterized() const
virtual RooAbsReal * highBoundFunc() const
virtual RooAbsReal * lowBoundFunc() const
void setInterpolationOrder(Int_t order)
Set interpolation order of RooHistFunct representing cache histogram.
RooAbsCategoryLValue is the common abstract base class for objects that represent a discrete value th...
virtual bool setIndex(value_type index, bool printError=true)=0
Change category state by specifying the index code of the desired state.
virtual value_type getCurrentIndex() const
Return index number of current state.
Bool_t isSignType(Bool_t mustHaveZero=kFALSE) const
Determine if category has 2 or 3 states with index values -1,0,1.
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...
virtual Bool_t replace(const RooAbsArg &var1, const RooAbsArg &var2)
Replace var1 with var2 and return kTRUE for success.
void sort(Bool_t reverse=false)
Sort collection using std::sort and name comparison.
virtual TObject * Clone(const char *newname=0) const
Make a clone of an object using the Streamer facility.
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
RooAbsCollection * selectByName(const char *nameList, Bool_t verbose=kFALSE) const
Create a subset of the current collection, consisting only of those elements with names matching the ...
Bool_t overlaps(const RooAbsCollection &otherColl) const
Check if this and other collection have common entries.
std::string contentsString() const
Return comma separated list of contained object names as STL string.
virtual Bool_t remove(const RooAbsArg &var, Bool_t silent=kFALSE, Bool_t matchByNameOnly=kFALSE)
Remove the specified argument from our list.
TIterator * createIterator(Bool_t dir=kIterForward) const
TIterator-style iteration over contained elements.
RooAbsArg * find(const char *name) const
Find object with given name in list.
RooAbsData is the common abstract base class for binned and unbinned datasets.
virtual const RooArgSet * get() const
RooAbsData * reduce(const RooCmdArg &arg1, const RooCmdArg &arg2=RooCmdArg(), const RooCmdArg &arg3=RooCmdArg(), const RooCmdArg &arg4=RooCmdArg(), const RooCmdArg &arg5=RooCmdArg(), const RooCmdArg &arg6=RooCmdArg(), const RooCmdArg &arg7=RooCmdArg(), const RooCmdArg &arg8=RooCmdArg())
Create a reduced copy of this dataset.
Abstract interface for evaluating a real-valued function of one real variable and performing numerica...
Abstract base class for objects that are lvalues, i.e.
RooAbsMoment represents the first, second, or third order derivative of any RooAbsReal as calculated ...
void constOptimizeTestStatistic(ConstOpCode opcode, Bool_t doAlsoTrackingOpt=kTRUE)
Driver function to propagate constant term optimizations in test statistic.
const char * normRange() const
virtual Double_t expectedEvents(const RooArgSet *nset) const
Return expected number of events from this p.d.f for use in extended likelihood calculations.
virtual ExtendMode extendMode() const
RooAbsProxy is the abstact interface for proxy classes.
virtual void print(std::ostream &os, Bool_t addContents=kFALSE) const
Print proxy name.
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 const RooAbsBinning & getBinning(const char *name=0, Bool_t verbose=kTRUE, Bool_t createOnTheFly=kFALSE) const =0
Retrive binning configuration with given name or default binning.
TH1 * createHistogram(const char *name, const RooCmdArg &arg1=RooCmdArg::none(), 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
virtual Double_t getMin(const char *name=0) const
Get miniminum of currently defined range.
void setServerValues(const char *tmp)
void setMessage(const char *tmp)
RooAbsReal is the common abstract base class for objects that represent a real value and implements f...
void plotOnCompSelect(RooArgSet *selNodes) const
Helper function for plotting of composite p.d.fs.
std::unique_ptr< TreeReadBuffer > _treeReadBuffer
RooGenFunction * iGenFunction(RooRealVar &x, const RooArgSet &nset=RooArgSet())
virtual RooPlot * plotSliceOn(RooPlot *frame, const RooArgSet &sliceSet, Option_t *drawOptions="L", Double_t scaleFactor=1.0, ScaleType stype=Relative, const RooAbsData *projData=0) const
virtual Bool_t readFromStream(std::istream &is, Bool_t compact, Bool_t verbose=kFALSE)
Read object contents from stream (dummy for now)
TH1 * fillHistogram(TH1 *hist, const RooArgList &plotVars, Double_t scaleFactor=1, const RooArgSet *projectedVars=0, Bool_t scaling=kTRUE, const RooArgSet *condObs=0, Bool_t setError=kTRUE) const
Fill the ROOT histogram 'hist' with values sampled from this function at the bin centers.
virtual void selectNormalizationRange(const char *rangeName=0, Bool_t force=kFALSE)
Interface function to force use of a given normalization range to interpret function value.
TString getTitle(Bool_t appendUnit=kFALSE) const
Return this variable's title string.
static Int_t numEvalErrorItems()
static void setHideOffset(Bool_t flag)
RooAbsReal()
coverity[UNINIT_CTOR] Default constructor
virtual Int_t getAnalyticalIntegralWN(RooArgSet &allVars, RooArgSet &analVars, const RooArgSet *normSet, const char *rangeName=0) const
Variant of getAnalyticalIntegral that is also passed the normalization set that should be applied to ...
friend class RooRealBinding
RooAbsReal * createScanRI(const RooArgSet &iset, const RooArgSet &nset, Int_t numScanBins, Int_t intOrder)
Utility function for createRunningIntegral that construct an object implementing the numeric scanning...
virtual Bool_t isIdentical(const RooAbsArg &other, Bool_t assumeSameType=kFALSE) const
Bool_t isSelectedComp() const
If true, the current pdf is a selected component (for use in plotting)
virtual RooAbsReal * createChi2(RooDataHist &data, const RooLinkedList &cmdList)
Create a variable from a histogram and this function.
RooDerivative * derivative(RooRealVar &obs, Int_t order=1, Double_t eps=0.001)
Return function representing first, second or third order derivative of this function.
virtual Bool_t isValidReal(Double_t value, Bool_t printError=kFALSE) const
Interface function to check if given value is a valid value for this object.
virtual RooPlot * plotOn(RooPlot *frame, const RooCmdArg &arg1=RooCmdArg(), const RooCmdArg &arg2=RooCmdArg(), const RooCmdArg &arg3=RooCmdArg(), const RooCmdArg &arg4=RooCmdArg(), const RooCmdArg &arg5=RooCmdArg(), const RooCmdArg &arg6=RooCmdArg(), const RooCmdArg &arg7=RooCmdArg(), const RooCmdArg &arg8=RooCmdArg(), const RooCmdArg &arg9=RooCmdArg(), const RooCmdArg &arg10=RooCmdArg()) const
Plot (project) PDF on specified frame.
void setParameterizeIntegral(const RooArgSet ¶mVars)
static Bool_t hideOffset()
virtual Double_t evaluate() const =0
Evaluate this PDF / function / constant. Needs to be overridden by all derived classes.
virtual void setTreeBranchStatus(TTree &t, Bool_t active)
(De)Activate associated tree branch
virtual RooAbsReal * createProfile(const RooArgSet ¶msOfInterest)
Create a RooProfileLL object that eliminates all nuisance parameters in the present function.
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< double > evaluateBatch(std::size_t begin, std::size_t maxSize) const
Evaluate function for a batch of input data points.
Double_t _DEBUG_getVal(const RooArgSet *normalisationSet) const
Debug version of getVal(), which is slow and does error checking.
Double_t traceEval(const RooArgSet *set) const
Calculate current value of object, with error tracing wrapper.
virtual void copyCache(const RooAbsArg *source, Bool_t valueOnly=kFALSE, Bool_t setValDirty=kTRUE)
Copy the cached value of another RooAbsArg to our cache.
virtual void fixAddCoefRange(const char *rangeName=0, Bool_t force=kTRUE)
Fix the interpretation of the coefficient of any RooAddPdf component in the expression tree headed by...
TH1 * createHistogram(const char *varNameList, Int_t xbins=0, Int_t ybins=0, Int_t zbins=0) const
Create and fill a ROOT histogram TH1, TH2 or TH3 with the values of this function for the variables w...
static EvalErrorIter evalErrorIter()
RooAbsReal & operator=(const RooAbsReal &other)
Assign values, name and configs from another RooAbsReal.
Double_t findRoot(RooRealVar &x, Double_t xmin, Double_t xmax, Double_t yval)
Return value of x (in range xmin,xmax) at which function equals yval.
RooAbsReal * createIntRI(const RooArgSet &iset, const RooArgSet &nset=RooArgSet())
Utility function for createRunningIntegral.
TF1 * asTF(const RooArgList &obs, const RooArgList &pars=RooArgList(), const RooArgSet &nset=RooArgSet()) const
Return a ROOT TF1,2,3 object bound to this RooAbsReal with given definition of observables and parame...
virtual Double_t analyticalIntegralWN(Int_t code, const RooArgSet *normSet, const char *rangeName=0) const
Implements the actual analytical integral(s) advertised by getAnalyticalIntegral.
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 void attachToVStore(RooVectorDataStore &vstore)
TString _unit
Value storage for batches of events.
static RooNumIntConfig * defaultIntegratorConfig()
Returns the default numeric integration configuration for all RooAbsReals.
RooFitResult * chi2FitDriver(RooAbsReal &fcn, RooLinkedList &cmdList)
Internal driver function for chi2 fits.
virtual std::list< Double_t > * binBoundaries(RooAbsRealLValue &, Double_t, Double_t) const
virtual void attachToTree(TTree &t, Int_t bufSize=32000)
Attach object to a branch of given TTree.
Bool_t matchArgs(const RooArgSet &allDeps, RooArgSet &numDeps, const RooArgProxy &a) const
Utility function for use in getAnalyticalIntegral().
RooNumIntConfig * _specIntegratorConfig
virtual RooPlot * plotAsymOn(RooPlot *frame, const RooAbsCategoryLValue &asymCat, PlotOpt o) const
static ErrorLoggingMode evalErrorLoggingMode()
Return current evaluation error logging mode.
void setIntegratorConfig()
Remove the specialized numeric integration configuration associated with this object.
virtual void selectNormalization(const RooArgSet *depSet=0, Bool_t force=kFALSE)
Interface function to force use of a given set of observables to interpret function value.
virtual void printMultiline(std::ostream &os, Int_t contents, Bool_t verbose=kFALSE, TString indent="") const
Structure printing.
friend class RooRealIntegral
Double_t getVal(const RooArgSet *normalisationSet=nullptr) const
Evaluate object.
RooAbsMoment * moment(RooRealVar &obs, Int_t order, Bool_t central, Bool_t takeRoot)
Return function representing moment of function of given order.
static ErrorLoggingMode _evalErrorMode
Double_t getPropagatedError(const RooFitResult &fr, const RooArgSet &nset=RooArgSet()) const
Calculate error on self by propagated errors on parameters with correlations as given by fit result T...
static Int_t _evalErrorCount
virtual Bool_t isValid() const
Check if current value is valid.
static void globalSelectComp(Bool_t flag)
Global switch controlling the activation of the selectComp() functionality.
RooAbsArg * createFundamental(const char *newname=0) const
Create a RooRealVar fundamental object with our properties.
std::map< constRooAbsArg *, std::pair< std::string, std::list< EvalError > > >::const_iterator EvalErrorIter
virtual Double_t maxVal(Int_t code) const
Return maximum value for set of observables identified by code assigned in getMaxVal.
static Bool_t _globalSelectComp
Component selection flag for RooAbsPdf::plotCompOn.
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 ...
RooNumIntConfig * specialIntegratorConfig() const
Returns the specialized integrator configuration for this RooAbsReal.
RooAbsReal * createIntObj(const RooArgSet &iset, const RooArgSet *nset, const RooNumIntConfig *cfg, const char *rangeName) const
Internal utility function for createIntegral() that creates the actual integral object.
BatchHelpers::BatchData _batchData
static std::map< const RooAbsArg *, std::pair< std::string, std::list< EvalError > > > _evalErrorList
Bool_t plotSanityChecks(RooPlot *frame) const
Utility function for plotOn(), perform general sanity check on frame to ensure safe plotting operatio...
const char * getPlotLabel() const
Get the label associated with the variable.
virtual RooFitResult * chi2FitTo(RooDataHist &data, const RooCmdArg &arg1=RooCmdArg::none(), 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())
Perform a fit to given histogram.
virtual void fillTreeBranch(TTree &t)
Fill the tree branch that associated with this object with its current value.
Bool_t operator==(Double_t value) const
Equality operator comparing to a Double_t.
RooAbsFunc * bindVars(const RooArgSet &vars, const RooArgSet *nset=0, Bool_t clipInvalid=kFALSE) const
Create an interface adaptor f(vars) that binds us to the specified variables (in arbitrary order).
static Int_t numEvalErrors()
Return the number of logged evaluation errors since the last clearing.
static void setEvalErrorLoggingMode(ErrorLoggingMode m)
Set evaluation error logging mode.
Bool_t matchArgsByName(const RooArgSet &allArgs, RooArgSet &matchedArgs, const TList &nameList) const
Check if allArgs contains matching elements for each name in nameList.
virtual Double_t analyticalIntegral(Int_t code, const char *rangeName=0) const
Implements the actual analytical integral(s) advertised by getAnalyticalIntegral.
virtual Int_t getAnalyticalIntegral(RooArgSet &allVars, RooArgSet &analVars, const char *rangeName=0) const
Interface function getAnalyticalIntergral advertises the analytical integrals that are supported.
virtual void preferredObservableScanOrder(const RooArgSet &obs, RooArgSet &orderedObs) const
Interface method for function objects to indicate their preferred order of observables for scanning t...
void findInnerMostIntegration(const RooArgSet &allObs, RooArgSet &innerObs, const char *rangeName) const
Utility function for createIntObj() that aids in the construct of recursive integrals over functions ...
void selectComp(Bool_t flag)
static Bool_t _hideOffset
const Text_t * getUnit() const
virtual void printValue(std::ostream &os) const
Print object value.
const RooNumIntConfig * getIntegratorConfig() const
Return the numeric integration configuration used for this object.
RooPlot * plotOnWithErrorBand(RooPlot *frame, const RooFitResult &fr, Double_t Z, const RooArgSet *params, const RooLinkedList &argList, Bool_t method1) const
Plot function or PDF on frame with support for visualization of the uncertainty encoded in the given ...
virtual Double_t getValV(const RooArgSet *normalisationSet=nullptr) const
Return value of object.
static void printEvalErrors(std::ostream &os=std::cout, Int_t maxPerNode=10000000)
Print all outstanding logged evaluation error on the given ostream.
virtual void fixAddCoefNormalization(const RooArgSet &addNormSet=RooArgSet(), Bool_t force=kTRUE)
Fix the interpretation of the coefficient of any RooAddPdf component in the expression tree headed by...
void logEvalError(const char *message, const char *serverValueString=0) const
Log evaluation error message.
static void clearEvalErrorLog()
Clear the stack of evaluation error messages.
void checkBatchComputation(std::size_t evtNo, const RooArgSet *normSet=nullptr, double relAccuracy=1.E-13) const
RooFunctor * functor(const RooArgList &obs, const RooArgList &pars=RooArgList(), const RooArgSet &nset=RooArgSet()) const
Return a RooFunctor object bound to this RooAbsReal with given definition of observables and paramete...
void makeProjectionSet(const RooAbsArg *plotVar, const RooArgSet *allVars, RooArgSet &projectedVars, Bool_t silent) const
Utility function for plotOn() that constructs the set of observables to project when plotting ourselv...
void setPlotLabel(const char *label)
Set the label associated with this variable.
virtual ~RooAbsReal()
Destructor.
virtual Int_t getMaxVal(const RooArgSet &vars) const
Advertise capability to determine maximum value of function for given set of observables.
virtual void writeToStream(std::ostream &os, Bool_t compact) const
Write object contents to stream (dummy for now)
const RooAbsReal * createPlotProjection(const RooArgSet &depVars, const RooArgSet &projVars) const
Utility function for plotOn() that creates a projection of a function or p.d.f to be plotted on a Roo...
RooDataHist * fillDataHist(RooDataHist *hist, const RooArgSet *nset, Double_t scaleFactor, Bool_t correctForBinVolume=kFALSE, Bool_t showProgress=kFALSE) const
Fill a RooDataHist with values sampled from this function at the bin centers.
virtual Double_t offset() const
RooAbsReal * createRunningIntegral(const RooArgSet &iset, const RooArgSet &nset=RooArgSet())
Calls createRunningIntegral(const RooArgSet&, const RooCmdArg&, const RooCmdArg&, const RooCmdArg&,...
RooAddition calculates the sum of a set of RooAbsReal terms, or when constructed with two sets,...
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.
RooArgProxy is the abstract interface for RooAbsArg proxy classes.
RooArgSet is a container object that can hold multiple RooAbsArg objects.
virtual Bool_t addOwned(const RooAbsCollection &col, Bool_t silent=kFALSE)
Add a collection of arguments to this collection by calling addOwned() for each element in the source...
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...
virtual void addClone(const RooAbsCollection &col, Bool_t silent=kFALSE)
Add a collection of arguments to this collection by calling addOwned() for each element in the source...
Class RooBinning is an implements RooAbsBinning in terms of an array of boundary values,...
Implement the abstract 1-dimensional root finding interface using the Brent-Decker method.
virtual Bool_t findRoot(Double_t &result, Double_t xlo, Double_t xhi, Double_t value=0) const
Do the root finding using the Brent-Decker method.
RooCachedReal is an implementation of RooAbsCachedReal that can cache any external RooAbsReal input f...
void setCacheSource(Bool_t flag)
RooCategory is an object to represent discrete states.
virtual Bool_t setLabel(const char *label, bool printError=true) override
Set value by specifying the name of the desired state.
Class RooChi2Var implements a simple chi^2 calculation from a binned dataset and a PDF.
RooCmdArg is a named container for two doubles, two integers two object points and three string point...
static const RooCmdArg & none()
Return reference to null argument.
const char * getString(Int_t idx) const
void setString(Int_t idx, const char *value)
virtual TObject * Clone(const char *newName=0) const
Make a clone of an object using the Streamer facility.
Class RooCmdConfig is a configurable parser for RooCmdArg named arguments.
A RooCurve is a one-dimensional graphical representation of a real-valued function.
RooCurve * makeErrorBand(const std::vector< RooCurve * > &variations, Double_t Z=1) const
Construct filled RooCurve represented error band that captures alpha% of the variations of the curves...
RooCustomizer is a factory class to produce clones of a prototype composite PDF object with the same ...
void setCloneBranchSet(RooArgSet &cloneBranchSet)
Install the input RooArgSet as container in which all cloned branches will be stored.
void replaceArg(const RooAbsArg &orig, const RooAbsArg &subst)
Replace any occurence of arg 'orig' with arg 'subst'.
void setOwning(Bool_t flag)
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...
The RooDataHist is a container class to hold N-dimensional binned data.
void set(Double_t weight, Double_t wgtErr=-1)
Set the weight and weight error of the bin enclosing the current (i.e.
virtual Int_t numEntries() const
Return the number of bins.
virtual const RooArgSet * get() const
Double_t binVolume() const
RooDataSet is a container class to hold unbinned data.
Class RooDataWeightedAverage calculate a weighted average of a function or p.d.f given a dataset with...
RooDerivative represents the first, second, or third order derivative of any RooAbsReal as calculated...
A one-time forward iterator working on RooLinkedList or RooAbsCollection.
RooAbsArg * next()
Return next element or nullptr if at end.
RooFirstMoment represents the first, second, or third order derivative of any RooAbsReal as calculate...
RooFitResult is a container class to hold the input and output of a PDF fit to a dataset.
const TMatrixDSym & covarianceMatrix() const
Return covariance matrix.
TMatrixDSym reducedCovarianceMatrix(const RooArgList ¶ms) const
Return a reduced covariance matrix (Note that Vred is a simple sub-matrix of V, row/columns are order...
const RooArgList & floatParsFinal() const
RooAbsPdf * createHessePdf(const RooArgSet ¶ms) const
Return a p.d.f that represents the fit result as a multi-variate probability densisty function on the...
Lightweight interface adaptor that exports a RooAbsPdf as a functor.
Lightweight interface adaptor that exports a RooAbsReal as a ROOT::Math::IGenFunction.
A RooHist is a graphical representation of binned data based on the TGraphAsymmErrors class.
RooLinkedList is an collection class for internal use, storing a collection of RooAbsArg pointers in ...
TObject * FindObject(const char *name) const
Return pointer to obejct with given name.
RooFIter fwdIterator() const
Create a one-time-use forward iterator for this list.
TIterator * MakeIterator(Bool_t forward=kTRUE) const
Create a TIterator for this list.
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)
RooMinimizer is a wrapper class around ROOT::Fit:Fitter that provides a seamless interface between th...
RooMinuit is a wrapper class around TFitter/TMinuit that provides a seamless interface between the MI...
RooMoment represents the first, second, or third order derivative of any RooAbsReal as calculated (nu...
static RooMsgService & instance()
Return reference to singleton instance.
Lightweight interface adaptor that exports a RooAbsReal as a ROOT::Math::IMultiGenFunction.
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...
void setNameList(const char *givenList)
RooNumIntConfig holds the configuration parameters of the various numeric integrators used by RooReal...
static RooNumIntConfig & defaultConfig()
Return reference to instance of default numeric integrator configuration object.
Class RooNumRunningInt is an implementation of RooAbsCachedReal that represents a running integral.
Class RooParamBinning is an implementation of RooAbsBinning that constructs a binning with a range de...
A RooPlot is a plot frame and a container for graphics objects within that frame.
Bool_t drawBefore(const char *before, const char *target)
Change the order in which our contained objects are drawn so that the target object is drawn just bef...
TAttMarker * getAttMarker(const char *name=0) const
Return a pointer to the marker attributes of the named object in this plot, or zero if the named obje...
Double_t GetMinimum(Double_t minval=-FLT_MAX) const
const RooArgSet * getNormVars() const
TAttFill * getAttFill(const char *name=0) const
Return a pointer to the fill attributes of the named object in this plot, or zero if the named object...
virtual void SetMinimum(Double_t minimum=-1111)
Set minimum value of Y axis.
TObject * getObject(Int_t idx) const
Return the name of the object at slot 'idx' in this RooPlot.
RooAbsRealLValue * getPlotVar() const
TAttLine * getAttLine(const char *name=0) const
Return a pointer to the line attributes of the named object in this plot, or zero if the named object...
Double_t GetMaximum(Double_t maxval=FLT_MAX) const
void addPlotable(RooPlotable *plotable, Option_t *drawOptions="", Bool_t invisible=kFALSE, Bool_t refreshNorm=kFALSE)
Add the specified plotable object to our plot.
virtual void SetMaximum(Double_t maximum=-1111)
Set maximum value of Y axis.
TObject * findObject(const char *name, const TClass *clas=0) const
Find the named object in our list of items and return a pointer to it.
RooCurve * getCurve(const char *name=0) const
Return a RooCurve pointer of the named object in this plot, or zero if the named object does not exis...
void remove(const char *name=0, Bool_t deleteToo=kTRUE)
Remove object with given name, or last object added if no name is given.
virtual void printStream(std::ostream &os, Int_t contents, StyleOption style, TString indent="") const
Print description of object on ostream, printing contents set by contents integer,...
Class RooProfileLL implements the profile likelihood estimator for a given likelihood and set of para...
Lightweight interface adaptor that binds a RooAbsReal object to a subset of its servers and present i...
RooRealIntegral performs hybrid numerical/analytical integrals of RooAbsReal objects.
const RooArgSet & numIntRealVars() const
RooRealVar represents a variable that can be changed from the outside.
RooSpan< const double > getValBatch(std::size_t begin, std::size_t batchSize, const RooArgSet *=nullptr) const final
Return batch of data starting at begin.
void setBins(Int_t nBins, const char *name=0)
Create a uniform binning under name 'name' for this variable.
void setRange(const char *name, Double_t min, Double_t max)
Set a fit or plotting range.
void removeRange(const char *name=0)
Remove range limits for binning with given name. Empty name means default range.
Double_t getError() const
void setBinning(const RooAbsBinning &binning, const char *name=0)
Add given binning under name 'name' with this variable.
virtual TObject * clone(const char *newname) const
virtual void setVal(Double_t value)
Set value of variable to 'value'.
Lightweight RooAbsFunction implementation that applies a constant scale factor to another RooAbsFunc.
RooSecondMoment represents the first, second, or third order derivative of any RooAbsReal as calculat...
A simple container to hold a batch of data values.
constexpr bool empty() const noexcept
void setBuffer(RooAbsReal *real, Double_t *newBuf)
const std::vector< double > & data() const
RooVectorDataStore is the abstract base class for data collection that use a TTree as internal storag...
RealVector * addReal(RooAbsReal *real)
virtual void SetFillColor(Color_t fcolor)
Set the fill area color.
virtual void SetFillStyle(Style_t fstyle)
Set the fill area style.
virtual void SetLineStyle(Style_t lstyle)
Set the line style.
virtual void SetLineWidth(Width_t lwidth)
Set the line width.
virtual void SetLineColor(Color_t lcolor)
Set the line color.
virtual void SetMarkerColor(Color_t mcolor=1)
Set the marker color.
virtual void SetMarkerStyle(Style_t mstyle=1)
Set the marker style.
virtual void SetMarkerSize(Size_t msize=1)
Set the marker size.
Class to manage histogram axis.
virtual Double_t GetBinCenter(Int_t bin) const
Return center of bin.
A TTree is a list of TBranches.
TObjArray * GetListOfLeaves()
virtual Int_t GetSize() const
Return the capacity of the collection, i.e.
virtual void SetParName(Int_t ipar, const char *name)
Set name of parameter number ipar.
virtual void SetParameter(Int_t param, Double_t value)
A 2-Dim function with parameters.
A 3-Dim function with parameters.
virtual void SetName(const char *name="")
Set graph name.
virtual Int_t GetNbinsY() const
virtual Int_t GetNbinsZ() const
virtual Int_t GetDimension() const
TAxis * GetXaxis()
Get the behaviour adopted by the object about the statoverflows. See EStatOverflows for more informat...
virtual Int_t GetBin(Int_t binx, Int_t biny=0, Int_t binz=0) const
Return Global bin number corresponding to binx,y,z.
virtual Int_t GetNbinsX() const
virtual void SetBinError(Int_t bin, Double_t error)
Set the bin Error Note that this resets the bin eror option to be of Normal Type and for the non-empt...
virtual void SetBinContent(Int_t bin, Double_t content)
Set bin content see convention for numbering bins in TH1::GetBin In case the bin number is greater th...
Iterator abstract base class.
virtual TObject * Next()=0
A TLeaf describes individual elements of a TBranch See TBranch structure in TTree.
virtual const char * GetTypeName() const
virtual TLeaf * GetLeafCounter(Int_t &countval) const
Return a pointer to the counter of this leaf (if any) or store the number of elements that the leaf c...
virtual void Add(TObject *obj)
virtual TIterator * MakeIterator(Bool_t dir=kIterForward) const
Return a list iterator.
virtual void Delete(Option_t *option="")
Remove all objects from the list AND delete all heap based objects.
virtual void SetTitle(const char *title="")
Set the title of the TNamed.
virtual const char * GetTitle() const
Returns title of object.
virtual const char * GetName() const
Returns name of object.
TObject * At(Int_t idx) const
Collectable string class.
Mother of all ROOT objects.
virtual const char * GetName() const
Returns name of object.
virtual const char * ClassName() const
Returns name of class to which the object belongs.
virtual Bool_t InheritsFrom(const char *classname) const
Returns kTRUE if object inherits from class "classname".
int CompareTo(const char *cs, ECaseCompare cmp=kExact) const
Compare a string to char *cs2.
const char * Data() const
TString & Prepend(const char *cs)
TString & Append(const char *cs)
Bool_t Contains(const char *pat, ECaseCompare cmp=kExact) const
virtual void StackTrace()
Print a stack trace.
A TTree represents a columnar dataset.
virtual TBranch * GetBranch(const char *name)
Return pointer to the branch with the given name in this tree or its friends.
virtual Int_t SetBranchAddress(const char *bname, void *add, TBranch **ptr=0)
Change branch address, dealing with clone trees properly.
TBranch * Branch(const char *name, T *obj, Int_t bufsize=32000, Int_t splitlevel=99)
Add a new branch, and infer the data type from the type of obj being passed.
virtual void SetBranchStatus(const char *bname, Bool_t status=1, UInt_t *found=0)
Set branch status to Process or DoNotProcess.
RooCmdArg YVar(const RooAbsRealLValue &var, const RooCmdArg &arg=RooCmdArg::none())
RooCmdArg ZVar(const RooAbsRealLValue &var, const RooCmdArg &arg=RooCmdArg::none())
RooCmdArg SupNormSet(const RooArgSet &nset)
RooCmdArg Binning(const RooAbsBinning &binning)
RooCmdArg NormRange(const char *rangeNameList)
VecExpr< UnaryOp< Fabs< T >, VecExpr< A, T, D >, T >, T, D > fabs(const VecExpr< A, T, D > &rhs)
std::vector< std::string > tokenise(const std::string &str, const std::string &delims, bool returnEmptyToken=true)
Tokenise the string by splitting at the characters in delims.
static constexpr double pc
static constexpr double second
static constexpr double ps
Double_t Erfc(Double_t x)
Compute the complementary error function erfc(x).
Bool_t postRangeFracScale
const char * normRangeName
RooFit::MPSplit interleave
const char * projectionRangeName
const RooArgSet * projDataSet
const char * curveNameSuffix
const char * addToCurveName
const RooFitResult * errorFR
const RooArgSet * projSet
const RooAbsData * projData
static long int sum(long int i)
static void output(int code)