Thanks very much for the feedback and advice! That was also my interpretation, that the distribution here is actually a Lorentzian and that each bin or x-value of the plot is Poisson distributed, so I’m glad that the way I’ve implemented it seems correct.