diff options
| -rw-r--r-- | fg21sim/configs/config.spec | 4 | ||||
| -rw-r--r-- | fg21sim/extragalactic/clusters/main.py | 10 |
2 files changed, 13 insertions, 1 deletions
diff --git a/fg21sim/configs/config.spec b/fg21sim/configs/config.spec index 2faac91..1bfb27a 100644 --- a/fg21sim/configs/config.spec +++ b/fg21sim/configs/config.spec @@ -445,3 +445,7 @@ stream = option("stderr", "stdout", "", default="stderr") bolo_freq_min = float(default=10) bolo_freq_max = float(default=1e5) bolo_freq_num = integer(default=256) + + # The threshold for the bolometric emissivity acceleration factor to + # determine whether the halo is genuine/formed. + genuine_threshold = float(default=10) diff --git a/fg21sim/extragalactic/clusters/main.py b/fg21sim/extragalactic/clusters/main.py index 11aeaf6..c84ca01 100644 --- a/fg21sim/extragalactic/clusters/main.py +++ b/fg21sim/extragalactic/clusters/main.py @@ -92,6 +92,7 @@ class GalaxyClusters: sec = "extragalactic/halos" self.eta_b = configs.getn(sec+"/x_cr") + self.genuine_threshold = configs.getn(sec+"/genuine_threshold") if self.use_dump_halos_data and (not self.use_dump_catalog_data): self.use_dump_catalog_data = True @@ -303,6 +304,7 @@ class GalaxyClusters: spectrum = halo.calc_electron_spectrum() spectrum_fiducial = halo.calc_electron_spectrum(fiducial=True) factor_acc = halo.calc_acc_factor(spectrum, spectrum_fiducial) + genuine = factor_acc >= self.genuine_threshold theta = halo.radius / (clinfo["DA"]*1e3) * AUC.rad2arcsec # [arcsec] haloinfo = OrderedDict( @@ -314,6 +316,7 @@ class GalaxyClusters: gamma=halo.gamma, # Lorentz factors Ke=halo.injection_rate, # [cm^-3 Gyr^-1] factor_acc=factor_acc, + genuine=genuine, # bool ) return haloinfo @@ -330,14 +333,19 @@ class GalaxyClusters: if cdict["merger_num"] > 0] num = len(idx_hasmerger) logger.info("Simulating halos for %d clusters with mergers ..." % num) + + nhalo = 0 self.halos = [] for i, idx in enumerate(idx_hasmerger): ii = i + 1 if ii % 50 == 0: logger.info("[%d/%d] %.1f%% ..." % (ii, num, 100*ii/num)) haloinfo = self._simulate_halo1(self.catalog[idx]) + nhalo += haloinfo["genuine"] self.halos.append(haloinfo) - logger.info("Simulated radio halos.") + + logger.info("Simulated radio halos: %d (%.1f%%) genuine." % + (nhalo, 100*nhalo/num)) def _calc_halos_emission(self): """ |