gstlal_inspiral_create_p_of_ifos_given_horizon: program to figure out what...

gstlal_inspiral_create_p_of_ifos_given_horizon: program to figure out what probability you have of forming coincs in an instrument network

gstlal-inspiral/bin/gstlal_inspiral_create_p_of_ifos_given_horizon

+#!/usr/bin/env python
+#
+# Copyright (C) 2018 Chad Hanna
+#
+# This program is free software; you can redistribute it and/or modify it
+# under the terms of the GNU General Public License as published by the
+# Free Software Foundation; either version 2 of the License, or (at your
+# option) any later version.
+#
+# This program is distributed in the hope that it will be useful, but
+# WITHOUT ANY WARRANTY; without even the implied warranty of
+# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General
+# Public License for more details.
+#
+# You should have received a copy of the GNU General Public License along
+# with this program; if not, write to the Free Software Foundation, Inc.,
+# 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA.
+
+import sys
+from gstlal.stats import inspiral_extrinsics
+
+parser = argparse.ArgumentParser(description='Generate PDFs of extrinsic parameters')
+parser.add_argument('--snr-thresh', type=float, default = 4.0, help = 'set the snr minimum to define found')
+parser.add_argument('--output-file', default = 'extparams.h5', help = 'set the output hdf5 file. Default extparams.h5')
+parser.add_argument('--instrument', action = "append", help='add instrument. Can be given multiple times')
+args = parser.parse_args()
+
+PofI = inspiral_extrinsics.p_of_instruments_given_horizons(instruments = args.instrument, snr_thresh = args.snr_thresh)
+PofI.to_hdf5(args.output_file)