StdHist.h

/*
 * Hal::StdHist.h
 *
 *  Created on: 30-04-2022
 *      Author: Daniel Wielanek
 *      E-mail: daniel.wielanek@gmail.com
 *      Warsaw University of Technology, Faculty of Physics
 */
#ifndef HALSTDHIST_H_
#define HALSTDHIST_H_
 
#include <TString.h>
#include <TVector3.h>
 
class TH1;
class TH1D;
class TH2;
class TH2D;
class TH3;
class TH3D;
class TAxis;
class TVirtualPad;
 
namespace Hal {
  namespace Std {
    namespace {
      int anonymCounter = {0};
    }
    void RemoveNan(TH1* h, Double_t fill = 0, Double_t fill_e = 1);
    void CopyAxisProp(const TAxis* from, TAxis* to, TString opt = "");
    void CopyHistProp(const TH1& from, TH1& to, TString opt = "");
    TH1D* GetDiagonalProjection1D(TH3* h, TString dir, Double_t start, Double_t start2 = 0.);
    TH1D* GetProjection1D(const TH2* histo, Double_t min, Double_t max, Option_t* opt);
    TH1D* GetProjection1D(const TH3* histo, Double_t min1, Double_t max1, Double_t min2, Double_t max2, Option_t* opt);
    TH1D* SmearHistogram(TH1D* input_histogram, TH2D* smear_matrix, Option_t* opt = "");
    TH2D* GetProjection2D(const TH3* histo, Double_t min, Double_t max, Option_t* opt);
    void HistogramEdges(TH1* h, TString option = "x+uv", Double_t value = 0);
    void HistogramExtend(TH1* h, Char_t axis = 'x', Double_t factor = 1);
    Bool_t AreSimilar(const TH1* x, const TH1* y, Bool_t classes = kTRUE);
    void GetAxisPar(const TH1& obj, Int_t& nbins, Double_t& min, Double_t& xmax, Option_t* opt = "x");
    std::vector<double> GetAxisCenters(const TH1& obj, Option_t* opt = "x");
    TH1* MakeHisto1D(TString name, TString title, TVector3 Xaxis, Char_t type = 'D');
    TH2* MakeHisto2D(TString name, TString title, TVector3 Xaxis, TVector3 Yaxis, Char_t type = 'D');
    TH3* MakeHisto3D(TString name, TString title, TVector3 Xaxis, TVector3 Yaxis, TVector3 Zaxis, Char_t type = 'D');
    void MakeBeautiful();
    void Fold1D(Double_t val, TH1& h);
    void Fold2D(Double_t val, TH2& h, TString opt = "x");
    void Fold3D(Double_t val, TH3& h, TString opt = "x");
    void SetColor(TH1& h, Color_t color);
    void SetColorAndMarker(TH1& h, Color_t color, Marker_t m, Size_t s = 1);
    void SetRainbow(TH2& h, Double_t x1, Double_t y1, Double_t x2, Double_t y2);
    Int_t GetAntiColor(Int_t col);
    Int_t GetListOfSubPads(TVirtualPad* pad);
    TH1D* Crop1D(const TH1& h, Double_t min, Double_t max, TString opt = "vals");
    TH2D* Crop2D(const TH2& h, Double_t minX, Double_t maxX, Double_t minY, Double_t maxY, TString opt = "vals");
    TH3D* Crop3D(const TH3& h,
                 Double_t minX,
                 Double_t maxX,
                 Double_t minY,
                 Double_t maxY,
                 Double_t minZ,
                 Double_t maxZ,
                 TString opt = "vals");
    std::vector<TObject*> GetPadChildren(TString objName, TString className = "", TVirtualPad* pad = nullptr);
    Bool_t CheckHistogramData(const TH1& h1, const TH1& h2, Double_t thres = 0, Option_t* opt = "");
  }  // namespace Std
}  // namespace Hal
#endif /* HALSTDHIST_H_ */