Histo2D.cpp

Go to the documentation of this file.

//$Id: Histo2D.cpp,v 1.3 2006/03/16 17:08:23 fmachefe Exp $
//------------------------------------------------------------------------------
//
//  Package    : Histo2D
//
//  Description:
//
//  Author(s)  : F. Machefert -
//  Date       : 5 March 2004
//
//------------------------------------------------------------------------------

// include files
#include <math.h>

// CATCore include Files
#include "Tools.h"

// local include file
#include "Histo2D.h"

Histo2D::Histo2D( TH2D *histo ) : m_zero ( 1.e-30 ) , m_infini(1.e30) {
  //m_histo = histo; 
  
  m_nxbins = histo->GetNbinsX();
  m_nybins = histo->GetNbinsY();

  setTitle(std::string(histo -> GetTitle()));
  
  m_minX = histo->GetXaxis()->GetXmin();
  m_maxX = histo->GetXaxis()->GetXmax();
  m_minY = histo->GetYaxis()->GetXmin();
  m_maxY = histo->GetYaxis()->GetXmax();

  m_minZ = histo->GetMinimum();
  m_maxZ = histo->GetMaximum();

  unsigned long nxbin=0;
  unsigned long nybin=0;
  for (nybin=0;nxbin<m_nybins; ++nybin){
    /*
      std::vector<double> tmpbins;
      std::vector<double> tmpxbin;
      std::vector<double> tmpybin;
    */
    for (nxbin=0;nxbin<m_nxbins; ++nxbin){
      m_bins.push_back(histo->GetBinContent(nybin*nxbin+1+2*(nybin-1)));
      m_xcenters.push_back(histo->GetXaxis()->GetBinCenter(nxbin+1));
      m_ycenters.push_back(histo->GetYaxis()->GetBinCenter(nybin+1));
    }
    /*
      m_bins.push_back(tmpbins);
      m_xcenters.push_back(tmpxbin);
      m_ycenters.push_back(tmpybin);    
    */
  }
  
  /*  
  m_underflow = histo->GetBinContent(0);
  m_overflow = histo->GetBinContent(histo->GetNbinsX()+1);  
  */

  m_content = histo->GetEntries();
  m_mean    = histo->GetMean();
  m_rms     = histo->GetRMS();
}

/*
double Histo2D::bin(double x){
    if (x>=m_maxX) {return m_overflow;}
    else if(x<m_minX) {return m_underflow;}
    unsigned long bin=(unsigned long)(floor(m_nbins*(x-m_minX)/m_delta));
    return m_bins[bin];
}   
*/

void Histo2D::statistics(double &content, 
                         double &mean, 
                         double &rms){
  content=m_content;
  mean   =m_mean;
  rms    =m_rms;
}