#
# Run coincidence algorithm.
#
burca.burca(xmldoc=xmldoc,
process_id=process.process_id,
EventListType=EventListType,
CoincTables=CoincTables,
coinc_definer_row=CoincDef,
thresholds=options.threshold,
ntuple_comparefunc=ntuple_comparefunc,
verbose=options.verbose)
#
# Close out the process table.
#
ligolw_process.set_process_end_time(process)
#
# Write back to disk, and clean up.
#
ligolw_utils.write_filename(xmldoc,
filename,
verbose=options.verbose,
gz=(filename or "stdout").endswith(".gz"))
xmldoc.unlink()
lsctables.reset_next_ids(lsctables.TableByName.values())
h5file.create_dataset("/sim_inspiral",
data=ligolw_table_to_array(sims),
compression='gzip',
compression_opts=1)
h5file.create_dataset("/sngl_inspiral",
data=ligolw_table_to_array(tmplts),
compression='gzip',
compression_opts=1)
h5file.create_dataset("/match_map",
data=match_map,
compression='gzip',
compression_opts=1)
h5file.flush()
# merge process and process_params tables, then complete ourselves
lsctables.reset_next_ids(
(lsctables.ProcessTable, lsctables.ProcessParamsTable))
ligolw_add.reassign_ids(fake_xmldoc)
ligolw_add.merge_ligolws(fake_xmldoc)
ligolw_add.merge_compatible_tables(fake_xmldoc)
ligolw_process.set_process_end_time(process)
# output process
proc = lsctables.ProcessTable.get_table(fake_xmldoc)
for p in proc:
p.cvs_entry_time = 0
p.end_time = 0
h5file.create_dataset("/process", data=ligolw_table_to_array(proc))
pp = lsctables.ProcessParamsTable.get_table(fake_xmldoc)
h5file.create_dataset("/process_params", data=ligolw_table_to_array(pp))
h5file.close()
print >>sys.stderr, "skipping (use --force to force)"
continue
if options.verbose:
print >>sys.stderr, "continuing by --force"
#
# add process metadata to document
#
process = binjfind.append_process(xmldoc, match_algorithm = options.match_algorithm, comment = options.comment)
#
# run binjfind algorithm
#
binjfind.binjfind(xmldoc, process, search, snglcomparefunc, nearcoinccomparefunc, verbose = options.verbose)
#
# close out the process metadata
#
ligolw_process.set_process_end_time(process)
#
# done
#
utils.write_filename(xmldoc, filename, verbose = options.verbose, gz = (filename or "stdout").endswith(".gz"))
xmldoc.unlink()
lsctables.reset_next_ids(lsctables.TableByName.values())
match_map[j] = (match_tup[0], inj_wf.sigmasq)
tmplts.append(bank._templates[match_tup[1]].to_sngl())
if verbose:
print "total number of match calculations:", bank._nmatch
# write out results
h5file = H5File("%s.h5" % usertag, "w")
h5file.create_dataset("/sim_inspiral", data=ligolw_table_to_array(sims), compression='gzip', compression_opts=1)
h5file.create_dataset("/sngl_inspiral", data=ligolw_table_to_array(tmplts), compression='gzip', compression_opts=1)
h5file.create_dataset("/match_map", data=match_map, compression='gzip', compression_opts=1)
h5file.flush()
# merge process and process_params tables, then complete ourselves
lsctables.reset_next_ids((lsctables.ProcessTable, lsctables.ProcessParamsTable))
ligolw_add.reassign_ids(fake_xmldoc)
ligolw_add.merge_ligolws(fake_xmldoc)
ligolw_add.merge_compatible_tables(fake_xmldoc)
ligolw_process.set_process_end_time(process)
# output process
proc = lsctables.ProcessTable.get_table(fake_xmldoc)
for p in proc:
p.cvs_entry_time = 0
p.end_time = 0
h5file.create_dataset("/process", data=ligolw_table_to_array(proc))
pp = lsctables.ProcessParamsTable.get_table(fake_xmldoc)
h5file.create_dataset("/process_params", data=ligolw_table_to_array(pp))
h5file.close()