In the two-body scattering problem in general relativity, we study the final graviton particle distribution using a perturbative approach. We compute the mean, the variance and the factorial moments of the distribution from the expectation value of the graviton number operator in the KMOC formalism. For minimally coupled scalar particles, the leading deviation from the Poissonian distribution is given by the unitarity cut involving the six-point tree amplitude with the emission of two gravitons. We compute this amplitude in two independent ways. First, we use an extension of the Cheung-Remmen parametrization that includes minimally coupled scalars. We then repeat the calculation using on-shell BCFW-like techniques, finding complete agreement. In the classical limit, this amplitude gives a purely quantum contribution, proving that we can describe the final radiative state is coherent at least up to order O(G 4 ). Finally, we give general arguments about why we expect this to hold also at higher order in perturbation theory and what are the consequences for classical radiative observables.