Fix and deprecate histogram filling with `operator+=`

At some point the filling of histograms with h += was broken, probably with !1136 (merged) .

Here I propose a fix and also that we deprecate that interface. It is rather unintuitive and not (yet) used much.