def magfilt_cmesh(self, f0=5.35e9, alpha=0.45): Magcom=self.Magcom #(self.Magcom.transpose()-self.Magcom[:, 0]).transpose() fq_vec=self.frequency #array([sqrt(f*(f-2*qdt.call_func("Lamb_shift", f=f, f0=f0, couple_mult=alpha))) for f in self.frequency]) Magfilt=array([fft_filter(Magcom[:,n], self.filt_start_ind, self.filt_end_ind) for n in range(len(self.yoko))]).transpose() Magfilt2=array([fft_filter(Magcom[:,n], 0, 34) for n in range(len(self.yoko))]).transpose() pl=Plotter(fig_width=9.0, fig_height=6.0, name="magabs_{}".format(self.name)) pl, pf=colormesh(self.yoko, fq_vec/1e9, (absolute(Magfilt.transpose()-0.0*Magfilt[:,0])).transpose(), plotter=pl)
def magfilt_cmesh(self, f0=5.35e9, alpha=0.45): Magcom=self.Magcom #(self.Magcom.transpose()-self.Magcom[:, 0]).transpose() fq_vec=self.frequency #array([sqrt(f*(f-2*qdt.call_func("Lamb_shift", f=f, f0=f0, couple_mult=alpha))) for f in self.frequency]) Magfilt=array([fft_filter(Magcom[:,n], self.filt_start_ind, self.filt_end_ind) for n in range(len(self.yoko))]).transpose() Magfilt2=array([fft_filter(Magcom[:,n], 0, 34) for n in range(len(self.yoko))]).transpose() pl=Plotter(fig_width=9.0, fig_height=6.0, name="magabs_{}".format(self.name)) pl, pf=colormesh(self.yoko, fq_vec/1e9, (absolute(Magfilt.transpose()-0.0*Magfilt[:,0])).transpose(), plotter=pl)