Fix syst merging and shape systematic underflow
- Systematics that only affect some sub-processes were not properly expanded with "default" values (1 for lnN, nominal yields for shape) to obtain full arrays with all sub-processes.
- With many shape systematics, the interpolation resulted in an underflow from the product of normalization factors over all sources. Now moved the normalization before the product over sources, seems to be more stable.