ChiSqMap2D.h

/*
 * HalChiSqMap.h
 *
 *  Created on: 30-04-2022
 *      Author: Daniel Wielanek
 *      E-mail: daniel.wielanek@gmail.com
 *      Warsaw University of Technology, Faculty of Physics
 */
#ifndef HALCHISQMAP2D_H_
#define HALCHISQMAP2D_H_
 
#include "Object.h"
 
#include <TString.h>
 
 
class TH2D;
class TLine;
namespace Hal {
  class ChiSqMap2D : public Object {
    TH2D* fHist;
    TLine* fLineX;
    TLine* fLineY;
    TLine* fLineXlow;
    TLine* fLineXhigh;
    TLine* fLineYlow;
    TLine* fLineYhigh;
    TString fLabel;
    void DrawLineX(TLine* l);
    void DrawLineY(TLine* l);
    Bool_t CheckX(Double_t x) const;
    Bool_t CheckY(Double_t y) const;
    TLine* MakeXLine(Double_t x, Style_t style) const;
    TLine* MakeYLine(Double_t x, Style_t style) const;
 
  protected:
    ChiSqMap2D();
 
  public:
    ChiSqMap2D(TString name, Int_t xbins, Double_t xlow, Double_t xhigh, Int_t ybins, Double_t ylow, Double_t yhigh);
    void SetParNames(TString xpar, TString ypar);
    void SetLabel(TString label) { fLabel = label; };
    TString GetLabel() const { return fLabel; };
    static TString ChiLabel() { return "#chi^{2}"; };
    static TString ChiNDFLabel() { return "#chi^{2}/NDF"; };
    Double_t GetFitX() const;
    Double_t GetFitY() const;
    Double_t GetFitXErrorUp() const;
    Double_t GetFitXErrorLow() const;
    Double_t GetFitYErrorUp() const;
    Double_t GetFitYErrorLow() const;
    Double_t GetEstErrorX(Double_t thres = 0.1, Bool_t around_fit = kTRUE) const;
    Double_t GetEstErrorY(Double_t thres = 0.1, Bool_t around_fit = kTRUE) const;
    Double_t GetEstX() const;
    Double_t GetEstY() const;
 
    TH2D* GetHist() const { return fHist; };
    void SetXFit(Double_t Value, Double_t errLow = 0, Double_t errHigh = -1);
    void SetYFit(Double_t Value, Double_t errLow = 0, Double_t errHigh = -1);
    void GetMin(Double_t& x, Double_t& y) const;
    void Draw(Option_t* opt = "");
    void GetXYE(Double_t& x, Double_t& eX, Double_t& y, Double_t& eY) const;
    virtual ~ChiSqMap2D();
    ClassDef(ChiSqMap2D, 1)
  };
}  // namespace Hal
#endif /* HALCHISQMAP2D_H_ */