Commit a711b90e authored by Peter Jansweijer's avatar Peter Jansweijer

removed linear fit and 1st order approximation from "CRTT versus wavelength…

removed linear fit and 1st order approximation from "CRTT versus wavelength figure" and "fiber Delay Coefficient (alpha) as a function of wavelength"
parent 9bc9f91b
......@@ -580,11 +580,11 @@ def analyze_plot(insitu_file, analyse_single, x, y, name, tolerance, use_itu_cha
if use_itu_channels:
ax.set_title("CRTT versus ITU Channel number")
ax.set_xlabel('ITU Channel number')
ax.text(0.01, 0.85, "1st order linear tangent: " + '{:5.3f}'.format(tangent) + " +/- " + '{:5.3f}'.format(tangent_stdev) + " [ps/ITU channel]", transform=ax.transAxes)
#ax.text(0.01, 0.85, "1st order linear tangent: " + '{:5.3f}'.format(tangent) + " +/- " + '{:5.3f}'.format(tangent_stdev) + " [ps/ITU channel]", transform=ax.transAxes)
else:
ax.set_title("CRTT versus wavelength")
ax.set_xlabel('Wavelenth [nm]')
ax.text(0.01, 0.85, "1st order linear tangent: " + '{:5.3f}'.format(tangent) + " +/- " + '{:5.3f}'.format(tangent_stdev) + " [ps/nm]", transform=ax.transAxes)
#ax.text(0.01, 0.85, "1st order linear tangent: " + '{:5.3f}'.format(tangent) + " +/- " + '{:5.3f}'.format(tangent_stdev) + " [ps/nm]", transform=ax.transAxes)
with_tolerances = False
ax.text(0.01, 0.95, str(insitu_file), transform=ax.transAxes)
......@@ -602,7 +602,7 @@ def analyze_plot(insitu_file, analyse_single, x, y, name, tolerance, use_itu_cha
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]
lns2 = ax.plot(x_clean, ly_lin_clean, "-", color = 'blue', dashes=[5,5,5,5], label='linear fit') # final clean linear 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()
ax2.set_ylabel('data versus fit err [ps]')
......@@ -610,7 +610,8 @@ def analyze_plot(insitu_file, analyse_single, x, y, name, tolerance, use_itu_cha
# Calcualte standard deviation of data to fit error and format string
stdev_str = "{0:.3f}".format(data_vs_fit_err_clean.std(ddof=1))
lns3 = ax2.plot(x_clean, data_vs_fit_err_clean, color='grey', label='data vs. 5-term Sellmeier fit error\n StDev: ' + stdev_str + '[ps]')
lns=lns+lns1+lns2+lns3
#lns=lns+lns1+lns2+lns3
lns=lns+lns1+lns3 # Removed linear fit from legend
labels=[l.get_label() for l in lns]
......@@ -637,7 +638,7 @@ def analyze_plot(insitu_file, analyse_single, x, y, name, tolerance, use_itu_cha
ax.axhline(0, color='gray')
ax.axvline(fixed_lambda*1e9, color='red')
ax.plot(x_clean, alpha_sellmeier, color='green', label='dirived from 5-term Sellmeier fit')
ax.plot(x_clean, alpha_tangent, color='blue', dashes=[5,5,5,5], label='dirived from linear fit')
#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')
ax.legend(loc='lower right', fontsize='medium')
#pdb.set_trace()
......
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