Changeset 401 for trunk/src/SDMemTable.cc

Timestamp:

02/09/05 22:28:48 (19 years ago)

Author:

kil064

Message:

add rest freq specification by line name
add function spectraLines to list known lines

File:

• trunk/src/SDMemTable.cc (modified) (5 diffs)

trunk/src/SDMemTable.cc

@@ -742 +742 @@
 // Set reference frame conversion  (requires row)
 
-  MDirection direct = getDirection(whichRow);
+  MDirection dir = getDirection(whichRow);
   MEpoch epoch = getEpoch(whichRow);
   MPosition pos = getAntennaPosition();
-  if (!spec.setReferenceConversion(mtype,epoch,pos,direct)) {
+//
+  if (!spec.setReferenceConversion(mtype,epoch,pos,dir)) {
     throw(AipsError("Couldn't convert frequency frame."));
   }
@@ -1326 +1327 @@
 }
 
-casa::Bool SDMemTable::setRestFreqs (const casa::Vector<Double>& restFreqs, const String& sUnit,
-                                     const casa::String& source,
-                                     casa::Int whichIF)
+Bool SDMemTable::setRestFreqs (const Vector<Double>& restFreqsIn, 
+                                     const String& sUnit,
+                                     const vector<string>& lines,
+                                     const String& source,
+                                     Int whichIF)
 {
    const Int nIFs = nIF();
@@ -1335 +1338 @@
    }
 
+// FInd vector of restfrequencies
+// Double takes precedence over String
+
+   Unit unit;
+   Vector<Double> restFreqs;
+   if (restFreqsIn.nelements()>0) {
+      restFreqs.resize(restFreqsIn.nelements());
+      restFreqs = restFreqsIn;
+      unit = Unit(sUnit);
+   } else if (lines.size()>0) {
+      const uInt nLines = lines.size();
+      unit = Unit("Hz");
+      restFreqs.resize(nLines);
+      MFrequency lineFreq;
+      for (uInt i=0; i<nLines; i++) {
+         String tS(lines[i]);
+         tS.upcase();
+         if (MeasTable::Line(lineFreq, tS)) {
+            restFreqs[i] = lineFreq.getValue().getValue();          // Hz
+         } else {
+            String s = String(lines[i]) + String(" is an unrecognized spectral line");
+            throw(AipsError(s));
+         }
+      }
+   } else {
+      throw(AipsError("You have not specified any rest frequencies or lines"));
+   }
+
 // If multiple restfreqs, must be length nIF. In this
 // case we will just replace the rest frequencies
@@ -1340 +1371 @@
 
    const uInt nRestFreqs = restFreqs.nelements();
-   Unit unit(sUnit);
    Int idx = -1;
    SDFrequencyTable sdft = getSDFreqTable();
@@ -1408 +1438 @@
 }
 
+void SDMemTable::spectralLines () const
+{
+   Vector<String> lines = MeasTable::Lines();
+   MFrequency lineFreq;
+   Double freq;
+//
+   cerr.flags(std::ios_base::left);
+   cerr << "Line      Frequency (Hz)" << endl;
+   cerr << "-----------------------" << endl;
+   for (uInt i=0; i<lines.nelements(); i++) {
+     MeasTable::Line(lineFreq, lines[i]);
+     freq = lineFreq.getValue().getValue();          // Hz
+//
+     cerr << setw(11) << lines[i] << setprecision(10) << freq << endl;
+   }
+}
 
 void SDMemTable::renumber()