Evolution of the BS maximum significance, ${\sigma }_{\max }^{\mathrm{BS}}$, with E_th and ψ, for the baseline astrophysical scenario with source composition model A, 1 nG EGMF and the Sun+Planck GMF model with λ_c = 100 pc. Top left: evolution of ${\sigma }_{\max }^{\mathrm{BS}}$ with E_th for a smoothing angle ψ = 27°. The thick blue dashed line shows the evolution averaged over 300 datasets. The thin lines correspond to 10 individual datasets chosen randomly. The light blue shaded area shows the 90% interval of our simulated datasets. The brown shaded area shows the 99% interval for isotropic datasets. Top right: evolution of ${\sigma }_{\max }^{\mathrm{BS}}$ with ψ, for a threshold energy of E_th = 40 EeV. Bottom panels: same as the top panels, for a wider range of the E_th and ψ scans.