Fix merging shape systematics across channels
Shape systematics need to know the yield integrals for channels separately.
When systematics didn't affect all processes present, and if several channels were present, those integrals were computed across all channels instead of separately. Solved by propagating the indices of the channels when creating "empty" systematics.