diff options
Diffstat (limited to 'fg21sim')
-rw-r--r-- | fg21sim/extragalactic/clusters/halo.py | 10 | ||||
-rw-r--r-- | fg21sim/utils/units.py | 2 |
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) |