135 ccoutE(InputArguments) <<
"Parameter " << arg->
GetName() <<
"is not in the list of LikelihoodInterval parameters "
136 <<
" - do not use for plotting " << std::endl;
145 ccoutE(InputArguments) <<
"LikelihoodIntervalPlot::Draw(" <<
GetName()
146 <<
") ERROR: contours for more than 2 dimensions not implemented!" << std::endl;
224 TF1 *
f1 =
static_cast<TF1*
>(tmp->Clone());
344 for (std::size_t i = 0; i < params.
size(); ++i) {
373 int nBins =
int( std::sqrt(
double(
nPoints)) + 0.5 );
387 for (
int k = 1; k <
nLevels; ++k) {
394 hist2D->DrawClone(
"COLZ");
403 double confLevels[
nLevels] = { 0.1,0.3,0.5,0.683,0.95,0.9973,0.9999366575,0.9999994267};
404 for (
int k = 0; k <
nLevels; ++k) {
413 tmpOpt.ReplaceAll(
"same",
"");
414 if (
tmpOpt.Length() < 3) opt +=
"cont3";
430 h->Draw(
"CONT LIST");
445 if (
tmpOpt.Contains(
"cont4")) {
446 double bm =
gPad->GetBottomMargin();
447 double lm =
gPad->GetLeftMargin();
448 double rm =
gPad->GetRightMargin();
449 double tm =
gPad->GetTopMargin();
450 double x1 =
hist2D->GetXaxis()->GetXmin();
451 double y1 =
hist2D->GetYaxis()->GetXmin();
452 double x2 =
hist2D->GetXaxis()->GetXmax();
453 double y2 =
hist2D->GetYaxis()->GetXmax();
455 TPad *null=
new TPad(
"null",
"null",0,0,1,1);
456 null->SetFillStyle(0);
457 null->SetFrameFillStyle(0);
489 ccoutE(InputArguments) <<
"LikelihoodIntervalPlot::Draw(" <<
GetName()
490 <<
") ERROR: no contours found in ListOfSpecial" << std::endl;
505 std::cout <<
"Warning - Less points calculated in contours np = " <<
ncp <<
" / " <<
nPoints << std::endl;
546 gr0->SetMarkerStyle(33);
552 delete bestFitParams;
ROOT::Detail::TRangeCast< T, true > TRangeDynCast
TRangeDynCast is an adapter class that allows the typed iteration through a TCollection.
void GetParameters(TFitEditor::FuncParams_t &pars, TF1 *func)
Stores the parameters of the given function into pars.
Option_t Option_t TPoint TPoint const char x2
Option_t Option_t TPoint TPoint const char x1
Option_t Option_t TPoint TPoint const char y2
Option_t Option_t TPoint TPoint const char y1
char * Form(const char *fmt,...)
Formats a string in a circular formatting buffer.
double getRealValue(const char *name, double defVal=0.0, bool verbose=false) const
Get value of a RooAbsReal stored in set with given name.
virtual bool remove(const RooAbsArg &var, bool silent=false, bool matchByNameOnly=false)
Remove the specified argument from our list.
const char * GetName() const override
Returns name of object.
void assign(const RooAbsCollection &other) const
Sets the value, cache and constant attribute of any argument in our set that also appears in the othe...
Storage_t::size_type size() const
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.
TObject * clone(const char *newname) const override
Named container for two doubles, two integers two object points and three string pointers that can be...
Plot frame and a container for graphics objects within that frame.
void SetTitle(const char *name) override
Set the title of the RooPlot to 'title'.
void addObject(TObject *obj, Option_t *drawOptions="", bool invisible=false)
Add a generic object to this plot.
void Draw(Option_t *options=nullptr) override
Draw this plot and all of the elements it contains.
Implements the profile likelihood estimator for a given likelihood and set of parameters of interest.
Variable that can be changed from the outside.
void setVal(double value) override
Set value of variable to 'value'.
This class provides simple and straightforward utilities to plot a LikelihoodInterval object.
TObject * fPlotObject
plotted object
Color_t fColor
color for the contour (for 2D) or function (in 1D)
void SetPlotParameters(const RooArgSet *params)
Int_t fNPoints
number of points used to scan the PL, default depends if 1D or 2D
Style_t fFillStyle
fill style for contours, half transparent by default
void Draw(const Option_t *options=nullptr) override
draw the likelihood interval or contour for the 1D case a RooPlot is drawn by default of the profiled...
LikelihoodIntervalPlot()
LikelihoodIntervalPlot default constructor with default parameters.
double fPrecision
RooCurve precision, use default in case of -1.
double fMaximum
function maximum
LikelihoodInterval * fInterval
void SetLikelihoodInterval(LikelihoodInterval *theInterval)
Color_t fLineColor
line color for the interval (1D) or for other contours (2D)
LikelihoodInterval is a concrete implementation of the RooStats::ConfInterval interface.
double ConfidenceLevel() const override
return confidence level
double UpperLimit(const RooRealVar ¶m)
return the upper bound of the interval on a given parameter
Int_t GetContourPoints(const RooRealVar ¶mX, const RooRealVar ¶mY, double *x, double *y, Int_t npoints=30)
return the 2D-contour points for the given subset of parameters by default make the contour using 30 ...
RooArgSet * GetParameters() const override
return a cloned list of parameters of interest. User manages the return object
double LowerLimit(const RooRealVar ¶m)
return the lower bound of the interval on a given parameter
RooAbsReal * GetLikelihoodRatio()
return the profile log-likelihood ratio function
const RooArgSet * GetBestFitParameters() const
return a pointer to a snapshot with best fit parameter of interest
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 TH1 * GetHistogram() const
Return a pointer to the histogram used to visualise the function Note that this histogram is managed ...
TAxis * GetYaxis() const
Get y axis of the function.
virtual void SetMaximum(Double_t maximum=-1111)
Set the maximum value along Y for this function In case the function is already drawn,...
virtual void SetRange(Double_t xmin, Double_t xmax)
Initialize the upper and lower bounds to draw the function.
void SetTitle(const char *title="") override
Set function title if title has the form "fffffff;xxxx;yyyy", it is assumed that the function title i...
void Draw(Option_t *option="") override
Draw this function with its current attributes.
virtual Double_t GetX(Double_t y, Double_t xmin=0, Double_t xmax=0, Double_t epsilon=1.E-10, Int_t maxiter=100, Bool_t logx=false) const
Returns the X value corresponding to the function value fy for (xmin<x<xmax).
TAxis * GetXaxis() const
Get x axis of the function.
A TGraph is an object made of two arrays X and Y with npoints each.
virtual void SetPoint(Int_t i, Double_t x, Double_t y)
Set x and y values for point number i.
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.
void Draw(Option_t *chopt="") override
Draw this graph with its current attributes.
@ kNoStats
Don't draw stats box.
virtual void SetContour(Int_t nlevels, const Double_t *levels=nullptr)
Set the number and values of contour levels.
2-D histogram with a double per channel (see TH1 documentation)
2-D histogram with a float per channel (see TH1 documentation)
Service class for 2-D histogram classes.
Use the TLine constructor to create a simple line.
virtual void SetTitle(const char *title="")
Set the title of the TNamed.
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 TClass * IsA() const
The most important graphics class in the ROOT system.
void ToLower()
Change string to lower-case.
const char * Data() const
TString & ReplaceAll(const TString &s1, const TString &s2)
TString & Append(const char *cs)
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
TVirtualPad is an abstract base class for the Pad and Canvas classes.
RooCmdArg Precision(double prec)
RooCmdArg LineColor(Color_t color)
double chisquared_quantile(double z, double r)
Inverse ( ) of the cumulative distribution function of the lower tail of the distribution with degr...
Namespace for the RooStats classes.