Skip to content

W
White Rabbit Calibration
Commit c2affb46 authored by Peter Jansweijer's avatar Peter Jansweijer
Browse files
Options

optimized figures readability

parent 7af1f7aa
Expand all Hide whitespace changes
Inline Side-by-side
Showing with 61 additions and 61 deletions
sw/insitu_alpha/analyze_sellmeier.py
View file @ c2affb46
......@@ -783,7 +783,7 @@ def analyze_plot(insitu_file, analyse_single, x, y, name, tolerance, use_itu_cha
with_tolerances = False
ax.text(0.01, 0.95, str(insitu_file), transform=ax.transAxes)
ax.axvline(fixed_lambda*1e9, color='red')
lns = ax.plot(x, y, color = 'blue', label='data') # Original datapoints
lns = ax.plot(x, y, color = 'green', label='data') # Original datapoints
if with_tolerances:
ax.text(0.01, 0.90, "tolerance: " + str(tolerance) + "; outliers: " + str(outlier_cnt), transform=ax.transAxes)
#ax.text(0.01, 0.85, "fit y = a * x + b", transform=ax.transAxes)
......@@ -795,7 +795,7 @@ def analyze_plot(insitu_file, analyse_single, x, y, name, tolerance, use_itu_cha
lns3 = ax.plot(x_clean, y_clean, color = 'red', label='data, outliers removed') # cleaned array in [nm]
lns=lns+lns1+lns2+lns3
lns1 = ax.plot(x_clean, ly_clean, color = 'green', label='5-term Sellmeier fit') # final clean 5th order sellmeier fit in [nm]
lns1 = ax.plot(x_clean, ly_clean, color = 'blue', label='5-term Sellmeier fit') # final clean 5th order sellmeier fit in [nm]
#lns2 = ax.plot(x_clean, ly_lin_clean, "-", color = 'blue', dashes=[5,5,5,5], label='linear fit') # final clean linear fit in [nm]
ax2 = ax.twinx()
......@@ -834,7 +834,7 @@ def analyze_plot(insitu_file, analyse_single, x, y, name, tolerance, use_itu_cha
#ax.text(0.01, 0.90, 'linear and 5-term Sellmeier fit', transform=ax.transAxes)
ax.axhline(0, color='gray')
ax.axvline(fixed_lambda*1e9, color='red')
label = ax.plot(x_clean, alpha_sellmeier, color='green', label='dirived from 5-term Sellmeier fit')
label = ax.plot(x_clean, alpha_sellmeier, color='blue', label='dirived from 5-term Sellmeier fit')
lns.append(label[0])
#ax.plot(x_clean, alpha_tangent, color='blue', dashes=[5,5,5,5], label='dirived from linear fit')
#ax.plot(x_clean, alpha_linear, color='blue', dashes=[5,5,5,5], label='linear')
......@@ -975,25 +975,25 @@ def analyze_plot(insitu_file, analyse_single, x, y, name, tolerance, use_itu_cha
#PLOT_A = False
#if PLOT_A:
lns_a.append(ax.errorbar(selected_wavelength_0, selected_alpha_3wl_a_0, yerr = selected_alpha_3wl_a_0_stddev, fmt='d', capsize=6, color='red', label=r'$\alpha_{3-\lambda}$ A, $\lambda_{2}$ = 1568.77 [nm]'))
lns_a.append(ax.errorbar(selected_wavelength_16, selected_alpha_3wl_a_16, yerr = selected_alpha_3wl_a_16_stddev, fmt='d', capsize=6, color='blue', label=r'$\alpha_{3-\lambda}$ A, $\lambda_{2}$ = 1554.94 [nm]'))
lns_a.append(ax.errorbar(selected_wavelength_32, selected_alpha_3wl_a_32, yerr = selected_alpha_3wl_a_32_stddev, fmt='d', capsize=6, color='grey', label=r'$\alpha_{3-\lambda}$ A, $\lambda_{2}$ = 1542.14 [nm]'))
lns_a.append(ax.errorbar(selected_wavelength_49, selected_alpha_3wl_a_49, yerr = selected_alpha_3wl_a_49_stddev, fmt='d', capsize=6, color='cyan', label=r'$\alpha_{3-\lambda}$ A, $\lambda_{2}$ = 1528.77 [nm]'))
lns_a.append(ax.errorbar(selected_wavelength_7, selected_alpha_3wl_a_7, yerr = selected_alpha_3wl_a_7_stddev, fmt='d', capsize=6, color='orange', label=r'$\alpha_{3-\lambda}$ A, $\lambda_{2}$ = 1563.05 [nm]'))
lns_a.append(ax.errorbar(selected_wavelength_8, selected_alpha_3wl_a_8, yerr = selected_alpha_3wl_a_8_stddev, fmt='d', capsize=6, color='yellow', label=r'$\alpha_{3-\lambda}$ A, $\lambda_{2}$ = 1562.23 [nm]'))
lns_a.append(ax.errorbar(selected_wavelength_9, selected_alpha_3wl_a_9, yerr = selected_alpha_3wl_a_9_stddev, fmt='d', capsize=6, color='brown', label=r'$\alpha_{3-\lambda}$ A, $\lambda_{2}$ = 1560.61 [nm]'))
lns_a.append(ax.errorbar(selected_wavelength_10, selected_alpha_3wl_a_10, yerr = selected_alpha_3wl_a_10_stddev, fmt='d', capsize=6, color='pink', label=r'$\alpha_{3-\lambda}$ A, $\lambda_{2}$ = 1559.79 [nm]'))
lns_a.append(ax.errorbar(selected_wavelength_0, selected_alpha_3wl_a_0, yerr = selected_alpha_3wl_a_0_stddev, fmt='^', capsize=6, color='C1', label=r'$\alpha_{3-\lambda}$ A, $\lambda_{2}$ = 1568.77 [nm]'))
lns_a.append(ax.errorbar(selected_wavelength_16, selected_alpha_3wl_a_16, yerr = selected_alpha_3wl_a_16_stddev, fmt='^', capsize=6, color='C2', label=r'$\alpha_{3-\lambda}$ A, $\lambda_{2}$ = 1554.94 [nm]'))
lns_a.append(ax.errorbar(selected_wavelength_32, selected_alpha_3wl_a_32, yerr = selected_alpha_3wl_a_32_stddev, fmt='^', capsize=6, color='C3', label=r'$\alpha_{3-\lambda}$ A, $\lambda_{2}$ = 1542.14 [nm]'))
lns_a.append(ax.errorbar(selected_wavelength_49, selected_alpha_3wl_a_49, yerr = selected_alpha_3wl_a_49_stddev, fmt='^', capsize=6, color='C4', label=r'$\alpha_{3-\lambda}$ A, $\lambda_{2}$ = 1528.77 [nm]'))
lns_a.append(ax.errorbar(selected_wavelength_7, selected_alpha_3wl_a_7, yerr = selected_alpha_3wl_a_7_stddev, fmt='^', capsize=6, color='C5', label=r'$\alpha_{3-\lambda}$ A, $\lambda_{2}$ = 1563.05 [nm]'))
lns_a.append(ax.errorbar(selected_wavelength_8, selected_alpha_3wl_a_8, yerr = selected_alpha_3wl_a_8_stddev, fmt='^', capsize=6, color='C6', label=r'$\alpha_{3-\lambda}$ A, $\lambda_{2}$ = 1562.23 [nm]'))
lns_a.append(ax.errorbar(selected_wavelength_9, selected_alpha_3wl_a_9, yerr = selected_alpha_3wl_a_9_stddev, fmt='^', capsize=6, color='C7', label=r'$\alpha_{3-\lambda}$ A, $\lambda_{2}$ = 1560.61 [nm]'))
lns_a.append(ax.errorbar(selected_wavelength_10, selected_alpha_3wl_a_10, yerr = selected_alpha_3wl_a_10_stddev, fmt='^', capsize=6, color='C8', label=r'$\alpha_{3-\lambda}$ A, $\lambda_{2}$ = 1559.79 [nm]'))
#else:
lns_b.append(ax.errorbar(selected_wavelength_0, selected_alpha_3wl_b_0, yerr = selected_alpha_3wl_b_0_stddev, fmt='h', capsize=6, color='red', label=r'$\alpha_{3-\lambda}$ B, $\lambda_{2}$ = 1568.77 [nm]'))
lns_b.append(ax.errorbar(selected_wavelength_16, selected_alpha_3wl_b_16, yerr = selected_alpha_3wl_b_16_stddev, fmt='h', capsize=6, color='blue', label=r'$\alpha_{3-\lambda}$ B, $\lambda_{2}$ = 1554.94 [nm]'))
lns_b.append(ax.errorbar(selected_wavelength_32, selected_alpha_3wl_b_32, yerr = selected_alpha_3wl_b_32_stddev, fmt='h', capsize=6, color='grey', label=r'$\alpha_{3-\lambda}$ B, $\lambda_{2}$ = 1542.14 [nm]'))
lns_b.append(ax.errorbar(selected_wavelength_49, selected_alpha_3wl_b_49, yerr = selected_alpha_3wl_b_49_stddev, fmt='h', capsize=6, color='cyan', label=r'$\alpha_{3-\lambda}$ B, $\lambda_{2}$ = 1528.77 [nm]'))
lns_b.append(ax.errorbar(selected_wavelength_7, selected_alpha_3wl_b_7, yerr = selected_alpha_3wl_b_7_stddev, fmt='h', capsize=6, color='orange', label=r'$\alpha_{3-\lambda}$ B, $\lambda_{2}$ = 1563.05 [nm]'))
lns_b.append(ax.errorbar(selected_wavelength_8, selected_alpha_3wl_b_8, yerr = selected_alpha_3wl_b_8_stddev, fmt='h', capsize=6, color='yellow', label=r'$\alpha_{3-\lambda}$ B, $\lambda_{2}$ = 1562.23 [nm]'))
lns_b.append(ax.errorbar(selected_wavelength_9, selected_alpha_3wl_b_9, yerr = selected_alpha_3wl_b_9_stddev, fmt='h', capsize=6, color='brown', label=r'$\alpha_{3-\lambda}$ B, $\lambda_{2}$ = 1560.61 [nm]'))
lns_b.append(ax.errorbar(selected_wavelength_10, selected_alpha_3wl_b_10, yerr = selected_alpha_3wl_b_10_stddev, fmt='h', capsize=6, color='pink', label=r'$\alpha_{3-\lambda}$ B, $\lambda_{2}$ = 1559.79 [nm]'))
lns.append(ax.errorbar(selected_wavelength, selected_alpha_sel, yerr = selected_alpha_sel_stddev, fmt='*', capsize=6, color='green', label='Sellmeier'))
lns_b.append(ax.errorbar(selected_wavelength_0, selected_alpha_3wl_b_0, yerr = selected_alpha_3wl_b_0_stddev, fmt='v', capsize=6, color='C1', label=r'$\alpha_{3-\lambda}$ B, $\lambda_{2}$ = 1568.77 [nm]'))
lns_b.append(ax.errorbar(selected_wavelength_16, selected_alpha_3wl_b_16, yerr = selected_alpha_3wl_b_16_stddev, fmt='v', capsize=6, color='C2', label=r'$\alpha_{3-\lambda}$ B, $\lambda_{2}$ = 1554.94 [nm]'))
lns_b.append(ax.errorbar(selected_wavelength_32, selected_alpha_3wl_b_32, yerr = selected_alpha_3wl_b_32_stddev, fmt='v', capsize=6, color='C3', label=r'$\alpha_{3-\lambda}$ B, $\lambda_{2}$ = 1542.14 [nm]'))
lns_b.append(ax.errorbar(selected_wavelength_49, selected_alpha_3wl_b_49, yerr = selected_alpha_3wl_b_49_stddev, fmt='v', capsize=6, color='C4', label=r'$\alpha_{3-\lambda}$ B, $\lambda_{2}$ = 1528.77 [nm]'))
lns_b.append(ax.errorbar(selected_wavelength_7, selected_alpha_3wl_b_7, yerr = selected_alpha_3wl_b_7_stddev, fmt='v', capsize=6, color='C5', label=r'$\alpha_{3-\lambda}$ B, $\lambda_{2}$ = 1563.05 [nm]'))
lns_b.append(ax.errorbar(selected_wavelength_8, selected_alpha_3wl_b_8, yerr = selected_alpha_3wl_b_8_stddev, fmt='v', capsize=6, color='C6', label=r'$\alpha_{3-\lambda}$ B, $\lambda_{2}$ = 1562.23 [nm]'))
lns_b.append(ax.errorbar(selected_wavelength_9, selected_alpha_3wl_b_9, yerr = selected_alpha_3wl_b_9_stddev, fmt='v', capsize=6, color='C7', label=r'$\alpha_{3-\lambda}$ B, $\lambda_{2}$ = 1560.61 [nm]'))
lns_b.append(ax.errorbar(selected_wavelength_10, selected_alpha_3wl_b_10, yerr = selected_alpha_3wl_b_10_stddev, fmt='v', capsize=6, color='C8', label=r'$\alpha_{3-\lambda}$ B, $\lambda_{2}$ = 1559.79 [nm]'))
lns.append(ax.errorbar(selected_wavelength, selected_alpha_sel, yerr = selected_alpha_sel_stddev, fmt='*', capsize=6, color='blue', label='Sellmeier'))
#############################################################################################
labels = [l.get_label() for l in lns]
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment