BCH1D.h

Go to the documentation of this file.

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

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

#ifndef __BCH1D__H
#define __BCH1D__H

#include <vector>

#include <TH1.h>

#include "BCLog.h"

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

class BCH1D
{
  
 public:
  
   /* @{ */ 

   BCH1D(); 

   ~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);

   void 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

---