def maxandmin_intime():
plot(time*1e6, dB(S11c[:, 148]), label="Max refl, yoko=3.15 V".format(yoko[315]))
plot(time*1e6, dB(S11c[:, 200]), label="Min refl, yoko=2.21 V".format(yoko[221]))
xlabel("Time (us)")
ylabel("Reflection (dB)")
title("Time domain reflection at 4.46 GHz")
ylim(-65, -30)
legend()
def maxandmin_intime():
plot(time * 1e6,
dB(S11c[:, 148]),
label="Max refl, yoko=3.15 V".format(yoko[315]))
plot(time * 1e6,
dB(S11c[:, 200]),
label="Min refl, yoko=2.21 V".format(yoko[221]))
xlabel("Time (us)")
ylabel("Reflection (dB)")
title("Time domain reflection at 4.46 GHz")
ylim(-65, -30)
legend()
def VNA_twotone_colormesh():
pcolormesh(yoko2, pwr2[:,0]-fridge_att, dB(Magcom[0, :, :]))
ylim(-187, -87)
xlim(5.0, 6.0)
xlabel("Flux (V)")
ylabel("Control power (dBm)")
title("Reflection (dB) Probe at 4.403 GHz, -137 dBm,\n Control at 4.285 GHz")
def plotmapdBtime():
pcolormesh(yoko, time * 1e6, dB(S11c), vmin=-65, vmax=-30)
title("Reflection (dB) vs flux \n and time (1 us pulse) at 4.46 GHz")
xlabel("Flux (V)")
ylabel("Time (us)")
#ylim(0, 1.5)
colorbar()
def plotmapdBtime():
pcolormesh(yoko, time*1e6, dB(S11c), vmin=-65, vmax=-30)
title("Reflection (dB) vs flux \n and time (1 us pulse) at 4.46 GHz")
xlabel("Flux (V)")
ylabel("Time (us)")
#ylim(0, 1.5)
colorbar()
def VNA_twotonesat():
powsat = mean(dB(Magcom[0, :, 19:31]), axis=1)
print shape(pwr2)
print pwr2[:, 0]
print yoko2[0, :]
print shape(yoko2)
fridge_att = 87.0 + 15.0 * 0.0
pwrlin = 0.001 * 10.0**((pwr2[:, 0] - fridge_att) / 10.0)
powsat_theory = lorentzian(flux_parabola(flux_rescale(yoko2[0, :])),
4.285e9, pwrlin, 50.0e6)
print shape(powsat_theory)
plot(pwr2[:, 0] - fridge_att,
normalize(absolute(powsat_theory[:, 56])),
label="Theory 50 MHz coupling")
plot(pwr2[:, 0] - fridge_att, normalize(powsat),
label="Data") #mean(dB(Magcom[0, :, 19:31]), axis=1))
xlabel("Control power (dBm)")
ylabel("Reflection (normalized)")
title(
"Reflection (normalized) Probe at 4.403 GHz, -137 dBm, \n Control at 4.285 GHz"
)
legend(loc=3)
def cs_dB_meanyoko_refl():
plot(freq, dB(mean(Magcom[:, :], axis=1)))
#plot(freq, dB(Magcom[:, :]))
xlabel("Frequency (Hz)")
ylabel("Reflection (dB)")
title("Reflection")
ylim(-50, -20)
def cs_dB_meanyoko_refl():
plot(freq, dB(mean(Magcom[:, :], axis=1)))
#plot(freq, dB(Magcom[:, :]))
xlabel("Frequency (Hz)")
ylabel("Reflection (dB)")
title("Reflection")
ylim(-50, -20)
def plotdB_colormap(pwi):
pcolormesh(yoko, freq, dB(Magcom[:, pwi, :]))
xlabel("Flux (V)")
ylabel("Frequency (Hz)")
caption="{} dBm".format(pwr[pwi]-fridge_attn)
title("Reflection fluxmap at {}".format(caption))
ylim(amin(freq), amax(freq))
colorbar()
plotdB_colormap.caption=caption
def VNA_twotone_colormesh():
pcolormesh(yoko2, pwr2[:, 0] - fridge_att, dB(Magcom[0, :, :]))
ylim(-187, -87)
xlim(5.0, 6.0)
xlabel("Flux (V)")
ylabel("Control power (dBm)")
title(
"Reflection (dB) Probe at 4.403 GHz, -137 dBm,\n Control at 4.285 GHz"
)
def plotdB_colormap_refl():
pcolormesh(yoko, freq, dB(Magcom[:, :]))
xlabel("Flux (V)")
ylabel("Frequency (Hz)")
caption = "{} dBm".format(-30 - fridge_attn)
title("Reflection fluxmap at {}".format(caption))
ylim(amin(freq), amax(freq))
colorbar()
plotdB_colormap.caption = caption
def time_cuts():
plot(yoko, dB(S11c[35, :]), label="20 ns".format(time[35]-time[30]))
plot(yoko, dB(S11c[46, :]), label="64 ns".format(time[46]-time[30]))
plot(yoko, dB(S11c[83, :]), label="212 ns".format(time[83]-time[30]))
plot(yoko, dB(S11c[90, :]), label="240 ns".format(time[90]-time[30]))
plot(yoko, dB(S11c[100, :]), label="280 ns".format(time[100]-time[30]))
plot(yoko, dB(S11c[110, :]), label="320 ns".format(time[110]-time[30]))
plot(yoko, dB(S11c[153, :]), label="492 ns".format(time[153]-time[30]))
plot(yoko, dB(S11c[230, :]), label="800 ns".format(time[230]-time[30]))
plot(yoko, dB(S11c[294, :]), label="1056 ns".format(time[294]-time[30]))
title("Flux modulation at various times at 4.46 GHz")
xlabel("Flux (V)")
ylabel("Reflection (dB)")
legend()
def time_cuts():
plot(yoko, dB(S11c[35, :]), label="20 ns".format(time[35] - time[30]))
plot(yoko, dB(S11c[46, :]), label="64 ns".format(time[46] - time[30]))
plot(yoko, dB(S11c[83, :]), label="212 ns".format(time[83] - time[30]))
plot(yoko, dB(S11c[90, :]), label="240 ns".format(time[90] - time[30]))
plot(yoko, dB(S11c[100, :]), label="280 ns".format(time[100] - time[30]))
plot(yoko, dB(S11c[110, :]), label="320 ns".format(time[110] - time[30]))
plot(yoko, dB(S11c[153, :]), label="492 ns".format(time[153] - time[30]))
plot(yoko, dB(S11c[230, :]), label="800 ns".format(time[230] - time[30]))
plot(yoko, dB(S11c[294, :]), label="1056 ns".format(time[294] - time[30]))
title("Flux modulation at various times at 4.46 GHz")
xlabel("Flux (V)")
ylabel("Reflection (dB)")
legend()
def VNA_twotonesat():
powsat=mean(dB(Magcom[0, :, 19:31]), axis=1)
print shape(pwr2)
print pwr2[:,0]
print yoko2[0,:]
print shape(yoko2)
fridge_att=87.0+15.0*0.0
pwrlin=0.001*10.0**((pwr2[:,0]-fridge_att)/10.0)
powsat_theory=lorentzian(flux_parabola(flux_rescale(yoko2[0,:])), 4.285e9, pwrlin, 50.0e6)
print shape(powsat_theory)
plot(pwr2[:,0]-fridge_att, normalize(absolute(powsat_theory[:,56])), label="Theory 50 MHz coupling")
plot(pwr2[:,0]-fridge_att, normalize(powsat), label="Data")#mean(dB(Magcom[0, :, 19:31]), axis=1))
xlabel("Control power (dBm)")
ylabel("Reflection (normalized)")
title("Reflection (normalized) Probe at 4.403 GHz, -137 dBm, \n Control at 4.285 GHz")
legend(loc=3)
def cs_dB(fqi, pwi):
plot(yoko, dB(Magcom[fqi, pwi, :]))
print shape(Magcom)
Magcom=squeeze(Magcom)
# print shape(Magcom)
# print shape(yoko)
# print yoko
#print freq
#Magcom=mean(Magcom, axis=1)
powind=4
print pwr[powind]
Magabs=Magcom[:, :, :]-mean(Magcom[:, 197:200, :], axis=1, keepdims=True)
fridge_att=87.0+20.0+5.0
pwrlin=0*0.001*10.0**((pwr[powind]-fridge_att)/10.0)
if 0:
pcolormesh(dB(Magcom[:, :, powind]))
show()
def cs_refl_lowfrq(fqi, pwi):
plot(yoko, normalize(absolute(Magabs[fqi, :, pwi])))
RR=lorentzian(flux_parabola(flux_rescale(yoko, offset=-0.02)), freq[fqi], [pwrlin], g)
#RR=1/(1-1j*detuning(yoko*0.195)/(2.0*pi*10.0e6))
plot(yoko, normalize(absolute(RR)), label="50 MHz {}".format(g))
title("Cross section reflection (normalized) \n at {frq:.3f} GHz and {pwr} dBm".format(frq=freq[fqi]/1.0e9, pwr=pwr[pwi]))
xlabel("Flux (V)")
ylabel("Reflection (normalized)")
def cs_refl_lowfrq_dB(fqi, pwi):
fqi=103
plot(yoko, normalize(dB(Magabs[fqi, :, pwi])))
RR=lorentzian(flux_parabola(flux_rescale(yoko, offset=-0.07)), freq[fqi], [0*pwrlin], g)
def gatecolormesh(pwi):
pcolormesh(yoko, freq, dB(Magcom[:, 0, :]))
xlabel("Flux (V)")
ylabel("Frequency (Hz)")
title("Reflection fluxmap at -117 dBm")
colorbar()
def gatecolormesh_noindex(pwi):
pcolormesh(dB(Magcom[:, pwi, :]))
def cs_refl_lowfrq_dB(fqi, pwi):
fqi=103
plot(yoko, normalize(dB(Magabs[fqi, :, pwi])))
RR=lorentzian(flux_parabola(flux_rescale(yoko, offset=-0.07)), freq[fqi], [0*pwrlin], g)
#RR=1/(1-1j*detuning(yoko*0.195)/(2.0*pi*10.0e6))
plot(yoko, normalize(absolute(RR)), label="50 MHz {}".format(g))
def cm_refl_lowfrq_dB(pwi):
pcolormesh(yoko, freq, dB(Magcom[:, :, pwi]))
ylim(amin(freq), amax(freq))
def cs_dB(fqi, pwi):
plot(yoko, dB(Magcom[fqi, pwi, :]))
def cs_dB_meanyoko_trans():
plot(freq, dB(mean(Magcom[:, :], axis=1)))
xlabel("Frequency (Hz)")
ylabel("Transmission (dB)")
title("Transmission")
ylim(-60, -20)
def plotdBnoindex():
pcolormesh( dB(Magcom[:, 0, :]))
def cs_dB_meanyoko_trans():
plot(freq, dB(mean(Magcom[:, :], axis=1)))
xlabel("Frequency (Hz)")
ylabel("Transmission (dB)")
title("Transmission")
ylim(-60, -20)
def plotdBnoindex():
pcolormesh(dB(Magcom[:, :]))