BCH1D.h

Go to the documentation of this file.

#ifndef __BCH1D__H
#define __BCH1D__H

/*
 * Copyright (C) 2008, Daniel Kollar and Kevin Kroeninger.
 * All rights reserved.
 *
 * For the licensing terms see doc/COPYING.
 */

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

#include <vector>

#include <TH1.h>

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

class BCH1D
{
   public:

      /* @{ */

      BCH1D();

      BCH1D(TH1D * hist)
         { fHistogram = hist; };

      ~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 GetIntegral(double valuemin, double valuemax);

      double GetPValue(double probability);

      /* @} */

      /* @{ */

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

      void SetDefaultCLLimit(double limit);

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

      /* @} */

      /* @{ */

      void Print(const char * filename, int options=0, double ovalue=0., int ww=0, int wh=0);

      void Draw(int options=0, double ovalue=0.);

      void DrawShadedLimits(double mode, double min, double max, double limit=0);

      void DrawSmallest(double mode, double prob, bool drawmean=true);

      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;

};

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

#endif

---