Skip to content
Snippets Groups Projects
Commit 190b7e24 authored by Chad Hanna's avatar Chad Hanna
Browse files

gstlal_segments_operations: fix default output name and add useful verbosity

parent f9465eb7
No related branches found
No related tags found
No related merge requests found
......@@ -40,7 +40,7 @@ def parse_command_line():
parser.add_option("--intersection", action = "store_true", help = "intersection of --segment-name1 in segment-file1 and --segment-name2 in segment-file2")
parser.add_option("--union", action = "store_true", help = "union of --segment-name1 in segment-file1 and --segment-name2 in segment-file2")
parser.add_option("--diff", action = "store_true", help = "diff of --segment-name1 in segment-file1 and --segment-name2 in segment-file2")
parser.add_option("--output-file", metavar = "file", default = ".", help = "Set the name of the output file")
parser.add_option("--output-file", metavar = "file", help = "Set the name of the output file")
parser.add_option("--output-segment-name", metavar = "name", default = "datasegments", help = "Set the name of output segments default datasegments")
parser.add_option("-v", "--verbose", action = "store_true", help = "Be verbose.")
options, filenames = parser.parse_args()
......@@ -59,6 +59,12 @@ def extract_segs(fname, segname):
segs1 = extract_segs(options.segment_file1, options.segment_name1)
segs2 = extract_segs(options.segment_file2, options.segment_name2)
if options.verbose:
print "segs 1 live time: ", abs(segs1)
print "segs 2 live time: ", abs(segs2)
print "two or more detectors on segs1: ", abs(segmentsUtils.vote(segs1.values(), 2))
print "two or more detectors on segs2: ", abs(segmentsUtils.vote(segs2.values(), 2))
if options.diff:
segs = segs1 - segs2
if options.intersection:
......@@ -66,17 +72,21 @@ if options.intersection:
if options.union:
segs = segs1 | segs2
xmldoc = ligolw.Document()
xmldoc.appendChild(ligolw.LIGO_LW())
if options.verbose:
print "segs live time: ", abs(segs)
print "two or more detectors on segs: ", abs(segmentsUtils.vote(segs.values(), 2))
process = ligolw_process.register_to_xmldoc(xmldoc, sys.argv[0], options.__dict__)
ligolw_process.set_process_end_time(process)
if options.output_file:
xmldoc = ligolw.Document()
xmldoc.appendChild(ligolw.LIGO_LW())
lwseglists = ligolw_segments.LigolwSegments(xmldoc)
lwseglists.insert_from_segmentlistdict(segs, options.output_segment_name)
lwseglists.optimize()
lwseglists.finalize(process)
process = ligolw_process.register_to_xmldoc(xmldoc, sys.argv[0], options.__dict__)
ligolw_process.set_process_end_time(process)
lwseglists = ligolw_segments.LigolwSegments(xmldoc)
lwseglists.insert_from_segmentlistdict(segs, options.output_segment_name)
lwseglists.optimize()
lwseglists.finalize(process)
ligolw_utils.write_filename(xmldoc, options.output_file, verbose = options.verbose, gz = options.output_file.endswith('.gz'))
ligolw_utils.write_filename(xmldoc, options.output_file, verbose = options.verbose, gz = options.output_file.endswith('.gz'))
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment