Task #2655

Added function to set the reference type (of a moving object)

LCS/MSLofar/include/MSLofar/MSLofarFieldColumns.h

@@ -106,6 +106,11 @@ namespace LOFAR {
       { return rwTileBeamDirMeas_p; }
     // </group>
 
+    // Set the direction reference type for the various direction columns.
+    // This can only be done when the table has no rows.
+    // Trying to do so at other times will throw an exception.
+    void setDirectionRef(casa::MDirection::Types ref);
+
   protected:
     //# Default constructor creates a object that is not usable. Use the attach
     //# function correct this.

LCS/MSLofar/src/MSLofarFieldColumns.cc

@@ -72,4 +72,10 @@ namespace LOFAR {
     rwTileBeamDirMeas_p.attach (msLofarField, "LOFAR_TILE_BEAM_DIR");
   }
 
+  void MSLofarFieldColumns::setDirectionRef(MDirection::Types ref)
+  {
+    rwTileBeamDirMeas_p.setDescRefCode (ref);
+    MSFieldColumns::setDirectionRef (ref);
+  }
+
 } //# end namespace

LCS/MSLofar/test/tMSLofar.cc

@@ -40,13 +40,15 @@ using namespace LOFAR;
 using namespace casa;
 using namespace std;
 
-void createMS()
+void createMS (MDirection::Types dirRef)
 {
   Timer timer;
   TableDesc td(MSLofar::requiredTableDesc());
   SetupNewTable setup("tMSLofar_tmp.ms", td, Table::New);
   MSLofar ms(setup);
   ms.createDefaultSubtables(Table::New);
+  MSLofarFieldColumns fieldCols(ms.field());
+  fieldCols.setDirectionRef (dirRef);
   ms.flush();
   timer.show("create");
 }
@@ -171,7 +173,7 @@ void checkMS (const MSLofar& ms)
   ASSERT (fld.flagRow()(0) == False);
   MDirection tileBeamDir = fld.tileBeamDirMeasCol()(0);
   tileBeamDir.print(cout);
-  cout<<endl;
+  cout << ' ' << tileBeamDir.getRefString() << endl;
 
   ASSERT (obs.telescopeName()(0) == "LOFAR");
   ASSERT (allEQ(obs.timeRange()(0), Vector<double>(2,10.)));
@@ -220,7 +222,13 @@ void checkMS (const MSLofar& ms)
 int main()
 {
   try {
-    createMS();
+    // Create for a J2000 direction.
+    createMS (MDirection::J2000);
+    openMS();
+    fillMS();
+    checkMS (MSLofar("tMSLofar_tmp.ms"));
+    // Create for a SUN direction.
+    createMS (MDirection::SUN);
     openMS();
     fillMS();
     checkMS (MSLofar("tMSLofar_tmp.ms"));