42xRooNLLVar::xRooHypoSpace::xRooHypoSpace(
const char *
name,
const char *title)
45 if (
name ==
nullptr || strlen(
name) == 0) {
64 spaceSize *= (
v->getMax() -
v->getMin());
69 double ssize = spaceSize;
74 ssize /= (
v->getMax() -
v->getMin());
75 double remain = std::fmod(xVal, ssize);
76 v->setVal((xVal - remain) / ssize);
82 for (
auto &
p : *
this) {
83 for (
auto a : *
p.coords) {
84 if (!
fPars->find(
a->GetName()))
92 return pdf(toArgs(parValues));
101 std::shared_ptr<xRooNode> out =
nullptr;
103 for (
auto &[_range, _pdf] : fPdfs) {
106 bool collision =
true;
107 for (
auto &_lhs : *_range) {
108 auto _rhs = rhs.
find(*_lhs);
113 if (!(
v->getMin() <=
v2->getMax() &&
v2->getMin() <=
v->getMax())) {
118 if (!(
v->getMin() <=
c2->getVal() &&
c2->getVal() <=
v->getMax())) {
125 if (!(
c->getVal() <=
v2->getMax() &&
v2->getMin() <=
c->getVal())) {
130 if (!(
c->getVal() ==
c2->getVal())) {
139 throw std::runtime_error(
"Multiple pdf possibilities");
157 auto _idx = s.
Index(
'=');
166 _idx = _val.
Index(
',');
185 throw std::runtime_error(
"Unknown parameter");
186 _par->setAttribute(
"axis");
189 _par->setVal((high + low) * 0.5);
194 double step = (high - low) / (nPoints - 1);
196 throw std::runtime_error(
"Invalid steps");
198 for (
size_t i = 0; i < nPoints; i++) {
199 _par->setVal((i == nPoints - 1) ? high : (low + step * i));
207 if (axes().empty()) {
210 throw std::runtime_error(
"No POI to scan");
212 poi().first()->setAttribute(
"axis");
218 poi().setAttribAll(
"poi",
false);
219 axes().setAttribAll(
"poi");
226 const std::vector<double> &nSigmas,
double relUncert)
234 throw std::runtime_error(
"scan type must be equal to one of: plr, cls, ts, pnull");
238 if (axes().empty()) {
241 throw std::runtime_error(
"No POI to scan");
243 poi().first()->setAttribute(
"axis");
248 poi().setAttribAll(
"poi",
false);
249 axes().setAttribAll(
"poi");
251 auto p =
dynamic_cast<RooRealVar *
>(axes().first());
253 throw std::runtime_error(
TString::Format(
"%s not scannable", axes().first()->GetName()));
257 if (empty() && relUncert == std::numeric_limits<double>::infinity()) {
259 ::Info(
"xRooHypoSpace::scan",
"Using default precision of 10%% for auto-scan");
262 for (
auto a : axes()) {
263 if (!
a->hasRange(
"physical")) {
264 ::Info(
"xRooHypoSpace::scan",
"No physical range set for %s, setting to [0,inf]",
p->
GetName());
265 dynamic_cast<RooRealVar *
>(
a)->setRange(
"physical", 0, std::numeric_limits<double>::infinity());
267 if (!
a->getStringAttribute(
"altVal") || !strlen(
p->getStringAttribute(
"altVal"))) {
268 ::Info(
"xRooHypoSpace::scan",
"No altVal set for %s, setting to 0",
a->GetName());
269 a->setStringAttribute(
"altVal",
"0");
272 double altVal =
TString(
a->getStringAttribute(
"altVal")).
Atof();
274 if (
v->getMin() >= altVal) {
275 ::Info(
"xRooHypoSpace::scan",
"range of POI does not straddle alt value, adjusting minimum to %g",
277 v->setMin(altVal - 1
e-5);
279 if (
v->getMax() <= altVal) {
280 ::Info(
"xRooHypoSpace::scan",
"range of POI does not straddle alt value, adjusting maximum to %g",
282 v->setMax(altVal + 1
e-5);
284 for (
auto &[
pdf, nll] : fNlls) {
285 if (
auto _v =
dynamic_cast<RooRealVar *
>(nll->pars()->find(*
a))) {
286 _v->setRange(
v->getMin(),
v->getMax());
296 if (
p && high <= low) {
298 low =
p->getMin(
"scan");
299 high =
p->getMax(
"scan");
300 ::Info(
"xRooHypoSpace::scan",
"Using %s range: %g - %g",
p->
GetName(), low, high);
302 if (!std::isnan(low) && !std::isnan(high) && !(std::isinf(low) && std::isinf(high))) {
303 p->setRange(
"scan", low, high);
307 for (
auto nSigma : nSigmas) {
308 if (std::isnan(nSigma)) {
317 relUncert = std::numeric_limits<double>::infinity();
326#if ROOT_VERSION_CODE >= ROOT_VERSION(6, 26, 00)
346 for (
auto file : *
gROOT->GetListOfFiles()) {
347 if (
auto f =
dynamic_cast<TMemFile *
>(file)) {
353 new TMemFile(
"fitDatabase",
"RECREATE");
362 for (
double nSigma : nSigmas) {
364 if (std::isnan(nSigma)) {
372 if (std::isnan(res.first) || std::isnan(res.second)) {
374 }
else if (std::isinf(res.second)) {
379 throw std::runtime_error(
TString::Format(
"Automatic scanning not yet supported for %s",
type));
384 AddPoint(
TString::Format(
"%s=%g", poi().first()->GetName(), (high + low) / 2.));
387 double step = (high - low) / (nPoints - 1);
388 for (
size_t i = 0; i < nPoints; i++) {
389 AddPoint(
TString::Format(
"%s=%g", poi().first()->GetName(), low + step * i));
401std::map<std::string, std::pair<double, double>>
408 relUncert = std::numeric_limits<double>::infinity();
411 scan(opt, nSigmas, relUncert);
413 std::map<std::string, std::pair<double, double>> out;
414 for (
auto nSigma : nSigmas) {
415 auto lim = limit(opt, nSigma);
417 lim.second = -lim.second;
435 auto _idx = s.
Index(
'=');
440 auto _v =
dynamic_cast<RooRealVar *
>(fPars->find(_name));
446 _v->setVal(_val.
Atof());
453 throw std::runtime_error(
"no model at coordinates");
460 if (fNlls.find(_pdf) == fNlls.end()) {
461 fNlls[_pdf] = std::make_shared<xRooNLLVar>(_pdf->nll(
"" , {}));
466 out.nllVar = fNlls[_pdf];
467 out.fData = fNlls[_pdf]->getData();
468 out.isExpected =
dynamic_cast<RooDataSet *
>(out.fData.first.get()) &&
469 dynamic_cast<RooDataSet *
>(out.fData.first.get())->weightVar()->getAttribute(
"expected");
472 out.coords.reset(fPars->snapshot());
474 out.poi().setAttribAll(
"Constant");
481 ->remove(*std::unique_ptr<RooAbsCollection>(out.coords->selectByAttrib(
"Constant",
false)),
true,
true);
483 double value = out.fNullVal();
484 double alt_value = out.fAltVal();
486 auto _type = fTestStatType;
489 if (std::isnan(alt_value)) {
491 }
else if (
value >= alt_value) {
498 out.fPllType = _type;
501 for (
auto &
p : *
this) {
502 if (
p.nllVar != out.nllVar)
504 if (
p.fData != out.fData)
506 if (!
p.alt_poi().equals(out.alt_poi()))
509 for (
auto c :
p.alt_poi()) {
511 v && std::abs(
v->getVal() - out.alt_poi().getRealValue(
v->GetName())) > 1
e-12) {
518 if (!
p.coords->equals(*out.coords))
520 for (
auto c : *
p.coords) {
521 if (
c->getAttribute(
"poi")) {
525 v && std::abs(
v->getVal() - out.coords->getRealValue(
v->GetName())) > 1
e-12) {
534 if (
auto cfit =
p.cfit_alt(
true)) {
535 out.fAlt_cfit = cfit;
537 if (
p.asimov(
true) &&
p.asimov(
true)->fData.first && (!out.asimov(
true) || !out.asimov(
true)->fData.first)) {
538 out.asimov()->fData =
p.asimov(
true)->fData;
540 if (!
p.poi().equals(out.poi()))
542 for (
auto c :
p.poi()) {
544 v && std::abs(
v->getVal() - out.poi().getRealValue(
v->GetName())) > 1
e-12) {
555 std::string coordString;
556 for (
auto a : axes()) {
557 coordString +=
TString::Format(
"%s=%g",
a->GetName(), out.coords->getRealValue(
a->GetName()));
560 coordString.erase(coordString.end() - 1);
562 ::Info(
"xRooHypoSpace::AddPoint",
"Added new point @ %s", coordString.c_str());
563 return emplace_back(out);
570 throw std::runtime_error(
"Not a pdf");
576 auto vpars = toArgs(validity);
578 vpars.remove(
pars,
true,
true);
581 if (
auto existing =
pdf(
pars)) {
582 throw std::runtime_error(std::string(
"Clashing model: ") + existing->GetName());
585 auto myPars = std::shared_ptr<RooArgList>(
dynamic_cast<RooArgList *
>(
pars.snapshot()));
588 pars.remove(*fPars,
true,
true);
590 fPars->addClone(
pars);
592 fPdfs.insert(std::make_pair(myPars, std::make_shared<xRooNode>(_pdf)));
603 out.add(*std::unique_ptr<RooAbsCollection>(
604 fPars->selectByAttrib(
"axis",
true)));
610 std::set<std::vector<double>> coords;
611 for (
auto &
p : *
this) {
612 std::vector<double> p_coords;
614 auto _v =
dynamic_cast<RooRealVar *
>(
p.coords->find(o->GetName()));
616 (_v && _v->isConstant())
618 : std::numeric_limits<double>::infinity());
621 if (coords.find(p_coords) != coords.end()) {
625 coords.insert(p_coords);
630 std::map<std::string, std::unordered_set<double>> values;
631 for (
auto &par : *
pars()) {
634 for (
auto p : *
this) {
635 auto _v =
dynamic_cast<RooRealVar *
>(
p.coords->find(par->GetName()));
636 values[par->GetName()].insert(
637 (_v && _v->isConstant())
639 : std::numeric_limits<double>::infinity());
648 for (
auto &[k,
v] : values) {
649 if (
v.size() > maxDiff || (
v.size() == maxDiff && !isPOI &&
pars()->find(k.c_str())->getAttribute(
"poi"))) {
651 isPOI =
pars()->find(k.c_str())->getAttribute(
"poi");
652 maxDiff = std::max(maxDiff,
v.size());
655 if (bestVar.empty()) {
659 out.add(*
pars()->find(bestVar.c_str()));
664 std::unique_ptr<RooAbsCollection> poi(out.selectByAttrib(
"poi",
true));
675 out.add(*std::unique_ptr<RooAbsCollection>(
pars()->selectByAttrib(
"poi",
true)));
699 Error(
"LoadFits",
"Path not found %s", apath);
709 std::map<std::set<std::tuple<std::string, double, std::string>>, std::set<std::set<std::string>>> cfits;
710 std::set<std::string> allpois;
715 std::cout <<
"Processing " << _dir->GetName() << std::endl;
716 if (
auto keys = _dir->GetListOfKeys(); keys) {
719 bool isNllDir =
false;
720 for (
auto &&k : *keys) {
721 TKey *key =
dynamic_cast<TKey *
>(k);
722 if (strcmp(key->
GetClassName(),
"RooLinkedList") == 0) {
728 for (
auto &&k : *keys) {
729 if (
auto subdir = _dir->GetDirectory(k->GetName()); subdir) {
731 LoadFits(subdir->GetPath());
738 if (cl->InheritsFrom(
"RooFitResult")) {
739 if (
auto cachedFit = _dir->Get<
RooFitResult>(k->GetName()); cachedFit) {
743 std::unique_ptr<RooAbsCollection> snap(cachedFit->floatParsFinal().snapshot());
744 snap->remove(*fPars,
true,
true);
745 fPars->addClone(*snap);
747 for (
auto &
p : cachedFit->constPars()) {
748 if (
p->getAttribute(
"global"))
754 if (!fPars->contains(*
v))
759 std::set<std::string> floatPars;
760 for (
auto &
p : cachedFit->floatParsFinal())
765 std::set<std::tuple<std::string, double, std::string>> constPars;
766 for (
auto &
p : cachedFit->constPars()) {
767 if (
p->getAttribute(
"global"))
774 std::make_tuple(
v->GetName(),
v->getVal(),
775 v->getStringAttribute(
"altVal") ?
v->getStringAttribute(
"altVal") :
""));
778 for (
auto &&[key,
value] : cfits) {
779 if (constPars == key)
781 if (std::includes(constPars.begin(), constPars.end(), key.begin(), key.end())) {
785 std::set<std::string> pois;
786 for (
auto &&par : constPars) {
787 if (key.find(par) == key.end()) {
788 pois.insert(std::get<0>(par));
789 allpois.insert(std::get<0>(par));
793 cfits[constPars].insert(pois);
828 ::Info(
"xRooHypoSpace::xRooHypoSpace",
"%s - Loaded %d fits", apath, nFits);
830 if (allpois.size() == 1) {
831 ::Info(
"xRooHypoSpace::xRooHypoSpace",
"Detected POI: %s", allpois.begin()->c_str());
833 auto nll = std::make_shared<xRooNLLVar>(
nullptr,
nullptr);
834 auto dummyNll = std::make_shared<RooRealVar>(apath,
"Dummy NLL", std::numeric_limits<double>::quiet_NaN());
835 nll->std::shared_ptr<
RooAbsReal>::operator=(dummyNll);
839 for (
auto &&
p : *fPars) {
844 for (
auto &
p : allpois) {
845 fPars->find(
p.c_str())->setAttribute(
"poi",
true);
850 for (
auto &&[key,
value] : cfits) {
853 auto _coords = std::make_shared<RooArgSet>();
854 for (
auto &k : key) {
855 auto v = _coords->addClone(
RooRealVar(std::get<0>(k).c_str(), std::get<0>(k).c_str(), std::get<1>(k)));
856 v->setAttribute(
"poi", allpois.find(std::get<0>(k)) != allpois.end());
857 if (!std::get<2>(k).empty()) {
858 v->setStringAttribute(
"altVal", std::get<2>(k).c_str());
873 if (std::isnan(hp.
fAltVal())) {
884 }
else if (nFits > 0) {
885 std::cout <<
"possible POI: ";
886 for (
auto p : allpois)
887 std::cout <<
p <<
",";
888 std::cout << std::endl;
899 for (
size_t i = 0; i <
size(); i++) {
900 std::cout << i <<
") ";
901 for (
auto a : _axes) {
902 if (
a != _axes.first())
904 std::cout <<
a->GetName() <<
"="
905 << at(i).coords->getRealValue(
a->GetName(), std::numeric_limits<double>::quiet_NaN());
907 std::cout <<
" status=[ufit:";
912 std::cout << ufit->
status();
915 std::cout <<
",cfit_null:";
916 auto cfit =
const_cast<xRooHypoPoint &
>(at(i)).cfit_null(
true);
920 std::cout << cfit->
status();
923 std::cout <<
",cfit_alt:";
924 auto afit =
const_cast<xRooHypoPoint &
>(at(i)).cfit_alt(
true);
928 std::cout << afit->
status();
931 if (
auto asiPoint =
const_cast<xRooHypoPoint &
>(at(i)).asimov(
true)) {
932 std::cout <<
",asimov.ufit:";
933 auto asi_ufit = asiPoint->ufit(
true);
937 std::cout << asi_ufit->status();
940 std::cout <<
",asimov.cfit_null:";
941 auto asi_cfit = asiPoint->cfit_null(
true);
945 std::cout << asi_cfit->status();
950 auto sigma_mu =
const_cast<xRooHypoPoint &
>(at(i)).sigma_mu(
true);
951 if (!std::isnan(sigma_mu.first)) {
952 std::cout <<
" sigma_mu=" << sigma_mu.first;
954 std::cout <<
" +/- " << sigma_mu.second;
956 std::cout << std::endl;
958 std::cout <<
"--------------------------" << std::endl;
959 std::cout <<
"Number of bad fits: " << badFits << std::endl;
963 const char *opt )
const
970 bool readOnly = sOpt.
Contains(
"readonly");
971 bool visualize = sOpt.
Contains(
"visualize") && !readOnly;
978 : std::numeric_limits<double>::quiet_NaN();
980 !std::isnan(nSigma) && nSigma && !(sOpt(sOpt.
Index(
"exp") + 3) ==
'+' || sOpt(sOpt.
Index(
"exp") + 3) ==
'-');
983 if (_axes.size() != 1)
986 auto out = std::make_shared<TGraphErrors>();
987 out->SetName(GetName());
988 out->SetEditable(
false);
989 const char *sCL = (doCLs) ?
"CLs" :
"null";
995 : (!nSigma ?
"Expected"
997 expBand || !nSigma ?
"" : ((nSigma < 0) ?
"-" :
"+"),
int(nSigma))
999 _axes.at(0)->GetTitle(), sCL);
1001 if (std::isnan(nSigma)) {
1003 out->SetMarkerStyle(20);
1004 out->SetMarkerSize(0.5);
1006 out->SetNameTitle(
"obs_ts",
TString::Format(
"Observed;%s;%s", _axes.at(0)->GetTitle(),
1007 (empty() ?
"" : front().tsTitle(
true).Data())));
1010 out->SetNameTitle(
TString::Format(
"exp%d_p%s",
int(nSigma), sCL), title);
1011 out->SetMarkerStyle(0);
1012 out->SetMarkerSize(0);
1013 out->SetLineStyle(2 +
int(nSigma));
1014 if (expBand && nSigma) {
1017 out->SetLineStyle(0);
1018 out->SetLineWidth(0);
1019 auto x = out->Clone(
"up");
1023 out->GetListOfFunctions()->Add(
x,
"F");
1024 x = out->Clone(
"down");
1028 out->GetListOfFunctions()->Add(
x,
"F");
1032 TString::Format(
"Expected;%s;%s", _axes.at(0)->GetTitle(), front().tsTitle(
true).Data()));
1036 auto badPoints = [&]() {
1037 auto badPoints2 =
dynamic_cast<TGraph *
>(out->GetListOfFunctions()->FindObject(
"badPoints"));
1041 badPoints2->SetName(
"badPoints");
1042 badPoints2->SetMarkerStyle(5);
1043 badPoints2->SetMarkerColor(std::isnan(nSigma) ?
kRed :
kBlue);
1044 badPoints2->SetMarkerSize(1);
1045 out->GetListOfFunctions()->Add(badPoints2,
"P");
1049 int nPointsDown = 0;
1054 for (
auto &
p : *
this) {
1058 auto gra =
graph(sOpt +
" readOnly");
1059 if (gra && gra->GetN()) {
1061 gROOT->GetSelectedPad()->cd();
1064 gra->DrawClone(expBand ?
"AF" :
"ALP")->SetBit(
kCanDelete);
1065#if ROOT_VERSION_CODE >= ROOT_VERSION(6, 30, 00)
1066 if (
auto pad =
gROOT->GetSelectedPad()) {
1067 pad->GetCanvas()->ResetUpdated();
1073 ::Info(
"xRooHypoSpace::graph",
"Completed %d/%d points for %s",
int(nDone),
int(
size()), sOpt.
Data());
1079 double _x =
p.coords->getRealValue(_axes.at(0)->GetName(), std::numeric_limits<double>::quiet_NaN());
1081 auto idx = out->GetN() - nPointsDown;
1083 if (std::isnan(pval.first)) {
1084 if (
p.status() != 0) {
1085 badPoints()->SetPoint(badPoints()->GetN(), _x, 0);
1088 out->InsertPointBefore(idx, _x, pval.first);
1089 out->SetPointError(idx, 0, pval.second);
1092 if (expBand && nSigma) {
1096 if (std::isnan(pval.first)) {
1097 if (
p.status() != 0) {
1098 badPoints()->SetPoint(badPoints()->GetN(), _x, 0);
1101 out->InsertPointBefore(idx + 1, _x, pval.first);
1102 out->SetPointError(idx + 1, 0, pval.second);
1104 if (out->GetPointY(idx) < pval.first)
1111 if (out->GetN() == 0)
1116 if (out->GetListOfFunctions()->FindObject(
"badPoints")) {
1118 for (
int i = 0; i < badPoints()->GetN(); i++) {
1119 badPoints()->SetPointY(i, out->Eval(badPoints()->GetPointX(i)));
1124 out->Sort(&
TGraph::CompareX,
false, out->GetN() - nPointsDown, out->GetN() - 1);
1127 auto up =
dynamic_cast<TGraph *
>(out->GetListOfFunctions()->FindObject(
"up"));
1128 auto down =
dynamic_cast<TGraph *
>(out->GetListOfFunctions()->FindObject(
"down"));
1130 for (
int i = 0; i < out->GetN(); i++) {
1131 if (i < out->GetN() - nPointsDown) {
1132 up->
SetPoint(up->GetN(), out->GetPointX(i), out->GetPointY(i) + out->GetErrorY(i) * (above ? 1. : -1.));
1133 down->SetPoint(down->GetN(), out->GetPointX(i), out->GetPointY(i) - out->GetErrorY(i) * (above ? 1. : -1.));
1135 up->SetPoint(up->GetN(), out->GetPointX(i), out->GetPointY(i) - out->GetErrorY(i) * (above ? 1. : -1.));
1136 down->SetPoint(down->GetN(), out->GetPointX(i), out->GetPointY(i) + out->GetErrorY(i) * (above ? 1. : -1.));
1144 gROOT->GetSelectedPad()->cd();
1147 out->DrawClone(expBand ?
"AF" :
"ALP")->SetBit(
kCanDelete);
1148#if ROOT_VERSION_CODE >= ROOT_VERSION(6, 30, 00)
1149 if (
auto pad =
gROOT->GetSelectedPad()) {
1150 pad->GetCanvas()->ResetUpdated();
1163 std::shared_ptr<TMultiGraph> out;
1166 bool visualize = sOpt.
Contains(
"visualize");
1169 auto exp2 =
graph(sOpt +
" exp2");
1170 auto exp1 =
graph(sOpt +
" exp1");
1171 auto exp =
graph(sOpt +
" exp");
1174 auto obs = (doObs) ?
graph(sOpt) :
nullptr;
1176 out = std::make_shared<TMultiGraph>(GetName(), GetTitle());
1177 if (exp2 && exp2->GetN() > 1)
1178 out->Add(
static_cast<TGraph *
>(exp2->Clone()),
"FP");
1179 if (exp1 && exp1->GetN() > 1)
1180 out->Add(
static_cast<TGraph *
>(exp1->Clone()),
"FP");
1181 if (exp && exp->GetN() > 1)
1182 out->Add(
static_cast<TGraph *
>(exp->Clone()),
"LP");
1183 if (obs && obs->GetN() > 1)
1184 out->Add(
static_cast<TGraph *
>(obs->Clone()),
"LP");
1186 if (!out->GetListOfGraphs()) {
1190 TGraph *testedPoints =
nullptr;
1193 line->SetName(
"alpha");
1195 line->SetEditable(
false);
1196 line->SetPoint(
line->GetN(), out->GetHistogram()->GetXaxis()->GetXmin() - 10, 0.05);
1197 testedPoints =
new TGraph;
1198 testedPoints->
SetName(
"hypoPoints");
1203 for (
int i = 0; i < exp->GetN(); i++) {
1204 testedPoints->
SetPoint(testedPoints->
GetN(), exp->GetPointX(i), 0.05);
1207 line->SetPoint(
line->GetN(), out->GetHistogram()->GetXaxis()->GetXmax() + 10, 0.05);
1209 out->GetListOfFunctions()->Add(
line,
"L");
1212 out->GetHistogram()->GetXaxis()->SetTitle(exp->GetHistogram()->GetXaxis()->GetTitle());
1213 out->GetHistogram()->GetYaxis()->SetTitle(exp->GetHistogram()->GetYaxis()->GetTitle());
1217 leg->SetName(
"legend");
1220 out->GetListOfFunctions()->Add(
leg);
1223 for (
auto g : *out->GetListOfGraphs()) {
1225 leg->AddEntry(o,
"",
"F");
1227 leg->AddEntry(
g,
"",
"LPE");
1233 if (exp2 && exp2->GetN() > 1) {
1237 if (exp1 && exp1->GetN() > 1) {
1241 if (exp && exp->GetN() > 1) {
1245 if (exp1 && exp1->GetN() > 1) {
1249 if (exp2 && exp2->GetN() > 1) {
1253 if (obs && obs->GetN() > 1) {
1259 out->Add(testedPoints,
"P");
1263 gROOT->GetSelectedPad()->cd();
1266 auto gra2 =
static_cast<TMultiGraph *
>(out->DrawClone(
"A"));
1269 gra2->GetHistogram()->SetMinimum(1
e-6);
1273 gPad->GetCanvas()->Paint();
1274 gPad->GetCanvas()->Update();
1275#if ROOT_VERSION_CODE >= ROOT_VERSION(6, 30, 00)
1276 gPad->GetCanvas()->ResetUpdated();
1289 if (std::isnan(
target)) {
1293 auto gr = std::make_shared<TGraph>(pValues);
1296 std::set<double> existingX;
1297 while (i < gr->GetN()) {
1300 }
else if (existingX.find(
gr->
GetPointX(i)) != existingX.end()) {
1314 return std::pair(std::numeric_limits<double>::quiet_NaN(), 0);
1317 double alpha = log(
target);
1320 for (
int ii = 1; ii <
gr->
GetN(); ii++) {
1331 return std::pair(lim, err);
1346 return std::pair(std::numeric_limits<double>::infinity(), std::numeric_limits<double>::infinity());
1347 return std::pair((alpha -
y1) /
m +
x1, std::numeric_limits<double>::infinity());
1354 return std::pair(-std::numeric_limits<double>::infinity(), -std::numeric_limits<double>::infinity());
1355 return std::pair((alpha -
y1) /
m +
x1, -std::numeric_limits<double>::infinity());
1362 if (std::isnan(nSigma)) {
1367 return GetLimit(*
graph(sOpt +
" readonly"));
1374 bool visualize = sOpt.
Contains(
"visualize");
1376 std::shared_ptr<TGraphErrors>
gr =
graph(sOpt +
" readonly");
1378 auto gra = graphs(sOpt.
Contains(
"toys") ?
"pcls readonly toys" :
"pcls readonly");
1385 gra->GetHistogram()->SetMinimum(1
e-9);
1386 gra->GetHistogram()->GetYaxis()->SetRangeUser(1
e-9, 1);
1388#if ROOT_VERSION_CODE >= ROOT_VERSION(6, 30, 00)
1389 gPad->GetCanvas()->ResetUpdated();
1396 for (
auto p : axes()) {
1397 for (
auto &[
pdf, nll] : fNlls) {
1398 if (
auto _v =
dynamic_cast<RooRealVar *
>(nll->pars()->find(*
p))) {
1399 dynamic_cast<RooRealVar *
>(
p)->setRange(_v->getMin(), _v->getMax());
1405 auto v = (axes().empty()) ?
nullptr :
dynamic_cast<RooRealVar *
>(*axes().rbegin());
1407 return std::pair(std::numeric_limits<double>::quiet_NaN(), 0.);
1408 double muMax = std::min(std::min(
v->getMax(
"physical"),
v->getMax()),
v->getMax(
"scan"));
1409 double muMin = std::max(std::max(
v->getMin(
"physical"),
v->getMin()),
v->getMin(
"scan"));
1411 if (maxTries == 0 || std::isnan(AddPoint(
TString::Format(
"%s=%g",
v->GetName(), muMin)).getVal(sOpt).first)) {
1413 Error(
"findlimit",
"Problem evaluating %s @ %s=%g", sOpt.
Data(),
v->GetName(), muMin);
1414 return std::pair(std::numeric_limits<double>::quiet_NaN(), 0.);
1417 return findlimit(opt, relUncert, maxTries - 1);
1424 double nextPoint = muMin + (muMax - muMin) / 50;
1430 auto expLim = findlimit(sOpt2, std::numeric_limits<double>::infinity(), 0);
1431 if (!std::isnan(expLim.first) && expLim.first < nextPoint)
1432 nextPoint = expLim.first;
1436 (sOpt.
Contains(
"exp")) ? back().asimov() : std::shared_ptr<xRooHypoPoint>(&back(), [](
xRooHypoPoint *) {});
1438 if (point && point->ufit()) {
1439 double rough_sigma_mu = point->mu_hat().getError();
1442 nextPoint = another_estimate;
1443 ::Info(
"xRooHypoSpace::findlimit",
"Guessing %g based on rough sigma_mu = %g", nextPoint, rough_sigma_mu);
1447 if (maxTries == 0 || std::isnan(AddPoint(
TString::Format(
"%s=%g",
v->GetName(), nextPoint)).getVal(sOpt).first)) {
1449 Error(
"findlimit",
"Problem evaluating %s @ %s=%g", sOpt.
Data(),
v->GetName(), nextPoint);
1450 return std::pair(std::numeric_limits<double>::quiet_NaN(), 0.);
1453 return findlimit(opt, relUncert, maxTries - 1);
1456 auto lim = GetLimit(*
gr);
1458 if (std::isnan(lim.first)) {
1462 auto v =
dynamic_cast<RooRealVar *
>(*axes().rbegin());
1463 double maxMu = std::min(std::min(
v->getMax(
"physical"),
v->getMax()),
v->getMax(
"scan"));
1464 double minMu = std::max(std::max(
v->getMin(
"physical"),
v->getMin()),
v->getMin(
"scan"));
1467 if (lim.first > -std::numeric_limits<double>::infinity() && lim.first < std::numeric_limits<double>::infinity() &&
1468 (std::abs(lim.second) <= relUncert * std::abs(lim.first) ))
1473 if (lim.second == std::numeric_limits<double>::infinity()) {
1475 nextPoint = lim.first;
1476 if (nextPoint == std::numeric_limits<double>::infinity() || nextPoint > maxMu) {
1483 (sOpt.
Contains(
"exp")) ? back().asimov() : std::shared_ptr<xRooHypoPoint>(&back(), [](
xRooHypoPoint *) {});
1485 if (point && point->ufit()) {
1486 double rough_sigma_mu = point->mu_hat().getError();
1489 nextPoint = std::max(nextPoint, another_estimate);
1490 ::Info(
"xRooHypoSpace::findlimit",
"Guessing %g based on rough sigma_mu = %g", nextPoint, rough_sigma_mu);
1493 nextPoint = std::min(nextPoint + nextPoint * relUncert * 0.99, maxMu);
1495 if (nextPoint > maxMu)
1497 }
else if (lim.second == -std::numeric_limits<double>::infinity()) {
1499 nextPoint = lim.first;
1500 if (nextPoint < minMu)
1501 nextPoint =
gr->
GetPointX(0) - (maxMu - minMu) / 50;
1502 if (nextPoint < minMu)
1505 nextPoint = lim.first + lim.second * relUncert * 0.99;
1511 ::Info(
"xRooHypoSpace::findlimit",
"%s -- Testing new point @ %s=%g (delta=%g)", sOpt.
Data(),
v->GetName(),
1512 nextPoint, lim.second);
1513 if (maxTries == 0 || std::isnan(AddPoint(
TString::Format(
"%s=%g",
v->GetName(), nextPoint)).getVal(sOpt).first)) {
1514 if (maxTries == 0) {
1515 Warning(
"findlimit",
"Reached max number of point evaluations");
1517 Error(
"findlimit",
"Problem evaluating %s @ %s=%g", sOpt.
Data(),
v->GetName(), nextPoint);
1522 return findlimit(opt, relUncert, maxTries - 1);
1531 if ((sOpt ==
"" || sOpt ==
"same") && !empty()) {
1534 for (
auto &hp : *
this) {
1535 if (!hp.nullToys.empty() || !hp.altToys.empty()) {
1546 auto _axes = axes();
1550 if (sOpt ==
"status") {
1552 if (_axes.size() <= 2) {
1555 out->SetName(
"points");
1556 out->SetMarkerSize(0.5);
1567 badPoints->
SetName(
"bad_ufit");
1573 badPoints2->
SetName(
"bad_cfit_null");
1579 out->SetTitle(
TString::Format(
"%s;%s;%s", GetTitle(), _axes.at(0)->GetTitle(),
1580 (_axes.size() == 1) ?
"" : _axes.at(1)->GetTitle()));
1581 for (
auto &
p : *
this) {
1582 bool _readOnly =
p.nllVar ?
p.nllVar->get()->getAttribute(
"readOnly") :
false;
1584 p.nllVar->get()->setAttribute(
"readOnly",
true);
1585 double x =
p.coords->getRealValue(_axes.at(0)->GetName());
1586 double y = _axes.size() == 1 ?
p.ts_asymp().first :
p.coords->getRealValue(_axes.at(1)->GetName());
1587 out->SetPoint(out->GetN(),
x,
y);
1588 if (!std::isnan(
p.ts_asymp().first)) {
1589 if (_axes.size() == 1)
1590 out->SetPointError(out->GetN() - 1, 0,
p.ts_asymp().second);
1592 }
else if (
p.fUfit && (std::isnan(
p.fUfit->minNll()) ||
1596 }
else if (
p.fNull_cfit && (std::isnan(
p.fNull_cfit->minNll()) ||
1601 if (!std::isnan(
p.ts_asymp(0).first)) {
1603 }
else if (
p.asimov() &&
p.asimov()->fUfit &&
1604 (std::isnan(
p.asimov()->fUfit->minNll()) ||
1608 }
else if (
p.asimov() &&
p.asimov()->fNull_cfit &&
1609 (std::isnan(
p.asimov()->fNull_cfit->minNll()) ||
1614 p.nllVar->get()->setAttribute(
"readOnly", _readOnly);
1617 if (_axes.size() == 1) {
1619 for (
int i = 0; i < out->GetN(); i++) {
1620 if (!std::isnan(out->GetPointY(i)))
1621 tmp.SetPoint(tmp.GetN(), out->GetPointX(i), out->GetPointY(i));
1623 auto fixPoints = [&](
TGraph *
g) {
1624 for (
int i = 0; i <
g->GetN(); i++) {
1625 if (std::isnan(
g->GetPointY(i)))
1626 g->SetPointY(i, std::isnan(tmp.Eval(
g->GetPointX(i))) ? 0. : tmp.Eval(
g->GetPointX(i)));
1631 fixPoints(expAvail);
1632 fixPoints(badPoints);
1633 fixPoints(badPoints2);
1636 out->SetMarkerStyle(4);
1638 auto leg =
new TLegend(1. -
gPad->GetRightMargin() - 0.3, 1. -
gPad->GetTopMargin() - 0.35,
1639 1. -
gPad->GetRightMargin() - 0.05, 1. -
gPad->GetTopMargin() - 0.05);
1640 leg->SetName(
"legend");
1641 leg->AddEntry(out,
"Uncomputed",
"P");
1643 if (tsAvail->
GetN()) {
1644 out->GetListOfFunctions()->Add(tsAvail,
"P");
1645 leg->AddEntry(tsAvail,
"Computed",
"P");
1649 if (expAvail->
GetN()) {
1650 out->GetListOfFunctions()->Add(expAvail,
"P");
1651 leg->AddEntry(expAvail,
"Expected computed",
"P");
1655 if (badPoints->
GetN()) {
1656 out->GetListOfFunctions()->Add(badPoints,
"P");
1657 leg->AddEntry(badPoints,
"Bad ufit",
"P");
1661 if (badPoints2->
GetN()) {
1662 out->GetListOfFunctions()->Add(badPoints2,
"P");
1663 leg->AddEntry(badPoints2,
"Bad null cfit",
"P");
1670 gPad->SetGrid(
true, _axes.size() > 1);
1671 if (_axes.size() == 1)
1672 gPad->SetLogy(
false);
1681 auto gra = graphs(sOpt +
" readonly");
1683 gROOT->GetSelectedPad()->cd();
1691 gra2->GetHistogram()->SetMinimum(1
e-6);
1695 gPad->GetCanvas()->Paint();
1696 gPad->GetCanvas()->Update();
1697#if ROOT_VERSION_CODE >= ROOT_VERSION(6, 30, 00)
1698 gPad->GetCanvas()->ResetUpdated();
1719 out->SetName(GetName());
1724 if (!empty() && axes().
size() == 1) {
1725 for (
auto &
p : *
this) {
1727 pllType =
p.fPllType;
1731 title += front().tsTitle(
true);
1734 out->SetTitle(title);
1735 *
dynamic_cast<TAttFill *
>(out) = *
this;
1736 *
dynamic_cast<TAttLine *
>(out) = *
this;
1746 bool doFits =
false;
1752 auto mainPad =
gPad;
1754 out->SetEditable(
false);
1759 gPad->SetBottomMargin(
gPad->GetBottomMargin() * 2.);
1763 mainPad = basePad->
GetPad(1);
1769 std::pair<double, double> minMax(std::numeric_limits<double>::infinity(), -std::numeric_limits<double>::infinity());
1770 for (
auto &
p : *
this) {
1771 if (
p.fPllType != pllType)
1773 auto val =
p.pll(
true).first;
1774 if (std::isnan(val))
1776 minMax.first = std::min(minMax.first, val);
1777 minMax.second = std::max(minMax.second, val);
1779 if (minMax.first < std::numeric_limits<double>::infinity())
1780 out->GetHistogram()->SetMinimum(minMax.first);
1781 if (minMax.second > -std::numeric_limits<double>::infinity())
1782 out->GetHistogram()->SetMaximum(minMax.second);
1784 TGraph *badPoints =
nullptr;
1788 std::shared_ptr<const RooFitResult> ufr;
1789 for (
auto &
p : *
this) {
1790 if (
p.fPllType != pllType)
1792 auto val =
p.pll().first;
1795 if (out->GetN() == 0 && ufr && ufr->status() == 0) {
1796 out->SetPoint(out->GetN(),
1797 ufr->floatParsFinal().getRealValue(axes().first()->GetName(),
1798 ufr->constPars().getRealValue(axes().first()->GetName())),
1800 out->SetPointError(out->GetN() - 1, 0, ufr->edm());
1802 if (
auto fr =
p.fNull_cfit;
1807 new TPad(fr->GetName(),
TString::Format(
"%s = %g", poi().first()->GetTitle(),
p.fNullVal()), 0, 0, 1., 1);
1808 pad->SetNumber(out->GetN() + 1);
1815 if (std::isnan(val) &&
p.status() != 0) {
1819 badPoints->
SetName(
"badPoints");
1823 out->GetListOfFunctions()->Add(badPoints,
"P");
1825 badPoints->
SetPoint(badPoints->
GetN(),
p.fNullVal(), out->Eval(
p.fNullVal()));
1826 mainPad->Modified();
1827 }
else if (!std::isnan(val)) {
1828 out->SetPoint(out->GetN(),
p.coords->getRealValue(axes().first()->GetName()),
p.pll().first);
1829 out->SetPointError(out->GetN() - 1, 0,
p.pll().second);
1834 for (
int i = 0; i < badPoints->
GetN(); i++)
1838 mainPad->Modified();
1851#if ROOT_VERSION_CODE >= ROOT_VERSION(6, 30, 00)
1857 if (ufr && doFits) {
1859 auto pad =
new TPad(ufr->GetName(),
"unconditional fit", 0, 0, 1., 1.);
1867 pad =
new TPad(
"selected",
"selected", 0, 0, 1, 1);
1882 gPad->GetCanvas()->
Connect(
"Highlighted(TVirtualPad*,TObject*,Int_t,Int_t)",
"xRooNode::InteractiveObject",
1894 auto _axes = axes();
1899 out->SetTitle(GetTitle());
1901 for (
auto &
p : *
this) {
1902 double _x =
p.coords->getRealValue(_axes.at(0)->GetName(), std::numeric_limits<double>::quiet_NaN());
1903 out->Add(_x,
p.result());
size_t size(const MatrixT &matrix)
retrieve the size of a square matrix
void Info(const char *location, const char *msgfmt,...)
Use this function for informational messages.
void Error(const char *location, const char *msgfmt,...)
Use this function in case an error occurred.
void Warning(const char *location, const char *msgfmt,...)
Use this function in warning situations.
winID h TVirtualViewer3D TVirtualGLPainter p
Option_t Option_t SetFillStyle
Option_t Option_t TPoint TPoint const char GetTextMagnitude GetFillStyle GetLineColor GetLineWidth GetMarkerStyle GetTextAlign GetTextColor GetTextSize void char Point_t Rectangle_t WindowAttributes_t Float_t Float_t Float_t Int_t Int_t UInt_t UInt_t Rectangle_t Int_t Int_t Window_t TString Int_t GCValues_t GetPrimarySelectionOwner GetDisplay GetScreen GetColormap GetNativeEvent const char const char dpyName wid window const char font_name cursor keysym reg const char only_if_exist regb h Point_t winding char text const char depth char const char Int_t count const char ColorStruct_t color const char Pixmap_t Pixmap_t PictureAttributes_t attr const char char ret_data h unsigned char height h offset
Option_t Option_t TPoint TPoint const char GetTextMagnitude GetFillStyle GetLineColor GetLineWidth GetMarkerStyle GetTextAlign GetTextColor GetTextSize void char Point_t Rectangle_t WindowAttributes_t Float_t Float_t Float_t Int_t Int_t UInt_t UInt_t Rectangle_t Int_t Int_t Window_t TString Int_t GCValues_t GetPrimarySelectionOwner GetDisplay GetScreen GetColormap GetNativeEvent const char const char dpyName wid window const char font_name cursor keysym reg const char only_if_exist regb h Point_t winding char text const char depth char const char Int_t count const char ColorStruct_t color const char Pixmap_t Pixmap_t PictureAttributes_t attr const char char ret_data h unsigned char height h Atom_t Int_t ULong_t ULong_t unsigned char prop_list Atom_t Atom_t target
Option_t Option_t TPoint TPoint const char GetTextMagnitude GetFillStyle GetLineColor GetLineWidth GetMarkerStyle GetTextAlign GetTextColor GetTextSize void char Point_t Rectangle_t WindowAttributes_t Float_t Float_t Float_t Int_t Int_t UInt_t UInt_t Rectangle_t result
Option_t Option_t TPoint TPoint const char GetTextMagnitude GetFillStyle GetLineColor GetLineWidth GetMarkerStyle GetTextAlign GetTextColor GetTextSize void value
Option_t Option_t TPoint TPoint const char x1
Option_t Option_t SetFillColor
Option_t Option_t TPoint TPoint const char GetTextMagnitude GetFillStyle GetLineColor GetLineWidth GetMarkerStyle GetTextAlign GetTextColor GetTextSize void char Point_t Rectangle_t WindowAttributes_t Float_t Float_t Float_t Int_t Int_t UInt_t UInt_t Rectangle_t Int_t Int_t Window_t TString Int_t GCValues_t GetPrimarySelectionOwner GetDisplay GetScreen GetColormap GetNativeEvent const char const char dpyName wid window const char font_name cursor keysym reg const char only_if_exist regb h Point_t winding char text const char depth char const char Int_t count const char ColorStruct_t color const char Pixmap_t Pixmap_t PictureAttributes_t attr const char char ret_data h unsigned char height h Atom_t Int_t ULong_t ULong_t unsigned char prop_list Atom_t Atom_t Atom_t Time_t type
Option_t Option_t TPoint TPoint const char y1
R__EXTERN TStyle * gStyle
R__EXTERN TSystem * gSystem
static std::pair< double, double > matchPrecision(const std::pair< double, double > &in)
static std::set< int > allowedStatusCodes
std::shared_ptr< const RooAbsCollection > coords
std::shared_ptr< xRooNLLVar > nllVar
xRooFit::Asymptotics::PLLType fPllType
xValueWithError limit(const char *type="cls", double nSigma=std::numeric_limits< double >::quiet_NaN()) const
xRooHypoPoint & point(size_t i)
RooStats::HypoTestInverterResult * result()
bool AddModel(const xRooNode &pdf, const char *validity="")
std::shared_ptr< TMultiGraph > graphs(const char *opt)
static RooArgList toArgs(const char *str)
xValueWithError findlimit(const char *opt, double relUncert=std::numeric_limits< double >::infinity(), unsigned int maxTries=20)
int AddPoints(const char *parName, size_t nPoints, double low, double high)
static std::pair< double, double > GetLimit(const TGraph &pValues, double target=std::numeric_limits< double >::quiet_NaN())
std::shared_ptr< TGraphErrors > graph(const char *opt) const
std::shared_ptr< RooArgSet > fPars
void Print(Option_t *opt="") const override
Print TNamed name and title.
xRooHypoPoint & AddPoint(double value)
xRooHypoSpace(const char *name="", const char *title="")
int scan(const char *type, size_t nPoints, double low=std::numeric_limits< double >::quiet_NaN(), double high=std::numeric_limits< double >::quiet_NaN(), const std::vector< double > &nSigmas={0, 1, 2, -1, -2, std::numeric_limits< double >::quiet_NaN()}, double relUncert=0.1)
std::shared_ptr< xRooNode > pdf(const RooAbsCollection &parValues) const
void Draw(Option_t *opt="") override
Default Draw method for all objects.
std::map< std::string, std::pair< double, double > > limits(const char *opt="cls", const std::vector< double > &nSigmas={0, 1, 2, -1, -2, std::numeric_limits< double >::quiet_NaN()}, double relUncert=std::numeric_limits< double >::infinity())
void LoadFits(const char *apath)
std::shared_ptr< RooAbsPdf > pdf() const
std::shared_ptr< RooArgSet > pars(bool stripGlobalObs=true) const
The xRooNode class is designed to wrap over a TObject and provide functionality to aid with interacti...
static InteractiveObject * gIntObj
void Draw(Option_t *opt="") override
Default Draw method for all objects.
RooArgList argList() const
xRooNode pars() const
List of parameters (non-observables) of this node.
void setAttribute(const Text_t *name, bool value=true)
Set (default) or clear a named boolean attribute of this object.
Abstract container object that can hold multiple RooAbsArg objects.
virtual bool add(const RooAbsArg &var, bool silent=false)
Add the specified argument to list.
void sort(bool reverse=false)
Sort collection using std::sort and name comparison.
RooAbsArg * find(const char *name) const
Find object with given name in list.
Abstract interface for all probability density functions.
Abstract base class for objects that represent a real value that may appear on the left hand side of ...
void setConstant(bool value=true)
Abstract base class for objects that represent a real value and implements functionality common to al...
RooArgList is a container object that can hold multiple RooAbsArg objects.
RooArgSet is a container object that can hold multiple RooAbsArg objects.
Represents a constant real-valued object.
Container class to hold unbinned data.
RooFitResult is a container class to hold the input and output of a PDF fit to a dataset.
Variable that can be changed from the outside.
HypoTestInverterResult class holds the array of hypothesis test results and compute a confidence inte...
HypoTestResult * GetResult(int index) const
return a pointer to the i^th result object
int ArraySize() const
number of entries in the results array
double GetXValue(int index) const
function to return the value of the parameter of interest for the i^th entry in the results
RooArgSet * GetParameters() const override
return a cloned list with the parameter of interest
Fill Area Attributes class.
virtual void SetLineStyle(Style_t lstyle)
Set the line style.
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.
static TCanvas * MakeDefCanvas()
Static function to build a default canvas.
void Paint(Option_t *option="") override
Paint canvas.
void Update() override
Update canvas pad buffers.
static TClass * GetClass(const char *name, Bool_t load=kTRUE, Bool_t silent=kFALSE)
Static method returning pointer to TClass of the specified class name.
Describe directory structure in memory.
virtual Bool_t cd()
Change current directory to "this" directory.
virtual Int_t GetValue(const char *name, Int_t dflt) const
Returns the integer value for a resource.
static TFile * Open(const char *name, Option_t *option="", const char *ftitle="", Int_t compress=ROOT::RCompressionSetting::EDefaults::kUseCompiledDefault, Int_t netopt=0)
Create / open a file.
A TGraphErrors is a TGraph with error bars.
A TGraph is an object made of two arrays X and Y with npoints each.
virtual Double_t GetPointX(Int_t i) const
Get x value for point i.
virtual void SetPoint(Int_t i, Double_t x, Double_t y)
Set x and y values for point number i.
static Bool_t CompareX(const TGraph *gr, Int_t left, Int_t right)
Return kTRUE if fX[left] > fX[right]. Can be used by Sort.
virtual void Sort(Bool_t(*greater)(const TGraph *, Int_t, Int_t)=&TGraph::CompareX, Bool_t ascending=kTRUE, Int_t low=0, Int_t high=-1111)
Sorts the points of this TGraph using in-place quicksort (see e.g.
TList * GetListOfFunctions() const
virtual Int_t RemovePoint()
Delete point close to the mouse position Returns index of removed point (or -1 if nothing was changed...
void SetName(const char *name="") override
Set graph name.
virtual Double_t GetPointY(Int_t i) const
Get y value for point i.
virtual void SetPointY(Int_t i, Double_t y)
Set y value for point i.
virtual void SetEditable(Bool_t editable=kTRUE)
if editable=kFALSE, the graph cannot be modified with the mouse by default a TGraph is editable
Book space in a file, create I/O buffers, to fill them, (un)compress them.
virtual const char * GetClassName() const
This class displays a legend box (TPaveText) containing several legend entries.
TObject * FindObject(const char *name) const override
Find an object in this list using its name.
A TMemFile is like a normal TFile except that it reads and writes only from memory.
A TMultiGraph is a collection of TGraph (or derived) objects.
The TNamed class is the base class for all named ROOT classes.
const char * GetName() const override
Returns name of object.
const char * GetTitle() const override
Returns title of object.
virtual void SetName(const char *name)
Set the name of the TNamed.
virtual void SetNameTitle(const char *name, const char *title)
Set all the TNamed parameters (name and title).
Mother of all ROOT objects.
virtual const char * GetName() const
Returns name of object.
virtual void AppendPad(Option_t *option="")
Append graphics object to current pad.
void SetBit(UInt_t f, Bool_t set)
Set or unset the user status bits as specified in f.
The most important graphics class in the ROOT system.
Bool_t Connect(const char *signal, const char *receiver_class, void *receiver, const char *slot)
Non-static method is used to connect from the signal of this object to the receiver slot.
Double_t RealTime()
Stop the stopwatch (if it is running) and return the realtime (in seconds) passed between the start a...
void Start(Bool_t reset=kTRUE)
Start the stopwatch.
void Continue()
Resume a stopped stopwatch.
Provides iteration through tokens of a given string.
Bool_t NextToken()
Get the next token, it is stored in this TString.
void ToLower()
Change string to lower-case.
Double_t Atof() const
Return floating-point value contained in string.
Bool_t IsFloat() const
Returns kTRUE if string contains a floating point or integer number.
const char * Data() const
TString & ReplaceAll(const TString &s1, const TString &s2)
Bool_t BeginsWith(const char *s, ECaseCompare cmp=kExact) const
static TString Format(const char *fmt,...)
Static method which formats a string using a printf style format descriptor and return a TString.
Bool_t Contains(const char *pat, ECaseCompare cmp=kExact) const
Ssiz_t Index(const char *pat, Ssiz_t i=0, ECaseCompare cmp=kExact) const
Float_t GetPadRightMargin() const
Float_t GetPadTopMargin() const
virtual Bool_t ProcessEvents()
Process pending events (GUI, timers, sockets).
This class defines a UUID (Universally Unique IDentifier), also known as GUIDs (Globally Unique IDent...
TVirtualPad is an abstract base class for the Pad and Canvas classes.
virtual void SetGrid(Int_t valuex=1, Int_t valuey=1)=0
virtual TVirtualPad * cd(Int_t subpadnumber=0)=0
virtual TVirtualPad * GetPad(Int_t subpadnumber) const =0
virtual TCanvas * GetCanvas() const =0
void Clear(Option_t *option="") override=0
double gaussian_quantile(double z, double sigma)
Inverse ( ) of the cumulative distribution function of the lower tail of the normal (Gaussian) distri...
#define END_XROOFIT_NAMESPACE