Histo1D.cpp

Go to the documentation of this file.

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

// include files
#include <math.h>

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

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

Histo1D::Histo1D( TH1D *histo ) : m_zero ( 1.e-30 ) , m_infini(1.e30) {
  //m_histo = histo; 
  
  m_nbins = histo->GetNbinsX();

  setTitle(std::string(histo -> GetTitle()));

  m_bins.resize (m_nbins,0.);
  m_centers.resize (m_nbins,0.);
  m_edges.resize (m_nbins+1,0.);

  unsigned long nbin=0;
  for (nbin=0;nbin<m_nbins; ++nbin){
    m_bins[nbin]=histo->GetBinContent(nbin+1);
    m_centers[nbin]=histo->GetBinCenter(nbin+1);
  }

  m_delta = histo->GetBinWidth(0);

  for (nbin=0;nbin<m_nbins+1; ++nbin){
    m_edges[nbin]=histo->GetBinLowEdge(1)+nbin*m_delta;
  }

  m_underflow = histo->GetBinContent(0);
  m_overflow = histo->GetBinContent(histo->GetNbinsX()+1);  

  m_minX = histo->GetBinLowEdge(1);
  m_maxX = histo->GetBinLowEdge(histo->GetNbinsX()+1);

  m_minY = histo->GetMinimum();
  m_maxY = histo->GetMaximum();

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

double Histo1D::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];
}   

double Histo1D::minY(){
    return m_minY;
}

double Histo1D::maxY(){
  return m_maxY;
}

double Histo1D::binEdgeMin( unsigned int bin ){
  return m_minX+bin*(m_delta);
}

double Histo1D::binEdgeMax ( unsigned int bin ){
  return m_minX+(bin+1)*(m_delta);
}

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