casa6

/*

 * PointingDirectionCache.cc

 *

 *  Created on: Dec 1, 2016

 *      Author: jjacobs

 */

​

#include <msvis/MSVis/PointingDirectionCache.h>

​

#include <casacore/casa/Exceptions/Error.h>

#include <casacore/casa/BasicSL/String.h>

#include <casacore/measures/Measures/MDirection.h>

#include <casacore/ms/MeasurementSets/MSPointingColumns.h>

​

#include <memory>

​

using namespace casacore;

​

namespace casa {

namespace vi {

​

PointingColumns::PointingColumns (const MSPointingColumns & pointingColumns)

: pointingColumns_p (pointingColumns)

{}

​

Pointing

PointingColumns::getPointingRow (int row, double /*targetTime*/, bool asMeasure) const

{

    Pointing pointing;

​

    // See if the row's time is no earlier than 5 minutes before target

    // An optimization to avoid unnecessary direction measure conversions

    // for rows unlikely to be used (the conversions are fairly expensive)

​

    pointing.time = pointingColumns_p.time () (row);

    pointing.interval = pointingColumns_p.interval () (row);

    pointing.antennaId = pointingColumns_p.antennaId() (row);

    pointing.row = row;

    pointing.source = this;

​

    if (asMeasure){

​

        // Pointing is likely to be useful so get all the data.

​

        pointing.direction.reset (new MDirection (pointingColumns_p.directionMeas (row)));

    }

​

    return pointing;

}

​

int

PointingColumns::nRows () const

{

    return pointingColumns_p.nrow();

}

​

​

PointingDirectionCache::PointingDirectionCache (int nAntennas,

                                                const PointingSource & pointingSource)

: antennaLevelCache_p (MinTimeEntries, MaxTimeEntries, nAntennas, this),

  antennaEarliestTime_p (nAntennas, -1.0),

  lastRowRead_p (-1),

  nFallbacks_p (0),