rarNLL.hh

Go to the documentation of this file.

/*****************************************************************************
* Project: BaBar detector at the SLAC PEP-II B-factory
* Package: RooRarFit
 *    File: $Id: rarNLL.rdl,v 1.7 2011/06/16 13:18:50 fwilson Exp $
 * Authors: Lei Zhang
 * History:
 *
 * Copyright (C) 2005-2012, University of California, Riverside
 *****************************************************************************/
#ifndef RAR_NLL
#define RAR_NLL

#include "TList.h"
#include "TMatrixD.h"
#include "TString.h"
#include "TObject.h"

#include "RooFitCore/RooCurve.hh"

class rarNLL : public TNamed {
  
public:
  rarNLL();
  rarNLL(RooCurve *curve,
         const char *name="theRARNLL", const char *title="the RARNLL",
         const Bool_t verbose=kFALSE);
  virtual ~rarNLL();
  
  virtual void setVerbose(Bool_t verbose=kTRUE) {_verbose=verbose;}
  
  virtual void init(RooCurve *curve=0);
  Double_t getNLL(Double_t x);
  Double_t getY(Double_t x);
  TArrayD getX(Double_t y);
  void getMin(Double_t &x, Double_t &y);
  TArrayD getMin();
  Double_t getLIntegral(Double_t x);
  Double_t getLIntegral(Double_t xl, Double_t xh);
  Double_t getLIntegral();
  Double_t getLIntegralInverse(Double_t xl, Double_t iVal);
  Double_t getLIntegralInverse(Double_t iVal);
  Double_t lIntegralFunc(Double_t x, TMatrixD &A, TMatrixD &lA);
  Double_t lIntegralFuncInverse(Double_t x0, Double_t x2, Int_t iter,
                                Double_t &thisXI,
                                Double_t &la, Double_t &lb, Double_t &lc);

protected:
  void getMin(TArrayD &xy, Double_t x, Double_t a, Double_t b, Double_t c);
  TArrayD getMin(Double_t x0, Double_t x1, Double_t x2,
                 Double_t a, Double_t b, Double_t c);
  
  Int_t _nPoints; 
  TArrayD _xs; 
  TArrayD _ys; 
  Int_t _mIdx; 
  Int_t _nSteps; 
  TList _coeffMList; 
  TList _lCoefMList; 
  TArrayD _x0s; 
  TArrayD _x1s; 
  TArrayD _x2s; 
  TArrayD _iLIntegrals; 
  TArrayD _tLIntegrals; 
  
  Bool_t _verbose; 
  
private:
  rarNLL(const rarNLL&);
  ClassDef(rarNLL, 0) // RooRarFit NLL class
    ;
};

#endif

---