clusters/halo.py: Add property "angular_radius" (unit: [arcsec])

Signed-off-by: Aaron LI <aly@aaronly.me>

1 files changed, 10 insertions, 0 deletions

diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py
index b65cd08..adc2bc5 100644
--- a/fg21sim/extragalactic/clusters/halo.py
+++ b/fg21sim/extragalactic/clusters/halo.py
@@ -159,6 +159,16 @@ class RadioHalo:
         return r_halo
 
     @property
+    def angular_radius(self):
+        """
+        The angular radius of the radio halo.
+        Unit: [arcsec]
+        """
+        DA = COSMO.DA(self.z) * 1e3  # [Mpc] -> [kpc]
+        theta = self.radius / DA  # [rad]
+        return theta * AUC.rad2arcsec
+
+    @property
     def volume(self):
         """
         The halo volume, calculated from the above radius.