ctapipe is not stable yet, so expect large and rapid changes to structure and functionality as we explore various design choices before the 1.0 release.

neg_log_likelihood_numeric#

ctapipe.image.pixel_likelihood.neg_log_likelihood_numeric(image, prediction, spe_width, pedestal, confidence=0.999)[source]#

Calculate likelihood of prediction given the measured signal, full numerical integration from [DNR09].

Parameters:
image: ndarray

Pixel amplitudes from image (\(s\)).

prediction: ndarray

Predicted pixel amplitudes from model (\(μ\)).

spe_width: ndarray

Width of single p.e. peak (\(σ_γ\)).

pedestal: ndarray

Width of pedestal (\(σ_p\)).

confidence: float, 0 < x < 1

Upper end of Poisson confidence interval of maximum prediction. Determines upper end of poisson integration.

Returns:
ndarray