CGridder.hpp

/* -*- c++ -*- */

#ifndef __CGRIDDER_H__
#define __CGRIDDER_H__

#include "Matrix.h"
#include "Cvals.h"

//______________________________________________________________________________
class CGridder
{
public:
  CGridder(Cvals * xval, Cvals * yval, Cvals * binval,
           double xmin, double xmax, double ymin, double ymax,
           int istart, int iend);
  virtual ~CGridder(){}; 
  virtual int doGridding()=0; 
  Matrix<float> * getDataPtr(){ return _data; };

protected:
  void allocateMemory(int nx, int ny);

  int  getStartIndex () const { return _istart; }; 
  int  getEndIndex   () const { return _iend; };   
  void setNx(const int & nx) { _nx = nx; };
  void setNy(const int & ny) { _ny = ny; };

  Matrix<float> * _data;  
  int _nx;                
  int _ny;                
  Cvals * _xval;          
  Cvals * _yval;          
  Cvals * _binval;        
  double _xmin;
  double _xmax;
  double _ymin;
  double _ymax;

  int _ngridded;           
private:

  int _istart;             
  int _iend;               
};


//______________________________________________________________________________
class CProj3DKernel : public CGridder
{
public:
  CProj3DKernel(int nx, int ny, Cvals * xval, Cvals * yval, Cvals * binval,
                double xmin, double xmax, double ymin, double ymax,
                int istart, int iend);

  int doGridding();
  
private:
  static void  setupIntegratedKernel();
  static float wFromKernTab(const float & R2);

  static bool   _setup_int_kernel;          
  static int    _int_kernel_size;           
  static vector<float> _intkerntab;         
  static double _radkernel;                 /*<! max radius of the kernel (in h) (default: 2, r >= 2h -> kernel = 0) */
  static double _radkernel2;                /*<! radkernel * _radkernel */
  static double _dR2table;                  /*<! convert index to position */
  static double _ddR2table;                 /*<! inverse of _dR2table */

};

//______________________________________________________________________________
class CProj3DGridder : public CGridder
{
public:
  CProj3DGridder(int nx, int ny, Cvals * xval, Cvals * yval, Cvals * binval,
                 double xmin, double xmax, double ymin, double ymax,
                 int istart, int iend);
  ~CProj3DGridder();

  int doGridding();
  
  Matrix<int> * getNumberPtr(){ return _number; };

private:
  Matrix<int> * _number;
};

//______________________________________________________________________________
class CParticleIDFinder : public CGridder
{
public:
  CParticleIDFinder(Cvals * xval, Cvals * yval, Cvals * binval,
                    double xmin, double xmax, double ymin, double ymax,
                    int istart, int iend);

  int doGridding();
};


#endif // __CGRIDDER_H__

---