UvwCoords.h

Go to the documentation of this file.

#include <vector>
#include <map>
using namespace std;

#include "Enum.hpp"
#include "CorrelationMode.h"

#include <measures/Measures.h>
#include <measures/Measures/MEpoch.h>
#include <measures/Measures/MPosition.h>
#include <measures/Measures/MDirection.h>
#include <measures/Measures/MCBaseline.h>
#include <measures/Measures/MBaseline.h>
#include <measures/Measures/MCuvw.h>
#include <measures/Measures/Muvw.h>
#include <casa/Quanta.h>

using namespace casa;

#include "ASDMEntities.h"
using namespace asdm;

class UvwCoords{
public:

  UvwCoords(ASDM* const datasetPtr);

  ~UvwCoords();

  void uvw_bl( Tag configDescriptionId, 
               const vector<vector<Angle> >& phaseDir,
               double timeCentroid,
               Enum<CorrelationMode> correlationMode,
               bool reverse, bool autoTrailing,
               vector<Vector<casa::Double> >& v_uvw );

  void uvw_bl( Tag configDescriptionId, 
               const vector<vector<Angle> >& phaseDir,
               const vector<double>& v_timeCentroid,
               Enum<CorrelationMode> correlationMode,
               bool reverse, bool autoTrailing, 
               vector<Vector<casa::Double> >& v_uvw);

  void uvw_bl( asdm::MainRow* mainRow, vector<pair<unsigned int,double> > v_tci, 
               Enum<CorrelationMode> correlationMode,
               pair<bool,bool> dataOrder, vector<Vector<casa::Double> >& v_uvw);


  void uvw_bl( asdm::MainRow* mainRow, vector<double> v_timeCentroid, 
               Enum<CorrelationMode> correlationMode,
               pair<bool,bool> dataOrder,
               vector<casa::Vector<casa::Double> >& v_uvw );

private:

  struct ArrayParam{
  public:
    Tag                   subarrayId;        //<! (sub)array identifier
    vector<Tag>           v_ant;             //<! sequence of antennas
    unsigned int          nrepeat;           //<! number of MS main table rows per antenna baseline
    Enum<CorrelationMode> e_correlationMode; //<! correlation mode (original mode passed through the user filter)
    bool operator== ( ArrayParam & rhs){
      set<Tag> a, b;
      for(unsigned int na=0; na<v_ant.size(); na++){
        a.insert(v_ant[na]);
      }

      for (unsigned int na=0; na<rhs.v_ant.size(); na++) {
        b.insert(rhs.v_ant[na]);
      }

      if( a != b ) return false;
      return true;
    }
    string show(){
      ostringstream os;
      for(unsigned int n=0; n<v_ant.size(); n++)os << v_ant[n].toString() << " ";
      os << " nrepeat " << nrepeat;
      return os.str();
    }
  };

  map<Tag,ArrayParam>            m_array_;
  map<Tag,MPosition>             m_antPos_;
  map<Tag,Vector<casa::Double> > m_antUVW_;

  // the 3 fundamental attributes of the state machine
  double                         timeCentroid_;
  vector<vector<Angle> >         phaseDir_; 
  Tag                            subarrayId_;

  Vector<casa::Double>           sduvw_;

  vector<double> antPos(const vector<double>& stationPos, const vector<double>& antOffset);

  void uvw_an(double timeCentroid, const vector<vector<Angle> >& phaseDir, const vector<Tag>& v_antId);


  void uvw_bl( const vector<Tag>& v_antennaId, unsigned int nrep, bool reverse, vector<Vector<casa::Double> >& v_uvw);

};