From d94bd5da14b9b27bad8df4c0eb4cdbd15205da76 Mon Sep 17 00:00:00 2001
From: Chad Hanna <chad.hanna@ligo.org>
Date: Sat, 2 Jun 2018 14:47:37 -0400
Subject: [PATCH] gstlal_inspiral_create_dt_dphi_snr_ratio_pdfs: optionally
only create a subset of the marginalized dataset
---
...gstlal_inspiral_create_dt_dphi_snr_ratio_pdfs | 16 ++++++++++++----
1 file changed, 12 insertions(+), 4 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 d67c9b6b21..4e4477a10c 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
@@ -19,9 +19,17 @@
import sys
from gstlal.stats import inspiral_extrinsics
-if len(sys.argv) > 1:
- print "USAGE: %s" % sys.argv[0]
+if len(sys.argv) != 1 and len(sys.argv) != 3:
+ print "USAGE: %s <marg start> <marg stop>" % sys.argv[0]
sys.exit()
-IE = inspiral_extrinsics.TimePhaseSNR(verbose = True)
-IE.to_hdf5("inspiral_dtdphi_pdf.h5")
+if len(sys.argv) == 3:
+ start = int(sys.argv[1])
+ stop = int(sys.argv[2])
+ if start % 50 or stop % 50:
+ raise ValueError("start and stop must be multiple of 50")
+ IE = inspiral_extrinsics.TimePhaseSNR(verbose = True, margstart = start, margstop = stop)
+ IE.to_hdf5("inspiral_dtdphi_%d_%d_pdf.h5" % (start, stop))
+else:
+ IE = inspiral_extrinsics.TimePhaseSNR(verbose = True)
+ IE.to_hdf5("inspiral_dtdphi_pdf.h5")
--
GitLab