GloriaSkyCatalog.h

// -*- mode: c++; mode: auto-fill; mode: flyspell-prig; -*-


#ifndef GloriaSkyCatalog_h
#define GloriaSkyCatalog_h 1


/*
 *   This source code is part of the LUIZA software package for GLORIA.
 *   You are free to use this source files for your own development as
 *   long as it stays in a public research context. You are not
 *   allowed to use it for commercial purpose. You must put this
 *   header with author names in all development based on this file.
 *
 */


#include "GloriaFitsTable.h"

using namespace std;

namespace gloria{


  class GloriaSkyCatalog : public GloriaFitsTable {


  public:


    GloriaSkyCatalog();


    // Data access methods
    // ===================

    inline void AddObject(double ra, double dec) { AddRow(ra,dec); };


    inline void AddObject() { AddRow(); };


    inline int GetObjectNumber() {return _nRow;};

    inline double GetObjectRA( int iObject) 
          { return ((GloriaTableColumnOf<double> *)_columns.at(0))->At(iObject); };

    inline double GetObjectDEC( int iObject) 
             { return ((GloriaTableColumnOf<double> *)_columns.at(1))->At(iObject); };

    inline vector<double>& GetRAVector() 
             { return ((GloriaTableColumnOf<double> *)_columns.at(0))->Data(); };

    inline vector<double>& GetDECVector() 
             { return ((GloriaTableColumnOf<double> *)_columns.at(1))->Data(); };


    inline void SetObjectRA( int iObject, double ra) 
             { ((GloriaTableColumnOf<double> *)_columns.at(0))->At(iObject) = ra; };

    inline void SetObjectDEC( int iObject, double dec) 
             { ((GloriaTableColumnOf<double> *)_columns.at(1))->At(iObject) = dec; };

    bool  CheckCatalog();


    inline GloriaFitsTable& GetTable() { return *(( GloriaFitsTable *) this ); };


    int MatchObject(double ra, double dec, double cut);



    vector<long> MatchAllObjects(double ra, double dec, double distMax, double distMin=0.);



    long* MatchObjectArray(int nObject, double* RaArr, double* DecArr, double cut);

    

    vector<long> MatchObjectVector(vector<double>* RaVec, vector<double>* DecVec, double cut);

    

    vector<long> MatchObjectTable(GloriaFitsTable *table, double cut,
 vector<string> RaDecNames, vector<string> CopyNames, vector<string> NewNames );

    vector<long> MatchObjectTable(GloriaFitsTable *table, double cut);

    vector<long> MatchObjectTable(GloriaFitsTable *table, double cut, vector<string> RaDecNames);

  private:

  double  SkyDistance(double ra1, double dec1, double ra2, double dec2);
    
  } ;


}

#endif