Loading src/lsst_inaf_agile/catalog_agn.py +20 −9 Original line number Diff line number Diff line Loading @@ -165,8 +165,19 @@ class CatalogAGN: ("magabs_" + b.strip() + "_point", np.float64, f"AGN {b} absolute magnitude", "ABmag") for b in bands ] + [(b.strip() + "_tau", np.float64, f"AGN damped random walk timescale {b}", "days") for b in bands] + [(b.strip() + "_sf_inf", np.float64, f"AGN damped random walk structure function at infinity", "ABmag") for b in bands] + [ (b.strip() + "_tau", np.float64, f"AGN damped random walk timescale {b}", "days") for b in bands ] + [ ( b.strip() + "_sf_inf", np.float64, "AGN damped random walk structure function at infinity", "ABmag", ) for b in bands ] ) def get_dtype(self) -> np.dtype: Loading Loading @@ -726,14 +737,14 @@ class CatalogAGN: Return whether the galaxy contains a SMBH according to the occupation fraction. """ U = np.random.rand(self["M"].size) return self["occupation_fraction"] > U return self.get_occupation_fraction() > U def get_is_agn_focc(self): """ Return the joint probability that a galaxy has a BH and the BH is active. """ return self["is_agn"] & self["has_bh"] return self["is_agn"] & self.get_has_bh() def get_lightcurve(self, i, band, *args, **kwargs): """ Loading Loading
src/lsst_inaf_agile/catalog_agn.py +20 −9 Original line number Diff line number Diff line Loading @@ -165,8 +165,19 @@ class CatalogAGN: ("magabs_" + b.strip() + "_point", np.float64, f"AGN {b} absolute magnitude", "ABmag") for b in bands ] + [(b.strip() + "_tau", np.float64, f"AGN damped random walk timescale {b}", "days") for b in bands] + [(b.strip() + "_sf_inf", np.float64, f"AGN damped random walk structure function at infinity", "ABmag") for b in bands] + [ (b.strip() + "_tau", np.float64, f"AGN damped random walk timescale {b}", "days") for b in bands ] + [ ( b.strip() + "_sf_inf", np.float64, "AGN damped random walk structure function at infinity", "ABmag", ) for b in bands ] ) def get_dtype(self) -> np.dtype: Loading Loading @@ -726,14 +737,14 @@ class CatalogAGN: Return whether the galaxy contains a SMBH according to the occupation fraction. """ U = np.random.rand(self["M"].size) return self["occupation_fraction"] > U return self.get_occupation_fraction() > U def get_is_agn_focc(self): """ Return the joint probability that a galaxy has a BH and the BH is active. """ return self["is_agn"] & self["has_bh"] return self["is_agn"] & self.get_has_bh() def get_lightcurve(self, i, band, *args, **kwargs): """ Loading