Changeset 1007

Timestamp:

04/07/06 13:29:08 (18 years ago)

Author:

mar637

Message:

Fixed computiaons of stokes and back after discussions with Jim Caswell.

Location:

trunk/src

Files:

• STPolLinear.cpp (modified) (3 diffs)
• STPolLinear.h (modified) (1 diff)
• STPolStokes.cpp (modified) (2 diffs)

trunk/src/STPolLinear.cpp

@@ -40 +40 @@
         break;
       case 1:
-        out = Vector<Float>(q * cos(phase) - getSpectrum(2) * sin(phase));
+        out = Vector<Float>(q * cos(phase) - Float(2.0)*getSpectrum(2) * sin(phase));
         break;
       case 2:
-        out = Vector<Float>(q * sin(phase) + getSpectrum(2) * cos(phase));
+        out = Vector<Float>(q * sin(phase) + Float(2.0)*getSpectrum(2) * cos(phase));
         break;
       case 3:
-        out = getFeedHand() * Vector<Float>(getSpectrum(3));
+        out = getFeedHand() * Float(2.0) * Vector<Float>(getSpectrum(3));
         break;
     }
@@ -94 +94 @@
   switch(index) {
   case 0:
-    out = (I + V);
+    out = (I + V)/Float(2.0);
     break;
   case 1:
-    out = (I - V);
+    out = (I - V)/Float(2.0);
     break;
   default:
@@ -146 +146 @@
   Q = Q2;
   Matrix<Float>& specs = getSpectra();
-  specs.column(0) = I+Q;
-  specs.column(1) = I-Q;
-  specs.column(2) = U;
+  specs.column(0) = (I+Q)/Float(2.0);
+  specs.column(1) = (I-Q)/Float(2.0);
+  specs.column(2) = U/Float(2.0);
 
 }

trunk/src/STPolLinear.h

@@ -19 +19 @@
 
 /**
-The linear representation of polarisation
+  The linear representation of polarisation.
+  We are using the following convention:
+  @li I = XX + YY
+  @li Q = XX - YY
+  @li U = 2*Real(XY)
+  @li V = 2*Imag(XY)
 
-@author Malte Marquarding
+  resulting in:
+  @li I' = I
+  @li Q' = Q * cos(theta) - V *sin(theta)
+  @li U' = Q * sin(theta) + U * cos(theta)
+  @li V' = V
+  @author Malte Marquarding
+
 */
 class STPolLinear : public STPol

trunk/src/STPolStokes.cpp

@@ -60 +60 @@
   case 0:
     out = (getSpectrum(0)+getSpectrum(1))/Float(2.0);
+    break;
   case 1:
     out = (getSpectrum(0)-getSpectrum(1))/Float(2.0);
+    break;
+  case 2:
+    out = getSpectrum(2)/Float(2.0);
+    break;
+  case 3:
+    out = getSpectrum(3)/Float(2.0);
     break;
   default:
@@ -81 +88 @@
   switch(index) {
   case 0:
-    out = (getSpectrum(0) + getSpectrum(3));
+    out = (getSpectrum(0) + getSpectrum(3))/Float(2.0);
     break;
   case 1:
-    out = (getSpectrum(0) - getSpectrum(3));
+    out = (getSpectrum(0) - getSpectrum(3))/Float(2.0);
     break;
   default: