clusters/main: Add 'genuine' property for halos

Configuration option 'genuine_threshold' is added together.

1 files changed, 9 insertions, 1 deletions

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):
         """