Use pol_axis directly

AlecThomson · Oct 9, 2023 · 71f84e8 · 71f84e8
1 parent 5d2c17c
commit 71f84e8
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/fixms/fix_ms_corrs.py b/fixms/fix_ms_corrs.py
@@ -39,6 +39,8 @@ def get_pol_axis(ms: Path, feed_idx: Optional[int] = None) -> u.Quantity:
     """
     with table((ms / "FEED").as_posix(), readonly=True, ack=False) as tf:
         ms_feed = tf.getcol("RECEPTOR_ANGLE") * u.rad
+        # PAF is at 45deg to feeds
+        # 45 - feed_angle = pol_angle
         pol_axes = -(ms_feed - 45.0 * u.deg)
 
     if feed_idx is None:
@@ -167,7 +169,7 @@ def convert_correlations(correlations: np.ndarray, pol_axis: u.Quantity) -> np.n
 
 
     """
-    theta = (pol_axis + 45.0 * u.deg).to(u.rad).value
+    theta = pol_axis.to(u.rad).value
     correction_matrix = np.matrix(
         [
             [