BCH2D.h

#ifndef __BCH2D__H
#define __BCH2D__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 <string>

// ROOT classes

class TH1D;
class TH2D;
class TGraph;
class TObject;

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

class BCH2D
{

public:

   BCH2D(TH2D* h = 0);

   ~BCH2D();

   TH2D* GetHistogram()
   { return fHistogram; };

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

   void SetColorScheme(int scheme);

   void SetHistogram(TH2D* hist);

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

   void Print(const char* filename, std::string options="BTfB3CS1meangmode", 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="BTfB3CS1meangmodelmode", std::vector<double> intervals=std::vector<double>(0));

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

   void CalculateIntegratedHistogram();

   void PrintIntegratedHistogram(const char* filename);

   double GetLevel(double p);

   double GetArea(double p);

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

   TGraph* CalculateProfileGraph(int axis, std::string options="mode");

   TGraph* DrawProfile(int axis, std::string options, std::string drawoptions="blacksolid");

   TGraph* DrawProfileX(std::string options, std::string drawoptions)
   { return DrawProfile(0, options, drawoptions); };

   TGraph* DrawProfileY(std::string options, std::string drawoptions)
   { return DrawProfile(1, options, drawoptions); };

private:

   TH2D* fHistogram;

   TH1D* fIntegratedHistogram;

   double fMode[2];

   int fModeFlag;

   std::vector<int> fColors;

   std::vector<TObject*> fROOTObjects;

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

   static unsigned int fHCounter;
};

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

#endif