Feature #2662
ctlike/ctulimit tolerance
Status: | Closed | Start date: | 08/15/2018 | |
---|---|---|---|---|
Priority: | Normal | Due date: | ||
Assigned To: | Knödlseder Jürgen | % Done: | 100% | |
Category: | - | |||
Target version: | 1.6.0 | |||
Duration: |
Description
We noticed that for a weak dark matter signal above CR background, ctlike and ctulimit return results too good to be true in comparison with a simple signal/noise estimation for 500 hours of simulated data.
Artificial reduction of time (to have no more than 1e4 photons in 10% broad energy bin) seems to solve the issue. Therefore, the issue connected to numerical accuracy.
It seems that in some cases, especially for the long observational time, default fitting/integration precision is not sufficient.
It will be good to be able to control the fitting/integration tolerance with e.g. the additional hidden parameter for ctlike/ctulimit even by the price of a potentially significant increase of the computational time.
Recurrence
No recurrence.
History
#1 Updated by Knödlseder Jürgen over 6 years ago
- Assigned To set to Knödlseder Jürgen
- Target version set to 1.6.0
We can add the fit tolerance as a hidden parameter. It would be good if you could provide a specific test case that reproduces your problem to see whether fit tolerance is actually the issue.
#2 Updated by Knödlseder Jürgen over 6 years ago
- Status changed from New to Pull request
- % Done changed from 0 to 100
I added a like_accuracy
parameter to all tools and scripts that do maximum likelihood fitting. I will now merge the code into devel
.
#3 Updated by Knödlseder Jürgen over 6 years ago
- Status changed from Pull request to Closed
Code is integrated in devel
branch.