ctools

I added the comskymap script that performs a Richardson-Lucy deconvolution of a COMPTEL observation. The script handles multiple energy bins and a spectral model that is underlying the reconstructed sky map.

On entry, the script takes a model definition XML file that may come from an earlier model fit. The script then replaces the spatial component of the specified source component, replacing it by an isotropic diffuse map. In an initial fit, the map is fitted together with all other model components to the data. Based on the result of this fitting, a Richardson-Lucy correction map is computed that is then fitted together with the initial map to the data. Both maps (initial map and Richardson-Lucy correction map) share the same spectral model. After model fitting, the map is updated with the fitted version of the Richardson-Lucy correction map and the procedure is repeated.

Out output, the script writes the initial sky map and the maps obtained after each iteration into a FITS file. The script also writes the resulting model definition file, which may be used on input for a new run of comskymap to continue the iterations.