Changeset 2662 for tags/asap-4.1.0/src/STFiller.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/STFiller.cpp (modified) (5 diffs)

tags/asap-4.1.0/src/STFiller.cpp

@@ -36 +36 @@
 #include <atnf/PKSIO/NROReader.h>
 #include <casa/Logging/LogIO.h>
+#include <fits/FITS/FITSSpectralUtil.h>
 
 #include <time.h>
@@ -137 +138 @@
   // Get basic parameters.
   if ( anyEQ(haveXPol_, True) ) {
-    pushLog("Cross polarization present");
     for (uInt i=0; i< npols.nelements();++i) {
       if (npols[i] < 3) npols[i] += 2;// Convert Complex -> 2 Floats
@@ -164 +164 @@
     throw(AipsError("Failed to get header."));
   }
+  LogIO os( casa::LogOrigin( "STFiller") );
+  
   if ((header_->obstype).matches("*SW*")) {
     // need robust way here - probably read ahead of next timestamp
-    pushLog("Header indicates frequency switched observation.\n"
-               "setting # of IFs = 1 ");
+    os << "Header indicates frequency switched observation.\n"
+       << "setting # of IFs = 1 " << LogIO::POST;
     nIF_ = 1;
     header_->obstype = String("fswitch");
@@ -226 +228 @@
   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 +252 @@
   //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);
 }