|
| | imrtgr_imr_consistency_test.preamble |
| |
| | imrtgr_imr_consistency_test.usetex |
| |
| | imrtgr_imr_consistency_test.family |
| |
| | imrtgr_imr_consistency_test.serif |
| |
| | imrtgr_imr_consistency_test.default |
| |
| | imrtgr_imr_consistency_test.markeredgewidth |
| |
| | imrtgr_imr_consistency_test.linewidth |
| |
| | imrtgr_imr_consistency_test.labelsize |
| |
| | imrtgr_imr_consistency_test.fontsize |
| |
| | imrtgr_imr_consistency_test.pad |
| |
| | imrtgr_imr_consistency_test.size |
| |
| | imrtgr_imr_consistency_test.start_time = time.time() |
| |
| | imrtgr_imr_consistency_test.parser = OptionParser() |
| |
| | imrtgr_imr_consistency_test.dest |
| |
| | imrtgr_imr_consistency_test.help |
| |
| | imrtgr_imr_consistency_test.type |
| |
| | imrtgr_imr_consistency_test.options |
| |
| | imrtgr_imr_consistency_test.args |
| |
| | imrtgr_imr_consistency_test.MfafKDE = options.MfafKDE |
| |
| | imrtgr_imr_consistency_test.insp_post = options.insp_post |
| |
| | imrtgr_imr_consistency_test.ring_post = options.ring_post |
| |
| | imrtgr_imr_consistency_test.imr_post = options.imr_post |
| |
| | imrtgr_imr_consistency_test.prior_Mfchif_file = options.prior_Mfchif_file |
| |
| | imrtgr_imr_consistency_test.out_dir = options.out_dir |
| |
| | imrtgr_imr_consistency_test.fit_formula = options.fit_formula |
| |
| | imrtgr_imr_consistency_test.debug_plots = options.debug_plots |
| |
| | imrtgr_imr_consistency_test.pepostproc = options.pepostproc |
| |
| | imrtgr_imr_consistency_test.insp_fhigh = float(options.insp_fhigh) |
| |
| | imrtgr_imr_consistency_test.ring_flow = float(options.ring_flow) |
| |
| | imrtgr_imr_consistency_test.waveform = options.waveform |
| |
| | imrtgr_imr_consistency_test.N_bins = int(options.N_bins) |
| |
| | imrtgr_imr_consistency_test.dMfbyMf_lim = float(options.dMfbyMf_lim) |
| |
| | imrtgr_imr_consistency_test.dchifbychif_lim = float(options.dchifbychif_lim) |
| |
| | imrtgr_imr_consistency_test.lalinference_datadir = os.getenv('LALINFERENCE_DATADIR') |
| |
| | imrtgr_imr_consistency_test.run_command = open('%s/command.txt'%(out_dir),'w') |
| |
| | imrtgr_imr_consistency_test.insp_post_link = insp_post |
| |
| | imrtgr_imr_consistency_test.ring_post_link = ring_post |
| |
| | imrtgr_imr_consistency_test.imr_post_link = imr_post |
| |
| | imrtgr_imr_consistency_test.insp_posplots = os.path.realpath(os.path.dirname(insp_post_link)) |
| |
| | imrtgr_imr_consistency_test.ring_posplots = os.path.realpath(os.path.dirname(ring_post_link)) |
| |
| | imrtgr_imr_consistency_test.imr_posplots = os.path.realpath(os.path.dirname(imr_post_link)) |
| |
| | imrtgr_imr_consistency_test.insp_target = os.path.join(out_dir, 'pe_insp') |
| |
| | imrtgr_imr_consistency_test.ring_target = os.path.join(out_dir, 'pe_ring') |
| |
| | imrtgr_imr_consistency_test.imr_target = os.path.join(out_dir, 'pe_imr') |
| |
| | imrtgr_imr_consistency_test.m1_inj = options.m1_inj |
| |
| | imrtgr_imr_consistency_test.m2_inj = options.m2_inj |
| |
| | imrtgr_imr_consistency_test.chi1_inj = options.chi1_inj |
| |
| | imrtgr_imr_consistency_test.chi2_inj = options.chi2_inj |
| |
| | imrtgr_imr_consistency_test.chi1z_inj = options.chi1z_inj |
| |
| | imrtgr_imr_consistency_test.chi2z_inj = options.chi2z_inj |
| |
| | imrtgr_imr_consistency_test.phi12_inj = options.phi12_inj |
| |
| bool | imrtgr_imr_consistency_test.plot_injection_lines = False |
| |
| | imrtgr_imr_consistency_test.q_inj = m1_inj/m2_inj |
| |
| int | imrtgr_imr_consistency_test.eta_inj = q_inj/(1.+q_inj)**2. |
| |
| | imrtgr_imr_consistency_test.Mf_inj |
| |
| | imrtgr_imr_consistency_test.chif_inj |
| |
| | imrtgr_imr_consistency_test.m1_id |
| | Read the posteriors from the inspiral, ringdown and imr PE runs (after post-processing) More...
|
| |
| | imrtgr_imr_consistency_test.m2_id |
| |
| | imrtgr_imr_consistency_test.a1_id |
| |
| | imrtgr_imr_consistency_test.a2_id |
| |
| | imrtgr_imr_consistency_test.a1z_id |
| |
| | imrtgr_imr_consistency_test.a2z_id |
| |
| | imrtgr_imr_consistency_test.phi12_id |
| |
| | imrtgr_imr_consistency_test.insp_data = np.genfromtxt(insp_post, dtype=None, names=True) |
| |
| | imrtgr_imr_consistency_test.m1_i |
| |
| | imrtgr_imr_consistency_test.m2_i |
| |
| | imrtgr_imr_consistency_test.chi1_i |
| |
| | imrtgr_imr_consistency_test.chi2_i |
| |
| | imrtgr_imr_consistency_test.chi1z_i |
| |
| | imrtgr_imr_consistency_test.chi2z_i |
| |
| | imrtgr_imr_consistency_test.phi12_i = insp_data[phi12_id] |
| |
| | imrtgr_imr_consistency_test.Mf_i |
| |
| | imrtgr_imr_consistency_test.chif_i |
| |
| | imrtgr_imr_consistency_test.ring_data = np.genfromtxt(ring_post, dtype=None, names=True) |
| |
| | imrtgr_imr_consistency_test.m1_r |
| |
| | imrtgr_imr_consistency_test.m2_r |
| |
| | imrtgr_imr_consistency_test.chi1_r |
| |
| | imrtgr_imr_consistency_test.chi2_r |
| |
| | imrtgr_imr_consistency_test.chi1z_r |
| |
| | imrtgr_imr_consistency_test.chi2z_r |
| |
| | imrtgr_imr_consistency_test.phi12_r = ring_data[phi12_id] |
| |
| | imrtgr_imr_consistency_test.Mf_r |
| |
| | imrtgr_imr_consistency_test.chif_r |
| |
| | imrtgr_imr_consistency_test.imr_data = np.genfromtxt(imr_post, dtype=None, names=True) |
| |
| | imrtgr_imr_consistency_test.m1_imr |
| |
| | imrtgr_imr_consistency_test.m2_imr |
| |
| | imrtgr_imr_consistency_test.chi1_imr |
| |
| | imrtgr_imr_consistency_test.chi2_imr |
| |
| | imrtgr_imr_consistency_test.chi1z_imr |
| |
| | imrtgr_imr_consistency_test.chi2z_imr |
| |
| | imrtgr_imr_consistency_test.phi12_imr = imr_data[phi12_id] |
| |
| | imrtgr_imr_consistency_test.Mf_imr |
| |
| | imrtgr_imr_consistency_test.chif_imr |
| |
| | imrtgr_imr_consistency_test.Mf_lim = max(abs(np.append(np.append(Mf_i, Mf_r), Mf_imr))) |
| | compute the limits of integration for computing delta_Mf and delta_chif More...
|
| |
| | imrtgr_imr_consistency_test.chif_lim = max(abs(np.append(np.append(chif_i, chif_r), chif_imr))) |
| |
| | imrtgr_imr_consistency_test.Mf_bins = np.linspace(-Mf_lim, Mf_lim, N_bins) |
| |
| | imrtgr_imr_consistency_test.chif_bins = np.linspace(-chif_lim, chif_lim, N_bins) |
| |
| | imrtgr_imr_consistency_test.dMf = np.mean(np.diff(Mf_bins)) |
| |
| | imrtgr_imr_consistency_test.dchif = np.mean(np.diff(chif_bins)) |
| |
| tuple | imrtgr_imr_consistency_test.Mf_intp = (Mf_bins[:-1] + Mf_bins[1:])/2. |
| |
| tuple | imrtgr_imr_consistency_test.chif_intp = (chif_bins[:-1] + chif_bins[1:])/2. |
| |
| | imrtgr_imr_consistency_test.M_i |
| |
| | imrtgr_imr_consistency_test.C_i |
| |
| | imrtgr_imr_consistency_test.joint_data = gaussian_kde(joint_data) |
| |
| | imrtgr_imr_consistency_test.f_i = lambda x,y:kernel.evaluate([x,y]) |
| |
| | imrtgr_imr_consistency_test.P_Mfchif_i = np.vectorize(f_i)(M_i,C_i)/kernel.integrate_box([-Mf_lim,-chif_lim],[Mf_lim,chif_lim]) |
| |
| | imrtgr_imr_consistency_test.f_r = lambda x,y:kernel.evaluate([x,y]) |
| |
| | imrtgr_imr_consistency_test.P_Mfchif_r = np.vectorize(f_r)(M_i,C_i)/kernel.integrate_box([-Mf_lim,-chif_lim],[Mf_lim,chif_lim]) |
| |
| | imrtgr_imr_consistency_test.f_imr = lambda x,y:kernel.evaluate([x,y]) |
| |
| | imrtgr_imr_consistency_test.P_Mfchif_imr = np.vectorize(f_imr)(M_i,C_i)/kernel.integrate_box([-Mf_lim,-chif_lim],[Mf_lim,chif_lim]) |
| |
| | imrtgr_imr_consistency_test.bins |
| |
| | imrtgr_imr_consistency_test.density |
| |
| | imrtgr_imr_consistency_test.f = gzip.open(prior_Mfchif_file,'rb') |
| | Undo the effect of the prior from the lalinference posterior. More...
|
| |
| | imrtgr_imr_consistency_test.P_Mfchif_pr_interp_obj = pickle.load(f) |
| |
| | imrtgr_imr_consistency_test.P_Mfchif_pr = P_Mfchif_pr_interp_obj(Mf_intp, chif_intp) |
| |
| | imrtgr_imr_consistency_test.P_Mfchif_i_interp_object = scipy.interpolate.interp2d(Mf_intp, chif_intp, P_Mfchif_i, fill_value=0., bounds_error=False) |
| | compute the posterior of (delta_Mf/Mf, delta_chif/chif) More...
|
| |
| | imrtgr_imr_consistency_test.P_Mfchif_r_interp_object = scipy.interpolate.interp2d(Mf_intp, chif_intp, P_Mfchif_r, fill_value=0., bounds_error=False) |
| |
| | imrtgr_imr_consistency_test.dMfbyMf_vec = np.linspace(-dMfbyMf_lim, dMfbyMf_lim, N_bins) |
| |
| | imrtgr_imr_consistency_test.dchifbychif_vec = np.linspace(-dchifbychif_lim, dchifbychif_lim, N_bins) |
| |
| | imrtgr_imr_consistency_test.diff_dMfbyMf = np.mean(np.diff(dMfbyMf_vec)) |
| |
| | imrtgr_imr_consistency_test.diff_dchifbychif = np.mean(np.diff(dchifbychif_vec)) |
| |
| | imrtgr_imr_consistency_test.P_dMfbyMf_dchifbychif = np.zeros(shape=(N_bins,N_bins)) |
| |
| | imrtgr_imr_consistency_test.P_dMfbyMf = np.sum(P_dMfbyMf_dchifbychif, axis=0) * diff_dchifbychif |
| |
| | imrtgr_imr_consistency_test.P_dchifbychif = np.sum(P_dMfbyMf_dchifbychif, axis=1) * diff_dMfbyMf |
| |
| | imrtgr_imr_consistency_test.conf_v1v2 = confidence(P_dMfbyMf_dchifbychif) |
| |
| | imrtgr_imr_consistency_test.gr_height = P_dMfbyMf_dchifbychif[np.argmin(abs(dMfbyMf_vec)), np.argmin(abs(dchifbychif_vec))] |
| |
| | imrtgr_imr_consistency_test.gr_conf_level = conf_v1v2.level_from_height(gr_height) |
| |
| list | imrtgr_imr_consistency_test.param_table |
| |
| | imrtgr_imr_consistency_test.delimiter |
| |
| | imrtgr_imr_consistency_test.fmt |
| |
| | imrtgr_imr_consistency_test.P_m1m2_i = P_m1m2_i.T |
| | plotting More...
|
| |
| | imrtgr_imr_consistency_test.m1_bins_i |
| |
| | imrtgr_imr_consistency_test.m2_bins_i |
| |
| | imrtgr_imr_consistency_test.P_chi1chi2_i = P_chi1chi2_i.T |
| |
| | imrtgr_imr_consistency_test.chi1_bins_i |
| |
| | imrtgr_imr_consistency_test.chi2_bins_i |
| |
| | imrtgr_imr_consistency_test.conf_m1m2_i = confidence(P_m1m2_i) |
| |
| | imrtgr_imr_consistency_test.s1_m1m2_i = conf_m1m2_i.height_from_level(0.68) |
| |
| | imrtgr_imr_consistency_test.s2_m1m2_i = conf_m1m2_i.height_from_level(0.95) |
| |
| | imrtgr_imr_consistency_test.conf_chi1chi2_i = confidence(P_chi1chi2_i) |
| |
| | imrtgr_imr_consistency_test.s1_chi1chi2_i = conf_chi1chi2_i.height_from_level(0.68) |
| |
| | imrtgr_imr_consistency_test.s2_chi1chi2_i = conf_chi1chi2_i.height_from_level(0.95) |
| |
| | imrtgr_imr_consistency_test.conf_Mfchif_i = confidence(P_Mfchif_i) |
| |
| | imrtgr_imr_consistency_test.s1_Mfchif_i = conf_Mfchif_i.height_from_level(0.68) |
| |
| | imrtgr_imr_consistency_test.s2_Mfchif_i = conf_Mfchif_i.height_from_level(0.95) |
| |
| | imrtgr_imr_consistency_test.figsize |
| |
| | imrtgr_imr_consistency_test.cmap |
| |
| | imrtgr_imr_consistency_test.levels |
| |
| | imrtgr_imr_consistency_test.linewidths |
| |
| | imrtgr_imr_consistency_test.x |
| |
| | imrtgr_imr_consistency_test.ls |
| |
| | imrtgr_imr_consistency_test.color |
| |
| | imrtgr_imr_consistency_test.y |
| |
| | imrtgr_imr_consistency_test.dpi |
| |
| | imrtgr_imr_consistency_test.ms |
| |
| | imrtgr_imr_consistency_test.P_m1m2_r = P_m1m2_r.T |
| |
| | imrtgr_imr_consistency_test.m1_bins_r |
| |
| | imrtgr_imr_consistency_test.m2_bins_r |
| |
| | imrtgr_imr_consistency_test.P_chi1chi2_r = P_chi1chi2_r.T |
| |
| | imrtgr_imr_consistency_test.chi1_bins_r |
| |
| | imrtgr_imr_consistency_test.chi2_bins_r |
| |
| | imrtgr_imr_consistency_test.conf_m1m2_r = confidence(P_m1m2_r) |
| |
| | imrtgr_imr_consistency_test.s1_m1m2_r = conf_m1m2_r.height_from_level(0.68) |
| |
| | imrtgr_imr_consistency_test.s2_m1m2_r = conf_m1m2_r.height_from_level(0.95) |
| |
| | imrtgr_imr_consistency_test.conf_chi1chi2_r = confidence(P_chi1chi2_r) |
| |
| | imrtgr_imr_consistency_test.s1_chi1chi2_r = conf_chi1chi2_r.height_from_level(0.68) |
| |
| | imrtgr_imr_consistency_test.s2_chi1chi2_r = conf_chi1chi2_r.height_from_level(0.95) |
| |
| | imrtgr_imr_consistency_test.conf_Mfchif_r = confidence(P_Mfchif_r) |
| |
| | imrtgr_imr_consistency_test.s1_Mfchif_r = conf_Mfchif_r.height_from_level(0.68) |
| |
| | imrtgr_imr_consistency_test.s2_Mfchif_r = conf_Mfchif_r.height_from_level(0.95) |
| |
| | imrtgr_imr_consistency_test.P_m1m2_imr = P_m1m2_imr.T |
| |
| | imrtgr_imr_consistency_test.m1_bins_imr |
| |
| | imrtgr_imr_consistency_test.m2_bins_imr |
| |
| | imrtgr_imr_consistency_test.P_chi1chi2_imr = P_chi1chi2_imr.T |
| |
| | imrtgr_imr_consistency_test.chi1_bins_imr |
| |
| | imrtgr_imr_consistency_test.chi2_bins_imr |
| |
| | imrtgr_imr_consistency_test.conf_m1m2_imr = confidence(P_m1m2_imr) |
| |
| | imrtgr_imr_consistency_test.s1_m1m2_imr = conf_m1m2_imr.height_from_level(0.68) |
| |
| | imrtgr_imr_consistency_test.s2_m1m2_imr = conf_m1m2_imr.height_from_level(0.95) |
| |
| | imrtgr_imr_consistency_test.conf_chi1chi2_imr = confidence(P_chi1chi2_imr) |
| |
| | imrtgr_imr_consistency_test.s1_chi1chi2_imr = conf_chi1chi2_imr.height_from_level(0.68) |
| |
| | imrtgr_imr_consistency_test.s2_chi1chi2_imr = conf_chi1chi2_imr.height_from_level(0.95) |
| |
| | imrtgr_imr_consistency_test.conf_Mfchif_imr = confidence(P_Mfchif_imr) |
| |
| | imrtgr_imr_consistency_test.s1_Mfchif_imr = conf_Mfchif_imr.height_from_level(0.68) |
| |
| | imrtgr_imr_consistency_test.s2_Mfchif_imr = conf_Mfchif_imr.height_from_level(0.95) |
| |
| | imrtgr_imr_consistency_test.CSi = plt.contour(Mf_bins[:-1], chif_bins[:-1], tgr.gf(P_Mfchif_i), levels=(s2_Mfchif_i,s1_Mfchif_i), linewidths=(1,1.5), colors='orange') |
| |
| | imrtgr_imr_consistency_test.CSr = plt.contour(Mf_bins[:-1], chif_bins[:-1], tgr.gf(P_Mfchif_r), levels=(s2_Mfchif_r,s1_Mfchif_r), linewidths=(1,1.5), colors='red') |
| |
| | imrtgr_imr_consistency_test.CSimr = plt.contour(Mf_bins[:-1], chif_bins[:-1], tgr.gf(P_Mfchif_imr), levels=(s2_Mfchif_imr,s1_Mfchif_imr), linewidths=(1,1.5), colors='k') |
| |
| list | imrtgr_imr_consistency_test.strs_i = [ 'inspiral', 'inspiral' ] |
| |
| list | imrtgr_imr_consistency_test.strs_r = [ 'ringdown', 'ringdown' ] |
| |
| list | imrtgr_imr_consistency_test.strs_imr = [ 'IMR', 'IMR' ] |
| |
| dictionary | imrtgr_imr_consistency_test.fmt_i = {} |
| |
| dictionary | imrtgr_imr_consistency_test.fmt_r = {} |
| |
| dictionary | imrtgr_imr_consistency_test.fmt_imr = {} |
| |
| | imrtgr_imr_consistency_test.inline |
| | Label every other level using strings. More...
|
| |
| | imrtgr_imr_consistency_test.True |
| |
| | imrtgr_imr_consistency_test.use_clabeltext |
| |
| | imrtgr_imr_consistency_test.s1_v1v2 = conf_v1v2.height_from_level(0.68) |
| |
| | imrtgr_imr_consistency_test.s2_v1v2 = conf_v1v2.height_from_level(0.95) |
| |
| | imrtgr_imr_consistency_test.conf_v1 = confidence(P_dMfbyMf) |
| |
| | imrtgr_imr_consistency_test.s1_v1 = conf_v1.height_from_level(0.68) |
| |
| | imrtgr_imr_consistency_test.s2_v1 = conf_v1.height_from_level(0.95) |
| |
| | imrtgr_imr_consistency_test.conf_v2 = confidence(P_dchifbychif) |
| |
| | imrtgr_imr_consistency_test.s1_v2 = conf_v2.height_from_level(0.68) |
| |
| | imrtgr_imr_consistency_test.s2_v2 = conf_v2.height_from_level(0.95) |
| |
| | imrtgr_imr_consistency_test.left1_v1 = min(dMfbyMf_vec[np.where(P_dMfbyMf>=s1_v1)[0]]) |
| |
| | imrtgr_imr_consistency_test.right1_v1 = max(dMfbyMf_vec[np.where(P_dMfbyMf>=s1_v1)[0]]) |
| |
| | imrtgr_imr_consistency_test.left2_v1 = min(dMfbyMf_vec[np.where(P_dMfbyMf>=s2_v1)[0]]) |
| |
| | imrtgr_imr_consistency_test.right2_v1 = max(dMfbyMf_vec[np.where(P_dMfbyMf>=s2_v1)[0]]) |
| |
| | imrtgr_imr_consistency_test.left1_v2 = min(dchifbychif_vec[np.where(P_dchifbychif>s1_v2)[0]]) |
| |
| | imrtgr_imr_consistency_test.right1_v2 = max(dchifbychif_vec[np.where(P_dchifbychif>s1_v2)[0]]) |
| |
| | imrtgr_imr_consistency_test.left2_v2 = min(dchifbychif_vec[np.where(P_dchifbychif>s2_v2)[0]]) |
| |
| | imrtgr_imr_consistency_test.right2_v2 = max(dchifbychif_vec[np.where(P_dchifbychif>s2_v2)[0]]) |
| |
| | imrtgr_imr_consistency_test.colspan |
| |
| | imrtgr_imr_consistency_test.lw |
| |
| | imrtgr_imr_consistency_test.rowspan |
| |
| | imrtgr_imr_consistency_test.mew |
| |
1.9.4