diff --git a/final_plots.py b/final_plots.py index 1d79a032f14785c296aa79cb8bf3da0eac73bf84..3cec91f8de04ca78396f9010a84396442786c398 100644 --- a/final_plots.py +++ b/final_plots.py @@ -1514,6 +1514,7 @@ def export_bruhat18_pgfplots( bg1 = np.array([quad_vec((lambda t: splev((vdc1+vac*np.cos(t))*omega1/omega2, g_exp2_bg_spl, ext=3)), 0, np.pi)[0]/np.pi for vac in vac_omega1]) bg1_pat = np.array([quad_vec((lambda t: splev((vdc1_pat+vac*np.cos(t))*omega1/omega2, g_exp2_bg_spl, ext=3)), 0, np.pi)[0]/np.pi for vac in vac_omega1]) + np.savetxt("figdata/bruhat18_19GHz.dat", np.array([vdc1, vdc1*omega1_muV, *g1, *bg1, *(g1 + bg1)]).T, header = " ".join(( @@ -1558,6 +1559,43 @@ def export_bruhat18_pgfplots( fmt = "%.6g", comments = "") + vdc_bg_min = vdc_omega_exp2.min() + vdc_bg_max = vdc_omega_exp2.max() + gfull1_good = [b+g for b,g in zip(bg1, g1)] + gfull1_bad = [b+g for b,g in zip(bg1, g1)] + for i, v in enumerate(vac_omega1): + idx_min = vdc1.searchsorted(vdc_bg_min * 12/19 + v, side="right") + idx_max = vdc1.searchsorted(vdc_bg_max * 12/19 - v, side="left") + gfull1_good[i][:idx_min] = np.nan + gfull1_good[i][idx_max:] = np.nan + gfull1_bad[i][idx_min+1:idx_max-1] = np.nan + gfull2_good = [b+g for b,g in zip(bg2, g2)] + gfull2_bad = [b+g for b,g in zip(bg2, g2)] + for i, v in enumerate(vac_omega2): + idx_min = vdc2.searchsorted(vdc_bg_min + v, side="right") + idx_max = vdc2.searchsorted(vdc_bg_max - v, side="left") + gfull2_good[i][:idx_min] = np.nan + gfull2_good[i][idx_max:] = np.nan + gfull2_bad[i][idx_min+1:idx_max-1] = np.nan + np.savetxt("figdata/bruhat18_19GHz_filtered.dat", + np.array([vdc1, vdc1*omega1_muV, *gfull1_good, *gfull1_bad]).T, + header = " ".join(( + "vdc vdc_muV", + " ".join(f"ggood{v:.3g}" for v in vac_omega1), + " ".join(f"gbad{v:.3g}" for v in vac_omega1), + )), + fmt = "%.6g", + comments = "") + np.savetxt("figdata/bruhat18_12GHz_filtered.dat", + np.array([vdc2, vdc2*omega2_muV, *gfull2_good, *gfull2_bad]).T, + header = " ".join(( + "vdc vdc_muV", + " ".join(f"ggood{v:.3g}" for v in vac_omega2), + " ".join(f"gbad{v:.3g}" for v in vac_omega2), + )), + fmt = "%.6g", + comments = "") + def export_pulse_current_pgfplots(omega=1.5, pulse_duration=0.01): dm = DataManager()