diff --git a/icecube_tools/neutrino_calculator.py b/icecube_tools/neutrino_calculator.py
index 22ba2a8..d3e9a80 100644
--- a/icecube_tools/neutrino_calculator.py
+++ b/icecube_tools/neutrino_calculator.py
@@ -168,17 +168,17 @@ def __call__(self, time=1, min_energy=1e2, max_energy=1e9, min_cosz=-1, max_cosz
 
         # @TODO: Add contribution from bins on boundary.
         self._selected_effective_area_values[
-            self.effective_area.true_energy_bins[1:] < min_energy
+            self.effective_area.true_energy_bins[1:] <= min_energy
         ] = 0
         self._selected_effective_area_values[
-            self.effective_area.true_energy_bins[:-1] > max_energy
+            self.effective_area.true_energy_bins[:-1] >= max_energy
         ] = 0
 
         self._selected_effective_area_values.T[
-            self.effective_area.cos_zenith_bins[1:] < min_cosz
+            self.effective_area.cos_zenith_bins[1:] <= min_cosz
         ] = 0
         self._selected_effective_area_values.T[
-            self.effective_area.cos_zenith_bins[:-1] > max_cosz
+            self.effective_area.cos_zenith_bins[:-1] >= max_cosz
         ] = 0
 
         N = []
@@ -190,10 +190,10 @@ def __call__(self, time=1, min_energy=1e2, max_energy=1e9, min_cosz=-1, max_cosz
 
             self._selected_aeff = self._selected_effective_area_values.copy()
             self._selected_aeff[
-                self.effective_area.true_energy_bins[1:] < src_min_energy
+                self.effective_area.true_energy_bins[1:] <= src_min_energy
             ] = 0
             self._selected_aeff[
-                self.effective_area.true_energy_bins[:-1] > src_max_energy
+                self.effective_area.true_energy_bins[:-1] >= src_max_energy
             ] = 0
 
             if source.source_type == DIFFUSE: