From f6331e3f8548eed6579ad3ed869003a4cb4f8f7a Mon Sep 17 00:00:00 2001
From: Aaron LI <aly@aaronly.me>
Date: Sat, 26 Jan 2019 15:25:54 +0800
Subject: clusters/halo: Some cleanups

---
 fg21sim/extragalactic/clusters/halo.py | 43 ++++++++++++----------------------
 1 file changed, 15 insertions(+), 28 deletions(-)

(limited to 'fg21sim')

diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py
index 706426e..ed05edf 100644
--- a/fg21sim/extragalactic/clusters/halo.py
+++ b/fg21sim/extragalactic/clusters/halo.py
@@ -256,24 +256,6 @@ class RadioHalo1M:
         return helper.radius_stripping(M_main, M_sub, z,
                                        f_rc=self.f_rc, beta=self.beta)
 
-    def calc_radius(self):
-        """
-        The estimated radius of the simulated radio halo.
-        Unit: [kpc]
-        """
-        return self.radius_turbulence(self.age_merger) * self.f_radius
-
-    @lru_cache()
-    def kT(self, t):
-        """
-        The ICM mean temperature of the main cluster.
-        Unit: [keV]
-        """
-        kT_out = self.configs.getn("extragalactic/clusters/kT_out")
-        M_main = self.mass_main(t)
-        z = COSMO.redshift(t)
-        return helper.kT_cluster(mass=M_main, z=z, kT_out=kT_out)
-
     def tau_acceleration(self, t):
         """
         Calculate the electron acceleration timescale due to turbulent
@@ -636,20 +618,26 @@ class RadioHalo1M:
     def mass_main(self, t):
         """
         Calculate the main cluster mass at the given (cosmic) time.
-        Unit: [Msun]
-
-        NOTE
-        ----
-        Since we currently only consider the last major merger event,
-        there may be a long time between ``z_merger`` and ``z_obs``.
-        So we assume that the main cluster grows linearly in time from
+        The main cluster is assumed to grow linearly in time from
         (M_main, z_merger) to (M_obs, z_obs).
+
+        Unit: [Msun]
         """
         t0 = self.age_begin
         rate = (self.M_obs - self.M_main) / (self.age_obs - t0)
         mass = rate * (t - t0) + self.M_main  # [Msun]
         return mass
 
+    def kT(self, t):
+        """
+        The ICM mean temperature of the main cluster.
+        Unit: [keV]
+        """
+        kT_out = self.configs.getn("extragalactic/clusters/kT_out")
+        M_main = self.mass_main(t)
+        z = COSMO.redshift(t)
+        return helper.kT_cluster(mass=M_main, z=z, kT_out=kT_out)
+
     def magnetic_field(self, t):
         """
         Calculate the mean magnetic field strength of the main cluster mass
@@ -713,7 +701,6 @@ class RadioHalo1M:
         rho_main = helper.density_number_thermal(M_main, z)  # [cm^-3]
         rho_main *= AC.mu*AC.u * AUC.g2Msun * AUC.kpc2cm**3  # [Msun/kpc^3]
         R_vir = helper.radius_virial(M_main, z)  # [kpc]
-        r_s = self.radius_stripping(t)  # [kpc]
 
         V_turb = np.pi * r_s**2 * (R_vir+r_s)  # [kpc^3]
         E_turb = self.eta_turb * rho_main * v_i**2 * V_turb
@@ -1020,8 +1007,8 @@ class RadioHalo:
         halos.sort(key=lambda h: h["radius_turb"], reverse=True)
         for hdict in halos:
             halo = hdict["halo"]
-            logger.info("Checking merger: %.2e & %.2e @ %.3f ..." %
-                        (halo.M_main, halo.M_sub, halo.z_merger))
+            logger.info("Checking merger: %.2e & %.2e @ %.3f -> %.3f ..." %
+                        (halo.M_main, halo.M_sub, halo.z_merger, halo.z_obs))
             n_e = halo.calc_electron_spectrum()
             genuine, em_factor = halo.is_genuine(n_e)
             hdict["n_e"] = n_e
-- 
cgit v1.2.2