Convolve exposure in ExcessMapEstimator
@@ -230,6 +230,7 @@ class ExcessMapEstimator(Estimator):
if dataset.exposure:
reco_exposure = estimate_exposure_reco_energy(dataset, self.spectral_model)
with np.errstate(invalid="ignore", divide="ignore"):
+ reco_exposure = reco_exposure.convolve(kernel.array) / mask.convolve(kernel.array)
flux = excess / reco_exposure
flux.quantity = flux.quantity.to("1 / (cm2 s)")
else:
