Changeset 2662 for tags/asap-4.1.0/src/STWriter.cpp

Timestamp:

10/10/12 14:53:53 (12 years ago)

Author:

Malte Marquarding

Message:

merge from trunk into release candidate

Location:

tags/asap-4.1.0

Files:

• . (modified) (1 prop)
• src (modified) (1 prop)
• src/STWriter.cpp (modified) (9 diffs)

tags/asap-4.1.0/src/STWriter.cpp

@@ -38 +38 @@
 #include <casa/Utilities/CountedPtr.h>
 #include <casa/Utilities/Assert.h>
+#include <casa/Logging/LogIO.h>
 
 #include <atnf/PKSIO/PKSrecord.h>
-#ifndef NOPKSMS
-#include <atnf/PKSIO/PKSMS2writer.h>
-#endif
 #include <atnf/PKSIO/PKSSDwriter.h>
 #include <atnf/PKSIO/SrcType.h>
@@ -51 +49 @@
 #include <tables/Tables/ArrayColumn.h>
 
+#include <fits/FITS/FITSSpectralUtil.h>
+
 #include "STFITSImageWriter.h"
 #include "STAsciiWriter.h"
@@ -68 +68 @@
   t.upcase();
   if (t == "MS2") {
-    #ifdef NOPKSMS
     throw (AipsError("MS2 OUTPUT FORMAT IS NO LONGER SUPPORTED"));
-    #else
-    writer_ = new PKSMS2writer();
-    #endif
   } else if (t == "SDFITS") {
     writer_ = new PKSSDwriter();
@@ -99 +95 @@
   t.upcase();
   if (t== "MS2") {
-    #ifdef NOPKSMS
     throw (AipsError("MS2 OUTPUT FORMAT IS NO LONGER SUPPORTED"));
-    #else
-    writer_ = new PKSMS2writer();
-    #endif
   } else if (t== "SDFITS") {
     writer_ = new PKSSDwriter();
@@ -121 +113 @@
   // in asap on-the-fly.
 
+  String freqframe;
+  FITSSpectralUtil::specsysFromFrame(freqframe, 
+                                     in->frequencies().getFrame(true));
   CountedPtr<Scantable> inst = in;
   if (in->frequencies().getFrame(true) != in->frequencies().getFrame(false)) {
+    FITSSpectralUtil::specsysFromFrame(freqframe, 
+                                     in->frequencies().getFrame(false));
+  
     STMath stm(false);
     inst = stm.frequencyAlign(in);
@@ -151 +149 @@
   // before. Need to check with the Offline User Test data
   STHeader hdr = inst->getHeader();
-  //const Int nPol  = hdr.npol;
-  //const Int nChan = hdr.nchan;
   std::vector<uint> ifs = inst->getIFNos();
   int nIF = inst->nif();//ifs.size();
@@ -177 +173 @@
   // Create the output file and write static data.
   Int status;
-//   status = writer_->create(String(filename), hdr.observer, hdr.project,
-//                            hdr.antennaname, hdr.antennaposition,
-//                            hdr.obstype, hdr.fluxunit,
-//                            hdr.equinox, hdr.freqref,
-//                            nChan, nPol, havexpol, False);
   status = writer_->create(String(filename), hdr.observer, hdr.project,
                            inst->getAntennaName(), hdr.antennaposition,
                            hdr.obstype, hdr.fluxunit,
-                           hdr.equinox, hdr.freqref,
+                           hdr.equinox, freqframe,
                            nChan, nPol, havexpol, False);
   if ( status ) {
@@ -331 +322 @@
     ++scanit;
   }
-  ostringstream oss;
-  oss << "STWriter: wrote " << count << " rows to " << filename;
-  pushLog(String(oss));
-  
+  LogIO os( casa::LogOrigin("STWriter"));
+  os << "STWriter: wrote " << count << " rows to " << filename
+     << casa::LogIO::POST;
+
   writer_->close();
-  //if MS2 delete POINTING table exists and copy the one in the keyword
-  if ( format_ == "MS2" ) {
-    replacePtTab(table, filename);
-  }
+
   return 0;
 }
@@ -390 +378 @@
 }
 
-// For writing MS data, if there is the reference to
-// original pointing table it replace it by it.
-void STWriter::replacePtTab (const Table& tab, const std::string& fname)
-{
-  String oldPtTabName = fname;
-  oldPtTabName.append("/POINTING");
-  if ( tab.keywordSet().isDefined("POINTING") ) {
-    String PointingTab = tab.keywordSet().asString("POINTING");
-    if ( Table::isReadable(PointingTab) ) {
-      Table newPtTab(PointingTab, Table::Old);
-      newPtTab.copy(oldPtTabName, Table::New);
-      ostringstream oss;
-      oss << "STWriter: copied  " <<PointingTab  << " to " << fname;
-      pushLog(String(oss));
-    }
-  }
-}
 
 // get obsType string from SRCTYPE value