Changeset 847 for trunk/src/Scantable.cpp

Timestamp:

02/24/06 10:03:21 (18 years ago)

Author:

mar637

Message:

numerous changes before move to new svn repository sourcecode.atnf.csiro.au

File:

• trunk/src/Scantable.cpp (modified) (8 diffs)

trunk/src/Scantable.cpp

@@ -224 +224 @@
   scanCol_.attach(table_, "SCANNO");
   beamCol_.attach(table_, "BEAMNO");
+  ifCol_.attach(table_, "IFNO");
+  polCol_.attach(table_, "POLNO");
   integrCol_.attach(table_, "INTERVAL");
   azCol_.attach(table_, "AZIMUTH");
@@ -292 +294 @@
 
 
-int Scantable::rowToScanIndex( int therow )
-{
-  int therealrow = -1;
-
-  return therealrow;
-}
-
 int Scantable::nscan() const {
   int n = 0;
@@ -344 +339 @@
 std::string Scantable::getFluxUnit() const
 {
-  String tmp;
-  table_.keywordSet().get("FluxUnit", tmp);
-  return tmp;
+  return table_.keywordSet().asString("FluxUnit");
 }
 
@@ -485 +478 @@
   }
   return 0;
+}
+
+
+int Scantable::getBeam(int whichrow) const
+{
+  return beamCol_(whichrow);
+}
+
+int Scantable::getIF(int whichrow) const
+{
+  return ifCol_(whichrow);
+}
+
+int Scantable::getPol(int whichrow) const
+{
+  return polCol_(whichrow);
 }
 
@@ -563 +572 @@
 }
 
-double Scantable::getInterval(int whichrow) const
-{
-  if (whichrow < 0) return 0.0;
-  Double intval;
-  integrCol_.get(Int(whichrow), intval);
-  return intval;
-}
-
 std::vector<bool> Scantable::getMask(int whichrow) const
 {
@@ -619 +620 @@
 {
   table_ = originalTable_;
+  attach();
   selector_.reset();
 }
@@ -629 +631 @@
   }
   table_ = tab;
+  attach();
   selector_ = selection;
 }
@@ -742 +745 @@
 }
 
+std::string Scantable::getAbcissaLabel( int whichrow ) const
+{
+  if ( whichrow > table_.nrow() ) throw(AipsError("Illegal ro number"));
+  const MPosition& mp = getAntennaPosition();
+  const MDirection& md = dirCol_(whichrow);
+  const MEpoch& me = timeCol_(whichrow);
+  const Double& rf = mmolidCol_(whichrow);
+  SpectralCoordinate spc =
+    freqTable_.getSpectralCoordinate(md, mp, me, rf, mfreqidCol_(whichrow));
+
+  String s = "Channel";
+  Unit u = Unit(freqTable_.getUnitString());
+  if (u == Unit("km/s")) {
+    s = CoordinateUtil::axisLabel(spc,0,True,True,True);
+  } else if (u == Unit("Hz")) {
+    Vector<String> wau(1);wau = u.getName();
+    spc.setWorldAxisUnits(wau);
+
+    s = CoordinateUtil::axisLabel(spc,0,True,True,False);
+  }
+  return s;
+
+}
+
+void asap::Scantable::setRestFrequencies( double rf, const std::string& unit )
+{
+  ///@todo lookup in line table
+  Unit u(unit);
+  Quantum<Double> urf(rf, u);
+  uInt id = moleculeTable_.addEntry(urf.getValue("Hz"), "", "");
+  TableVector<uInt> tabvec(table_, "MOLECULE_ID");
+  tabvec = id;
+}
+
+void asap::Scantable::setRestFrequencies( const std::string& name )
+{
+  throw(AipsError("setRestFrequencies( const std::string& name ) NYI"));
+  ///@todo implement
+}
+
 }//namespace asap