def compared_nikhef_senstech():
plot = Plot(height=r".67\linewidth")
plot_pi_ph(plot, data_nikhef_final_integral, 'square')
plot_pi_ph(plot, data_nikhef_integral, 'x')
plot_pi_ph(plot, data_nikhef_senstech_integral, 'triangle')
plot_pi_ph(plot, data_senstech_integral, 'o')
plot.plot([0, 7 * data_nikhef_final_integral.RATIO], [0, 7],
mark=None,
linestyle='gray')
plot.plot([0, 7 * data_nikhef_integral.RATIO], [0, 7],
mark=None,
linestyle='gray')
plot.plot([0, 7 * data_nikhef_senstech_integral.RATIO], [0, 7],
mark=None,
linestyle='gray')
plot.set_xlabel(r'Multiple-LED pulseintegrals [\si{\nano\volt\second}]')
plot.set_ylabel(r'Multiple-LED pulseheights [\si{\volt}]')
plot.set_xlimits(0, 7 * data_nikhef_final_integral.RATIO)
plot.set_ylimits(0, 7)
plot.save_as_pdf('plots/ph_pi_compared_nikhef_senstech')
LED_PI = [1, 2.43, 3.15, 2.48, 3.18, 3.60, 3.70, 2.72, 5.09, 2.49, 3.66,
3.57, 4.51, 3.15, 3.49, 3.81, 2.64, 4.93, 2.72, 3.13, 3.27,
2.55, 2.07, 3.40, 2.93]
LED_PI_ERR = [0] + [0.1] * 24
MULTI_LED = [((1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24), 5050, 86.50),
((1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22), 4723, 73.61),
((1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20), 4396, 68.58),
((1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18), 3950, 61.43),
((1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16), 3560, 55.46),
((1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14), 3150, 49.00),
((1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12), 2686, 42.50),
((1, 2, 3, 4, 5, 6, 7, 8, 9, 10), 2140, 34.34),
((1, 2, 3, 4, 5, 6, 7, 8), 1746, 28.13),
((1, 2, 3, 4, 5, 6), 1224, 20.36),
((1, 2, 3, 4), 716, 11.43),
((1, 2), 318, 5.00)]
M_PH, M_PH_ERR, M_PI, M_PI_ERR, E_PH, E_PH_ERR, E_PI, E_PI_ERR = get_measured_expected(
LED_PH, LED_PH_ERR, LED_PI, LED_PI_ERR, MULTI_LED)
RATIO = determine_pi_ph_ratio(LED_PI, LED_PH)
if __name__ == '__main__':
name = 'nikhef_final'
plot_ph(E_PH, M_PH, name, E_PH_ERR, M_PH_ERR)
plot_pi(E_PI, M_PI, name, E_PI_ERR, M_PI_ERR)
plot_pi_ph(M_PI, M_PH, name, RATIO, M_PI_ERR, M_PH_ERR)