aboutsummaryrefslogtreecommitdiffstats
diff options
context:
space:
mode:
authorAaron LI <aly@aaronly.me>2019-01-26 21:31:17 +0800
committerAaron LI <aly@aaronly.me>2019-01-26 21:31:17 +0800
commit55128a7f3bf5b01f5480b2e29bd2021353cb9991 (patch)
tree7a08d6ea54c89947e8780c28ec8044f07ce798e6
parent896e68f1173210c9b62a6f0d649ddf243d5ce30b (diff)
downloadfg21sim-55128a7f3bf5b01f5480b2e29bd2021353cb9991.tar.bz2
clusters/halo: Reimplement the 'radius' property
-rw-r--r--fg21sim/extragalactic/clusters/halo.py17
1 files changed, 11 insertions, 6 deletions
diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py
index 3ddfe5a..33998de 100644
--- a/fg21sim/extragalactic/clusters/halo.py
+++ b/fg21sim/extragalactic/clusters/halo.py
@@ -197,6 +197,14 @@ class RadioHalo1M:
"""
return self.t_merger
+ @property
+ def radius(self):
+ """
+ The estimated radius of the simulated radio halo.
+ Unit: [kpc]
+ """
+ return self.f_radius * self.radius_turb(self.t_merger)
+
@lru_cache()
def radius_strip(self, t_merger):
"""
@@ -743,9 +751,6 @@ class RadioHaloAM(RadioHalo1M):
The redshifts at each merger event, from small to large.
merger_num : int
Number of merger events traced for the cluster.
- radius : float
- The radius of the radio halo.
- Unit: [kpc]
"""
def __init__(self, M_obs, z_obs, M_main, M_sub, z_merger,
merger_num, radius, configs=CONFIGS):
@@ -756,15 +761,15 @@ class RadioHaloAM(RadioHalo1M):
M_main=M_main, M_sub=M_sub,
z_merger=z_merger, configs=configs)
self.merger_num = merger_num
- self.radius_ = radius # [kpc]
@property
def radius(self):
"""
- The radius of the final radio halo.
+ The halo radius estimated by using the maximum turbulence radius.
Unit: [kpc]
"""
- return self.radius_
+ r_turb = [self.radius_turb(tm) for tm in self.t_merger]
+ return self.f_radius * max(r_turb)
@property
def t_begin(self):