A complete characterization of quantum fluctuations in many-body systems is accessible through the full counting statistics. We present an exact computation of statistical properties of light in a basic model of light-matter interaction: a multimode
photonic field coupled to a single two-level emitter. We mostly consider an initial coherent state in a given mode and demonstrate how the original Poissonian statistics gets modified because of quantum many-body scattering effects leading to non-Poissonian distributions. We argue that measuring this statistics in a simple quantum optical setup provides an insight into many-body correlation effects with photons.
