diff options
author | Aaron LI <aly@aaronly.me> | 2019-02-21 17:57:18 +0800 |
---|---|---|
committer | Aaron LI <aly@aaronly.me> | 2019-02-21 17:57:18 +0800 |
commit | 528d71723f5c58261a6907528ab814f8f590b13f (patch) | |
tree | 1857f05672da6551f4ab102813ce278a5391ddc3 | |
parent | b4bd6c891ad3d03e0b3c59af19b1fca4f76c2015 (diff) | |
download | fg21sim-528d71723f5c58261a6907528ab814f8f590b13f.tar.bz2 |
clusters/halo: Update _validate_t_merger() to also check end time
Rename '_validate_t_merger()' to '_validate_time()' and update this
method to also check the end time of a merger.
-rw-r--r-- | fg21sim/extragalactic/clusters/halo.py | 23 |
1 files changed, 12 insertions, 11 deletions
diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py index 151dbab..d42d236 100644 --- a/fg21sim/extragalactic/clusters/halo.py +++ b/fg21sim/extragalactic/clusters/halo.py @@ -240,7 +240,7 @@ class RadioHalo1M: The stripping radius of the in-falling sub-cluster at time t. Unit: [kpc] """ - self._validate_t_merger(t_merger) + self._validate_time(t_merger) z = COSMO.redshift(t_merger) M_main = self.mass_main(t_merger) M_sub = self.mass_sub(t_merger) @@ -256,7 +256,7 @@ class RadioHalo1M: Unit: [kpc] """ - self._validate_t_merger(t_merger) + self._validate_time(t_merger) z = COSMO.redshift(t_merger) M_main = self.mass_main(t_merger) r_s = self.radius_strip(t_merger) @@ -274,7 +274,7 @@ class RadioHalo1M: Unit: [Gyr] """ - self._validate_t_merger(t_merger) + self._validate_time(t_merger) z_merger = COSMO.redshift(t_merger) M_main = self.mass_main(t=t_merger) M_sub = self.mass_sub(t=t_merger) @@ -315,7 +315,7 @@ class RadioHalo1M: The turbulence velocity dispersion Unit: [km/s] """ - self._validate_t_merger(t_merger) + self._validate_time(t_merger) z = COSMO.redshift(t_merger) M_main = self.mass_main(t_merger) M_sub = self.mass_sub(t_merger) @@ -343,7 +343,7 @@ class RadioHalo1M: """ The turbulence Mach number determined from its velocity dispersion. """ - self._validate_t_merger(t_merger) + self._validate_time(t_merger) cs = helper.speed_sound(self.kT(t_merger)) # [km/s] v_turb = self.velocity_turb(t_merger) # [km/s] return v_turb / cs @@ -387,7 +387,7 @@ class RadioHalo1M: * Ref.[pinzke2017],Eq.(37) * Ref.[miniati2015],Eq.(29) """ - self._validate_t_merger(t_merger) + self._validate_time(t_merger) k_L = 2 * np.pi / self.radius_turb(t_merger) # [kpc^-1] cs = helper.speed_sound(self.kT(t_merger)) # [km/s] @@ -657,13 +657,14 @@ class RadioHalo1M: """ return self.t_merger - def _validate_t_merger(self, t_merger): + def _validate_time(self, t): """ - Validate that the given time ``t_merger`` is the time when the - merger begins, otherwise raise an error. + Validate that the given time ``t`` is the time when a merger begins + or ends, otherwise raise an error. """ - if not np.any(np.isclose(t_merger, self.t_merger)): - raise ValueError("Not a merger time: %s" % t_merger) + if (not np.any(np.isclose(t, self.t_merger)) and + not np.any(np.isclose(t, self.t_merger_end))): + raise ValueError("Not a merger begin/end time: %f" % t) def mass_merged(self, t=None): """ |