GslHistogram.hh

/*
$Id: GslHistogram.hh 2001 2009-03-16 06:28:16Z jisuzuki $
 */

#ifndef GSLHISTOGRAM
#define GSLHISTOGRAM

#include "Header.hh"
#include <gsl/gsl_histogram.h>
#include "CppToPython.hh"


class GslHistogram {
private:
  gsl_histogram **h;
  void Set( Double *Bin, UInt4 Size );
  PyObject *VtoL( vector<Double> v );
  UInt4 NumOfMulTh;
public:

  GslHistogram( Double *Bin, UInt4 Size, UInt4 _NumOfMulTh=1 );
 ~GslHistogram();
  UInt4 PutHistSize(){ return gsl_histogram_bins( h[0] ); }
  PyObject *PutHistogramList(){ return VtoL( PutHistogram() );}
  vector<Double> PutHistogram();
  vector<Double> PutHistogramErr();
  PyObject *PutHistogramErrList(){ return VtoL( PutHistogramErr() ); }
  void Increment( Double x, UInt4 ThNum=0, Double Weight=1.0 );
  void Increment( vector<Double> x, UInt4 ThNum=0 );
  void Increment( vector<Int4> x, UInt4 ThNum=0 );
  void Dump(){ gsl_histogram_fprintf( stdout, h[0], "%g", "%g"); }
  //void Reset(){ gsl_histogram_reset( h ); }
  void Reset();
  gsl_histogram *PutHistPointer( UInt4 ThNum=0 ){ return h[ThNum]; }
  vector<Double> PutBin();
  PyObject *PutBinList(){ return VtoL( PutBin() ); }
  Double PutHistogramSum();
};

#endif