def plot(self,save_path=r'D:\analysis\output'): ZZ_corr_err=sscorr.get_correlation_errors(self.ZZ_corr) XX_corr_err=sscorr.get_correlation_errors(self.XX_corr) XmX_corr_err=sscorr.get_correlation_errors(self.XmX_corr) fig=plt.figure() plt.subplot(231) sscorr.plot_uncorrected(self.ZZ_corr, ZZ_corr_err) plt.subplot(232) sscorr.plot_uncorrected(self.XX_corr, XX_corr_err) plt.subplot(233) sscorr.plot_uncorrected(self.XmX_corr, XmX_corr_err) ZZ_c_corr, ZZ_c_corr_err=sscorr.ssro_correct_twoqubit_state_photon_numbers(self.ZZ_corr, self.F0a, self.F0b, self.F1a, self.F1b, dF0a=self.dF0a, dF0b=self.dF0b, dF1a=self.dF1a, dF1b=self.dF1b, verbose=False, return_error_bars=True) XX_c_corr, XX_c_corr_err=sscorr.ssro_correct_twoqubit_state_photon_numbers(self.XX_corr, self.F0a, self.F0b, self.F1a, self.F1b, dF0a=self.dF0a, dF0b=self.dF0b, dF1a=self.dF1a, dF1b=self.dF1b, verbose=False, return_error_bars=True) XmX_c_corr, XmX_c_corr_err=sscorr.ssro_correct_twoqubit_state_photon_numbers(self.XmX_corr, self.F0a, self.F0b, self.F1a, self.F1b, dF0a=self.dF0a, dF0b=self.dF0b, dF1a=self.dF1a, dF1b=self.dF1b, verbose=False, return_error_bars=True) plt.subplot(234) sscorr.plot_corrected(ZZ_c_corr, ZZ_c_corr_err) plt.subplot(235) sscorr.plot_corrected(XX_c_corr, XX_c_corr_err) plt.subplot(236) sscorr.plot_corrected(XmX_c_corr, XmX_c_corr_err) plt.suptitle('Correlations for state '+self.state+' with w_start ' + str(self.w_start) +\ ', w_length ' +str(self.w_length)+ ', w_dt ' + str(self.w_dt) + \ '\n and Fidelity F='+str(self.F)+'+/-'+str(self.dF)) fig.set_size_inches(12,12) fig.savefig(os.path.join(save_path,'fidelity'+'_'+self.state+'.pdf')) return fig
if estimate[0:3]=='mle': c_zz=MLE_psi2zz u_c_zz=u_MLE_psi2zz c_xx=MLE_psi2xx u_c_xx=u_MLE_psi2xx c_xmx=MLE_psi2xmx u_c_xmx=u_MLE_psi2xmx if estimate[-3:]=='yes': uub_c_zz=uub_MLE_psi2zz uub_c_xx=uub_MLE_psi2xx uub_c_xmx=uub_MLE_psi2xmx ulb_c_zz=ulb_MLE_psi2zz ulb_c_xx=ulb_MLE_psi2xx ulb_c_xmx=ulb_MLE_psi2xmx u_zz=sscorr.get_correlation_errors(zz) u_xx=sscorr.get_correlation_errors(xx) u_xmx=sscorr.get_correlation_errors(xmx) if estimate=='ro_correct': c_zz, u_c_zz=sscorr.ssro_correct_twoqubit_state_photon_numbers(zz[::-1], F0a, F0b, F1a, F1b, dF0a=dF0a, dF0b=dF0b, dF1a=dF1a, dF1b=dF1b, verbose=False, return_error_bars=True) c_xx,u_c_xx=sscorr.ssro_correct_twoqubit_state_photon_numbers(xx[::-1], F0a, F0b, F1a, F1b, dF0a=dF0a, dF0b=dF0b, dF1a=dF1a, dF1b=dF1b, verbose=False, return_error_bars=True) c_xmx, u_c_xmx=sscorr.ssro_correct_twoqubit_state_photon_numbers(xmx[::-1], F0a, F0b, F1a, F1b, dF0a=dF0a, dF0b=dF0b, dF1a=dF1a, dF1b=dF1b,