From b19a62c68a748b03390067b0b640621175981785 Mon Sep 17 00:00:00 2001
From: Daichi Tsuna <daichi.tsuna@ligo.org>
Date: Tue, 11 Dec 2018 21:12:29 -0800
Subject: [PATCH] cs_triggergen: fix bug at firbank stage
Fixed bug in assertion error due to the fir_matrix changing its size. Made it to not change.
A bug in framecppchanneldemux (assertion error due to numsrcpads being 0) still persists, which appears when running the DAG but not when running in the head nodes.
---
gstlal-burst/bin/gstlal_cs_triggergen | 14 ++++++++++----
1 file changed, 10 insertions(+), 4 deletions(-)
diff --git a/gstlal-burst/bin/gstlal_cs_triggergen b/gstlal-burst/bin/gstlal_cs_triggergen
index 928ba095ae..f22233b8bf 100755
--- a/gstlal-burst/bin/gstlal_cs_triggergen
+++ b/gstlal-burst/bin/gstlal_cs_triggergen
@@ -135,9 +135,15 @@ class PSDHandler(simplehandler.Handler):
firbank.set_property("fir_matrix", template_t)
self.firbank = firbank
else:
- template, _ = lalsimulation.GenerateStringCusp(1.0,30,1.0/options.sample_rate)
- firbank.set_property("fir_matrix",numpy.zeros((template.data.length,len(template_bank_table)), dtype=float))
- firbank.set_property("latency",-(template.data.length - 1) // 2)
+ # use templates with all zeros during burn-in period, that way we won't get any triggers.
+ template = [None] * len(template_bank_table)
+ for i, row in enumerate(template_bank_table):
+ template[i], _ = lalsimulation.GenerateStringCusp(1.0,30,1.0/options.sample_rate)
+ template[i] = lal.ResizeREAL8TimeSeries(template[i], -int(32*options.sample_rate - template[i].data.length)//2 ,int(32*options.sample_rate))
+ template[i] = template[i].data.data
+ template[i] *= 0.0
+ firbank.set_property("latency",-(len(template[0]) - 1) // 2)
+ firbank.set_property("fir_matrix", template)
self.firbank = firbank
return True
return False
@@ -203,7 +209,7 @@ template_bank_table = lsctables.SnglBurstTable.get_table(xmldoc)
# filter bank
#
-head = firbank = pipeparts.mkfirbank(pipeline, head, time_domain = False, block_stride = 4 * options.sample_rate)
+head = firbank = pipeparts.mkfirbank(pipeline, head, fir_matrix = numpy.zeros((len(template_bank_table),int(32*options.sample_rate)),dtype=numpy.float64), block_stride = 4 * options.sample_rate)
#
--
GitLab