source: trunk/src/histogram.cc @ 13

#include "histogram.h"


using namespace std;
using namespace thep_cpp_tools;


// Constructors and destructors
////////////////////////////////
histogram::histogram( const size_t& nbins, const double& min, const double& max )
{
  h_ = gsl_histogram_alloc( nbins );
  gsl_histogram_set_ranges_uniform( h_, min, max );
}


histogram::~histogram()
{
  gsl_histogram_free( h_ );
  if( p_ != NULL ) gsl_histogram_pdf_free( p_ );
}


// Functions (A-Z)
//////////////////
void histogram::init_probabilities()
{
  p_ = gsl_histogram_pdf_alloc( get_bins() );
  gsl_histogram_pdf_init( p_, h_ );
}


void histogram::load_from_text_file( const string& filename )
{
  ifstream fin( filename.c_str() );
  if( !fin ) { cerr << "No such file: " << filename << endl; exit( -1 ); }
  while( fin.peek() != EOF )
    {
     double val, intensity;
     fin >> val >> intensity;
     assert( intensity >= 0 );
     gsl_histogram_accumulate( h_, val, intensity );//intensity );
    }
}


// Operators
////////////
std::ostream& 
thep_cpp_tools::operator<< ( std::ostream& s_out, 
           const thep_cpp_tools::histogram& hist )
{
  for( size_t i = 0; i < hist.get_bins(); ++i )
    {
     double low, high;
     cout << i << endl;
     hist.get_range( i, low, high );
     s_out << "[ " << low << ", " << high << " ]\t=\t" 
     << hist.get_intensity( i ) << endl;
    }
  return s_out;
}