source: trunk/src/Histogram.cc @ 231

// $Id: Histogram.cc 231 2005-02-21 14:50:17Z peter $

#include "Histogram.h"

#include <cmath>
#include <fstream>


namespace theplu {
namespace cpptools {


Histogram::Histogram(void)
  : xmax_(0), xmin_(0), sum_all_(), sum_histogram_()
{
}



Histogram::Histogram(const Histogram& b)
{
  *this=b;
}



Histogram::Histogram(const double min, const double max, const size_t n)
  : histogram_(std::vector<double>(n,0.0)),
    xmax_(max), xmin_(min),
    sum_all_(), sum_histogram_()
{
}



Histogram::~Histogram(void)
{
}



int Histogram::add(const double x, const double w)
{
  sum_all_.add(x,w);
  if (x<xmin_)
    return -1;
  else if (x>=xmax_)
    return 1;
  
  sum_histogram_.add(x,w);
  histogram_[bin(x)] += w;
  return 0;
}



void Histogram::normalize(double scale_factor)
{
  for (size_t i=0; i<histogram_.size(); i++)
    histogram_[i]/=scale_factor;
}


void Histogram::reset(void)
{
  for (u_int i=0; i<histogram_.size(); i++)
    histogram_[i]=0;
  sum_all_.reset();
  sum_histogram_.reset();
}



const Histogram& Histogram::operator=(const Histogram& b)
{
  if (this==&b)
    return *this;
  histogram_=b.histogram_;
  xmax_=b.xmax_;
  xmin_=b.xmin_;
  sum_all_=b.sum_all_;
  sum_histogram_=b.sum_histogram_;
  return *this;
}



std::ostream& operator<<(std::ostream& s,const Histogram& histogram)
{
  s << "# histogram min : " << histogram.xmin() << '\n';
  s << "# histogram max : " << histogram.xmax() << '\n';
  s << "# number of bins: " << histogram.nof_bins() << '\n';
  s << "# nof points in histogram : " 
    << histogram.averager_histogram().sum_w() << '\n';
  s << "# nof points in total:      " 
    << histogram.averager_all().sum_w() << '\n';
  s << "# column 1: center of observation bin\n"
    << "# column 2: frequency\n";

  for (u_int i=0; i<histogram.nof_bins(); i++) {
    s.width(12);
    s << histogram.observation_value(i);
    s.width(12);
    s << histogram[i] << '\n';
  }

  return s;
}

}} // of namespace cpptools and namespace theplu