From cf8758de1255d12b94c0c88d9ece502c5f9dc8c1 Mon Sep 17 00:00:00 2001
From: Chad Hanna <chad.hanna@comp-hd-002.gwave.ics.psu.edu>
Date: Mon, 10 Jun 2019 15:46:15 -0400
Subject: [PATCH] gstlal_inspiral_create_dt_dphi_snr_ratio_pdfs: require
covarianc matrix .h5 as input
---
.../bin/gstlal_inspiral_create_dt_dphi_snr_ratio_pdfs | 5 +++--
1 file changed, 3 insertions(+), 2 deletions(-)
diff --git a/gstlal-inspiral/bin/gstlal_inspiral_create_dt_dphi_snr_ratio_pdfs b/gstlal-inspiral/bin/gstlal_inspiral_create_dt_dphi_snr_ratio_pdfs
index 89952cf653..9fffdfb349 100755
--- a/gstlal-inspiral/bin/gstlal_inspiral_create_dt_dphi_snr_ratio_pdfs
+++ b/gstlal-inspiral/bin/gstlal_inspiral_create_dt_dphi_snr_ratio_pdfs
@@ -21,6 +21,7 @@ from gstlal.stats import inspiral_extrinsics
import argparse
parser = argparse.ArgumentParser(description = "create dt dphi snr pdfs")
+parser.add_argument("--cov-mat-h5", help = "hdf5 file that contains the covariance matrix")
parser.add_argument("--start", metavar='N', type=int, help='The first index to marginalize')
parser.add_argument("--stop", metavar='N', type=int, help='The last index to marginalize')
parser.add_argument("--output", metavar='name', type=str, help='The output file name', default = "inspiral_dtdphi_pdf.h5")
@@ -29,8 +30,8 @@ options = parser.parse_args()
if options.start is not None and options.stop is not None:
if options.start % inspiral_extrinsics.TimePhaseSNR.numchunks or options.stop % inspiral_extrinsics.TimePhaseSNR.numchunks:
raise ValueError("start and stop must be multiple of %d" % TimePhaseSNR.numchunks)
- IE = inspiral_extrinsics.TimePhaseSNR(verbose = True, margstart = options.start, margstop = options.stop)
+ IE = inspiral_extrinsics.TimePhaseSNR.from_hdf5(fname = options.cov_mat_h5, verbose = True, margstart = options.start, margstop = options.stop)
else:
- IE = inspiral_extrinsics.TimePhaseSNR(verbose = True)
+ IE = inspiral_extrinsics.TimePhaseSNR.from_hdf5(fname = options.cov_mat_h5, verbose = True)
IE.to_hdf5(options.output)
--
GitLab