tested all the binning function for dN/dtau, dN/dx1, dN/dx2, and dN/detas
add a small random number on the spatial distribution of the Monte-Carlo samples to make binning more stable