angle = np.pi * float(sangle.text) /180
detuning = 2 * np.pi * float(sdetuning.text) * 1e6
transmission = float(stransmission.text)
omega_min = 2 * np.pi * float(sfreqmin.text) * 1e6 #s-1
omega_max = 2 * np.pi * float(sfreqmax.text) * 1e6 #s-1 <- range of sideband frequencies observed
# nb_freq = 1000 # <- freq resolution
omega = np.linspace(omega_min, omega_max, nb_freq)
for k in range(len(phase_mm)):
var_plot = np.zeros(nb_freq)
var_plot_vac = np.zeros(nb_freq)
for j in range(nb_freq):
if dB:
var_plot[j] = variance_dB(cov_ref_ro(omega[j], detuning, transmission, phase_mm[k]), angle)
var_plot_vac[j] = variance_dB(cov_ref_ro_vac(omega[j], detuning, transmission, phase_mm[k]), angle)
else:
var_plot[j] = Sxx(cov_ref_ro(omega[j], detuning, transmission, phase_mm[k]), angle, omega[j])
var_plot_vac[j] = Sxx(cov_ref_ro_vac(omega[j], detuning, transmission, phase_mm[k]), angle, omega[j])
l[k][0].set_ydata(var_plot)
v[k][0].set_ydata(var_plot_vac)
fig.canvas.draw()
sangle.on_submit(update)
sdetuning.on_submit(update)
stransmission.on_submit(update)
sfreqmin.on_submit(update)
sfreqmax.on_submit(update)
plt.show()
