code: def _get_aeff_corrections(intensity_ratio, mask):
          
          nebins = len(intensity_ratio.data)
          aeff_corrections = np.zeros((nebins))
          for i in range(nebins):
              bright_pixels_intensity = intensity_ratio.data[i][mask.data[i]]
              mean_bright_pixel = bright_pixels_intensity.mean()
              aeff_corrections[i] = 1. / mean_bright_pixel
  
          print("Aeff correction: ", aeff_corrections)
          return aeff_corrections
docs: Compute a correction for the effective area from the brighter pixesl