Changeset 472 for trunk/src/SDMemTable.cc

Timestamp:

02/18/05 11:16:27 (19 years ago)

Author:

kil064

Message:

add a new getStokesSPectrum function (for the writer) which gets all
Stokes or a beam/IF slice

File:

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

trunk/src/SDMemTable.cc

@@ -366 +366 @@
 std::vector<float> SDMemTable::getStokesSpectrum(Int whichRow, Bool doPol, 
                                                  Float paOffset) const 
-  // Gets
-  //  doPol=False  : I,Q,U,V
-  //  doPol=True   : I,P,PA,V   ; P = sqrt(Q**2+U**2), PA = 0.5*atan2(Q,U)
-  //
+//
+// Gets one STokes parameter 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);
@@ -450 +451 @@
 
   return stlout;
+}
+
+Array<Float> SDMemTable::getStokesSpectrum(Int whichRow, Int iBeam, Int iIF) const
+{
+
+// Get data
+
+  Array<Float> arr;
+  stokesCol_.get(whichRow, arr);
+
+// Set current cursor location and overwrite polarization axis
+
+  const IPosition& shape = arr.shape();
+  IPosition start(shape.nelements(),0);  
+  IPosition end(shape-1);
+  if (iBeam!=-1) {
+     start(asap::BeamAxis) = iBeam;
+     end(asap::BeamAxis) = iBeam;
+  }
+  if (iIF!=-1) {
+     start(asap::IFAxis) = iIF;
+     end(asap::IFAxis) = iIF;
+  }
+
+// Get slice
+
+  return arr(start,end);
 }
 
@@ -1172 +1200 @@
   histCol_.put(rno, sdc.getHistory());
   fitCol_.put(rno, sdc.getFitMap());
-
   return true;
 }