From ac379959f577fd9fbcbe1f2259926937b8b40d8c Mon Sep 17 00:00:00 2001
From: Aaron Viets <aaron.viets@ligo.org>
Date: Fri, 2 Mar 2018 11:21:48 -0800
Subject: [PATCH] calibration_parts.py:  Trying a new method to remove DC
 component from h(t)

---
 gstlal-calibration/python/calibration_parts.py | 14 ++++++++++++++
 1 file changed, 14 insertions(+)

diff --git a/gstlal-calibration/python/calibration_parts.py b/gstlal-calibration/python/calibration_parts.py
index f670fa1050..9c3134ad3d 100644
--- a/gstlal-calibration/python/calibration_parts.py
+++ b/gstlal-calibration/python/calibration_parts.py
@@ -161,6 +161,20 @@ def list_srcs(pipeline, *args):
 		out.append(src)
 	return tuple(out)
 
+def removeDC(pipeline, head, caps):
+	head = pipeparts.mktee(pipeline, head)
+	pipeparts.mknxydumpsink(pipeline, head, "head.txt")
+	DC = mkresample(pipeline, head, 3, True, "audio/x-raw, rate=16")
+	DC = pipeparts.mktee(pipeline, DC)
+	pipeparts.mknxydumpsink(pipeline, DC, "DC1.txt")
+	DC = pipeparts.mkgeneric(pipeline, DC, "lal_smoothkappas", default_kappa_re = 0, array_size = 1, avg_array_size = 64)
+	DC = pipeparts.mktee(pipeline, DC)
+	pipeparts.mknxydumpsink(pipeline, DC, "DC.txt")
+	DC = mkresample(pipeline, DC, 5, True, caps)
+	DC = pipeparts.mkaudioamplify(pipeline, DC, -1)
+
+	return mkadder(pipeline, list_srcs(pipeline, mkqueue(pipeline, head, 0), mkqueue(pipeline, DC, 0)))
+
 #
 # Calibration factor related functions
 #
-- 
GitLab