Once identified the convergence region of, one can compute statistical
errors in half of it, following [1]}. The error
on the generating function can be naturally estimated as
\(\mathsf {\mathfrak {\textrm {err}}}[\hat{G}(k,\tau _{b},N_{b})]=\sqrt{\mathrm {var}(\hat{G}(k,\tau _{b},N_{b}))/N_{b}},\) 
