Changeset 2657 for trunk/src/STFiller.cpp

Timestamp:

10/10/12 12:17:12 (12 years ago)

Author:

Malte Marquarding

Message:

Ticket #280: Use FITSSpectralUtil::specsysFromFrame to write correct FITS output frame

File:

• trunk/src/STFiller.cpp (modified) (3 diffs)

trunk/src/STFiller.cpp

@@ -36 +36 @@
 #include <atnf/PKSIO/NROReader.h>
 #include <casa/Logging/LogIO.h>
+#include <fits/FITS/FITSSpectralUtil.h>
 
 #include <time.h>
@@ -226 +227 @@
   Vector<Int> start(nIF_, 1);
   Vector<Int> end(nIF_, 0);
-  reader_->select(beams, ifs, start, end, ref, True, haveXPol_[0], False, getPt);
+  reader_->select(beams, ifs, start, end, ref, True, haveXPol_[0], False, 
+                  getPt);
   table_->setHeader(*header_);
   //For MS, add the location of POINTING of the input MS so one get
@@ -249 +251 @@
   //MS style (as PKSMS2reader converts the original frame
   //in FITS standard style)
-  if (freqFrame == "TOPOCENT") {
-    freqFrame = "TOPO";
-  } else if (freqFrame == "GEOCENER") {
-    freqFrame = "GEO";
-  } else if (freqFrame == "BARYCENT") {
-    freqFrame = "BARY";
-  } else if (freqFrame == "GALACTOC") {
-    freqFrame = "GALACTO";
-  } else if (freqFrame == "LOCALGRP") {
-    freqFrame = "LGROUP";
-  } else if (freqFrame == "CMBDIPOL") {
-    freqFrame = "CMB";
-  } else if (freqFrame == "SOURCE") {
-    freqFrame = "REST";
-  }
+  MFrequency::Types fframe;
+  FITSSpectralUtil::frameFromSpecsys(fframe, freqFrame);
   // set both "FRAME" and "BASEFRAME"
-  table_->frequencies().setFrame(freqFrame, false);
-  table_->frequencies().setFrame(freqFrame,true);
+  table_->frequencies().setFrame(fframe, false);
+  table_->frequencies().setFrame(fframe,true);
+  
   //table_->focus().setParallactify(true);
 }