Feature #874

Test gammalib morphology fitting agains other tools

Added by Deil Christoph almost 11 years ago. Updated over 8 years ago.

Status:NewStart date:05/02/2013
Priority:HighDue date:
Assigned To:Deil Christoph% Done:

0%

Category:-
Target version:-
Duration:

Description

In HESS we have started to check the morphology fitting results from the tools in use (HAP, ParisAnalysis, Sherpa) against each other.
It would be nice if gammalib / ctools would be included in this comparison, and this could also evolve into the OnOffFitting unit tests for gammalib.

Attached find a first simple example, once we have agreement there we can create a series of more complex examples (PSF convolution, non-flat background and exposure maps, On-off-likelihood statistics). We should do the same thing for other parts (statistics, background estimation, spectra) using fake and real data.

HAP has basically no unit test coverage, but the GoogleTest unit testing framework is in place in the scons build system
This gammalib / ctools effort is also a very nice opportunity to define exchangeable maps (basically a bunch of images) and spectrum (PHA?) data formats between the tools and develop some test cases to cross-check the implementations.

Jürgen, Pierrick: maybe you can organize the test cases you’d like in the issue tracker and then I’ll try to run / provide what is needed for the HD / HAP side.

simple.tar (89.5 KB) Deil Christoph, 05/02/2013 06:00 PM


Recurrence

No recurrence.

History

#1 Updated by Knödlseder Jürgen over 9 years ago

  • Priority changed from Normal to High
  • Target version set to 1.0.0

It would be good if someone could do this as soon as possible to validate the code prior to the release of version 1.0.0.

#2 Updated by Knödlseder Jürgen about 9 years ago

Some volunteer to work on that issue?

#3 Updated by Deil Christoph about 9 years ago

  • Assigned To set to Deil Christoph

I’ll do some tests against Sherpa and Gammapy next week.

If I remember correctly it’s not simple to just create a Gaussian or Multi-Gauss or some other simple test PSF ... I have to read the PSF model from FITS files?

#4 Updated by Knödlseder Jürgen over 8 years ago

  • Target version deleted (1.0.0)

I put this now off the 1.0 release. I have done substantial verification of the pull distributions for the various spatial models and everything looks fine. So we should go for the release without having done these tests.

I keep the issue, however, so that we don’t forget about doing these tests one day.

Also available in: Atom PDF