Gaussian fit in beams added in comments -- to be adapted....