Changeset 490 for trunk

Timestamp:

02/21/05 13:56:17 (20 years ago)

Author:

kil064

Message:

add functions

stokesLength
getStokesSpectrumLabel
getCircularSpectrumLabel

Location:

trunk/src

Files:

• SDMemTable.cc (modified) (3 diffs)
• SDMemTable.h (modified) (2 diffs)

trunk/src/SDMemTable.cc

@@ -370 +370 @@
 }
 
+
+int SDMemTable::stokesLength() const
+{
+   return stokesCol_.shape(0).nelements();        // All rows same shape
+}
+
+
 std::vector<float> SDMemTable::getStokesSpectrum(Int whichRow, Bool doPol, 
                                                  Float paOffset) const 
@@ -421 +428 @@
 }
 
+std::string SDMemTable::getStokesSpectrumLabel (Bool doPol) const
+//
+// Gets STokes label depending on cursor polSel location
+//  doPol=False  : I,Q,U,V
+//  doPol=True   : I,P,PA,V   ; P = sqrt(Q**2+U**2), PA = 0.5*atan2(Q,U)
+//
+{
+  AlwaysAssert(asap::nAxes==4,AipsError);
+  if (nPol()!=1 && nPol()!=2 && nPol()!=4) {
+     throw (AipsError("You must have 1,2 or 4 polarizations to get the Stokes parameters"));
+  }
+//
+   Stokes::StokesTypes type = Stokes::Undefined;
+   switch (polSel_) {
+      case 0:
+        {
+           type = Stokes::I;
+        }
+        break;
+      case 1:
+        {
+           if (doPol) {
+              type = Stokes::Plinear;
+           } else {
+              type = Stokes::Q;
+           }
+        }
+      case 2:
+        {
+           if (doPol) {
+              type = Stokes::Pangle;
+           } else {
+              type = Stokes::U;
+           }
+        }
+        break;
+      case 3:
+        {
+           type = Stokes::V;
+        }
+        break;
+      default:
+        {
+            throw(AipsError("Unknown Stokes type"));
+        }
+   }
+//
+   return SDPolUtil::stokesString(type);
+}
+
+
+
 std::vector<float> SDMemTable::getCircularSpectrum(Int whichRow, 
                                                    Bool doRR) const 
@@ -458 +517 @@
   return stlout;
 }
+
+
+std::string SDMemTable::getCircularSpectrumLabel (Bool doRR) const
+//
+// Gets Circular label 
+//
+{
+   Stokes::StokesTypes type = Stokes::Undefined;
+   if (doRR) {
+      type = Stokes::RR;
+   } else {
+      type = Stokes::LL;
+   }
+//
+   return SDPolUtil::stokesString(type);
+}
+
+
+
+
+
 
 Array<Float> SDMemTable::getStokesSpectrum(Int whichRow, Int iBeam, Int iIF) const

trunk/src/SDMemTable.h

@@ -103 +103 @@
   // getCircularSpectrum is 'getCircular'
 
+  // Get length of STokes spectrum. XX & YY -> I
+  // Anything else the length stays the same.
+  int stokesLength() const;
+
   // Get specific Stokes at cursor location. One of either I,Q,U,V or I,P,PA,V (doPol=True)
   // (determined by the polSel cursor location 0->3)
@@ -110 +114 @@
                                                casa::Float paOffset=0.0) const;
 
+  // Returns String for selected Stokes (polSel_)
+  std::string getStokesSpectrumLabel (casa::Bool doPol) const;
+
   // Get RR or LL at cursor location (except not polSel_)
   virtual std::vector<float> getCircularSpectrum(casa::Int whichRow=0, 
                                                  casa::Bool rr=casa::True) const;
+
+  // Get circular label
+  std::string getCircularSpectrumLabel (casa::Bool rr=casa::True) const;
 
   // Get all Stokes at the specified Beam/IF cursor location (ignoring