|
|
@@ -103,7 +103,7 @@ popt_CB, pcov = curve_fit(CrystalBall, bin_centers, bin_values, sigma=bin_errors
|
|
|
f, (ax1, ax2) = plt.subplots(2, sharex=True, figsize=(16, 12), gridspec_kw={'height_ratios': [3, 1]})
|
|
|
f.suptitle('Fitting inflated simulated signal', fontsize=22)
|
|
|
|
|
|
-ax1.scatter(bin_centers, bin_values, color='k', label='Simulated 10x signal')
|
|
|
+ax1.scatter(bin_centers, bin_values, color='k', label='Simulated 100x signal')
|
|
|
ax1.plot(xs, CrystalBall(xs, *popt_CB), color='r', label='Fitted signal')
|
|
|
ax1.axvline(popt_CB[-2], color='k', ls='--', alpha=0.3)
|
|
|
ax1.annotate(r'$m_{{Higgs}}^{{fit}}= {:.2f}$ GeV'.format(popt_CB[-2]) + '\n' + r'$N_{{Higgs}}^{{exp}} = {:d}$'.format(int(popt_CB[0])), (popt_CB[-2], 10000), xytext=(148, 2000), fontsize=20, arrowprops=dict(facecolor='black', shrink=0.05), size=20, bbox=dict(facecolor='w', edgecolor='gray', boxstyle='round,pad=0.5'))
|
Blaz Leban