def plot_eodform(self, ax, filtr): fu_freq = self.fund_freq low_pass = fu_freq*filtr if self.type_detector() == 'pulse': low_pass = 3000. f_eod = aux.butter_lowpass_filter(self.w_eod, low_pass, self._sample_rate) pt, tt, pe, te = self.w_pt plot_inxs = (self.w_time >= tt[1]) & (self.w_time <= tt[3]) ax.plot(1000. * (self.w_time[plot_inxs] - self.w_time[0]), f_eod[plot_inxs], color='green', lw=2) ax.set_ylabel('Amplitude [au]') # "au" stands for arbitrary unit ax.set_xlabel('Time [ms]') pass
def filtr_data(self, filter_fac=5.5): fund_freq = self.fund_freq() filtered_data = aux.butter_lowpass_filter(self._eod, fund_freq * filter_fac, self._sample_rate) return filtered_data