// This may look like C code, but it's really -*- C++ -*-

#if !defined(_Bolos2RING_H_SEEN_)
#define _Bolos2RING_H_SEEN_

#include <iostream.h>
#include <string>
#include "toiprocessor.h"
#include "spherehealpix.h"
#include "localmap.h"
#include "flagtoidef.h"

class Bolos2ring : public TOIProcessor
{
public:
  Bolos2ring(SphereHEALPix<r_8>* ringQ, 
	     SphereHEALPix<r_8>* ringU,
	     SphereHEALPix<r_8>* ringQW,
	     SphereHEALPix<r_8>* ringUW,
	     const vector<r_8>& table_angle,
	     int_4 *Bolos_OK,
	     int_4 wsz=8192);

  Bolos2ring(SphereHEALPix<r_8>* ringQ, 
	     SphereHEALPix<r_8>* ringU,
	     SphereHEALPix<r_8>* ringQW,
	     SphereHEALPix<r_8>* ringUW,
	     r_8 ang0, r_8 ang1, r_8 ang2, r_8 ang3,
	     int_4 *Bolos_OK,
	     int_4 wsz=8192);

  void init(int_4 wsz=8192);

  virtual ~Bolos2ring();
  inline void SetWSize(int_4 wsz) { wsize_ = (wsz<4) ? 4 : wsz; }
  inline int_8 ProcessedSampleCount() const { return totsncount_; }
  void SetNSide(int_4 nside);
  virtual void init();  virtual void run();
  
  virtual void PrintStatus(ostream& os);

protected:
  bool Bolos2ringdone;

  SphereHEALPix<r_8>* ringq;   // carte de Q
  SphereHEALPix<r_8>* ringu;   // carte de U
  

  SphereHEALPix<r_8>* ringqw;   // carte du nb de hit
  SphereHEALPix<r_8>* ringuw;
  vector<r_8> TableFP_;
  int_4* Bolos_OK;

  int wsize_;
  int_8 totsncount_;
  int_8 Npix_; 
  
};

#endif