From 9dd5e30488de653d78794210ae40979c9d2bdc01 Mon Sep 17 00:00:00 2001 From: Eric Teunis de Boone Date: Thu, 22 Dec 2022 12:06:34 +0100 Subject: [PATCH] ZH: rotate antenna clock_phase figures --- .../bd_antenna_phase_deltas.py | 29 ++++++++++--------- 1 file changed, 15 insertions(+), 14 deletions(-) diff --git a/simulations/airshower_beacon_simulation/bd_antenna_phase_deltas.py b/simulations/airshower_beacon_simulation/bd_antenna_phase_deltas.py index e45c24c..17cf4f1 100755 --- a/simulations/airshower_beacon_simulation/bd_antenna_phase_deltas.py +++ b/simulations/airshower_beacon_simulation/bd_antenna_phase_deltas.py @@ -106,42 +106,43 @@ if __name__ == "__main__": plot_residuals = i == 1 colors = ['blue', 'orange'] - fig, axs = plt.subplots(1,2, sharey=True) + fig, axs = plt.subplots(2, 1, sharex=True) if True: forward = lambda x: x/(2*np.pi*f_beacon) inverse = lambda x: 2*np.pi*x*f_beacon - secax = axs[-1].secondary_yaxis('right', functions=(forward, inverse)) - secax.set_ylabel('Time $\\Delta\\varphi/(2\\pi f_{beac})$ [ns]') + secax = axs[0].secondary_xaxis('top', functions=(forward, inverse)) + secax.set_xlabel('Time $\\Delta\\varphi/(2\\pi f_{beac})$ [ns]') if plot_residuals: phase_residuals = lib.phase_mod(mean_sigma_phase - actual_phase_shifts) fig.suptitle("Difference between Measured and Actual phases\n for Antenna $i$") - axs[0].set_ylabel("Antenna Phase Residual $\\Delta_\\varphi$") + axs[-1].set_xlabel("Antenna Phase Residual $\\Delta_\\varphi$") else: fig.suptitle("Comparison Measured and Actual phases\n for Antenna $i$") - axs[0].set_ylabel("Antenna Phase $\\Delta_\\varphi$") + axs[-1].set_xlabel("Antenna Phase $\\Delta_\\varphi$") i=0 - axs[i].set_xlabel("Antenna no.") + axs[i].set_ylabel("#") if plot_residuals: - axs[i].plot(np.arange(N_ant), phase_residuals, alpha=0.6, ls='none', marker='x', color=colors[0]) + axs[i].hist(phase_residuals, bins='sqrt', alpha=0.8, color=colors[0]) else: - axs[i].errorbar(np.arange(N_ant), mean_sigma_phase, yerr=std_sigma_phase, marker='4', alpha=0.7, ls='none', color=colors[0], label='Measured') - axs[i].plot(antenna_names, actual_phase_shifts, ls='none', marker='3', alpha=0.8, color=colors[1], label='Actual') - axs[i].legend() + axs[i].hist(mean_sigma_phase, bins='sqrt', density=False, alpha=0.8, color=colors[0], ls='solid' , histtype='step', label='Measured') + axs[i].hist(actual_phase_shifts, bins='sqrt', density=False, alpha=0.8, color=colors[1], ls='dashed', histtype='step', label='Actual') i=1 - axs[i].set_xlabel("#") + axs[i].set_ylabel("Antenna no.") if plot_residuals: - axs[i].hist(phase_residuals, bins='sqrt', alpha=0.8, color=colors[0], orientation='horizontal') + axs[i].plot(phase_residuals, np.arange(N_ant), alpha=0.6, ls='none', marker='x', color=colors[0]) else: - axs[i].hist(mean_sigma_phase, bins='sqrt', density=False, orientation='horizontal', color=colors[0], histtype='step', label='Measured') - axs[i].hist(actual_phase_shifts, bins='sqrt', density=False, orientation='horizontal', ls='dashed', color=colors[1], histtype='step', label='Actual') + axs[i].errorbar(mean_sigma_phase, np.arange(N_ant), yerr=std_sigma_phase, marker='4', alpha=0.7, ls='none', color=colors[0], label='Measured') + axs[i].plot(actual_phase_shifts, antenna_names, ls='none', marker='3', alpha=0.8, color=colors[1], label='Actual') + axs[i].legend() fig.tight_layout() + if fig_dir: extra_name = "measured" if plot_residuals: