#include <casa/BasicSL/Complex.h>
#include <synthesis/TransformMachines/PBMath.h>
#include <synthesis/TransformMachines2/VPSkyJones.h>
#include <ms/MeasurementSets/MeasurementSet.h>
#include <tables/Tables.h>
#include <tables/Tables/TableRecord.h>
#include <ms/MeasurementSets/MSColumns.h>
#include <measures/Measures/Stokes.h>
#include <casa/BasicSL/Constants.h>
#include <components/ComponentModels/Flux.h>
#include <components/ComponentModels/ComponentShape.h>
#include <images/Images/ImageInterface.h>
#include <lattices/Lattices/LatticeStepper.h>
#include <lattices/Lattices/LatticeIterator.h>
#include <coordinates/Coordinates/CoordinateSystem.h>
#include <coordinates/Coordinates/DirectionCoordinate.h>
#include <coordinates/Coordinates/SpectralCoordinate.h>
#include <coordinates/Coordinates/StokesCoordinate.h>
#include <coordinates/Coordinates/Projection.h>
#include <casa/Utilities/Assert.h>
using namespace casacore;
using namespace casacore;
using namespace casacore;
using namespace casa::refim;
VPSkyJones::VPSkyJones(const MSColumns& msc, Table& tab,
const Quantity &parAngleInc,
BeamSquint::SquintType doSquint,
const Quantity &skyPositionThreshold)
: BeamSkyJones(parAngleInc, doSquint,skyPositionThreshold)
LogIO os(LogOrigin("VPSkyJones", "VPSkyJones"));
const uInt nrow = tab.nrow();
ScalarColumn<String> telCol(tab, "telescope");
ScalarColumn<Int> antCol(tab, "antenna");
ScalarColumn<TableRecord> recCol(tab, "pbdescription");
for (uInt i=0; i < nrow; ++i) {
if (indexTelescope(telCol(i))<0) {
if(recCol(i).isDefined("name") && (recCol(i).asString("name") != "REFERENCE")){