diff --git a/simulations/airshower_beacon_simulation/bd_antenna_phase_deltas.py b/simulations/airshower_beacon_simulation/bd_antenna_phase_deltas.py
index d484c8c..6e088be 100755
--- a/simulations/airshower_beacon_simulation/bd_antenna_phase_deltas.py
+++ b/simulations/airshower_beacon_simulation/bd_antenna_phase_deltas.py
@@ -102,7 +102,7 @@ if __name__ == "__main__":
     if True:
         # show means and std
         fig, axs = plt.subplots(1,2, sharey=True)
-        axs[0].set_title("")
+        fig.suptitle("Comparison Measured and Actual Phases")
         axs[0].set_xlabel("Antenna no.")
         axs[0].set_ylabel("Antenna Phase $\\Delta_\\varphi$")
         axs[1].set_xlabel("#")
@@ -112,17 +112,17 @@ if __name__ == "__main__":
             secax = axs[1].secondary_yaxis('right', functions=(forward, inverse))
             secax.set_ylabel('Time $\\Delta\\varphi/(2\\pi f_{beac})$ [ns]')
 
-        l = axs[0].errorbar(np.arange(N_ant), mean_sigma_phase, yerr=std_sigma_phase, marker='.', alpha=0.7)
+        l = axs[0].errorbar(np.arange(N_ant), mean_sigma_phase, yerr=std_sigma_phase, marker='4', alpha=0.7, ls='none')
 
         axs[1].hist(mean_sigma_phase, bins='sqrt', density=False, orientation='horizontal', color=l[0].get_color(), histtype='step')
 
         # Actual time shifts
         if True:
             actual_phase_shifts = [ -1*lib.phase_mod(2*np.pi*f_beacon*v) for k,v in antenna_time_shifts.items() ]
-            antenna_names = [int(k) for k,v in antenna_time_shifts.items() ]
+            antenna_names = [int(k)-1 for k,v in antenna_time_shifts.items() ]
 
             # make sure to keep the same offset
-            if True:
+            if False:
                 phase_offset = mean_sigma_phase[0] - actual_phase_shifts[0]
 
                 actual_phase_shifts += phase_offset