import scipp as sc
import numpy as np

N = 1000
counts = sc.ones(sizes={'event': N}, dtype=int, unit='counts')
tof = sc.array(dims=['event'], values=np.random.normal(size=N))
da = sc.DataArray(data=counts, coords={'tof': tof})
da.hist(tof=100)