UvwCoords.h

Go to the documentation of this file.

#include <vector>
#include <map>

#include <alma/Enumtcl/Enum.hpp>
#include <alma/Enumtcl/CorrelationMode.h>

#include <ms/MeasurementSets/MSFieldColumns.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>

#include <alma/ASDM/ASDMEntities.h>

class UvwCoords{
public:

    UvwCoords(asdm::ASDM* const datasetPtr);

    ~UvwCoords();

    void uvw_bl( asdm::Tag configDescriptionId, 
                 const std::vector<std::vector<asdm::Angle> >& phaseDir,
                 double timeCentroid,
                 Enum<CorrelationModeMod::CorrelationMode> correlationMode,
                 bool reverse, bool autoTrailing,
                 std::vector<casacore::Vector<casacore::Double> >& v_uvw,
                 casacore::MSFieldColumns* msfc_p=0);

    void uvw_bl( asdm::Tag configDescriptionId, 
                 const std::vector<std::vector<asdm::Angle> >& phaseDir,
                 const std::vector<double>& v_timeCentroid,
                 Enum<CorrelationModeMod::CorrelationMode> correlationMode,
                 bool reverse, bool autoTrailing, 
                 std::vector<casacore::Vector<casacore::Double> >& v_uvw,
                 casacore::MSFieldColumns* msfc_p=0);

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


    void uvw_bl( asdm::MainRow* mainRow, std::vector<double> v_timeCentroid, 
                 Enum<CorrelationModeMod::CorrelationMode> correlationMode,
                 std::pair<bool,bool> dataOrder,
                 std::vector<casacore::Vector<casacore::Double> >& v_uvw,
                 casacore::MSFieldColumns* msfc_p=0 );

private:

    struct ArrayParam{
    public:
        asdm::Tag                   subarrayId;        //<! (sub)array identifier
        std::vector<asdm::Tag>           v_ant;        //<! sequence of antennas
        unsigned int                   nrepeat;        //<! number of casacore::MS main table rows per antenna baseline
        Enum<CorrelationModeMod::CorrelationMode> e_correlationMode; //<! correlation mode (original mode passed through the user filter)
        bool operator== ( ArrayParam & rhs){
            std::set<asdm::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(){
            std::ostringstream os;
            for(unsigned int n=0; n<v_ant.size(); n++)os << v_ant[n].toString() << " ";
            os << " nrepeat " << nrepeat;
            return os.str();
        }
    };

    std::map<asdm::Tag,ArrayParam>            m_array_;
    std::map<asdm::Tag,casacore::MPosition>             m_antPos_;
    std::map<asdm::Tag,casacore::Vector<casacore::Double> > m_antUVW_;

    // the 3 fundamental attributes of the state machine
    double                         timeCentroid_;
    std::vector<std::vector<asdm::Angle> >         phaseDir_; 
    asdm::Tag                            subarrayId_;
    
    casacore::Vector<casacore::Double>           sduvw_;

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

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


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

};