Most of histogram have poisson error. would be nice to have something like
hist( xxxxx, errorbar=True) Right now we have to do something similar to
h, e, p = hist(data, bins=100) errorbar((e[1:]+e[:-1])/2., h, sqrt(h))
Few additional point
- a plus if it does asymmetric error bar
- another plus if it can take into account the weight correctly for both using sum(weight^2) or just sum(weight).