From 59bd4eada4ecfc2a35fc035229093aeb4f95fb80 Mon Sep 17 00:00:00 2001
From: Patrick Godwin <patrick.godwin@ligo.org>
Date: Tue, 25 Apr 2017 07:53:59 -0700
Subject: [PATCH] gstlal_idq_trigger_gen: added latency column to trigger
output
---
gstlal-ugly/bin/gstlal_idq_trigger_gen | 7 +++++--
1 file changed, 5 insertions(+), 2 deletions(-)
diff --git a/gstlal-ugly/bin/gstlal_idq_trigger_gen b/gstlal-ugly/bin/gstlal_idq_trigger_gen
index 22666bad31..714091d28c 100755
--- a/gstlal-ugly/bin/gstlal_idq_trigger_gen
+++ b/gstlal-ugly/bin/gstlal_idq_trigger_gen
@@ -47,6 +47,7 @@ from gstlal import sngltriggertable
from gstlal import pipeparts
from gstlal import simplehandler
from glue.ligolw import utils as ligolw_utils
+from lal import gpstime
#
# Make sure we have sufficient resources
@@ -155,7 +156,7 @@ class MultiChannelHandler(simplehandler.Handler):
self.cadence = options.cadence
# create header for trigger file
#self.header = "# start_time stop_time time frequency unnormalized_energy normalized_energy chisqdof significance channel\n"
- self.header = "# %18s\t%20s\t%20s\t%6s\t%8s\t%8s\t%8s\t%8s\t%s\n" % ("start_time", "stop_time", "trigger_time", "phase", "snr", "chisq", "sigmasq", "rate", "channel")
+ self.header = "# %18s\t%20s\t%20s\t%6s\t%8s\t%8s\t%8s\t%9s\t%8s\t%s\n" % ("start_time", "stop_time", "trigger_time", "phase", "snr", "chisq", "sigmasq", "latency", "rate", "channel")
super(MultiChannelHandler, self).__init__(*args, **kwargs)
def do_on_message(self, bus, message):
@@ -186,9 +187,11 @@ class MultiChannelHandler(simplehandler.Handler):
if mapinfo.data:
for row in sngltriggertable.GSTLALSnglTrigger.from_buffer(mapinfo.data):
trigger_time = row.end_time + row.end_time_ns * 1e-9
+ current_time = gpstime.gps_time_now().gpsSeconds + gpstime.gps_time_now().gpsNanoSeconds * 1e-9
+ latency = numpy.ceil(current_time - buftime)
start_time = trigger_time - duration/2
stop_time = trigger_time + duration/2
- fdata += "%20.9f\t%20.9f\t%20.9f\t%6.3f\t%8.3f\t%8.3f\t%8.3f\t%8.1f\t%s\n" % (start_time, stop_time, trigger_time, row.phase, row.snr, row.chisq, row.sigmasq, float(rate), channel)
+ fdata += "%20.9f\t%20.9f\t%20.9f\t%6.3f\t%8.3f\t%8.3f\t%8.3f\t%9d\t%8.1f\t%s\n" % (start_time, stop_time, trigger_time, row.phase, row.snr, row.chisq, row.sigmasq, latency, float(rate), channel)
memory.unmap(mapinfo)
# Save a "latest"
--
GitLab