ctools

When computing a residual map using csresmap for a counts cube and a model containing a diffuse map or cube the script exists with the following error:

Traceback (most recent call last):
  File "/usr/local/gamma/bin/csresmap", line 367, in <module>
    app.execute()
  File "/usr/local/gamma/bin/csresmap", line 287, in execute
    self.run()
  File "/usr/local/gamma/bin/csresmap", line 214, in run
    model.run()
  File "/usr/local/gamma/lib/python2.7/site-packages/ctools/tools.py", line 1532, in run
    return _tools.ctmodel_run(self)
IndexError: *** ERROR in GSkyMap::operator(int&, int&): Sky map pixel index -1 is outside the valid range [0,89999].

The error arises due to the new model computation in ctmodel. The diffuse map or cube model computation for a stacked analysis makes use of the GCTAResponseCube class which cashes the diffuse model using the GCTAEventBin::ipix() and GCTAEventBin::ieng() method. Since the values of these indices is not set by ctmodel, but left at -1, an error occurs.

The problem should be solved by correctly setting the indices in ctmodel. This requires adding index setter methods to GCTAEventBin.