Open ayshih opened 9 years ago
The first-order (and possibly most important?) correction would be to add energy resolution. We could do this by multiplying the effective area array by a matrix containing Gaussian energy profiles with a FWHM corresponding to our detector energy resolution (~0.5 keV). This would be a relatively straightforward implementation, but the question of exactly where to implement in the code is maybe not so straightforward. Does it get its own routine? Does it go within the detector efficiency function?
All of the detector efficiency code would need to be refactored to produce matrices instead of just 1D arrays which I think is what Albert was saying was non-trivial to implement.
The current response matrixes consider only the diagonal elements. A full response matrix with off-diagonal elements is needed – how needed ought to be estimated – but this is non-trivial to implement (and to do right).