Shujie Li
--- a/benchmarks/tracking/scripts/tracking_performance_fit.py 0 → 100644
+++ b/benchmarks/tracking/scripts/tracking_performance_fit.py 0 → 100644
+
+    # ax.set_ylabel('Normalized Density / {:d} Bins'.format(nbins))#, fontsize=20)
+    ax.set_xlabel(r'$d\theta$ (rad)')#, fontsize=20)
+
+    # -----------------phi resolution----------------------
+    ax = axs.flat[5]
+    rec_th = df.groupby('event')['phi'].first().values
+    sim_th = dfc.groupby('event')['phi'].first().values
+    dth_th = (rec_th - sim_th)
+
+    # print(np.vstack([rec_th, sim_th]).T)
+#     hval, hbins, _ = ax.hist(dth_th, bins=np.linspace(-0.002, 0.002, 101),
+#                              weights=np.repeat(1./float(rec_th.shape[0]), rec_th.shape[0]), ec='k')
+    sig_ph,err_ph=hist_gaus(dth_th,ax,np.linspace(-ph_lim, ph_lim, nbins+1),klog=0,header=None)
+    # ax.set_ylabel('Normalized Counts / {:d} Bins'.format(nbins))#, fontsize=20)
+    ax.set_xlabel(r'$d\phi$ (rad)')#, fontsize=20)