GalacticCoords.hxx

Go to the documentation of this file.

//-----------------------------------------------------------------------------
// CVS: $Revision: 7503 $, $Date: 2010-09-07 16:34:15 +0100 (Tue, 07 Sep 2010) $, $Author: EckhardSutorius $
//-----------------------------------------------------------------------------
#ifndef GALACTICCOORDS_H
#define GALACTICCOORDS_H

extern "C"
{
  #include "slalib.h"
}
#include <limits>

#include "Constants.h"
#include "Converter.h"
#include "DataMethod.hxx"
#include "Logger.h"
#include "StringOps.h"
//-----------------------------------------------------------------------------
template<typename DataType>
class GalacticCoords : public DataMethod<DataType>
{
public:
  GalacticCoords() { DataMethod<DataType>::methodName = "GalacticCoords"; }

  virtual ~GalacticCoords() { }

  //---------------------------------------------------------------------------
  void setSource(IntMap& tbls, IntMap& cols, StringMap units)
  {
    raColno  = cols[0];
    decColno = cols[1];

    // Require radians for the calculation
    mDecFact = converter.getFactor(units[0], "RADIANS");
    mRAFact  = converter.getFactor(units[1], "RADIANS");
  }
  //---------------------------------------------------------------------------
  void setTarget(IntMap& cols, StringMap units)
  {
    lColno = cols[0];
    bColno = cols[1];

    // Check required output units for (l,b)
    mlFact = converter.getFactor("RADIANS", units[0]);
    mbFact = converter.getFactor("RADIANS", units[1]);
  }
  //---------------------------------------------------------------------------
  void doit(TableData<DataType>& data, int row1, int row2)
  {
    Logger log;
    for (int rowNo = row1; rowNo < row2; ++rowNo)
    {
      // Input RA and DEC
      double ra, dec;
      data.value(raColno, rowNo, ra);
      data.value(decColno, rowNo, dec);

      double dl, db;
      // Get l and b (2000J) from slalib
      cslaEqgal(ra*mRAFact, dec*mDecFact, &dl, &db);

      // Convert to required units
      dl *= mlFact;
      db *= mbFact;

      // Assign to correct columns
      if (dl != dl  or std::isinf(dl) or
          db != db or std::isinf(db))
      {
        string mess = "WCS error in l, b ";
        mess += StringOps::NumToString(lColno) + ",";
        mess += StringOps::NumToString(bColno) + " ";
        mess += StringOps::NumToString(rowNo) + " : ";
        mess += "value replaced with dbl_max";
        log.addError(mess);
        data.assign(lColno, rowNo, std::numeric_limits<double>::max( )-1);
        data.assign(bColno, rowNo, std::numeric_limits<double>::max( )-1);
      } else {
        data.assign(lColno, rowNo, dl);
        data.assign(bColno, rowNo, db);
      }
    }
  }

  // No data from raw data included
  void doit2(TableData<DataType>& rawData, TableData<DataType>& data, int row1, int row2) { }
  //---------------------------------------------------------------------------
private:
  int lColno;
  int bColno;
  int raColno;
  int decColno;
  double mDecFact;
  double mRAFact;
  double mlFact;
  double mbFact;
  Converter converter;
};
//-----------------------------------------------------------------------------
#endif
//-----------------------------------------------------------------------------
// Change log:
//
// 27-May-2004,  IAB: Original version.
// 22-Mar-2007,  RSC: Robustified against dodgy WCS values.
//  2-May-2007,  RSC: Reverted to original design of allowing dodgy WCS values.
//  7-Apr-2008, ETWS: Upgraded to use new detection table layout.