diff options
author | Aaron LI <aly@aaronly.me> | 2019-01-26 21:31:17 +0800 |
---|---|---|
committer | Aaron LI <aly@aaronly.me> | 2019-01-26 21:31:17 +0800 |
commit | 55128a7f3bf5b01f5480b2e29bd2021353cb9991 (patch) | |
tree | 7a08d6ea54c89947e8780c28ec8044f07ce798e6 | |
parent | 896e68f1173210c9b62a6f0d649ddf243d5ce30b (diff) | |
download | fg21sim-55128a7f3bf5b01f5480b2e29bd2021353cb9991.tar.bz2 |
clusters/halo: Reimplement the 'radius' property
-rw-r--r-- | fg21sim/extragalactic/clusters/halo.py | 17 |
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): |