Changeset 1779 for branches/mergetest/external/atnf/PKSIO/PKSFITSreader.cc

Timestamp:

07/29/10 19:13:46 (14 years ago)

Author:

Kana Sugimoto

Message:

New Development: Yes

JIRA Issue: No (test merging alma branch)

Ready for Test: Yes

Interface Changes: Yes

What Interface Changed:

Test Programs:

Put in Release Notes: No

Module(s):

Description:

Location:

branches/mergetest

Files:

• . (modified) (1 prop)
• external/atnf/PKSIO/PKSFITSreader.cc (modified) (18 diffs)

branches/mergetest/external/atnf/PKSIO/PKSFITSreader.cc

@@ -34 +34 @@
 //#---------------------------------------------------------------------------
 
-#include <atnf/PKSIO/PKSmsg.h>
 #include <atnf/PKSIO/MBFITSreader.h>
 #include <atnf/PKSIO/SDFITSreader.h>
@@ -44 +43 @@
 #include <casa/BasicMath/Math.h>
 #include <casa/Quanta/MVTime.h>
+#include <casa/Logging/LogIO.h>
 
 //----------------------------------------------- PKSFITSreader::PKSFITSreader
@@ -61 +61 @@
     cReader = new MBFITSreader(retry, interpolate ? 1 : 0);
   }
-
-  // By default, messages are written to stderr.
-  initMsg();
 }
 
@@ -76 +73 @@
 }
 
-//------------------------------------------------------ PKSFITSreader::setMsg
-
-// Set message disposition.  If fd is non-zero messages will be written
-// to that file descriptor, else stored for retrieval by getMsg().
-
-Int PKSFITSreader::setMsg(FILE *fd)
-{
-  PKSmsg::setMsg(fd);
-  cReader->setMsg(fd);
-
-  return 0;
-}
-
 //-------------------------------------------------------- PKSFITSreader::open
 
@@ -95 +79 @@
 Int PKSFITSreader::open(
         const String fitsName,
+        const String antenna,
         Vector<Bool> &beams,
         Vector<Bool> &IFs,
@@ -103 +88 @@
         Bool   &haveSpectra)
 {
-  clearMsg();
-
   int    extraSysCal, haveBase_, *haveXPol_, haveSpectra_, nBeam, *nChan_,
          nIF, *nPol_, status;
@@ -110 +93 @@
                          nChan_, nPol_, haveXPol_, haveBase_, haveSpectra_,
                          extraSysCal);
-  logMsg(cReader->getMsg());
-  cReader->clearMsg();
+  //logMsg(cReader->getMsg());
+  //cReader->clearMsg();
   if (status) {
     return status;
@@ -178 +161 @@
                               obsType_, bunit_, equinox_, radecsys,
                               dopplerFrame_, datobs, utc, refFreq, bandwidth);
-  logMsg(cReader->getMsg());
-  cReader->clearMsg();
+  //logMsg(cReader->getMsg());
+  //cReader->clearMsg();
   if (status) {
     return 1;
@@ -218 +201 @@
   Int status = cReader->getFreqInfo(nIF, startfreq, endfreq);
 
-  logMsg(cReader->getMsg());
-  cReader->clearMsg();
+  //logMsg(cReader->getMsg());
+  //cReader->clearMsg();
   if (!status) {
     startFreq.takeStorage(IPosition(1,nIF), startfreq, TAKE_OVER);
@@ -240 +223 @@
         const Bool getSpectra,
         const Bool getXPol,
+        const Bool getFeedPos,
+        const Bool getPointing,
         const Int  coordSys)
 {
@@ -307 +292 @@
   cGetSpectra = getSpectra;
   cGetXPol    = getXPol;
+  cGetFeedPos = getFeedPos;
+  cGetPointing = getPointing;
   cCoordSys   = coordSys;
 
   uInt maxNChan = cReader->select(start, end, ref, cGetSpectra, cGetXPol,
-                                  cCoordSys);
-  logMsg(cReader->getMsg());
-  cReader->clearMsg();
+                                  cGetFeedPos, cGetPointing, cCoordSys);
+  //logMsg(cReader->getMsg());
+  //cReader->clearMsg();
 
   delete [] end;
@@ -336 +323 @@
 
   Int status = cReader->findRange(nRow, nSel, dateSpan, utcSpan, posns);
-  logMsg(cReader->getMsg());
-  cReader->clearMsg();
+  //logMsg(cReader->getMsg());
+  //cReader->clearMsg();
 
   if (!status) {
@@ -361 +348 @@
 {
   Int status = cReader->read(cMBrec);
-  logMsg(cReader->getMsg());
-  cReader->clearMsg();
+  //logMsg(cReader->getMsg());
+  //cReader->clearMsg();
 
   if (status) {
@@ -378 +365 @@
   pksrec.scanNo  = cMBrec.scanNo;
   pksrec.cycleNo = cMBrec.cycleNo;
+  pksrec.polNo = cMBrec.polNo ;
 
   // Extract MJD.
   Int day, month, year;
-  sscanf(cMBrec.datobs, "%4d-%2d-%2d", &year, &month, &day);
-  pksrec.mjd = MVTime(year, month, Double(day)).day() + cMBrec.utc/86400.0;
+  if ( strstr( cMBrec.datobs, "T" ) == NULL ) {
+    sscanf(cMBrec.datobs, "%4d-%2d-%2d", &year, &month, &day);
+    pksrec.mjd = MVTime(year, month, Double(day)).day() + cMBrec.utc/86400.0;
+  }
+  else {
+    Double dd, hour, min, sec ;
+    sscanf( cMBrec.datobs, "%4d-%2d-%2lfT%lf:%lf:%lf", &year, &month, &dd, &hour, &min, &sec ) ;
+    dd = dd + ( hour * 3600.0 + min * 60.0 + sec ) / 86400.0 ;
+    pksrec.mjd = MVTime(year, month, dd).day() ;
+  }
 
   pksrec.interval  = cMBrec.exposure;
@@ -388 +384 @@
   pksrec.fieldName = trim(cMBrec.srcName);
   pksrec.srcName   = pksrec.fieldName;
+
+  int namelen = pksrec.srcName.length() ;
+  if ( namelen > 4 ) {
+    String srcsub = pksrec.srcName.substr( namelen-4, 4 ) ;
+    if ( srcsub.find( "_psc" ) != string::npos ) {
+      pksrec.fieldName = pksrec.srcName.substr( 0, namelen-4 ) ;
+      pksrec.srcName = pksrec.fieldName + "_ps_calon" ;
+    }
+    else if ( srcsub.find( "_pso" ) != string::npos ) {
+      pksrec.fieldName = pksrec.srcName.substr( 0, namelen-4 ) ;
+      pksrec.srcName = pksrec.fieldName + "_ps" ;
+    }
+    else if ( srcsub.find( "_prc" ) != string::npos ) {
+      pksrec.fieldName = pksrec.srcName.substr( 0, namelen-4 ) ;
+      pksrec.srcName = pksrec.fieldName + "_psr_calon" ;
+    }
+    else if ( srcsub.find( "_pro" ) != string::npos ) {
+      pksrec.fieldName = pksrec.srcName.substr( 0, namelen-4 ) ;
+      pksrec.srcName = pksrec.fieldName + "_psr" ;
+    }
+    else if ( srcsub.find( "_fsc" ) != string::npos ) {
+      pksrec.fieldName = pksrec.srcName.substr( 0, namelen-4 ) ;
+      pksrec.srcName = pksrec.fieldName + "_fs_calon" ;
+    }
+    else if ( srcsub.find( "_fso" ) != string::npos ) {
+      pksrec.fieldName = pksrec.srcName.substr( 0, namelen-4 ) ;
+      pksrec.srcName = pksrec.fieldName + "_fs" ;
+    }
+    else if ( srcsub.find( "_frc" ) != string::npos ) {
+      pksrec.fieldName = pksrec.srcName.substr( 0, namelen-4 ) ;
+      pksrec.srcName = pksrec.fieldName + "_fsr_calon" ;
+    }
+    else if ( srcsub.find( "_fro" ) != string::npos ) {
+      pksrec.fieldName = pksrec.srcName.substr( 0, namelen-4 ) ;
+      pksrec.srcName = pksrec.fieldName + "_fsr" ;
+    }
+    else if ( srcsub.find( "_nsc" ) != string::npos || srcsub.find( "_nrc" ) != string::npos ) {
+      pksrec.fieldName = pksrec.srcName.substr( 0, namelen-4 ) ;
+      pksrec.srcName = pksrec.fieldName + "_nod_calon" ;
+    }
+    else if ( srcsub.find( "_nso" ) != string::npos || srcsub.find( "_nro" ) != string::npos ) {
+      pksrec.fieldName = pksrec.srcName.substr( 0, namelen-4 ) ;
+      pksrec.srcName = pksrec.fieldName + "_nod" ;
+    }
+  }
 
   pksrec.srcDir.resize(2);
@@ -396 +437 @@
   pksrec.srcPM(0)  = 0.0;
   pksrec.srcPM(1)  = 0.0;
-  pksrec.srcVel    = 0.0;
+  pksrec.srcVel    = cMBrec.srcVelocity;
   pksrec.obsType   = trim(cMBrec.obsType);
 
@@ -404 +445 @@
   pksrec.bandwidth = chanWidth * nChan;
   pksrec.freqInc   = cMBrec.fqDelt[0];
-  pksrec.restFreq  = cMBrec.restFreq;
+  pksrec.restFreq.resize(1) ;
+  pksrec.restFreq(0)  = cMBrec.restFreq;
 
   pksrec.tcal.resize(nPol);
@@ -498 +540 @@
 {
   cReader->close();
-  logMsg(cReader->getMsg());
-  cReader->clearMsg();
+  //logMsg(cReader->getMsg());
+  //cReader->clearMsg();
 }