pl=pl88,
auto_ylim=False,
y_min=-40,
y_max=0,
auto_xlim=False,
x_min=4.2,
x_max=4.7,
xlabel="Frequency (GHz)",
ylabel="Transmission (dB)",
linewidth=1.0,
) #.show()
pl88.axes.set_xticks(linspace(4.2, 4.7, 3))
pl88.axes.set_yticks(linspace(-40, 0, 5))
d0317.read_data()
cdata = 10.0**(
(20.0 * log10(mean(absolute(d0317.MagcomData[64:76, :, 0]), axis=0)) -
d0317.probe_pwr) / 20.0)
pl88.axes.plot(d0317.frequency / 1e9,
20 * log10(cdata) + 10 - 1.0 - bg_A1(d0317.frequency),
"x",
10,
color="green")
pl88.nplot = 2
c.read_data()
c.filter_type = "FFT"
print c.MagcomFilt.shape
cdata = (20 * log10(absolute(c.MagcomFilt[:, :])).transpose() -
y_min=-40,
y_max=-10,
auto_xlim=False,
x_min=4e9,
x_max=5e9,
xlabel="Frequency (Hz)",
ylabel="Transmission (dB)",
title="Pickup IDT (37 fingers)") #.show()
line(a.frequency, angle(magfilt))
ifft_plot(a) #.show()
from D0317_S4A1_frq_pulse_flux import a as c
from numpy import mean
c.read_data()
pl = scatter(
c.frequency,
10 * log10(mean(absolute(c.MagcomData[64:76, :, 0]), axis=0)) + 8.5,
pl=pl,
color="purple",
marker="x")
#a.save_plots([pl,])
pl.show()
if __name__ == "__main__":
#b.filt_compare(a.on_res_ind)
print a.net_loss, a.rt_atten
#pl=a.magabs_colormesh()#magabs_colormesh3(s3a4_wg)
auto_ylim=False, y_min=-40, y_max=-10*0,
auto_xlim=False, x_min=4e9, x_max=5e9, xlabel="Frequency (Hz)", ylabel="Transmission (dB)",
title="Pickup IDT (37 fingers)")#.show()
line(a.frequency, angle(magfilt))
X=idt.Np*pi*(a.frequency-idt.f0)/idt.f0
#line(a.frequency, 20*log10(0.5*(sin(X)/(X))**2), pl=pl)
ifft_plot(a)#.show()
from D0317_S4A1_frq_pulse_flux import a as c
from numpy import mean
c.read_data()
print c.comment
cdata=10.0**((20.0*log10(mean(absolute(c.MagcomData[64:76, :, 0]), axis=0))-c.probe_pwr)/20.0)
pl=scatter(c.frequency, 20*log10(cdata)+10-1.0-bg_A1(c.frequency),
pl=pl, color="purple", marker="x")
#pl=scatter(c.frequency, 20*log10(mean(absolute(c.MagcomData[64:76, :, 0]), axis=0))+0*8.5-bg(c.frequency),
# pl=pl, color="purple", marker="x")
#a.save_plots([pl,])
pl.show()
if __name__=="__main__":
#b.filt_compare(a.on_res_ind)
print a.net_loss, a.rt_atten
print idt88.Np, idt88.K2
print idt88.f0
print d0527.comment
print -d0527.fridge_atten+d0527.fridge_gain-d0527.rt_atten+d0527.rt_gain-10
line(d0527.frequency/1e9, 20*log10(absolute(S13xS31))-4, color="red", pl=pl88,
auto_ylim=False, y_min=-40, y_max=0,
auto_xlim=False, x_min=4.2, x_max=4.7, xlabel="Frequency (GHz)",
ylabel="Transmission (dB)", linewidth=1.0,
)#.show()
pl88.axes.set_xticks(linspace(4.2, 4.7, 3))
pl88.axes.set_yticks(linspace(-40, 0, 5))
d0317.read_data()
cdata=10.0**((20.0*log10(mean(absolute(d0317.MagcomData[64:76, :, 0]), axis=0))-d0317.probe_pwr)/20.0)
pl88.axes.plot(d0317.frequency/1e9,
20*log10(cdata)+10-1.0-bg_A1(d0317.frequency), "x", 10,
color="green")
pl88.nplot=2
c.read_data()
c.filter_type="FFT"
print c.MagcomFilt.shape
cdata=(20*log10(absolute(c.MagcomFilt[:,:])).transpose()-bg_A4(c.frequency)).transpose()
flux_axis=c.flux_axis[c.flat_flux_indices]
freq_axis=c.freq_axis[c.indices]
#pl, pf=colormesh(flux_axis, freq_axis, (cdata[d1118.end_skip:-d1118.end_skip, :]/1.0),