Merge pull request #1164 from cta-observatory/fix_edisp_v0.9

Properly normalize edisp in lstchain v0.9

lstchain/tools/lstchain_create_irf_files.py

@@ -567,6 +567,12 @@ def start(self):
             fov_offset_bins,
             migration_bins,
         )
+        # Different normalization of EDISP had been assumed so far (to the sum of 1).
+        # According to GADF definition, it should be normalized to the integral of 1.
+        # This normalization definition is assumed in pyirf >= v0.10.0.
+        # See https://github.com/cta-observatory/pyirf/pull/250
+        bin_width = np.diff(migration_bins)
+        self.edisp /= bin_width[np.newaxis, :, np.newaxis]
         self.hdus.append(
             create_energy_dispersion_hdu(
                 self.edisp,
@@ -576,6 +582,7 @@ def start(self):
                 point_like=self.point_like,
                 extname="ENERGY DISPERSION",
                 **extra_headers,
+                FIXEDNRM=True,
             )
         )
         self.log.info("Energy Dispersion HDU created")