  We derive the Edgeworth streaming model (ESM) for the redshift space
correlation function starting from an arbitrary distribution function for
biased tracers of dark matter by considering its two-point statistics and show
that it reduces to the Gaussian streaming model (GSM) when neglecting
non-Gaussianities. We test the accuracy of the GSM and ESM independent of
perturbation theory using the Horizon Run 2 N-body halo catalog. While the
monopole of the redshift space halo correlation function is well described by
the GSM, higher multipoles improve upon including the leading order
non-Gaussian correction in the ESM: the GSM quadrupole breaks down on scales
below 30 Mpc/h whereas the ESM stays accurate to 2% within statistical errors
down to 10 Mpc/h. To predict the scale dependent functions entering the
streaming model we employ Convolution Lagrangian perturbation theory (CLPT)
based on the dust model and local Lagrangian bias. Since dark matter halos
carry an intrinsic length scale given by their Lagrangian radius, we extend
CLPT to the coarse-grained dust model and consider two different smoothing
approaches operating in Eulerian and Lagrangian space, respectively. The
coarse-graining in Eulerian space features modified fluid dynamics different
from dust while the coarse-graining in Lagrangian space is performed in the
initial conditions with subsequent single streaming dust dynamics, implemented
by smoothing the initial power spectrum in the spirit of the truncated
Zel'dovich approximation. Finally, we compare the predictions of the different
coarse-grained models for the streaming model ingredients to N-body
measurements and comment on the proper choice of both the tracer distribution
function and the smoothing scale. Since the perturbative methods we considered
are not yet accurate enough on small scales, the GSM is sufficient when applied
to perturbation theory.
