From 1abb6bfddb6858fd0140accaf7eece61b3764ae4 Mon Sep 17 00:00:00 2001
From: Aaron LI <aly@aaronly.me>
Date: Sun, 30 Jul 2017 10:19:18 +0800
Subject: clusters/halo.py: Add property "angular_radius" (unit: [arcsec])
Signed-off-by: Aaron LI <aly@aaronly.me>
---
fg21sim/extragalactic/clusters/halo.py | 10 ++++++++++
1 file changed, 10 insertions(+)
(limited to 'fg21sim/extragalactic/clusters')
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):
"""
--
cgit v1.2.2