Ok, this is starting to make sense now. For the parameterized case, would there be "smearing" kernels as well, it's just that they're built from parameterizations, rather than simulations? So it would be the exact same code in ```EventRateReco.py``` once you are given either kernels, correct?

What I did/am doing is: create the (4D) kernels from the input hdf5 file in the constructor. Then I have a ```get_kernels``` function that just returns these kernels as a dictionary. Then I'll have the ```EventRateReco``` module calculate the actual reconstucted maps from the kernals in the service. Is this what you had in mind?