Add particle read functionality to dataio. Add histogram for velocities