BCH1D.h

#ifndef __BCH1D__H
#define __BCH1D__H

/*
 * Copyright (C) 2007-2014, the BAT core developer team
 * All rights reserved.
 *
 * For the licensing terms see doc/COPYING.
 * For documentation see http://mpp.mpg.de/bat
 */

// ---------------------------------------------------------

#include <vector>

#include <TH1.h>

// ---------------------------------------------------------

class BCH1D
{
 public:

  BCH1D(TH1D * hist = 0);

  ~BCH1D();

  TH1D* GetHistogram()
    { return fHistogram; };

  double GetMean()
  { return fHistogram -> GetMean(); };

  double GetMode();

  double GetMedian()
  { return this -> GetQuantile(0.5); };

  double GetQuantile(double probablity);

  double GetLimit(double probability)
  { return this->GetQuantile(probability); };

  double GetRMS()
  { return fHistogram->GetRMS(); };

  double GetSTD()
  { return fHistogram->GetRMS(); };

  double GetVariance()
  { return (GetSTD()*GetSTD()); };

  double GetSkew()
  { return fHistogram->GetSkewness(); };

  double GetKurtosis()
  { return fHistogram->GetKurtosis(); };

  double GetIntegral(double valuemin, double valuemax);

  double GetPValue(double probability);

  int GetColor(int index)
  { return fColors.at(index); };

  void SetColorScheme(int scheme);

  void SetHistogram(TH1D * hist)
  { fHistogram = hist; };

  void SetDefaultCLLimit(double limit);

  void SetGlobalMode(double mode)
  { fMode=mode;
    fModeFlag=1; };

  void Print(const char * filename, std::string options="BTsiB3CS1D0Lmeanmode", std::vector<double> intervals=std::vector<double>(0), int ww=0, int wh=0);

  void Print(const char * filename, std::string options, double interval, int ww=0, int wh=0);

  void Draw(std::string options="BTsiB3CS1D0Lmeanmode", std::vector<double> intervals=std::vector<double>(0));

  void Draw(std::string options, double interval);

  void DrawDelta(double value) const;

  double GetSmallestInterval(double & min, double & max, double content=.68);

  TH1D* GetSmallestIntervalHistogram(double level);

  std::vector<double> GetSmallestIntervals(double content = 0.68);

  double IntegralWidth(double min, double max);

  TH1D* GetSubHisto(double min, double max, const char * name);

 private:

  TH1D* fHistogram;

  double fDefaultCLLimit;

  double fMode;

  int fModeFlag;

  std::vector<int> fColors;

  mutable std::vector<TObject*> fROOTObjects;

  static unsigned int getNextIndex()
    { return ++fHCounter; }

  static unsigned int fHCounter;
};

// ---------------------------------------------------------

#endif