prbi, npu = np.broadcast_arrays(q, numevents_per_unit)