Logo ROOT  
Reference Guide
RooCurve.h
Go to the documentation of this file.
1/*****************************************************************************
2 * Project: RooFit *
3 * Package: RooFitCore *
4 * File: $Id: RooCurve.h,v 1.24 2007/05/11 09:11:30 verkerke Exp $
5 * Authors: *
6 * WV, Wouter Verkerke, UC Santa Barbara, verkerke@slac.stanford.edu *
7 * DK, David Kirkby, UC Irvine, dkirkby@uci.edu *
8 * *
9 * Copyright (c) 2000-2005, Regents of the University of California *
10 * and Stanford University. All rights reserved. *
11 * *
12 * Redistribution and use in source and binary forms, *
13 * with or without modification, are permitted according to the terms *
14 * listed in LICENSE (http://roofit.sourceforge.net/license.txt) *
15 *****************************************************************************/
16#ifndef ROO_CURVE
17#define ROO_CURVE
18
19#include "TGraph.h"
20#include "RooPlotable.h"
21#include <list>
22#include <vector>
23#include "TMatrixDfwd.h"
24
25class RooAbsReal;
26class RooRealVar;
27class RooAbsFunc;
28class RooArgSet;
29class RooAbsRealLValue ;
30class RooHist ;
31
32class RooCurve : public TGraph, public RooPlotable {
33public:
34 RooCurve();
35 enum WingMode { NoWings=0 ,Straight=1, Extended=2 } ;
36 RooCurve(const RooAbsReal &func, RooAbsRealLValue &x, double xlo, double xhi, Int_t xbins,
37 double scaleFactor= 1, const RooArgSet *normVars= 0, double prec= 1e-3, double resolution= 1e-3,
38 bool shiftToZero=false, WingMode wmode=Extended, Int_t nEvalError=-1, Int_t doEEVal=false, double eeVal=0,
39 bool showProgress=false);
40 RooCurve(const char *name, const char *title, const RooAbsFunc &func, double xlo,
41 double xhi, UInt_t minPoints, double prec= 1e-3, double resolution= 1e-3,
42 bool shiftToZero=false, WingMode wmode=Extended, Int_t nEvalError=-1, Int_t doEEVal=false, double eeVal=0);
43 ~RooCurve() override;
44
45 RooCurve(const char* name, const char* title, const RooCurve& c1, const RooCurve& c2, double scale1=1., double scale2=1.) ;
46
47 void addPoint(double x, double y);
48
49 double getFitRangeBinW() const override;
50 double getFitRangeNEvt(double xlo, double xhi) const override ;
51 double getFitRangeNEvt() const override;
52
53
54 void printName(std::ostream& os) const override ;
55 void printTitle(std::ostream& os) const override ;
56 void printClassName(std::ostream& os) const override ;
57 void printMultiline(std::ostream& os, Int_t contents, bool verbose=false, TString indent="") const override;
58
59 inline void Print(Option_t *options= 0) const override {
60 // Printing interface
62 }
63
64 double chiSquare(const RooHist& hist, int nFitParam) const ;
65 Int_t findPoint(double value, double tolerance=1e-10) const ;
66 double average(double lo, double hi) const ;
67 double interpolate(double x, double tolerance=1e-10) const ;
68
69 bool isIdentical(const RooCurve& other, double tol=1e-6, bool verbose=true) const ;
70
71 RooCurve* makeErrorBand(const std::vector<RooCurve*>& variations, double Z=1) const ;
72 RooCurve* makeErrorBand(const std::vector<RooCurve*>& plusVar, const std::vector<RooCurve*>& minusVar, const TMatrixD& V, double Z=1) const ;
73
74protected:
75
76 void calcBandInterval(const std::vector<RooCurve*>& variations,Int_t i,double Z,double& lo, double& hi, bool approxGauss) const ;
77 void calcBandInterval(const std::vector<RooCurve*>& plusVar, const std::vector<RooCurve*>& minusVar, Int_t i, const TMatrixD& V,
78 double Z,double& lo, double& hi) const ;
79
80 void initialize();
81 void addPoints(const RooAbsFunc &func, double xlo, double xhi,
82 Int_t minPoints, double prec, double resolution, WingMode wmode,
83 Int_t numee=0, bool doEEVal=false, double eeVal=0.,std::list<double>* samplingHint=0) ;
84 void addRange(const RooAbsFunc& func, double x1, double x2, double y1,
85 double y2, double minDy, double minDx,
86 Int_t numee=0, bool doEEVal=false, double eeVal=0.) ;
87
88
89 void shiftCurveToZero(double prevYMax) ;
90
91 bool _showProgress ; ///<! Show progress indication when adding points
92
93 ClassDefOverride(RooCurve,1) // 1-dimensional smooth curve for use in RooPlots
94};
95
96#endif
#define e(i)
Definition: RSha256.hxx:103
const char Option_t
Definition: RtypesCore.h:66
#define ClassDefOverride(name, id)
Definition: Rtypes.h:339
static void indent(ostringstream &buf, int indent_level)
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 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 name[80]
Definition: TGX11.cxx:110
#define hi
Definition: THbookFile.cxx:128
Abstract interface for evaluating a real-valued function of one real variable and performing numerica...
Definition: RooAbsFunc.h:27
RooAbsRealLValue is the common abstract base class for objects that represent a real value that may a...
RooAbsReal is the common abstract base class for objects that represent a real value and implements f...
Definition: RooAbsReal.h:64
RooArgSet is a container object that can hold multiple RooAbsArg objects.
Definition: RooArgSet.h:57
A RooCurve is a one-dimensional graphical representation of a real-valued function.
Definition: RooCurve.h:32
void printTitle(std::ostream &os) const override
Print the title of this curve.
Definition: RooCurve.cxx:519
void initialize()
Perform initialization that is common to all curves.
Definition: RooCurve.cxx:257
~RooCurve() override
Destructor.
Definition: RooCurve.cxx:248
double getFitRangeBinW() const override
Get the bin width associated with this plotable object.
Definition: RooCurve.cxx:500
void Print(Option_t *options=0) const override
Print graph values.
Definition: RooCurve.h:59
void shiftCurveToZero(double prevYMax)
Find lowest point in curve and move all points in curve so that lowest point will go exactly through ...
Definition: RooCurve.cxx:271
void printName(std::ostream &os) const override
Print name of object.
Definition: RooCurve.cxx:508
void printMultiline(std::ostream &os, Int_t contents, bool verbose=false, TString indent="") const override
Print the details of this curve.
Definition: RooCurve.cxx:538
double interpolate(double x, double tolerance=1e-10) const
Return linearly interpolated value of curve at xvalue.
Definition: RooCurve.cxx:690
void addRange(const RooAbsFunc &func, double x1, double x2, double y1, double y2, double minDy, double minDx, Int_t numee=0, bool doEEVal=false, double eeVal=0.)
Fill the range (x1,x2) with points calculated using func(&x).
Definition: RooCurve.cxx:423
void printClassName(std::ostream &os) const override
Print the class name of this curve.
Definition: RooCurve.cxx:528
void addPoints(const RooAbsFunc &func, double xlo, double xhi, Int_t minPoints, double prec, double resolution, WingMode wmode, Int_t numee=0, bool doEEVal=false, double eeVal=0., std::list< double > *samplingHint=0)
Add points calculated with the specified function, over the range (xlo,xhi).
Definition: RooCurve.cxx:307
@ Extended
Definition: RooCurve.h:35
@ NoWings
Definition: RooCurve.h:35
@ Straight
Definition: RooCurve.h:35
RooCurve()
Default constructor.
Definition: RooCurve.cxx:68
bool _showProgress
! Show progress indication when adding points
Definition: RooCurve.h:91
RooCurve * makeErrorBand(const std::vector< RooCurve * > &variations, double Z=1) const
Construct filled RooCurve represented error band that captures alpha% of the variations of the curves...
Definition: RooCurve.cxx:737
double chiSquare(const RooHist &hist, int nFitParam) const
Calculate the chi^2/NDOF of this curve with respect to the histogram 'hist' accounting nFitParam floa...
Definition: RooCurve.cxx:556
void calcBandInterval(const std::vector< RooCurve * > &variations, Int_t i, double Z, double &lo, double &hi, bool approxGauss) const
Definition: RooCurve.cxx:847
double getFitRangeNEvt() const override
Return the number of events associated with the plotable object, it is always 1 for curves.
Definition: RooCurve.cxx:481
void addPoint(double x, double y)
Add a point with the specified coordinates. Update our y-axis limits.
Definition: RooCurve.cxx:468
bool isIdentical(const RooCurve &other, double tol=1e-6, bool verbose=true) const
Return true if curve is identical to other curve allowing for given absolute tolerance on each point ...
Definition: RooCurve.cxx:885
Int_t findPoint(double value, double tolerance=1e-10) const
Find the nearest point to xvalue.
Definition: RooCurve.cxx:668
double average(double lo, double hi) const
Return average curve value in [xFirst,xLast] by integrating curve between points and dividing by xLas...
Definition: RooCurve.cxx:604
A RooHist is a graphical representation of binned data based on the TGraphAsymmErrors class.
Definition: RooHist.h:27
Class RooPotable is a base class for objects that can be inserted into RooPlots and take advantage of...
Definition: RooPlotable.h:26
virtual StyleOption defaultPrintStyle(Option_t *opt) const
virtual Int_t defaultPrintContents(Option_t *opt) const
Default choice of contents to be printed (name and value)
static std::ostream & defaultPrintStream(std::ostream *os=0)
Return a reference to the current default stream to use in Print().
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,...
RooRealVar represents a variable that can be changed from the outside.
Definition: RooRealVar.h:40
A TGraph is an object made of two arrays X and Y with npoints each.
Definition: TGraph.h:41
Basic string class.
Definition: TString.h:136
Double_t y[n]
Definition: legend1.C:17
return c1
Definition: legend1.C:41
Double_t x[n]
Definition: legend1.C:17
return c2
Definition: legend2.C:14
const double xbins[xbins_n]