Logo ROOT   6.10/09
Reference Guide
RooFitResult.h
Go to the documentation of this file.
1 /*****************************************************************************
2  * Project: RooFit *
3  * Package: RooFitCore *
4  * File: $Id: RooFitResult.h,v 1.28 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_FIT_RESULT
17 #define ROO_FIT_RESULT
18 
19 #include "RooAbsArg.h"
20 #include "RooPrintable.h"
21 #include "RooDirItem.h"
22 #include "RooArgList.h"
23 
24 #include "RVersion.h"
25 #include "TMatrixFfwd.h"
26 #include "TMatrixDSym.h"
27 #include "TRootIOCtor.h"
28 
29 #include <vector>
30 #include <string>
31 #include <map>
32 
33 class RooArgSet ;
34 class RooAbsPdf ;
35 class RooPlot;
36 class TObject ;
37 class TH2 ;
38 typedef RooArgSet* pRooArgSet ;
39 
40 class RooFitResult : public TNamed, public RooPrintable, public RooDirItem {
41 public:
42 
43  // Constructors, assignment etc.
44  RooFitResult(const char* name=0, const char* title=0) ;
45  RooFitResult(const RooFitResult& other) ;
46  virtual TObject* Clone(const char* newname = 0) const {
47  RooFitResult* r = new RooFitResult(*this) ;
48  if (newname && *newname) r->SetName(newname) ;
49  return r ;
50  }
51  virtual TObject* clone() const { return new RooFitResult(*this); }
52  virtual ~RooFitResult() ;
53 
54  static RooFitResult* lastMinuitFit(const RooArgList& varList=RooArgList()) ;
55 
56  // Printing interface (human readable)
57  virtual void printValue(std::ostream& os) const ;
58  virtual void printName(std::ostream& os) const ;
59  virtual void printTitle(std::ostream& os) const ;
60  virtual void printClassName(std::ostream& os) const ;
61  virtual void printArgs(std::ostream& os) const ;
62  void printMultiline(std::ostream& os, Int_t contents, Bool_t verbose=kFALSE, TString indent="") const ;
63 
64  inline virtual void Print(Option_t *options= 0) const {
65  // Printing interface
67  }
68 
69  virtual Int_t defaultPrintContents(Option_t* opt) const ;
70  virtual StyleOption defaultPrintStyle(Option_t* opt) const ;
71 
72  RooAbsPdf* createHessePdf(const RooArgSet& params) const ;
73 
74  // Accessors
75  inline Int_t status() const {
76  // Return MINUIT status code
77  return _status ;
78  }
79 
80  inline UInt_t numStatusHistory() const { return _statusHistory.size() ; }
81  Int_t statusCodeHistory(UInt_t icycle) const ;
82  const char* statusLabelHistory(UInt_t icycle) const ;
83 
84  inline Int_t covQual() const {
85  // Return MINUIT quality code of covariance matrix
86  return _covQual ;
87  }
88  inline Int_t numInvalidNLL() const {
89  // Return number of NLL evaluations with problems
90  return _numBadNLL ;
91  }
92  inline Double_t edm() const {
93  // Return estimated distance to minimum
94  return _edm ;
95  }
96  inline Double_t minNll() const {
97  // Return minimized -log(L) value
98  return _minNLL ;
99  }
100  inline const RooArgList& constPars() const {
101  // Return list of constant parameters
102  return *_constPars ;
103  }
104  inline const RooArgList& floatParsInit() const {
105  // Return list of floating parameters before fit
106  return *_initPars ;
107  }
108  inline const RooArgList& floatParsFinal() const {
109  // Return list of floarting parameters after fit
110  return *_finalPars ;
111  }
112 
113  TH2* correlationHist(const char* name = "correlation_matrix") const ;
114 
115  Double_t correlation(const RooAbsArg& par1, const RooAbsArg& par2) const {
116  // Return correlation between par1 and par2
117  return correlation(par1.GetName(),par2.GetName()) ;
118  }
119  const RooArgList* correlation(const RooAbsArg& par) const {
120  // Return pointer to list of correlations of all parameters with par
121  return correlation(par.GetName()) ;
122  }
123 
124  Double_t correlation(const char* parname1, const char* parname2) const ;
125  const RooArgList* correlation(const char* parname) const ;
126 
127 
128  const TMatrixDSym& covarianceMatrix() const ;
129  const TMatrixDSym& correlationMatrix() const ;
130  TMatrixDSym reducedCovarianceMatrix(const RooArgList& params) const ;
131  TMatrixDSym conditionalCovarianceMatrix(const RooArgList& params) const ;
132 
133 
134  // Global correlation accessors
135  Double_t globalCorr(const RooAbsArg& par) { return globalCorr(par.GetName()) ; }
136  Double_t globalCorr(const char* parname) ;
137  const RooArgList* globalCorr() ;
138 
139 
140  // Add objects to a 2D plot
141  inline RooPlot *plotOn(RooPlot *frame, const RooAbsArg &par1, const RooAbsArg &par2,
142  const char *options= "ME") const {
143  // Plot error ellipse in par1 and par2 on frame
144  return plotOn(frame,par1.GetName(),par2.GetName(),options);
145  }
146  RooPlot *plotOn(RooPlot *plot, const char *parName1, const char *parName2,
147  const char *options= "ME") const;
148 
149  // Generate random perturbations of the final parameters using the covariance matrix
150  const RooArgList& randomizePars() const;
151 
152  Bool_t isIdentical(const RooFitResult& other, Double_t tol=5e-5, Double_t tolCorr=1e-4, Bool_t verbose=kTRUE) const ;
153 
154  void SetName(const char *name) ;
155  void SetNameTitle(const char *name, const char* title) ;
156 
157 protected:
158 
159  friend class RooMinuit ;
160  friend class RooMinimizer ;
161  void setCovarianceMatrix(TMatrixDSym& V) ;
162  void setConstParList(const RooArgList& list) ;
163  void setInitParList(const RooArgList& list) ;
164  void setFinalParList(const RooArgList& list) ;
165  inline void setMinNLL(Double_t val) { _minNLL = val ; }
166  inline void setEDM(Double_t val) { _edm = val ; }
167  inline void setStatus(Int_t val) { _status = val ; }
168  inline void setCovQual(Int_t val) { _covQual = val ; }
169  inline void setNumInvalidNLL(Int_t val) { _numBadNLL=val ; }
170  void fillCorrMatrix() ;
171  void fillCorrMatrix(const std::vector<double>& globalCC, const TMatrixDSym& corrs, const TMatrixDSym& covs) ;
172  void fillLegacyCorrMatrix() const ;
173  void setStatusHistory(std::vector<std::pair<std::string,int> >& hist) { _statusHistory = hist ; }
174 
175  Double_t correlation(Int_t row, Int_t col) const;
176  Double_t covariance(Int_t row, Int_t col) const;
177 
178  Int_t _status ; // MINUIT status code
179  Int_t _covQual ; // MINUIT quality code of covariance matrix
180  Int_t _numBadNLL ; // Number calls with bad (zero,negative) likelihood
181  Double_t _minNLL ; // NLL at minimum
182  Double_t _edm ; // Estimated distance to minimum
183  RooArgList* _constPars ; // List of constant parameters
184  RooArgList* _initPars ; // List of floating parameters with initial values
185  RooArgList* _finalPars ; // List of floating parameters with final values
186 
187  mutable RooArgList* _globalCorr ; //! List of global correlation coefficients
188  mutable TList _corrMatrix ; //! Correlation matrix (list of RooArgLists)
189 
190  mutable RooArgList *_randomPars; //! List of floating parameters with most recent random perturbation applied
191  mutable TMatrixF* _Lt; //! triangular matrix used for generate random perturbations
192 
193  TMatrixDSym* _CM ; // Correlation matrix
194  TMatrixDSym* _VM ; // Covariance matrix
195  TVectorD* _GC ; // Global correlation coefficients
196 
197  std::vector<std::pair<std::string,int> > _statusHistory ; // History of status codes
198 
199  ClassDef(RooFitResult,5) // Container class for fit result
200 };
201 
202 #endif
const TMatrixDSym & correlationMatrix() const
Return correlation matrix ;.
Double_t globalCorr(const RooAbsArg &par)
Definition: RooFitResult.h:135
const RooArgList * correlation(const RooAbsArg &par) const
Definition: RooFitResult.h:119
virtual const char * GetName() const
Returns name of object.
Definition: TNamed.h:47
double par[1]
Definition: unuranDistr.cxx:38
RooPlot * plotOn(RooPlot *frame, const RooAbsArg &par1, const RooAbsArg &par2, const char *options="ME") const
Definition: RooFitResult.h:141
Int_t numInvalidNLL() const
Definition: RooFitResult.h:88
Double_t _minNLL
Definition: RooFitResult.h:181
Double_t edm() const
Definition: RooFitResult.h:92
virtual void printArgs(std::ostream &os) const
Print arguments of fit result, i.e. the parameters of the fit.
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, which is interpreted as an OR of &#39;enum ContentsOptions&#39; values and in the style given by &#39;enum StyleOption&#39;.
TMatrixDSym conditionalCovarianceMatrix(const RooArgList &params) const
Return a reduced covariance matrix, which is calculated as ___ -1 Vred = V22 = V11 - V12 * V22 * V21...
const RooArgList & floatParsFinal() const
Definition: RooFitResult.h:108
const RooArgList & randomizePars() const
Return a list of floating parameter values that are perturbed from the final fit values by random amo...
virtual void printClassName(std::ostream &os) const
Print class name of fit result.
TMatrixDSym reducedCovarianceMatrix(const RooArgList &params) const
Return a reduced covariance matrix (Note that Vred is a simple sub-matrix of V, row/columns are order...
void setStatusHistory(std::vector< std::pair< std::string, int > > &hist)
Definition: RooFitResult.h:173
const char Option_t
Definition: RtypesCore.h:62
Double_t _edm
Definition: RooFitResult.h:182
const TMatrixDSym & covarianceMatrix() const
Return covariance matrix.
void fillLegacyCorrMatrix() const
Sanity check.
TMatrixDSym * _VM
Definition: RooFitResult.h:194
void SetNameTitle(const char *name, const char *title)
Change name and title of RooFitResult object.
TVectorT.
Definition: TMatrixTBase.h:77
void setNumInvalidNLL(Int_t val)
Definition: RooFitResult.h:169
void setStatus(Int_t val)
Definition: RooFitResult.h:167
Basic string class.
Definition: TString.h:129
RooArgList * _globalCorr
Definition: RooFitResult.h:187
int Int_t
Definition: RtypesCore.h:41
bool Bool_t
Definition: RtypesCore.h:59
virtual Int_t defaultPrintContents(Option_t *opt) const
Configure default contents to be printed.
void setInitParList(const RooArgList &list)
Fill the list of initial values of the floating parameters.
TVectorD * _GC
Definition: RooFitResult.h:195
virtual void Print(Option_t *options=0) const
Print TNamed name and title.
Definition: RooFitResult.h:64
virtual void printValue(std::ostream &os) const
Print the value of the fit result, i.e.g the status, minimized FCN, edm and covariance quality code...
TMatrixT.
Definition: TMatrixDfwd.h:22
void setEDM(Double_t val)
Definition: RooFitResult.h:166
virtual StyleOption defaultPrintStyle(Option_t *opt) const
Configure mapping of Print() arguments to RooPrintable print styles.
RooAbsPdf * createHessePdf(const RooArgSet &params) const
Return a p.d.f that represents the fit result as a multi-variate probability densisty function on the...
#define ClassDef(name, id)
Definition: Rtypes.h:297
RooPlotable is a &#39;mix-in&#39; base class that define the standard RooFit plotting and printing methods...
Definition: RooPrintable.h:25
const char * statusLabelHistory(UInt_t icycle) const
The TNamed class is the base class for all named ROOT classes.
Definition: TNamed.h:29
virtual TObject * Clone(const char *newname=0) const
Make a clone of an object using the Streamer facility.
Definition: RooFitResult.h:46
virtual TObject * clone() const
Definition: RooFitResult.h:51
TMatrixTSym.
RooArgList * _finalPars
Definition: RooFitResult.h:185
A doubly linked list.
Definition: TList.h:43
TH2 * correlationHist(const char *name="correlation_matrix") const
Return TH2D of correlation matrix.
void setMinNLL(Double_t val)
Definition: RooFitResult.h:165
TList _corrMatrix
List of global correlation coefficients.
Definition: RooFitResult.h:188
TMatrixDSym * _CM
triangular matrix used for generate random perturbations
Definition: RooFitResult.h:193
const double tol
TRandom2 r(17)
Service class for 2-Dim histogram classes.
Definition: TH2.h:30
Double_t covariance(Int_t row, Int_t col) const
Return the covariance matrix element addressed with numeric indices.
static RooFitResult * lastMinuitFit(const RooArgList &varList=RooArgList())
Import the results of the last fit performed by gMinuit, interpreting the fit parameters as the given...
Int_t statusCodeHistory(UInt_t icycle) const
unsigned int UInt_t
Definition: RtypesCore.h:42
Double_t minNll() const
Definition: RooFitResult.h:96
bool verbose
virtual ~RooFitResult()
Destructor.
void setFinalParList(const RooArgList &list)
Fill the list of final values of the floating parameters.
void printMultiline(std::ostream &os, Int_t contents, Bool_t verbose=kFALSE, TString indent="") const
Print fit result to stream &#39;os&#39;.
RooFitResult(const char *name=0, const char *title=0)
Constructor with name and title coverity[UNINIT_CTOR].
TMatrixF * _Lt
List of floating parameters with most recent random perturbation applied.
Definition: RooFitResult.h:191
void SetName(const char *name)
Change name of RooFitResult object.
RooArgList * _initPars
Definition: RooFitResult.h:184
A RooPlot is a plot frame and a container for graphics objects within that frame. ...
Definition: RooPlot.h:41
const Bool_t kFALSE
Definition: RtypesCore.h:92
static std::ostream & defaultPrintStream(std::ostream *os=0)
Return a reference to the current default stream to use in Print().
RooArgList * _randomPars
Correlation matrix (list of RooArgLists)
Definition: RooFitResult.h:190
void setConstParList(const RooArgList &list)
Fill the list of constant parameters.
double Double_t
Definition: RtypesCore.h:55
RooDirItem is a utility base class for RooFit objects that are to be attached to ROOT directories...
Definition: RooDirItem.h:22
virtual void printName(std::ostream &os) const
Print name of fit result.
you should not use this method at all Int_t Int_t Double_t Double_t Double_t e
Definition: TRolke.cxx:630
Int_t covQual() const
Definition: RooFitResult.h:84
Mother of all ROOT objects.
Definition: TObject.h:37
const RooArgList & constPars() const
Definition: RooFitResult.h:100
RooAbsPdf is the abstract interface for all probability density functions The class provides hybrid a...
Definition: RooAbsPdf.h:41
RooArgList * _constPars
Definition: RooFitResult.h:183
RooMinimizer is a wrapper class around ROOT::Fit:Fitter that provides a seamless interface between th...
Definition: RooMinimizer.h:38
UInt_t numStatusHistory() const
Definition: RooFitResult.h:80
const RooArgList * globalCorr()
Return the list of all global correlations.
const RooArgList & floatParsInit() const
Definition: RooFitResult.h:104
void setCovQual(Int_t val)
Definition: RooFitResult.h:168
RooArgSet * pRooArgSet
Definition: RooFitResult.h:37
RooAbsArg is the common abstract base class for objects that represent a value (of arbitrary type) an...
Definition: RooAbsArg.h:66
void fillCorrMatrix()
Internal utility method to extract the correlation matrix and the global correlation coefficients fro...
Int_t _numBadNLL
Definition: RooFitResult.h:180
const Bool_t kTRUE
Definition: RtypesCore.h:91
void setCovarianceMatrix(TMatrixDSym &V)
Store externally provided correlation matrix in this RooFitResult ;.
Double_t correlation(const RooAbsArg &par1, const RooAbsArg &par2) const
Definition: RooFitResult.h:115
Bool_t isIdentical(const RooFitResult &other, Double_t tol=5e-5, Double_t tolCorr=1e-4, Bool_t verbose=kTRUE) const
Return true if this fit result is identical to other within tolerance &#39;tol&#39; on fitted values and tole...
RooMinuit is a wrapper class around TFitter/TMinuit that provides a seamless interface between the MI...
Definition: RooMinuit.h:39
std::vector< std::pair< std::string, int > > _statusHistory
Definition: RooFitResult.h:197
Int_t status() const
Definition: RooFitResult.h:75
virtual void printTitle(std::ostream &os) const
Print title of fit result.