added fixed four lambda combinations as fixed probes

20d82320 · Peter Jansweijer · 54b0cd7b · 20d82320
Commit 20d82320 authored Jun 05, 2018 by Peter Jansweijer
Hide whitespace changes
Inline Side-by-side

Showing with 19 additions and 0 deletions

analyze_sellmeier.py sw/insitu_alpha/analyze_sellmeier.py +19 -0

No files found.
--- a/sw/insitu_alpha/analyze_sellmeier.py
+++ b/sw/insitu_alpha/analyze_sellmeier.py
@@ -794,6 +794,25 @@ def analyze_plot(insitu_file, analyse_single, x, y, name, tolerance, use_itu_cha
      alpha_err_file.write("l1[nm]\tl2[nm]\tcrtt_(l1)[ps]\tcrtt(l2)[ps]\talpha_sel(l1)\talpha_3wl(l1)\tSPEC_alpha_3wl(l1)\ttime_err[ps]\ttime_relative_err[t_err ps/link_length ps]\n")
      alpha_ddelay_file.write("# lambda-1[nm]\tITU-ch\tSPEC_alpha_3wl(l1)\tlambda-2[nm]\n")

+      for i in rng_sel:
+        lambda1 = l2[i[0]][0]
+        lambda2 = l2[i[1]][0]
+        lambda_idx = i[1],i[0]
+        alpha_spec_3wl=spec_alpha(alpha_3wl[lambda_idx])
+        itu_ch = itu_conv.wavelength_2_itu(lambda1/1e9)
+
+        alpha_err_writeln(alpha_err_file, lambda1, lambda2, fixed_lambda, lambda_idx, popt[0], popt[1], popt[2], popt[3], popt[4], crtt_fixed_lambda, alpha_3wl, time_err)
+        alpha_ddelay_file.write(str(lambda1)+", "+str(itu_ch)+", "+str(alpha_spec_3wl)+", "+str(lambda2)+"\n")
+
+      alpha_err_file.write("============================\n")
+      alpha_err_file.write("4 lambda combinations along the full band:\n")
+      alpha_ddelay_file.write("# ============================\n")
+      alpha_ddelay_file.write("# 4 lambda combinations along the full band:\n")
+
+      # fixed selection along the array:
+      rng_sel=([[0,16],[0,32],[0,49],[16,0],[16,32],[16,49],[32,0],[32,16],[32,49],[49,0],[49,16],[49,32]])
+      alpha_err_file.write("l1[nm]\tl2[nm]\tcrtt_(l1)[ps]\tcrtt(l2)[ps]\talpha_sel(l1)\talpha_3wl(l1)\tSPEC_alpha_3wl(l1)\ttime_err[ps]\ttime_relative_err[t_err ps/link_length ps]\n")
+      alpha_ddelay_file.write("# lambda-1[nm]\tITU-ch\tSPEC_alpha_3wl(l1)\tlambda-2[nm]\n")
      for i in rng_sel:
        lambda1 = l2[i[0]][0]
        lambda2 = l2[i[1]][0]