source: Sophya/trunk/Cosmo/RadioBeam/qhist.h@ 4050

//  Simple specific histo-like classes 
//  R. Ansari - Avril-Mai 2010 

#ifndef QHIST_SEEN
#define QHIST_SEEN

#include "histats.h"       // SOPHYA .h 

#include "machdefs.h"      // SOPHYA .h  
#include "sopnamsp.h"      // SOPHYA .h  
#include <iostream>

#include "histats.h"       // SOPHYA .h 

// -----------------------------------
// -- Classe ressemblant a un histo 1D 
class QHis1D {
  QHis1D();
  QHis1D(QHis1D const& a);
  QHis1D(r_8 xMin,r_8 xMax,int_4 nxBin);
  void Define(r_8 xMin,r_8 xMax,int_4 nxBin); 
  sa_size_t Add(r_8 x, r_8 w=1.);
  inline double WBinX() { return dxb_; }
  inline double NBinX() { return nx_; }

  Histo Convert();
  
  r_8 xmin_,xmax_;
  r_8 dxb_;
  sa_size_t nx_;
  TArray<r_8> aw_;
  sa_size_t nhis_, nunder_, nover_;
  r_8 swunder_, swover_ ; 
};

// -----------------------------------
// -- Classe ressemblant a un histo 2D 
class QHis2D {
public:
  QHis2D();
  QHis2D(r_8 xMin,r_8 xMax,int_4 nxBin,r_8 yMin,r_8 yMax,int_4 nyBin);
  void Define(r_8 xMin,r_8 xMax,int_4 nxBin,r_8 yMin,r_8 yMax,int_4 nyBin);
  double Add(r_8 x, r_8 y, r_8 w, bool fgfh);
  void SetZeroBin(r_8 x=0, r_8 y=0);
  inline r_8 WBinX() { return dxb; }
  inline r_8 WBinY() { return dyb; }
  inline sa_size_t NBinX() { return nx; }
  inline sa_size_t NBinY() { return ny; }
  inline r_8 X(sa_size_t i) { return (dxb*((double)i+0.5)+xmin); }
  inline r_8 Y(sa_size_t j) { return (dyb*((double)j+0.5)+ymin); }
  inline r_8 operator()(sa_size_t i, sa_size_t j) { return aw(i,j); }

  inline void GetMinMax(r_8& min, r_8& max) { return aw.MinMax(min,max); }

  inline double SumWBinZero() { return sumw0; }
  void Renormalize(r_8 mfac=1.) { aw*=mfac; }

  Histo2D Convert();
  
protected:
  r_8 xmin,xmax,ymin,ymax;
  r_8 dxb,dyb;
  sa_size_t nx,ny;
  sa_size_t ixb0, jyb0;
  TArray<r_8> aw;
  double sumw0;
};


#endif