started implementing Gaussian sampling of dn/dlogr
collisional system: use new readout routine