Cartesian.hxx

Go to the documentation of this file.

//-----------------------------------------------------------------------------
// CVS: $Revision: 4695 $, $Date: 2008-04-09 16:18:42 +0100 (Wed, 09 Apr 2008) $, $Author: EckhardSutorius $
//-----------------------------------------------------------------------------
#ifndef CARTESIAN_H
#define CARTESIAN_H

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

  virtual ~Cartesian() {}

  void setSource(IntMap& tbls, IntMap& cols, StringMap units)
  {
    raColno  = cols[0];
    decColno = cols[1];
    Converter converter;
    mDecFact = converter.getFactor(units[0], "RADIANS");
    mRAFact  = converter.getFactor(units[1], "RADIANS");
  }

  void setTarget(IntMap& cols, StringMap units)
  {
    cxColno = cols[0];
    cyColno = cols[1];
    czColno = cols[2];
  }

  void doit(TableData<DataType>& data, int row1, int row2)
  {
    for (int rowNo = row1; rowNo < row2; ++rowNo)
    {
      // Input RA and DEC
      double ra, dec;
      data.value(raColno, rowNo, ra);
      data.value(decColno, rowNo, dec);

      // Cartesian components
      double cx, cy, cz;
      DataOps::calcCartesian(mRAFact * ra, mDecFact * dec, cx, cy, cz);

      // Assign to correct columns
      data.assign(cxColno, rowNo, cx);
      data.assign(cyColno, rowNo, cy);
      data.assign(czColno, rowNo, cz);
    }
  }
  // No data from raw data included
  void doit2(TableData<DataType>& rawData, TableData<DataType>& data, int row1, int row2) { }

private:
  int cxColno;
  int cyColno;
  int czColno;
  int raColno;
  int decColno;
  double mDecFact;
  double mRAFact;
};
//-----------------------------------------------------------------------------
#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.

---