Data Structures
class	confidence

Functions
def	set_tick_sizes (ax, major, minor)

Variables
	preamble

	usetex

	family

	serif

	default

	markeredgewidth

	linewidth

	labelsize

	fontsize

	pad

	size

	start_time = time.time()

	parser = OptionParser()

	dest

	help

	type

	options

	args

	MfafKDE = options.MfafKDE

	insp_post = options.insp_post

	ring_post = options.ring_post

	imr_post = options.imr_post

	prior_Mfchif_file = options.prior_Mfchif_file

	out_dir = options.out_dir

	fit_formula = options.fit_formula

	debug_plots = options.debug_plots

	pepostproc = options.pepostproc

	insp_fhigh = float(options.insp_fhigh)

	ring_flow = float(options.ring_flow)

	waveform = options.waveform

	N_bins = int(options.N_bins)

	dMfbyMf_lim = float(options.dMfbyMf_lim)

	dchifbychif_lim = float(options.dchifbychif_lim)

	lalinference_datadir = os.getenv('LALINFERENCE_DATADIR')

	run_command = open('%s/command.txt'%(out_dir),'w')

	insp_post_link = insp_post

	ring_post_link = ring_post

	imr_post_link = imr_post

	insp_posplots = os.path.realpath(os.path.dirname(insp_post_link))

	ring_posplots = os.path.realpath(os.path.dirname(ring_post_link))

	imr_posplots = os.path.realpath(os.path.dirname(imr_post_link))

	insp_target = os.path.join(out_dir, 'pe_insp')

	ring_target = os.path.join(out_dir, 'pe_ring')

	imr_target = os.path.join(out_dir, 'pe_imr')

	m1_inj = options.m1_inj

	m2_inj = options.m2_inj

	chi1_inj = options.chi1_inj

	chi2_inj = options.chi2_inj

	chi1z_inj = options.chi1z_inj

	chi2z_inj = options.chi2z_inj

	phi12_inj = options.phi12_inj

bool	plot_injection_lines = False

	q_inj = m1_inj/m2_inj

int	eta_inj = q_inj/(1.+q_inj)**2.

	Mf_inj

	chif_inj

	m1_id
	Read the posteriors from the inspiral, ringdown and imr PE runs (after post-processing) More...

	m2_id

	a1_id

	a2_id

	a1z_id

	a2z_id

	phi12_id

	insp_data = np.genfromtxt(insp_post, dtype=None, names=True)

	m1_i

	m2_i

	chi1_i

	chi2_i

	chi1z_i

	chi2z_i

	phi12_i = insp_data[phi12_id]

	Mf_i

	chif_i

	ring_data = np.genfromtxt(ring_post, dtype=None, names=True)

	m1_r

	m2_r

	chi1_r

	chi2_r

	chi1z_r

	chi2z_r

	phi12_r = ring_data[phi12_id]

	Mf_r

	chif_r

	imr_data = np.genfromtxt(imr_post, dtype=None, names=True)

	m1_imr

	m2_imr

	chi1_imr

	chi2_imr

	chi1z_imr

	chi2z_imr

	phi12_imr = imr_data[phi12_id]

	Mf_imr

	chif_imr

	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...

	chif_lim = max(abs(np.append(np.append(chif_i, chif_r), chif_imr)))

	Mf_bins = np.linspace(-Mf_lim, Mf_lim, N_bins)

	chif_bins = np.linspace(-chif_lim, chif_lim, N_bins)

	dMf = np.mean(np.diff(Mf_bins))

	dchif = np.mean(np.diff(chif_bins))

tuple	Mf_intp = (Mf_bins[:-1] + Mf_bins[1:])/2.

tuple	chif_intp = (chif_bins[:-1] + chif_bins[1:])/2.

	M_i

	C_i

	joint_data = gaussian_kde(joint_data)

	f_i = lambda x,y:kernel.evaluate([x,y])

	P_Mfchif_i = np.vectorize(f_i)(M_i,C_i)/kernel.integrate_box([-Mf_lim,-chif_lim],[Mf_lim,chif_lim])

	f_r = lambda x,y:kernel.evaluate([x,y])

	P_Mfchif_r = np.vectorize(f_r)(M_i,C_i)/kernel.integrate_box([-Mf_lim,-chif_lim],[Mf_lim,chif_lim])

	f_imr = lambda x,y:kernel.evaluate([x,y])

	P_Mfchif_imr = np.vectorize(f_imr)(M_i,C_i)/kernel.integrate_box([-Mf_lim,-chif_lim],[Mf_lim,chif_lim])

	bins

	density

	f = gzip.open(prior_Mfchif_file,'rb')
	Undo the effect of the prior from the lalinference posterior. More...

	P_Mfchif_pr_interp_obj = pickle.load(f)

	P_Mfchif_pr = P_Mfchif_pr_interp_obj(Mf_intp, chif_intp)

	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...

	P_Mfchif_r_interp_object = scipy.interpolate.interp2d(Mf_intp, chif_intp, P_Mfchif_r, fill_value=0., bounds_error=False)

	dMfbyMf_vec = np.linspace(-dMfbyMf_lim, dMfbyMf_lim, N_bins)

	dchifbychif_vec = np.linspace(-dchifbychif_lim, dchifbychif_lim, N_bins)

	diff_dMfbyMf = np.mean(np.diff(dMfbyMf_vec))

	diff_dchifbychif = np.mean(np.diff(dchifbychif_vec))

	P_dMfbyMf_dchifbychif = np.zeros(shape=(N_bins,N_bins))

	P_dMfbyMf = np.sum(P_dMfbyMf_dchifbychif, axis=0) * diff_dchifbychif

	P_dchifbychif = np.sum(P_dMfbyMf_dchifbychif, axis=1) * diff_dMfbyMf

	conf_v1v2 = confidence(P_dMfbyMf_dchifbychif)

	gr_height = P_dMfbyMf_dchifbychif[np.argmin(abs(dMfbyMf_vec)), np.argmin(abs(dchifbychif_vec))]

	gr_conf_level = conf_v1v2.level_from_height(gr_height)

list	param_table

	delimiter

	fmt

	P_m1m2_i = P_m1m2_i.T
	plotting More...

	m1_bins_i

	m2_bins_i

	P_chi1chi2_i = P_chi1chi2_i.T

	chi1_bins_i

	chi2_bins_i

	conf_m1m2_i = confidence(P_m1m2_i)

	s1_m1m2_i = conf_m1m2_i.height_from_level(0.68)

	s2_m1m2_i = conf_m1m2_i.height_from_level(0.95)

	conf_chi1chi2_i = confidence(P_chi1chi2_i)

	s1_chi1chi2_i = conf_chi1chi2_i.height_from_level(0.68)

	s2_chi1chi2_i = conf_chi1chi2_i.height_from_level(0.95)

	conf_Mfchif_i = confidence(P_Mfchif_i)

	s1_Mfchif_i = conf_Mfchif_i.height_from_level(0.68)

	s2_Mfchif_i = conf_Mfchif_i.height_from_level(0.95)

	figsize

	cmap

	levels

	linewidths

	x

	ls

	color

	y

	dpi

	ms

	P_m1m2_r = P_m1m2_r.T

	m1_bins_r

	m2_bins_r

	P_chi1chi2_r = P_chi1chi2_r.T

	chi1_bins_r

	chi2_bins_r

	conf_m1m2_r = confidence(P_m1m2_r)

	s1_m1m2_r = conf_m1m2_r.height_from_level(0.68)

	s2_m1m2_r = conf_m1m2_r.height_from_level(0.95)

	conf_chi1chi2_r = confidence(P_chi1chi2_r)

	s1_chi1chi2_r = conf_chi1chi2_r.height_from_level(0.68)

	s2_chi1chi2_r = conf_chi1chi2_r.height_from_level(0.95)

	conf_Mfchif_r = confidence(P_Mfchif_r)

	s1_Mfchif_r = conf_Mfchif_r.height_from_level(0.68)

	s2_Mfchif_r = conf_Mfchif_r.height_from_level(0.95)

	P_m1m2_imr = P_m1m2_imr.T

	m1_bins_imr

	m2_bins_imr

	P_chi1chi2_imr = P_chi1chi2_imr.T

	chi1_bins_imr

	chi2_bins_imr

	conf_m1m2_imr = confidence(P_m1m2_imr)

	s1_m1m2_imr = conf_m1m2_imr.height_from_level(0.68)

	s2_m1m2_imr = conf_m1m2_imr.height_from_level(0.95)

	conf_chi1chi2_imr = confidence(P_chi1chi2_imr)

	s1_chi1chi2_imr = conf_chi1chi2_imr.height_from_level(0.68)

	s2_chi1chi2_imr = conf_chi1chi2_imr.height_from_level(0.95)

	conf_Mfchif_imr = confidence(P_Mfchif_imr)

	s1_Mfchif_imr = conf_Mfchif_imr.height_from_level(0.68)

	s2_Mfchif_imr = conf_Mfchif_imr.height_from_level(0.95)

	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')

	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')

	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	strs_i = [ 'inspiral', 'inspiral' ]

list	strs_r = [ 'ringdown', 'ringdown' ]

list	strs_imr = [ 'IMR', 'IMR' ]

dictionary	fmt_i = {}

dictionary	fmt_r = {}

dictionary	fmt_imr = {}

	inline
	Label every other level using strings. More...

	True

	use_clabeltext

	s1_v1v2 = conf_v1v2.height_from_level(0.68)

	s2_v1v2 = conf_v1v2.height_from_level(0.95)

	conf_v1 = confidence(P_dMfbyMf)

	s1_v1 = conf_v1.height_from_level(0.68)

	s2_v1 = conf_v1.height_from_level(0.95)

	conf_v2 = confidence(P_dchifbychif)

	s1_v2 = conf_v2.height_from_level(0.68)

	s2_v2 = conf_v2.height_from_level(0.95)

	left1_v1 = min(dMfbyMf_vec[np.where(P_dMfbyMf>=s1_v1)[0]])

	right1_v1 = max(dMfbyMf_vec[np.where(P_dMfbyMf>=s1_v1)[0]])

	left2_v1 = min(dMfbyMf_vec[np.where(P_dMfbyMf>=s2_v1)[0]])

	right2_v1 = max(dMfbyMf_vec[np.where(P_dMfbyMf>=s2_v1)[0]])

	left1_v2 = min(dchifbychif_vec[np.where(P_dchifbychif>s1_v2)[0]])

	right1_v2 = max(dchifbychif_vec[np.where(P_dchifbychif>s1_v2)[0]])

	left2_v2 = min(dchifbychif_vec[np.where(P_dchifbychif>s2_v2)[0]])

	right2_v2 = max(dchifbychif_vec[np.where(P_dchifbychif>s2_v2)[0]])

	colspan

	lw

	rowspan

	mew

Function Documentation

◆ set_tick_sizes()

def imrtgr_imr_consistency_test.set_tick_sizes	(	ax,
		major,
		minor
	)

Definition at line 45 of file imrtgr_imr_consistency_test.py.

Variable Documentation

◆ preamble

imrtgr_imr_consistency_test.preamble

Definition at line 27 of file imrtgr_imr_consistency_test.py.

◆ usetex

imrtgr_imr_consistency_test.usetex

Definition at line 29 of file imrtgr_imr_consistency_test.py.

◆ family

imrtgr_imr_consistency_test.family

Definition at line 30 of file imrtgr_imr_consistency_test.py.

◆ serif

imrtgr_imr_consistency_test.serif

Definition at line 31 of file imrtgr_imr_consistency_test.py.

◆ default

imrtgr_imr_consistency_test.default

Definition at line 32 of file imrtgr_imr_consistency_test.py.

◆ markeredgewidth

imrtgr_imr_consistency_test.markeredgewidth

Definition at line 33 of file imrtgr_imr_consistency_test.py.

◆ linewidth

imrtgr_imr_consistency_test.linewidth

Definition at line 34 of file imrtgr_imr_consistency_test.py.

◆ labelsize

imrtgr_imr_consistency_test.labelsize

Definition at line 35 of file imrtgr_imr_consistency_test.py.

◆ fontsize

imrtgr_imr_consistency_test.fontsize

Definition at line 39 of file imrtgr_imr_consistency_test.py.

◆ pad

imrtgr_imr_consistency_test.pad

Definition at line 40 of file imrtgr_imr_consistency_test.py.

◆ size

imrtgr_imr_consistency_test.size

Definition at line 42 of file imrtgr_imr_consistency_test.py.

◆ start_time

imrtgr_imr_consistency_test.start_time = time.time()

Definition at line 84 of file imrtgr_imr_consistency_test.py.

◆ parser

imrtgr_imr_consistency_test.parser = OptionParser()

Definition at line 87 of file imrtgr_imr_consistency_test.py.

◆ dest

imrtgr_imr_consistency_test.dest

Definition at line 88 of file imrtgr_imr_consistency_test.py.

◆ help

imrtgr_imr_consistency_test.help

Definition at line 88 of file imrtgr_imr_consistency_test.py.

◆ type

imrtgr_imr_consistency_test.type

Definition at line 105 of file imrtgr_imr_consistency_test.py.

◆ options

imrtgr_imr_consistency_test.options

Definition at line 111 of file imrtgr_imr_consistency_test.py.

◆ args

imrtgr_imr_consistency_test.args

Definition at line 111 of file imrtgr_imr_consistency_test.py.

◆ MfafKDE

imrtgr_imr_consistency_test.MfafKDE = options.MfafKDE

Definition at line 112 of file imrtgr_imr_consistency_test.py.

◆ insp_post

imrtgr_imr_consistency_test.insp_post = options.insp_post

Definition at line 113 of file imrtgr_imr_consistency_test.py.

◆ ring_post

imrtgr_imr_consistency_test.ring_post = options.ring_post

Definition at line 114 of file imrtgr_imr_consistency_test.py.

◆ imr_post

imrtgr_imr_consistency_test.imr_post = options.imr_post

Definition at line 115 of file imrtgr_imr_consistency_test.py.

◆ prior_Mfchif_file

imrtgr_imr_consistency_test.prior_Mfchif_file = options.prior_Mfchif_file

Definition at line 116 of file imrtgr_imr_consistency_test.py.

◆ out_dir

imrtgr_imr_consistency_test.out_dir = options.out_dir

Definition at line 117 of file imrtgr_imr_consistency_test.py.

◆ fit_formula

imrtgr_imr_consistency_test.fit_formula = options.fit_formula

Definition at line 118 of file imrtgr_imr_consistency_test.py.

◆ debug_plots

imrtgr_imr_consistency_test.debug_plots = options.debug_plots

Definition at line 119 of file imrtgr_imr_consistency_test.py.

◆ pepostproc

imrtgr_imr_consistency_test.pepostproc = options.pepostproc

Definition at line 120 of file imrtgr_imr_consistency_test.py.

◆ insp_fhigh

imrtgr_imr_consistency_test.insp_fhigh = float(options.insp_fhigh)

Definition at line 122 of file imrtgr_imr_consistency_test.py.

◆ ring_flow

imrtgr_imr_consistency_test.ring_flow = float(options.ring_flow)

Definition at line 127 of file imrtgr_imr_consistency_test.py.

◆ waveform

imrtgr_imr_consistency_test.waveform = options.waveform

Definition at line 131 of file imrtgr_imr_consistency_test.py.

◆ N_bins

imrtgr_imr_consistency_test.N_bins = int(options.N_bins)

Definition at line 135 of file imrtgr_imr_consistency_test.py.

◆ dMfbyMf_lim

imrtgr_imr_consistency_test.dMfbyMf_lim = float(options.dMfbyMf_lim)

Definition at line 136 of file imrtgr_imr_consistency_test.py.

◆ dchifbychif_lim

imrtgr_imr_consistency_test.dchifbychif_lim = float(options.dchifbychif_lim)

Definition at line 137 of file imrtgr_imr_consistency_test.py.

◆ lalinference_datadir

imrtgr_imr_consistency_test.lalinference_datadir = os.getenv('LALINFERENCE_DATADIR')

Definition at line 138 of file imrtgr_imr_consistency_test.py.

◆ run_command

imrtgr_imr_consistency_test.run_command = open('%s/command.txt'%(out_dir),'w')

Definition at line 148 of file imrtgr_imr_consistency_test.py.

◆ insp_post_link

imrtgr_imr_consistency_test.insp_post_link = insp_post

Definition at line 158 of file imrtgr_imr_consistency_test.py.

◆ ring_post_link

imrtgr_imr_consistency_test.ring_post_link = ring_post

Definition at line 159 of file imrtgr_imr_consistency_test.py.

◆ imr_post_link

imrtgr_imr_consistency_test.imr_post_link = imr_post

Definition at line 160 of file imrtgr_imr_consistency_test.py.

◆ insp_posplots

imrtgr_imr_consistency_test.insp_posplots = os.path.realpath(os.path.dirname(insp_post_link))

Definition at line 172 of file imrtgr_imr_consistency_test.py.

◆ ring_posplots

imrtgr_imr_consistency_test.ring_posplots = os.path.realpath(os.path.dirname(ring_post_link))

Definition at line 173 of file imrtgr_imr_consistency_test.py.

◆ imr_posplots

imrtgr_imr_consistency_test.imr_posplots = os.path.realpath(os.path.dirname(imr_post_link))

Definition at line 174 of file imrtgr_imr_consistency_test.py.

◆ insp_target

imrtgr_imr_consistency_test.insp_target = os.path.join(out_dir, 'pe_insp')

Definition at line 176 of file imrtgr_imr_consistency_test.py.

◆ ring_target

imrtgr_imr_consistency_test.ring_target = os.path.join(out_dir, 'pe_ring')

Definition at line 177 of file imrtgr_imr_consistency_test.py.

◆ imr_target

imrtgr_imr_consistency_test.imr_target = os.path.join(out_dir, 'pe_imr')

Definition at line 178 of file imrtgr_imr_consistency_test.py.

◆ m1_inj

imrtgr_imr_consistency_test.m1_inj = options.m1_inj

Definition at line 200 of file imrtgr_imr_consistency_test.py.

◆ m2_inj

imrtgr_imr_consistency_test.m2_inj = options.m2_inj

Definition at line 201 of file imrtgr_imr_consistency_test.py.

◆ chi1_inj

imrtgr_imr_consistency_test.chi1_inj = options.chi1_inj

Definition at line 202 of file imrtgr_imr_consistency_test.py.

◆ chi2_inj

imrtgr_imr_consistency_test.chi2_inj = options.chi2_inj

Definition at line 203 of file imrtgr_imr_consistency_test.py.

◆ chi1z_inj

imrtgr_imr_consistency_test.chi1z_inj = options.chi1z_inj

Definition at line 204 of file imrtgr_imr_consistency_test.py.

◆ chi2z_inj

imrtgr_imr_consistency_test.chi2z_inj = options.chi2z_inj

Definition at line 205 of file imrtgr_imr_consistency_test.py.

◆ phi12_inj

imrtgr_imr_consistency_test.phi12_inj = options.phi12_inj

Definition at line 206 of file imrtgr_imr_consistency_test.py.

◆ plot_injection_lines

bool imrtgr_imr_consistency_test.plot_injection_lines = False

Definition at line 209 of file imrtgr_imr_consistency_test.py.

◆ q_inj

imrtgr_imr_consistency_test.q_inj = m1_inj/m2_inj

Definition at line 219 of file imrtgr_imr_consistency_test.py.

◆ eta_inj

int imrtgr_imr_consistency_test.eta_inj = q_inj/(1.+q_inj)**2.

Definition at line 220 of file imrtgr_imr_consistency_test.py.

◆ Mf_inj

imrtgr_imr_consistency_test.Mf_inj

Definition at line 221 of file imrtgr_imr_consistency_test.py.

◆ chif_inj

imrtgr_imr_consistency_test.chif_inj

Definition at line 221 of file imrtgr_imr_consistency_test.py.

◆ m1_id

imrtgr_imr_consistency_test.m1_id

Read the posteriors from the inspiral, ringdown and imr PE runs (after post-processing)

Definition at line 229 of file imrtgr_imr_consistency_test.py.

◆ m2_id

imrtgr_imr_consistency_test.m2_id

Definition at line 229 of file imrtgr_imr_consistency_test.py.

◆ a1_id

imrtgr_imr_consistency_test.a1_id

Definition at line 229 of file imrtgr_imr_consistency_test.py.

◆ a2_id

imrtgr_imr_consistency_test.a2_id

Definition at line 229 of file imrtgr_imr_consistency_test.py.

◆ a1z_id

imrtgr_imr_consistency_test.a1z_id

Definition at line 229 of file imrtgr_imr_consistency_test.py.

◆ a2z_id

imrtgr_imr_consistency_test.a2z_id

Definition at line 229 of file imrtgr_imr_consistency_test.py.

◆ phi12_id

imrtgr_imr_consistency_test.phi12_id

Definition at line 229 of file imrtgr_imr_consistency_test.py.

◆ insp_data

imrtgr_imr_consistency_test.insp_data = np.genfromtxt(insp_post, dtype=None, names=True)

Definition at line 236 of file imrtgr_imr_consistency_test.py.

◆ m1_i

imrtgr_imr_consistency_test.m1_i

Definition at line 238 of file imrtgr_imr_consistency_test.py.

◆ m2_i

imrtgr_imr_consistency_test.m2_i

Definition at line 238 of file imrtgr_imr_consistency_test.py.

◆ chi1_i

imrtgr_imr_consistency_test.chi1_i

Definition at line 238 of file imrtgr_imr_consistency_test.py.

◆ chi2_i

imrtgr_imr_consistency_test.chi2_i

Definition at line 238 of file imrtgr_imr_consistency_test.py.

◆ chi1z_i

imrtgr_imr_consistency_test.chi1z_i

Definition at line 238 of file imrtgr_imr_consistency_test.py.

◆ chi2z_i

imrtgr_imr_consistency_test.chi2z_i

Definition at line 238 of file imrtgr_imr_consistency_test.py.

◆ phi12_i

imrtgr_imr_consistency_test.phi12_i = insp_data[phi12_id]

Definition at line 241 of file imrtgr_imr_consistency_test.py.

◆ Mf_i

imrtgr_imr_consistency_test.Mf_i

Definition at line 245 of file imrtgr_imr_consistency_test.py.

◆ chif_i

imrtgr_imr_consistency_test.chif_i

Definition at line 245 of file imrtgr_imr_consistency_test.py.

◆ ring_data

imrtgr_imr_consistency_test.ring_data = np.genfromtxt(ring_post, dtype=None, names=True)

Definition at line 248 of file imrtgr_imr_consistency_test.py.

◆ m1_r

imrtgr_imr_consistency_test.m1_r

Definition at line 250 of file imrtgr_imr_consistency_test.py.

◆ m2_r

imrtgr_imr_consistency_test.m2_r

Definition at line 250 of file imrtgr_imr_consistency_test.py.

◆ chi1_r

imrtgr_imr_consistency_test.chi1_r

Definition at line 250 of file imrtgr_imr_consistency_test.py.

◆ chi2_r

imrtgr_imr_consistency_test.chi2_r

Definition at line 250 of file imrtgr_imr_consistency_test.py.

◆ chi1z_r

imrtgr_imr_consistency_test.chi1z_r

Definition at line 250 of file imrtgr_imr_consistency_test.py.

◆ chi2z_r

imrtgr_imr_consistency_test.chi2z_r

Definition at line 250 of file imrtgr_imr_consistency_test.py.

◆ phi12_r

imrtgr_imr_consistency_test.phi12_r = ring_data[phi12_id]

Definition at line 253 of file imrtgr_imr_consistency_test.py.

◆ Mf_r

imrtgr_imr_consistency_test.Mf_r

Definition at line 257 of file imrtgr_imr_consistency_test.py.

◆ chif_r

imrtgr_imr_consistency_test.chif_r

Definition at line 257 of file imrtgr_imr_consistency_test.py.

◆ imr_data

imrtgr_imr_consistency_test.imr_data = np.genfromtxt(imr_post, dtype=None, names=True)

Definition at line 260 of file imrtgr_imr_consistency_test.py.

◆ m1_imr

imrtgr_imr_consistency_test.m1_imr

Definition at line 262 of file imrtgr_imr_consistency_test.py.

◆ m2_imr

imrtgr_imr_consistency_test.m2_imr

Definition at line 262 of file imrtgr_imr_consistency_test.py.

◆ chi1_imr

imrtgr_imr_consistency_test.chi1_imr

Definition at line 262 of file imrtgr_imr_consistency_test.py.

◆ chi2_imr

imrtgr_imr_consistency_test.chi2_imr

Definition at line 262 of file imrtgr_imr_consistency_test.py.

◆ chi1z_imr

imrtgr_imr_consistency_test.chi1z_imr

Definition at line 262 of file imrtgr_imr_consistency_test.py.

◆ chi2z_imr

imrtgr_imr_consistency_test.chi2z_imr

Definition at line 262 of file imrtgr_imr_consistency_test.py.

◆ phi12_imr

imrtgr_imr_consistency_test.phi12_imr = imr_data[phi12_id]

Definition at line 265 of file imrtgr_imr_consistency_test.py.

◆ Mf_imr

imrtgr_imr_consistency_test.Mf_imr

Definition at line 269 of file imrtgr_imr_consistency_test.py.

◆ chif_imr

imrtgr_imr_consistency_test.chif_imr

Definition at line 269 of file imrtgr_imr_consistency_test.py.

◆ Mf_lim

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

Definition at line 278 of file imrtgr_imr_consistency_test.py.

◆ chif_lim

imrtgr_imr_consistency_test.chif_lim = max(abs(np.append(np.append(chif_i, chif_r), chif_imr)))

Definition at line 279 of file imrtgr_imr_consistency_test.py.

◆ Mf_bins

imrtgr_imr_consistency_test.Mf_bins = np.linspace(-Mf_lim, Mf_lim, N_bins)

Definition at line 286 of file imrtgr_imr_consistency_test.py.

◆ chif_bins

imrtgr_imr_consistency_test.chif_bins = np.linspace(-chif_lim, chif_lim, N_bins)

Definition at line 287 of file imrtgr_imr_consistency_test.py.

◆ dMf

imrtgr_imr_consistency_test.dMf = np.mean(np.diff(Mf_bins))

Definition at line 289 of file imrtgr_imr_consistency_test.py.

◆ dchif

imrtgr_imr_consistency_test.dchif = np.mean(np.diff(chif_bins))

Definition at line 290 of file imrtgr_imr_consistency_test.py.

◆ Mf_intp

tuple imrtgr_imr_consistency_test.Mf_intp = (Mf_bins[:-1] + Mf_bins[1:])/2.

Definition at line 292 of file imrtgr_imr_consistency_test.py.

◆ chif_intp

tuple imrtgr_imr_consistency_test.chif_intp = (chif_bins[:-1] + chif_bins[1:])/2.

Definition at line 293 of file imrtgr_imr_consistency_test.py.

◆ M_i

imrtgr_imr_consistency_test.M_i

Definition at line 299 of file imrtgr_imr_consistency_test.py.

◆ C_i

imrtgr_imr_consistency_test.C_i

Definition at line 299 of file imrtgr_imr_consistency_test.py.

◆ joint_data

imrtgr_imr_consistency_test.joint_data = gaussian_kde(joint_data)

Definition at line 301 of file imrtgr_imr_consistency_test.py.

◆ f_i

imrtgr_imr_consistency_test.f_i = lambda x,y:kernel.evaluate([x,y])

Definition at line 302 of file imrtgr_imr_consistency_test.py.

◆ P_Mfchif_i

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])

Definition at line 304 of file imrtgr_imr_consistency_test.py.

◆ f_r

imrtgr_imr_consistency_test.f_r = lambda x,y:kernel.evaluate([x,y])

Definition at line 307 of file imrtgr_imr_consistency_test.py.

◆ P_Mfchif_r

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])

Definition at line 309 of file imrtgr_imr_consistency_test.py.

◆ f_imr

imrtgr_imr_consistency_test.f_imr = lambda x,y:kernel.evaluate([x,y])

Definition at line 312 of file imrtgr_imr_consistency_test.py.

◆ P_Mfchif_imr

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])

Definition at line 314 of file imrtgr_imr_consistency_test.py.

◆ bins

imrtgr_imr_consistency_test.bins

Definition at line 322 of file imrtgr_imr_consistency_test.py.

◆ density

imrtgr_imr_consistency_test.density

Definition at line 322 of file imrtgr_imr_consistency_test.py.

◆ f

imrtgr_imr_consistency_test.f = gzip.open(prior_Mfchif_file,'rb')

Undo the effect of the prior from the lalinference posterior.

Lalinference assumes a # uniform prior in component masses. We need to assume a uniform prior in Mf, chif #

Definition at line 348 of file imrtgr_imr_consistency_test.py.

◆ P_Mfchif_pr_interp_obj

imrtgr_imr_consistency_test.P_Mfchif_pr_interp_obj = pickle.load(f)

Definition at line 349 of file imrtgr_imr_consistency_test.py.

◆ P_Mfchif_pr

imrtgr_imr_consistency_test.P_Mfchif_pr = P_Mfchif_pr_interp_obj(Mf_intp, chif_intp)

Definition at line 350 of file imrtgr_imr_consistency_test.py.

◆ P_Mfchif_i_interp_object

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)

Definition at line 374 of file imrtgr_imr_consistency_test.py.

◆ P_Mfchif_r_interp_object

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)

Definition at line 375 of file imrtgr_imr_consistency_test.py.

◆ dMfbyMf_vec

imrtgr_imr_consistency_test.dMfbyMf_vec = np.linspace(-dMfbyMf_lim, dMfbyMf_lim, N_bins)

Definition at line 378 of file imrtgr_imr_consistency_test.py.

◆ dchifbychif_vec

imrtgr_imr_consistency_test.dchifbychif_vec = np.linspace(-dchifbychif_lim, dchifbychif_lim, N_bins)

Definition at line 379 of file imrtgr_imr_consistency_test.py.

◆ diff_dMfbyMf

imrtgr_imr_consistency_test.diff_dMfbyMf = np.mean(np.diff(dMfbyMf_vec))

Definition at line 382 of file imrtgr_imr_consistency_test.py.

◆ diff_dchifbychif

imrtgr_imr_consistency_test.diff_dchifbychif = np.mean(np.diff(dchifbychif_vec))

Definition at line 383 of file imrtgr_imr_consistency_test.py.

◆ P_dMfbyMf_dchifbychif

imrtgr_imr_consistency_test.P_dMfbyMf_dchifbychif = np.zeros(shape=(N_bins,N_bins))

Definition at line 384 of file imrtgr_imr_consistency_test.py.

◆ P_dMfbyMf

imrtgr_imr_consistency_test.P_dMfbyMf = np.sum(P_dMfbyMf_dchifbychif, axis=0) * diff_dchifbychif

Definition at line 396 of file imrtgr_imr_consistency_test.py.

◆ P_dchifbychif

imrtgr_imr_consistency_test.P_dchifbychif = np.sum(P_dMfbyMf_dchifbychif, axis=1) * diff_dMfbyMf

Definition at line 397 of file imrtgr_imr_consistency_test.py.

◆ conf_v1v2

imrtgr_imr_consistency_test.conf_v1v2 = confidence(P_dMfbyMf_dchifbychif)

Definition at line 401 of file imrtgr_imr_consistency_test.py.

◆ gr_height

imrtgr_imr_consistency_test.gr_height = P_dMfbyMf_dchifbychif[np.argmin(abs(dMfbyMf_vec)), np.argmin(abs(dchifbychif_vec))]

Definition at line 402 of file imrtgr_imr_consistency_test.py.

◆ gr_conf_level

imrtgr_imr_consistency_test.gr_conf_level = conf_v1v2.level_from_height(gr_height)

Definition at line 403 of file imrtgr_imr_consistency_test.py.

◆ param_table

list imrtgr_imr_consistency_test.param_table

Initial value:

 =  [['Upper cutoff freq for the inspiral analysis: %s Hz'%insp_fhigh],
                 ['Lower cutoff freq for the ringdown analysis: %s Hz'%ring_flow],
                 ['Waveform approximant: %s'%(waveform)],
                 ['Final mass/spin fitting formula: %s'%(fit_formula)],
                 ['No deviation from GR above %.1f%% confidence level'%(100.*gr_conf_level)]]

Definition at line 407 of file imrtgr_imr_consistency_test.py.

◆ delimiter

imrtgr_imr_consistency_test.delimiter

Definition at line 412 of file imrtgr_imr_consistency_test.py.

◆ fmt

imrtgr_imr_consistency_test.fmt

Definition at line 412 of file imrtgr_imr_consistency_test.py.

◆ P_m1m2_i

imrtgr_imr_consistency_test.P_m1m2_i = P_m1m2_i.T

plotting

inspiral

Definition at line 432 of file imrtgr_imr_consistency_test.py.

◆ m1_bins_i

imrtgr_imr_consistency_test.m1_bins_i

Definition at line 432 of file imrtgr_imr_consistency_test.py.

◆ m2_bins_i

imrtgr_imr_consistency_test.m2_bins_i

Definition at line 432 of file imrtgr_imr_consistency_test.py.

◆ P_chi1chi2_i

imrtgr_imr_consistency_test.P_chi1chi2_i = P_chi1chi2_i.T

Definition at line 433 of file imrtgr_imr_consistency_test.py.

◆ chi1_bins_i

imrtgr_imr_consistency_test.chi1_bins_i

Definition at line 433 of file imrtgr_imr_consistency_test.py.

◆ chi2_bins_i

imrtgr_imr_consistency_test.chi2_bins_i

Definition at line 433 of file imrtgr_imr_consistency_test.py.

◆ conf_m1m2_i

imrtgr_imr_consistency_test.conf_m1m2_i = confidence(P_m1m2_i)

Definition at line 438 of file imrtgr_imr_consistency_test.py.

◆ s1_m1m2_i

imrtgr_imr_consistency_test.s1_m1m2_i = conf_m1m2_i.height_from_level(0.68)

Definition at line 439 of file imrtgr_imr_consistency_test.py.

◆ s2_m1m2_i

imrtgr_imr_consistency_test.s2_m1m2_i = conf_m1m2_i.height_from_level(0.95)

Definition at line 440 of file imrtgr_imr_consistency_test.py.

◆ conf_chi1chi2_i

imrtgr_imr_consistency_test.conf_chi1chi2_i = confidence(P_chi1chi2_i)

Definition at line 442 of file imrtgr_imr_consistency_test.py.

◆ s1_chi1chi2_i

imrtgr_imr_consistency_test.s1_chi1chi2_i = conf_chi1chi2_i.height_from_level(0.68)

Definition at line 443 of file imrtgr_imr_consistency_test.py.

◆ s2_chi1chi2_i

imrtgr_imr_consistency_test.s2_chi1chi2_i = conf_chi1chi2_i.height_from_level(0.95)

Definition at line 444 of file imrtgr_imr_consistency_test.py.

◆ conf_Mfchif_i

imrtgr_imr_consistency_test.conf_Mfchif_i = confidence(P_Mfchif_i)

Definition at line 446 of file imrtgr_imr_consistency_test.py.

◆ s1_Mfchif_i

imrtgr_imr_consistency_test.s1_Mfchif_i = conf_Mfchif_i.height_from_level(0.68)

Definition at line 447 of file imrtgr_imr_consistency_test.py.

◆ s2_Mfchif_i

imrtgr_imr_consistency_test.s2_Mfchif_i = conf_Mfchif_i.height_from_level(0.95)

Definition at line 448 of file imrtgr_imr_consistency_test.py.

◆ figsize

imrtgr_imr_consistency_test.figsize

Definition at line 450 of file imrtgr_imr_consistency_test.py.

◆ cmap

imrtgr_imr_consistency_test.cmap

Definition at line 451 of file imrtgr_imr_consistency_test.py.

◆ levels

imrtgr_imr_consistency_test.levels

Definition at line 452 of file imrtgr_imr_consistency_test.py.

◆ linewidths

imrtgr_imr_consistency_test.linewidths

Definition at line 452 of file imrtgr_imr_consistency_test.py.

◆ x

imrtgr_imr_consistency_test.x

Definition at line 454 of file imrtgr_imr_consistency_test.py.

◆ ls

imrtgr_imr_consistency_test.ls

Definition at line 454 of file imrtgr_imr_consistency_test.py.

◆ color

imrtgr_imr_consistency_test.color

Definition at line 454 of file imrtgr_imr_consistency_test.py.

◆ y

imrtgr_imr_consistency_test.y

Definition at line 455 of file imrtgr_imr_consistency_test.py.

◆ dpi

imrtgr_imr_consistency_test.dpi

Definition at line 461 of file imrtgr_imr_consistency_test.py.

◆ ms

imrtgr_imr_consistency_test.ms

Definition at line 465 of file imrtgr_imr_consistency_test.py.

◆ P_m1m2_r

imrtgr_imr_consistency_test.P_m1m2_r = P_m1m2_r.T

Definition at line 520 of file imrtgr_imr_consistency_test.py.

◆ m1_bins_r

imrtgr_imr_consistency_test.m1_bins_r

Definition at line 520 of file imrtgr_imr_consistency_test.py.

◆ m2_bins_r

imrtgr_imr_consistency_test.m2_bins_r

Definition at line 520 of file imrtgr_imr_consistency_test.py.

◆ P_chi1chi2_r

imrtgr_imr_consistency_test.P_chi1chi2_r = P_chi1chi2_r.T

Definition at line 521 of file imrtgr_imr_consistency_test.py.

◆ chi1_bins_r

imrtgr_imr_consistency_test.chi1_bins_r

Definition at line 521 of file imrtgr_imr_consistency_test.py.

◆ chi2_bins_r

imrtgr_imr_consistency_test.chi2_bins_r

Definition at line 521 of file imrtgr_imr_consistency_test.py.

◆ conf_m1m2_r

imrtgr_imr_consistency_test.conf_m1m2_r = confidence(P_m1m2_r)

Definition at line 526 of file imrtgr_imr_consistency_test.py.

◆ s1_m1m2_r

imrtgr_imr_consistency_test.s1_m1m2_r = conf_m1m2_r.height_from_level(0.68)

Definition at line 527 of file imrtgr_imr_consistency_test.py.

◆ s2_m1m2_r

imrtgr_imr_consistency_test.s2_m1m2_r = conf_m1m2_r.height_from_level(0.95)

Definition at line 528 of file imrtgr_imr_consistency_test.py.

◆ conf_chi1chi2_r

imrtgr_imr_consistency_test.conf_chi1chi2_r = confidence(P_chi1chi2_r)

Definition at line 530 of file imrtgr_imr_consistency_test.py.

◆ s1_chi1chi2_r

imrtgr_imr_consistency_test.s1_chi1chi2_r = conf_chi1chi2_r.height_from_level(0.68)

Definition at line 531 of file imrtgr_imr_consistency_test.py.

◆ s2_chi1chi2_r

imrtgr_imr_consistency_test.s2_chi1chi2_r = conf_chi1chi2_r.height_from_level(0.95)

Definition at line 532 of file imrtgr_imr_consistency_test.py.

◆ conf_Mfchif_r

imrtgr_imr_consistency_test.conf_Mfchif_r = confidence(P_Mfchif_r)

Definition at line 534 of file imrtgr_imr_consistency_test.py.

◆ s1_Mfchif_r

imrtgr_imr_consistency_test.s1_Mfchif_r = conf_Mfchif_r.height_from_level(0.68)

Definition at line 535 of file imrtgr_imr_consistency_test.py.

◆ s2_Mfchif_r

imrtgr_imr_consistency_test.s2_Mfchif_r = conf_Mfchif_r.height_from_level(0.95)

Definition at line 536 of file imrtgr_imr_consistency_test.py.

◆ P_m1m2_imr

imrtgr_imr_consistency_test.P_m1m2_imr = P_m1m2_imr.T

Definition at line 607 of file imrtgr_imr_consistency_test.py.

◆ m1_bins_imr

imrtgr_imr_consistency_test.m1_bins_imr

Definition at line 607 of file imrtgr_imr_consistency_test.py.

◆ m2_bins_imr

imrtgr_imr_consistency_test.m2_bins_imr

Definition at line 607 of file imrtgr_imr_consistency_test.py.

◆ P_chi1chi2_imr

imrtgr_imr_consistency_test.P_chi1chi2_imr = P_chi1chi2_imr.T

Definition at line 608 of file imrtgr_imr_consistency_test.py.

◆ chi1_bins_imr

imrtgr_imr_consistency_test.chi1_bins_imr

Definition at line 608 of file imrtgr_imr_consistency_test.py.

◆ chi2_bins_imr

imrtgr_imr_consistency_test.chi2_bins_imr

Definition at line 608 of file imrtgr_imr_consistency_test.py.

◆ conf_m1m2_imr

imrtgr_imr_consistency_test.conf_m1m2_imr = confidence(P_m1m2_imr)

Definition at line 613 of file imrtgr_imr_consistency_test.py.

◆ s1_m1m2_imr

imrtgr_imr_consistency_test.s1_m1m2_imr = conf_m1m2_imr.height_from_level(0.68)

Definition at line 614 of file imrtgr_imr_consistency_test.py.

◆ s2_m1m2_imr

imrtgr_imr_consistency_test.s2_m1m2_imr = conf_m1m2_imr.height_from_level(0.95)

Definition at line 615 of file imrtgr_imr_consistency_test.py.

◆ conf_chi1chi2_imr

imrtgr_imr_consistency_test.conf_chi1chi2_imr = confidence(P_chi1chi2_imr)

Definition at line 617 of file imrtgr_imr_consistency_test.py.

◆ s1_chi1chi2_imr

imrtgr_imr_consistency_test.s1_chi1chi2_imr = conf_chi1chi2_imr.height_from_level(0.68)

Definition at line 618 of file imrtgr_imr_consistency_test.py.

◆ s2_chi1chi2_imr

imrtgr_imr_consistency_test.s2_chi1chi2_imr = conf_chi1chi2_imr.height_from_level(0.95)

Definition at line 619 of file imrtgr_imr_consistency_test.py.

◆ conf_Mfchif_imr

imrtgr_imr_consistency_test.conf_Mfchif_imr = confidence(P_Mfchif_imr)

Definition at line 621 of file imrtgr_imr_consistency_test.py.

◆ s1_Mfchif_imr

imrtgr_imr_consistency_test.s1_Mfchif_imr = conf_Mfchif_imr.height_from_level(0.68)

Definition at line 622 of file imrtgr_imr_consistency_test.py.

◆ s2_Mfchif_imr

imrtgr_imr_consistency_test.s2_Mfchif_imr = conf_Mfchif_imr.height_from_level(0.95)

Definition at line 623 of file imrtgr_imr_consistency_test.py.

◆ CSi

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')

Definition at line 695 of file imrtgr_imr_consistency_test.py.

◆ CSr

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')

Definition at line 696 of file imrtgr_imr_consistency_test.py.