What I wanted to do was to:
not add bins in order to find the +/- 34.1% errors to the red line, which can be binning dependent.
but instead fit a function on the distribution, and integrate that function around the red line to get the +/-34.1%.
Obviously the distribution is not Gaussian, but looks more like a Poisson...
I was thinking about rescaling the x-axis [since Poisson is accepting integer entries while I have floats], doing the fit, integrating, and then scale everything [together with the obtained variances] back to the original x-axis.
The point is that I have other distributions which look pretty much like Gaussians, and I wanted to make sure I could use a Poisson to fit them too [since the code should do that, I wouldn't want to check everytime the distribution and determine with what I could fit it with].