From 317f73b60da1562a414d98082f383766d8f1fd04 Mon Sep 17 00:00:00 2001 From: Eric Teunis de Boone Date: Wed, 11 Jan 2023 16:52:05 +0100 Subject: [PATCH] ZH: improve final time comparison plot --- .../bd_antenna_phase_deltas.py | 17 +++++++++++++---- 1 file changed, 13 insertions(+), 4 deletions(-) diff --git a/simulations/airshower_beacon_simulation/bd_antenna_phase_deltas.py b/simulations/airshower_beacon_simulation/bd_antenna_phase_deltas.py index 0d3bcf0..242e3fc 100755 --- a/simulations/airshower_beacon_simulation/bd_antenna_phase_deltas.py +++ b/simulations/airshower_beacon_simulation/bd_antenna_phase_deltas.py @@ -249,17 +249,26 @@ if __name__ == "__main__": if True: # indicate single beacon period span ax.plot((-1, -1), (-1/(2*f_beacon), 1/(2*f_beacon)), marker='3', ms=10, label='1/f_beacon') if i == 0: - ax.plot(np.arange(N_base), actual_baseline_time_shifts, marker='+', label='actual time shifts') + ax.plot(np.arange(N_base), actual_baseline_time_shifts, ls='none', marker='h', alpha=0.8, label='actual time shifts') else: - ax.plot(np.arange(N_base), (actual_baseline_time_shifts+1/(2*f_beacon))%(1/f_beacon) - 1/(2*f_beacon), marker='+', label='actual time shifts') - ax.plot(np.arange(N_base), measured_baseline_time_diffs, marker='x', label='calculated') + ax.plot(np.arange(N_base), (actual_baseline_time_shifts+1/(2*f_beacon))%(1/f_beacon) - 1/(2*f_beacon), ls='none', marker='h', label='actual time shifts', alpha=0.8) + ax.plot(np.arange(N_base), measured_baseline_time_diffs, ls='none', alpha=0.8, marker='x', label='calculated') ax.legend() if fig_dir: extra_name = '' if i == 1: extra_name = '.wrapped' - fig.savefig(path.join(fig_dir, path.basename(__file__) + f".time_comparison{extra_name}.pdf")) + + old_lims = (ax.get_xlim(), ax.get_ylim()) + for j in range(2): + if j == 1: + ax.set_xlim(-5, 50) + extra_name += '.n50' + fig.savefig(path.join(fig_dir, path.basename(__file__) + f".time_comparison{extra_name}.pdf")) + + ax.set_xlim(*old_lims[0]) + ax.set_ylim(*old_lims[1]) if show_plots: plt.show()