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_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_abs(fqi, pwi):
plot(yoko, normalize(absolute(Magabs[fqi, pwi, :])))
locpwr=pwr[pwi]-fridge_attn-extra_attn
pwrlin=0.001*10.0**((locpwr)/10.0)
RR=lorentzian(flux_parabola(flux_rescale(yoko, offset=-0.01)), 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=locpwr))
xlabel("Flux (V)")
ylabel("Reflection (normalized)")
def cs_abs(fqi, pwi):
plot(yoko, normalize(absolute(Magabs[fqi, pwi, :])))
locpwr = pwr[pwi] - fridge_attn - extra_attn
pwrlin = 0.001 * 10.0**((locpwr) / 10.0)
RR = lorentzian(flux_parabola(flux_rescale(yoko, offset=-0.01)), 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=locpwr))
xlabel("Flux (V)")
ylabel("Reflection (normalized)")
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_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))
yoko = data[:, 0, :].astype(float64)
lvlcorr = data[:, 1, :]
PXIIQ = data[:, 2, :] + 1j * data[:, 3, :]
print PXIIQ.dtype #=PXIIQ.astype(float)
print shape(PXIIQ)
plot(absolute(absolute(PXIIQ) - 0.02760) / 0.000319)
#plot(lvlcorr)
print shape(yoko)
yoko = squeeze(yoko)
fridge_att = 87.0 + 5.0
pwr = array([-50.0])
pwrlin = 0.001 * 10.0**((pwr - fridge_att) / 10.0)
plot(
normalize(
absolute(
lorentzian(flux_parabola(flux_rescale(yoko, offset=0.06)),
4.285e9, pwrlin, 50.0e6))))
show()
file_path = "/Users/thomasaref/Dropbox/Current stuff/Logbook/TA210715A46_cooldown1/Data_1010/TA46_VNA_refl_twotone_powsat.hdf5"
with File(file_path, 'r') as f:
print f["Traces"].keys()
print f.attrs["comment"]
print f["Instrument config"].keys()
ctl_frq = f["Instrument config"][
'PXI Aeroflex 302x Signal Generator - GPIB: PXI3::13::INSTR, PXI SigGen at localhost'].attrs[
"Frequency"] #["Power"]
probe_frq = f["Instrument config"][
'Rohde&Schwarz Network Analyzer - IP: 169.254.107.192, at localhost'].attrs[
"Start frequency"]
print shape(data[:])
#pwr=data[:,0,:]
yoko=data[:,0,:].astype(float64)
lvlcorr=data[:,1,:]
PXIIQ=data[:,2,:]+1j*data[:,3,:]
print PXIIQ.dtype #=PXIIQ.astype(float)
print shape(PXIIQ)
plot(absolute(absolute(PXIIQ)-0.02760)/0.000319)
#plot(lvlcorr)
print shape(yoko)
yoko=squeeze(yoko)
fridge_att=87.0+5.0
pwr=array([-50.0])
pwrlin=0.001*10.0**((pwr-fridge_att)/10.0)
plot(normalize(absolute(lorentzian(flux_parabola(flux_rescale(yoko, offset=0.06)), 4.285e9, pwrlin, 50.0e6))))
show()
file_path="/Users/thomasaref/Dropbox/Current stuff/Logbook/TA210715A46_cooldown1/Data_1010/TA46_VNA_refl_twotone_powsat.hdf5"
with File(file_path, 'r') as f:
print f["Traces"].keys()
print f.attrs["comment"]
print f["Instrument config"].keys()
ctl_frq=f["Instrument config"]['PXI Aeroflex 302x Signal Generator - GPIB: PXI3::13::INSTR, PXI SigGen at localhost'].attrs["Frequency"] #["Power"]
probe_frq=f["Instrument config"]['Rohde&Schwarz Network Analyzer - IP: 169.254.107.192, at localhost'].attrs["Start frequency"]
probe_pwr=f["Instrument config"]['Rohde&Schwarz Network Analyzer - IP: 169.254.107.192, at localhost'].attrs["Output power"]
print probe_frq, probe_pwr, ctl_frq
print f["Data"]["Channel names"][:]
Magvec=f["Traces"]["Rohde&Schwarz Network Analyzer - S12"]#[:]