correct slices.slice_constant_charge, center the binning around bunch.mean_dz, implement multiple longitudinal resonators (by passing the resonators as array)
prepare for fixed slices option, some fixes in bunch.slices, change BenchMark_BBresonator_longitudinal.py to precomupute wakefunction from BB resonator and use the interpolated wake table