aboutsummaryrefslogtreecommitdiffstats
path: root/fg21sim
diff options
context:
space:
mode:
Diffstat (limited to 'fg21sim')
-rw-r--r--fg21sim/extragalactic/clusters/halo.py10
-rw-r--r--fg21sim/utils/units.py2
2 files changed, 12 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.
diff --git a/fg21sim/utils/units.py b/fg21sim/utils/units.py
index e65f831..a6d7034 100644
--- a/fg21sim/utils/units.py
+++ b/fg21sim/utils/units.py
@@ -59,6 +59,8 @@ class UnitConversions:
# Angle
rad2deg = au.rad.to(au.deg)
deg2rad = au.deg.to(au.rad)
+ rad2arcsec = au.rad.to(au.arcsec)
+ arcsec2rad = au.arcsec.to(au.rad)
arcmin2deg = au.arcmin.to(au.deg)
deg2arcmin = au.deg.to(au.arcmin)
arcsec2deg = au.arcsec.to(au.deg)