calculate_sensitivity_lima

lstchain.mc.sensitivity.calculate_sensitivity_lima(n_signal, n_background, alpha)

Sensitivity calculation using the Li & Ma formula eq. 17 of Li & Ma (1983). https://ui.adsabs.harvard.edu/abs/1983ApJ…272..317L/abstract

We calculate the sensitivity in bins of energy and theta2

Parameters:

n_on_events: `numpy.ndarray` number of ON events in the signal region

n_background: `numpy.ndarray` number of events in the background region

alpha: `float` inverse of the number of off positions

n_bins_energy: `int` number of bins in energy

n_bins_theta2: `int` number of bins in theta2

Returns:

sensitivity: numpy.ndarray sensitivity in percentage of Crab units

n_excesses_5sigma: numpy.ndarray number of excesses corresponding to

a 5 sigma significance