aboutsummaryrefslogtreecommitdiffstats
diff options
context:
space:
mode:
authorAaron LI <aly@aaronly.me>2018-10-24 22:58:42 +0800
committerAaron LI <aly@aaronly.me>2018-10-24 22:58:42 +0800
commitf663daf132a803677262a4be774aa0b74e463902 (patch)
treeb64bb1423cc0d0e8e388821e3819c2a9356b987a
parent8d514dc79b37f5dec4000fff49cbe571ae3a8a69 (diff)
downloadfg21sim-f663daf132a803677262a4be774aa0b74e463902.tar.bz2
clusters/main: Minor improvements
-rw-r--r--fg21sim/extragalactic/clusters/main.py24
1 files changed, 11 insertions, 13 deletions
diff --git a/fg21sim/extragalactic/clusters/main.py b/fg21sim/extragalactic/clusters/main.py
index 0af41ab..0df271c 100644
--- a/fg21sim/extragalactic/clusters/main.py
+++ b/fg21sim/extragalactic/clusters/main.py
@@ -235,21 +235,19 @@ class GalaxyClusters:
def _simulate_halos(self):
"""
- Simulate the radio halo properties for each cluster with recent
- merger event.
+ Simulate the radio halo properties for each cluster that has mergers.
Attributes
----------
halos : list[dict]
- Simulated data for each cluster with recent merger.
+ Simulated data for each cluster with mergers.
"""
- # Select out the clusters with recent mergers
- idx_rmm = [idx for idx, cdict in enumerate(self.catalog)
- if cdict["merger_num"] > 0]
- num = len(idx_rmm)
+ idx_hasmerger = [idx for idx, cdict in enumerate(self.catalog)
+ if cdict["merger_num"] > 0]
+ num = len(idx_hasmerger)
logger.info("Simulating halos for %d clusters with mergers ..." % num)
self.halos = []
- for i, idx in enumerate(idx_rmm):
+ for i, idx in enumerate(idx_hasmerger):
ii = i + 1
if ii % 50 == 0:
logger.info("[%d/%d] %.1f%% ..." % (ii, num, 100*ii/num))
@@ -294,7 +292,7 @@ class GalaxyClusters:
("n_e", n_e), # [cm^-3]
])
self.halos.append(data)
- logger.info("Simulated radio halos for clusters with recent mergers.")
+ logger.info("Simulated radio halos.")
def _calc_halos_emission(self):
"""
@@ -322,18 +320,18 @@ class GalaxyClusters:
("flux", flux), # [Jy]
("Tb_mean", Tb_mean), # [K]
])
- logger.info("Done calculate the radio emissions.")
+ logger.info("Calculated the radio emissions.")
def _dropout_halos(self):
"""
Considering that the (very) massive galaxy clusters are very rare,
while the simulation sky area is rather small, therefore, once a
very massive cluster appears, its associated radio halo is also
- very powerful and (almost) dominate other intermediate/faint halos,
+ very bright and (almost) dominate other intermediate/faint halos,
causing the simulation results unstable and have large variation.
If ``halo_dropout`` is given, then select the specified number of
- most powerful radio halos from the catalog, and mark them with
+ most luminous radio halos from the catalog, and mark them with
a property ``drop=True``, which will then be excluded from the
following halo drawing step, in order to obtain more stable
simulation results.
@@ -380,7 +378,7 @@ class GalaxyClusters:
felong=hdict["felong"],
rotation=hdict["rotation"])
hdict["template"] = template
- logger.info("Done drawn halo template images.")
+ logger.info("Drew halo template images.")
def _save_catalog_data(self, outfile=None, dump=None, clobber=None):
"""