diff --git a/examples/reflected_region.yml b/examples/reflected_region.yml
index a6a4b92..aa89373 100644
--- a/examples/reflected_region.yml
+++ b/examples/reflected_region.yml
@@ -34,7 +34,7 @@ datasets:
   background:
     method: reflected
 fit:
-  model: pl
+  model: ecpl
   index: 2.5
   reference_energy: 1 TeV
 
diff --git a/v2dl5/light_curves/binary_plotting.py b/v2dl5/light_curves/binary_plotting.py
index 8f35a4a..b3e7660 100644
--- a/v2dl5/light_curves/binary_plotting.py
+++ b/v2dl5/light_curves/binary_plotting.py
@@ -131,6 +131,7 @@ def plot_flux_vs_time(
                     markersize=plotting_utilities.get_marker_size() * 0.5,
                 )
 
+        ax.set_ylim([self.config[0].get("flux_axis_min"), self.config[0].get("flux_axis_max")])
         ax.axhline(0, color="lightgray", linestyle="--")
         plt.xlabel(self._get_time_axis_label(time_axis), fontsize=fontsize)
         if mjd_min is not None and mjd_max is not None:
diff --git a/v2dl5/light_curves/data_reader.py b/v2dl5/light_curves/data_reader.py
index dfcd68a..c95bf81 100644
--- a/v2dl5/light_curves/data_reader.py
+++ b/v2dl5/light_curves/data_reader.py
@@ -192,6 +192,10 @@ def _read_fluxes_from_ecsv_file(
         table = Table.read(file_name)
         f = {}
 
+        if "time" not in table.colnames and "MJD" in table.colnames:
+            table.rename_column("MJD", "time")
+            time_min_max = False
+
         if not time_min_max:
             table["time_min"] = table["time"].data
             table["time_max"] = table["time"].data + 0.1