Changeset 272 for trunk/src/SDMath.cc

Timestamp:

01/24/05 01:00:14 (19 years ago)

Author:

kil064

Message:

add reference time arg. to function 'VelocityAlignment?'

File:

• trunk/src/SDMath.cc (modified) (8 diffs)

trunk/src/SDMath.cc

@@ -49 +49 @@
 #include <casa/Quanta/MVEpoch.h>
 #include <casa/Quanta/QC.h>
+#include <casa/Quanta/MVTime.h>
 #include <casa/Utilities/Assert.h>
 
@@ -107 +108 @@
 
 
-SDMemTable* SDMath::velocityAlignment (const SDMemTable& in) const
+SDMemTable* SDMath::velocityAlignment (const SDMemTable& in, const String& refTime) const
 {
 
@@ -139 +140 @@
 // Do it
 
-   return velocityAlign (in, velSystem, velUnit, doppler);
+   return velocityAlign (in, velSystem, velUnit, doppler, refTime);
 }
 
@@ -1103 +1104 @@
                                    MFrequency::Types velSystem,
                                    const String& velUnit,
-                                   MDoppler::Types doppler) const
+                                   MDoppler::Types doppler,
+                                   const String& refTime) const
 {
 // Get Header
@@ -1132 +1134 @@
    const uInt nSrcTab = srcTab.nelements();
    cerr << "Found " << srcTab.nelements() << " sources to align " << endl;
-/*
-   cerr << "source table = " << srcTab << endl;
-   cerr << "source idx = " << srcIdx << endl;
-   cerr << "first row = " << firstRow << endl;
-*/
-
-// Set reference Epoch to time of first row
-
-   MEpoch::Ref timeRef = MEpoch::Ref(in.getTimeReference());
+
+// Set reference Epoch to time of first row or given String
+
    Unit DAY(String("d"));
-   Quantum<Double> tQ(times[0], DAY);
-   MVEpoch mve(tQ);
-   MEpoch refTime(mve, timeRef);
+   MEpoch::Ref epochRef(in.getTimeReference());
+   MEpoch refEpoch;
+   if (refTime.length()>0) {
+      refEpoch = epochFromString(refTime, in.getTimeReference());
+   } else {
+      Quantum<Double> tQ(times[0], DAY);
+      MVEpoch mve(tQ);
+      refEpoch = MEpoch(mve, epochRef);
+   }
+   cerr << "Aligning at reference Epoch " << formatEpoch(refEpoch) << endl;
 
 // Set Reference Position
@@ -1166 +1169 @@
          MDirection refDir = in.getDirection(firstRow[iSrc]);
          uInt idx = (iSrc*nFreqIDs) + fqID;
-         vA[idx] = new VelocityAligner<Float>(sC, nChan, refTime, refDir, refPos,
+         vA[idx] = new VelocityAligner<Float>(sC, nChan, refEpoch, refDir, refPos,
                                               velUnit, doppler, velSystem);
       }
@@ -1196 +1199 @@
     Quantum<Double> tQ2(times[iRow],DAY);
     MVEpoch mv2(tQ2);
-    MEpoch epoch(mv2, timeRef);
+    MEpoch epoch(mv2, epochRef);
 
 // Get FreqID vector.  One freqID per IF
@@ -1610 +1613 @@
    }
 }
+
+MEpoch SDMath::epochFromString (const String& str, MEpoch::Types timeRef) const
+{
+   Quantum<Double> qt;
+   if (MVTime::read(qt,str)) {
+      MVEpoch mv(qt);
+      MEpoch me(mv, timeRef);
+      return me;
+   } else {
+      throw(AipsError("Invalid format for Epoch string"));
+   }
+}
+
+
+String SDMath::formatEpoch(const MEpoch& epoch)  const
+{
+   MVTime mvt(epoch.getValue());
+   return mvt.string(MVTime::YMD) + String(" (") + epoch.getRefString() + String(")");
+}
+