Better axes. Sort by XGM intensity.

pes_to_spec/test/offline_analysis.py

@@ -107,7 +107,7 @@ def plot_result(filename: str,
     ax.set(title=f"Beam intensity: {intensity*1e-3:.1f} mJ", #avg(stat unc) = {unc_stat}, avg(pca unc) = {unc_pca}",
            xlabel="Photon energy [eV]",
            ylabel="Intensity",
-           ylim=(0, 1.2*Y))
+           ylim=(0, 2*Y))
     if pes is not None:
         ax2 = plt.axes([0,0,1,1])
         # Manually set the position and relative size of the inset axes within ax1
@@ -471,7 +471,7 @@ def main():
     last = min(last-100, pes_raw_t["channel_1_D"].shape[1]-1)
     pes_to_show = 'channel_1_D'
     # plot
-    high_int_idx = np.argsort(xgm_flux_t[:,0])[::-1]
+    high_int_idx = np.argsort(xgm_flux_t[:,0])
     for idx in high_int_idx[-10:]:
         tid = test_tids[idx]
         plot_result(os.path.join(args.directory, f"test_{tid}.png"),