2020ApJ...896....6M

With 10 years of operation and an exceptional data set, the Fermi-Large Area Telescope allows us to unveil the detailed composition of the extragalactic γ-ray sky above 100 MeV. In this paper, we derive the intrinsic source-count distribution (logN-logS) of extragalactic sources (i.e., blazars) at | b| > 20^○ via the efficiency correction method. With this approach, we are able to measure the distribution down to a photon flux of ∼10^–10phcm^–2s^–1 and to an energy flux of ∼10^–12ergcm^–2s^–1. In both cases, the logN-logS becomes flatter at low fluxes. Moreover, we show that this logN-logS is representative of the blazar population (assuming the majority of unassociated sources are blazars) and allows us to constrain its evolution quite effectively. Among recently proposed evolutionary models, we find that the Pure Density Evolution model best describes the evolutionary properties of the blazar population and that their integrated emission accounts for ∼50_–5⁺¹⁰ % of the total extragalactic γ-ray background.