From 1703eb5e04f9e35a56e5a764a8b08be1350d5326 Mon Sep 17 00:00:00 2001
From: Aaron LI <aly@aaronly.me>
Date: Sat, 22 Dec 2018 21:15:07 +0800
Subject: clusters/main: Split out _simulate_halo1()
---
fg21sim/extragalactic/clusters/main.py | 98 +++++++++++++++++++---------------
1 file changed, 56 insertions(+), 42 deletions(-)
(limited to 'fg21sim')
diff --git a/fg21sim/extragalactic/clusters/main.py b/fg21sim/extragalactic/clusters/main.py
index f471d90..ba3ecf1 100644
--- a/fg21sim/extragalactic/clusters/main.py
+++ b/fg21sim/extragalactic/clusters/main.py
@@ -228,9 +228,62 @@ class GalaxyClusters:
nmax = max([cdict["merger_num"] for cdict in self.catalog])
logger.info("Maximum number of merger events: %d" % nmax)
+ def _simulate_halo1(self, clinfo):
+ """
+ Calculate the radio halo information for the given cluster.
+
+ Parameters
+ ----------
+ clinfo : dict
+ The cluster information from ``self._simulate_mergers()``.
+
+ Returns
+ -------
+ haloinfo : dict
+ The calculated radio halo information.
+ """
+ merger_num = clinfo["merger_num"]
+ M_obs = clinfo["mass"]
+ z_obs = clinfo["z"]
+ M1 = clinfo["merger_mass1"][merger_num-1]
+ z1 = clinfo["merger_z"][merger_num-1]
+ logger.info("M(%.2e)@z(%.3f) -> M(%.2e)@z(%.3f) with %d merger(s)" %
+ (M1, z1, M_obs, z_obs, merger_num))
+ halo = RadioHaloAM(M_obs=M_obs, z_obs=z_obs,
+ M_main=clinfo["merger_mass1"],
+ M_sub=clinfo["merger_mass2"],
+ z_merger=clinfo["merger_z"],
+ merger_num=merger_num,
+ configs=self.configs)
+ n_e = halo.calc_electron_spectrum()
+
+ return OrderedDict([
+ ("z0", z_obs),
+ ("M0", M_obs), # [Msun]
+ ("age0", halo.age_obs), # [Gyr]
+ ("merger_num", merger_num),
+ ("lon", clinfo["lon"]), # [deg] longitude
+ ("lat", clinfo["lat"]), # [deg] longitude
+ ("felong", clinfo["felong"]), # fraction of elongation
+ ("rotation", clinfo["rotation"]), # [deg] rotation angle
+ ("Rvir0", halo.radius_virial_obs), # [kpc]
+ ("kT0", halo.kT_obs), # [keV]
+ ("B0", halo.B_obs), # [uG] magnetic field @ z_obs
+ ("Rhalo", halo.radius), # [kpc]
+ ("Rhalo_angular", halo.angular_radius), # [arcsec]
+ ("volume", halo.volume), # [kpc^3]
+ ("Ke", halo.injection_rate), # [cm^-3 Gyr^-1]
+ ("time_turbulence", halo.time_turbulence_avg), # [Gyr]
+ ("Mach_turb", halo.mach_turbulence_avg), # Mach number
+ ("tau_acc", halo.tau_acceleration_avg), # [Gyr]
+ ("tfrac_acc", halo.time_acceleration_fraction),
+ ("gamma", halo.gamma), # Lorentz factors
+ ("n_e", n_e), # [cm^-3]
+ ])
+
def _simulate_halos(self):
"""
- Simulate the radio halo properties for each cluster that has mergers.
+ Calculate the radio halo information for each cluster that has mergers.
Attributes
----------
@@ -246,47 +299,8 @@ class GalaxyClusters:
ii = i + 1
if ii % 50 == 0:
logger.info("[%d/%d] %.1f%% ..." % (ii, num, 100*ii/num))
- cdict = self.catalog[idx]
- merger_num = cdict["merger_num"]
- M_obs = cdict["mass"]
- z_obs = cdict["z"]
- M1 = cdict["merger_mass1"][merger_num-1]
- z1 = cdict["merger_z"][merger_num-1]
- info = ("[%d/%d] " % (ii, num) +
- "M(%.2e)@z(%.3f) -> M(%.2e)@z(%.3f) with %d merger(s)" %
- (M1, z1, M_obs, z_obs, merger_num))
- logger.info(info)
- halo = RadioHaloAM(M_obs=M_obs, z_obs=z_obs,
- M_main=cdict["merger_mass1"],
- M_sub=cdict["merger_mass2"],
- z_merger=cdict["merger_z"],
- merger_num=merger_num,
- configs=self.configs)
- n_e = halo.calc_electron_spectrum()
- data = OrderedDict([
- ("z0", z_obs),
- ("M0", M_obs), # [Msun]
- ("age0", halo.age_obs), # [Gyr]
- ("merger_num", merger_num),
- ("lon", cdict["lon"]), # [deg] longitude
- ("lat", cdict["lat"]), # [deg] longitude
- ("felong", cdict["felong"]), # fraction of elongation
- ("rotation", cdict["rotation"]), # [deg] rotation angle
- ("Rvir0", halo.radius_virial_obs), # [kpc]
- ("kT0", halo.kT_obs), # [keV]
- ("B0", halo.B_obs), # [uG] magnetic field @ z_obs
- ("Rhalo", halo.radius), # [kpc]
- ("Rhalo_angular", halo.angular_radius), # [arcsec]
- ("volume", halo.volume), # [kpc^3]
- ("Ke", halo.injection_rate), # [cm^-3 Gyr^-1]
- ("time_turbulence", halo.time_turbulence_avg), # [Gyr]
- ("Mach_turb", halo.mach_turbulence_avg), # Mach number
- ("tau_acc", halo.tau_acceleration_avg), # [Gyr]
- ("tfrac_acc", halo.time_acceleration_fraction),
- ("gamma", halo.gamma), # Lorentz factors
- ("n_e", n_e), # [cm^-3]
- ])
- self.halos.append(data)
+ haloinfo = self._simulate_halo1(self.catalog[idx])
+ self.halos.append(haloinfo)
logger.info("Simulated radio halos.")
def _calc_halos_emission(self):
--
cgit v1.2.2