From 1abb6bfddb6858fd0140accaf7eece61b3764ae4 Mon Sep 17 00:00:00 2001 From: Aaron LI Date: Sun, 30 Jul 2017 10:19:18 +0800 Subject: clusters/halo.py: Add property "angular_radius" (unit: [arcsec]) Signed-off-by: Aaron LI --- fg21sim/extragalactic/clusters/halo.py | 10 ++++++++++ fg21sim/utils/units.py | 2 ++ 2 files changed, 12 insertions(+) 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 @@ -158,6 +158,16 @@ class RadioHalo: r_halo = helper.radius_halo(mass, self.z_merger) # [kpc] 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): """ 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) -- cgit v1.2.2