Changeset 426

Timestamp:

02/14/05 13:26:04 (20 years ago)

Author:

kil064

Message:

add paOffset to function SDmemTable::getSTokesSpectrum

Location:

trunk/src

Files:

• SDMemTable.cc (modified) (2 diffs)
• SDMemTable.h (modified) (1 diff)
• SDMemTableWrapper.h (modified) (1 diff)
• SDPol.cc (modified) (1 diff)
• python_SDMemTable.cc (modified) (1 diff)

trunk/src/SDMemTable.cc

@@ -346 +346 @@
 }
 
-std::vector<float> SDMemTable::getStokesSpectrum(Int whichRow, Bool doPol) const 
+std::vector<float> SDMemTable::getStokesSpectrum(Int whichRow, Bool doPol, Float paOffset) const 
 //
 // Gets
@@ -376 +376 @@
         out = SDPolUtil::polarizedIntensity(Q,U);
      } else if (polSel_==2) {                                 // P.A.
-        out = SDPolUtil::positionAngle(Q,U);
+        out = SDPolUtil::positionAngle(Q,U) + paOffset;
      }
 //

trunk/src/SDMemTable.h

@@ -94 +94 @@
   virtual std::vector<float> getSpectrum(casa::Int whichRow=0) const;
   virtual std::vector<bool> getMask(casa::Int whichRow=0) const;
+  // Get STokes at cursor location. One of either I,Q,U,V or I,P,PA,V (doPol=True)
+  // If the latter, you can add a PA offset (degrees)
   virtual std::vector<float> getStokesSpectrum(casa::Int whichRow=0, 
-                                               casa::Bool doPol=casa::False) const;
+                                               casa::Bool doPol=casa::False,
+                                               casa::Float paOffset=0.0) const;
 
   virtual casa::Float getTsys(casa::Int whichRow=0) const;

trunk/src/SDMemTableWrapper.h

@@ -81 +81 @@
   }
 
-  std::vector<float> getStokesSpectrum(int whichRow=0, bool doPol=false) const {
-    return table_->getStokesSpectrum(whichRow, doPol);
+  std::vector<float> getStokesSpectrum(int whichRow=0, bool doPol=false, float paOff=0.0) const {
+    return table_->getStokesSpectrum(whichRow, doPol, paOff);
   }
 

trunk/src/SDPol.cc

@@ -230 +230 @@
                                        const Array<Float>& U)
 {
-   return Float(180.0/C::pi/2.0)*atan2(Q,U);
+   return Float(180.0/C::pi/2.0)*atan2(Q,U);       // Degrees
 }
 

trunk/src/python_SDMemTable.cc

@@ -74 +74 @@
     .def("_getstokesspectrum", &SDMemTableWrapper::getStokesSpectrum,
          (boost::python::arg("whichRow")=0),
-         (boost::python::arg("pol")=false) )
+         (boost::python::arg("pol")=false),
+         (boost::python::arg("pa")=0.0) )
     .def("_setspectrum",&SDMemTableWrapper::setSpectrum,
          (boost::python::arg("whichRow")=0) )