The mass density field in simulated non-Gaussian scenarios