def flux_par4(self, offset=-0.08, flux_factor=0.16, Ejmax=h*44.0e9, f0=5.35e9, alpha=0.7, pl=None):
set_all_tags(qdt, log=False)
flux_o_flux0=flux_over_flux0(self.yoko, offset, flux_factor)
qEj=Ej(Ejmax, flux_o_flux0)
#flux_o_flux0=qdt.call_func("flux_over_flux0", voltage=self.yoko, offset=offset, flux_factor=flux_factor)
freq, frq2=flux_parabola(self.yoko, offset, flux_factor, Ejmax, qdt.Ec)
fq1=lamb_shifted_fq2(qEj/qdt.Ec, qdt.ft, qdt.Np, f0, qdt.epsinf, qdt.W, qdt.Dvv)
line(self.yoko, freq, plotter=pl, linewidth=1.0, alpha=0.5)
line(self.yoko, fq1/2, plotter=pl, linewidth=1.0, alpha=0.5)
def magabs_colormesh(self,
offset=-0.08,
flux_factor=0.52,
Ejmax=h * 44.0e9,
f0=5.35e9,
alpha=0.7,
pl=None):
fq_vec = array([
sqrt(f * (f + alpha * calc_freq_shift(f, qdt.ft, qdt.Np, f0,
qdt.epsinf, qdt.W, qdt.Dvv)))
for f in self.frequency
])
freq, frq2 = flux_parabola(self.yoko, offset, 0.16, Ejmax, qdt.Ec)
pl = Plotter(fig_width=9.0,
fig_height=6.0,
name="magabs_{}".format(self.name))
pl, pf = colormesh(freq,
fq_vec, (self.MagdB.transpose() -
self.MagdB[:, 0]).transpose(),
plotter=pl)
pf.set_clim(-0.3, 0.1)
line([min(freq), max(freq)], [min(freq), max(freq)], plotter=pl)
flux_o_flux0 = flux_over_flux0(self.yoko, offset, flux_factor)
qEj = Ej(Ejmax, flux_o_flux0)
EjdivEc = qEj / qdt.Ec
ls_fq = qdt.call_func("lamb_shifted_fq", EjdivEc=EjdivEc)
ls_fq2 = qdt.call_func("lamb_shifted_fq2", EjdivEc=EjdivEc)
frq2 = qdt.call_func("lamb_shifted_anharm", EjdivEc=EjdivEc) / h
line(ls_fq, ls_fq2, plotter=pl)
#pl.set_xlim(min(self.frequency/1e9), max(self.frequency/1e9))
#pl.set_ylim(min(self.yoko), max(self.yoko))
pl.ylabel = "Yoko (V)"
pl.xlabel = "Frequency (GHz)"
return pl