casa6

#include <synthesis/Utilities/FixVis.h>

#include <msvis/MSVis/MSUVWGenerator.h>

#include <msvis/MSVis/SubMS.h>

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

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

#include <casacore/casa/Logging/LogIO.h>

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

#include <casacore/casa/Quanta/MVAngle.h>

#include <casacore/coordinates/Coordinates/Projection.h>

#include <casacore/coordinates/Coordinates/DirectionCoordinate.h>

#include <casacore/coordinates/Coordinates/ObsInfo.h>

#include <casacore/images/Images/PagedImage.h>           // Used to pass coords

#include <casacore/images/Images/ImageInfo.h>            // to FTMachine.

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

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

#include <casacore/ms/MSSel/MSSelection.h>

#include <casacore/ms/MSSel/MSSelectionTools.h>

#include <msvis/MSVis/VisibilityIterator.h>

#include <msvis/MSVis/VisBuffer.h>

#include <casacore/casa/BasicSL/String.h>   // for parseColumnNames()

​

#include <iostream>

​

using namespace casacore;

namespace casa {

​

FixVis::FixVis(MeasurementSet& ms, const String& dataColName) :

  FTMachine(),

  ms_p(ms),

  msc_p(NULL),

  nsel_p(0),

  nAllFields_p(1),

  npix_p(2),

  cimageShape_p(4, npix_p, npix_p, 1, 1), // Can we get away with

  tileShape_p(4, npix_p, npix_p, 1, 1),   // (1, 1, 1, 1)?  Does it matter?

  tiledShape_p(cimageShape_p, tileShape_p),

  antennaSel_p(false),

  freqFrameValid_p(false)

  //  obsString_p("")

{

  logSink() << LogOrigin("FixVis", "") << LogIO::NORMAL3;

​

  antennaId_p.resize();

  antennaSelStr_p.resize();

  distances_p.resize();

  dataCols_p = SubMS::parseColumnNames(dataColName, ms);

  nDataCols_p = dataCols_p.nelements();

​

  nchan_p = 1; // imageNchan_p;

​

  spectralwindowids_p.resize(ms_p.spectralWindow().nrow());

  indgen(spectralwindowids_p);

​

  lockCounter_p = 0;

}

// Destructor

FixVis::~FixVis()

{

  if(!ms_p.isNull())

    ms_p.flush();

​

  delete msc_p;