From a9d43502ad01c24135d3771c26545424caa168bf Mon Sep 17 00:00:00 2001
From: Leo Tsukada <leo.tsukada@ligo.org>
Date: Mon, 3 Apr 2023 13:54:30 -0700
Subject: [PATCH] bin/gstlal_inspiral_bank_splitter : add extra field for
 manifest file

---
 gstlal-inspiral/bin/gstlal_inspiral_bank_splitter | 12 ++++++++++++
 1 file changed, 12 insertions(+)

diff --git a/gstlal-inspiral/bin/gstlal_inspiral_bank_splitter b/gstlal-inspiral/bin/gstlal_inspiral_bank_splitter
index aba17bf65c..b1c1c69291 100755
--- a/gstlal-inspiral/bin/gstlal_inspiral_bank_splitter
+++ b/gstlal-inspiral/bin/gstlal_inspiral_bank_splitter
@@ -297,7 +297,9 @@ svd_groups.sort(key=sort_svd, reverse = options.sort_by in ["template_duration"
 
 for n, svd in enumerate(svd_groups, start=svd_bin_start):
 	svd_mchirps = []
+	svd_mtotals = []
 	svd_etas = []
+	svd_mratios = []
 	svd_durs = []
 	svd_bw = []
 	svd_mu1 = []
@@ -330,6 +332,8 @@ for n, svd in enumerate(svd_groups, start=svd_bin_start):
 		else:
 			temp_id_output |= temp_id_splitbank
 		svd_mchirps.extend([r.mchirp for r in clipped_rows])
+		svd_mtotals.extend([r.mtotal for r in clipped_rows])
+		svd_mratios.extend([r.mass1 / r.mass2 for r in clipped_rows])
 		svd_etas.extend([r.eta for r in clipped_rows])
 		svd_durs.extend([r.template_duration for r in clipped_rows])
 		if options.sort_by == "mu":
@@ -346,10 +350,18 @@ for n, svd in enumerate(svd_groups, start=svd_bin_start):
 		"median_mchirp": numpy.median(svd_mchirps),
 		"min_mchirp": numpy.min(svd_mchirps),
 		"max_mchirp": numpy.max(svd_mchirps),
+		"mean_mtotal": numpy.mean(svd_mtotals),
+		"median_mtotal": numpy.median(svd_mtotals),
+		"min_mtotal": numpy.min(svd_mtotals),
+		"max_mtotal": numpy.max(svd_mtotals),
 		"mean_eta": numpy.mean(svd_etas),
 		"median_eta": numpy.median(svd_etas),
 		"min_eta": numpy.min(svd_etas),
 		"max_eta": numpy.max(svd_etas),
+		"mean_mratio": numpy.mean(svd_mratios),
+		"median_mratio": numpy.median(svd_mratios),
+		"min_mratio": numpy.min(svd_mratios),
+		"max_mratio": numpy.max(svd_mratios),
 		"mean_dur": numpy.mean(svd_durs),
 		"median_dur": numpy.median(svd_durs),
 		"min_dur": numpy.min(svd_durs),
-- 
GitLab