fix integrate_energy

VictorBarbosaMartins · Feb 6, 2024 · 9b489af · 9b489af
1 parent 0ecc198
commit 9b489af
Showing 1 changed file with 5 additions and 1 deletion.
diff --git a/ctao_cosmic_ray_spectra/spectral.py b/ctao_cosmic_ray_spectra/spectral.py
@@ -255,7 +255,11 @@ def integrate_energy(self, energy_min, energy_max):
         nominator = energy_max ** (self.index + 1) - energy_min ** (self.index + 1)
         denominator = (self.index + 1) * self.e_ref**self.index
 
-        return nominator/denominator * self.normalization
+        return PowerLaw(
+                normalization=(nominator/denominator * self.normalization),
+                index=self.index,
+                e_ref=self.e_ref,
+            )
 
     @u.quantity_input(inner=u.deg, outer=u.deg, area=u.cm**2, energy_min=u.TeV, energy_max=u.TeV)
     def derive_events_rate(self,