Action #1257

Implement upperlimit calculation

Added by Mayer Michael almost 10 years ago. Updated about 8 years ago.

Status:ClosedStart date:07/09/2014
Priority:NormalDue date:
Assigned To:-% Done:

0%

Category:-
Target version:-
Duration:

Description

To calculate spectral points or analyse weak sources, it is necessary to compute upper limits to properly handle the lack of statistics. Since ctlike provides full access to the likelihood function, we can calculate upper limits for any parameter for any source in the xml model.
In general, the upper limit value can be derived by inspecting the likelihood function and find the parameter value at which the likelihood changes by a certain value. For 95% confidence level, this value is, e.g., 1.92 (see Minuit manual)

The question remains, wether this should be implemented as a ctlike option, which could also add the UL value to the final xml model, or if we want to have a separate tool for that.

I attach a first approach using scipy to compute the upper limit value in a cs-script.

csupperlimit.py Magnifier (9.17 KB) Mayer Michael, 07/09/2014 02:43 PM

test_csupperlimit.py Magnifier (554 Bytes) Mayer Michael, 07/09/2014 02:43 PM


Recurrence

No recurrence.


Related issues

Duplicated by ctools - Feature #1363: Create ctulimit Closed 11/12/2014

History

#1 Updated by Mayer Michael almost 10 years ago

The script test_csupperlimit.py simulates a Crab observation and computes the upper limit of the spectral prefactor of that source.

#2 Updated by Mayer Michael about 9 years ago

done by related task #1363

#3 Updated by Knödlseder Jürgen about 8 years ago

  • Status changed from New to Closed
  • Remaining (hours) set to 0.0

Also available in: Atom PDF