BCH2D.h

Go to the documentation of this file.

#ifndef __BCH2D__H
#define __BCH2D__H

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

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

#include <vector>

// ROOT classes
class TH1D;
class TH2D;
class TGraph;

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

class BCH2D
{

   public:

      /* @{ */

      BCH2D();

      BCH2D(TH2D * h);

      ~BCH2D();

      /* @} */
      /* @{ */

      TH2D * GetHistogram()
         { return fHistogram; };

//    void GetMean(double& mean);

      void GetMode(double& mode);

      /* @} */
      /* @{ */

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

      void SetGlobalMode(double mode[2])
         { fMode[0] = mode[0]; fMode[1] = mode[1]; fModeFlag =1; };

      /* @} */
      /* @{ */

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

      void Draw(int options=0, bool drawmode=true);

      /*
       * Calculates the integral of the distribution as a function of the
       * height. */
      void CalculateIntegratedHistogram();

      /*
       * Calculates the height below which the integrated probability has
       * a certain value.
       * @param p The integrated probability in the region below the height to be estimated. */
      double GetLevel(double p);

      std::vector <int> GetNIntervalsY(TH2D * h, int &nfoundmax);

      TGraph * GetLowestBandGraph(TH2D * h, std::vector<int> nint);
      TGraph * GetLowestBandGraph(TH2D * h);


      std::vector <double> GetLevelBoundary(double level);
      std::vector <double> GetLevelBoundary(TH2D * h, double level);
      TGraph * GetBandGraph(double level1, double level2);
      TGraph * GetBandGraph(TH2D * h , double level1, double level2);

//    TGraph ** GetBandGraphs(TH2D * h);
      TGraph ** GetBandGraphs(TH2D * h, int &n);

      /* @} */

   private:

      TH2D * fHistogram;

      TH1D * fIntegratedHistogram;

      double fMode[2];

      int fModeFlag;

};

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

#endif

---