Commit 58d26226 authored by Leo P. Singer's avatar Leo P. Singer
Browse files

Add option to send output to file

Original: c94874fa5082dfa1318fa77db53824c03221978a
parent 740732c5
......@@ -17,6 +17,7 @@
# 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
from __future__ import division
from __future__ import print_function
Match sky maps with injections in an inspinjfind-style sqlite database and print
summary values for each sky map:
......@@ -38,7 +39,12 @@ from lalinference.bayestar import command
parser = OptionParser(
usage="%prog DATABASE.sqlite FILE1.fits[.gz] FILE2.fits[.gz] ...")
usage="%prog [-o OUTPUT] DATABASE.sqlite FILE1.fits[.gz] FILE2.fits[.gz] ...",
Option("-o", "--output", default="/dev/stdout",
help="Name of output file [default=%default]")
opts, args = parser.parse_args()
......@@ -47,6 +53,8 @@ try:
except IndexError:
parser.error("not enough command line arguments")
outfile = open(opts.output, "w")
# Imports.
import os
import numpy as np
......@@ -131,7 +139,8 @@ progress = ProgressBar()
progress.update(-1, 'opening database')
db = sqlite3.connect(dbfilename)
print 'objid,far,searched_area,searched_prob,offset,runtime'
print('objid', 'far', 'searched_area', 'searched_prob', 'offset', 'runtime',
sep=',', file=outfile)
for fitsfilename in progress.iterate(fitsfilenames):
sky_map, metadata = fits.read_sky_map(fitsfilename)
......@@ -145,4 +154,5 @@ for fitsfilename in progress.iterate(fitsfilenames):
true_ra, true_dec, far = db.execute(sql, (coinc_event_id,)).fetchone()
searched_area, searched_prob, offset = find_injection(sky_map, true_ra, true_dec)
print ','.join(str(item) for item in (coinc_event_id, far, searched_area, searched_prob, offset, runtime))
print(coinc_event_id, far, searched_area, searched_prob, offset, runtime,
sep=',', file=outfile)
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment