Merge pull request #390 from LCOGT/fix/aperture-phot-diameter

Fix to aperture photometry to correctly use diameter instead of radius
LCOGT · Apr 23, 2024 · afa292a · afa292a
2 parents fafdc5a + 55c3a50
commit afa292a
Show file tree

Hide file tree

Showing 2 changed files with 7 additions and 3 deletions.
diff --git a/CHANGES.md b/CHANGES.md
@@ -1,3 +1,7 @@
+1.16.1 (2023-04-23)
+-------------------
+- Correction to aperture photometry. We were incorrectly using the radius instead of the diameter
+
 1.16.0 (2023-04-18)
 -------------------
 - Calibration frames are now associated with output data products rather than frames

diff --git a/banzai/photometry.py b/banzai/photometry.py
@@ -125,9 +125,9 @@ def do_stage(self, image):
                              'xy': catalog.covar_sigxy.value,
                              'background': catalog.background_mean})
 
-            for r in range(1, 7):
-                radius_arcsec = r / image.pixel_scale
-                sources[f'fluxaper{r}'], sources[f'fluxerr{r}'] = catalog.circular_photometry(radius_arcsec)
+            for d in range(1, 7):
+                radius_arcsec = d / image.pixel_scale / 2.0
+                sources[f'fluxaper{d}'], sources[f'fluxerr{d}'] = catalog.circular_photometry(radius_arcsec)
 
             for r in [0.25, 0.5, 0.75]:
                 sources['fluxrad' + f'{r:.2f}'.lstrip("0.")] = catalog.fluxfrac_radius(r)