From 6d385885605f7ae77a39d2eab4bf078ff3bd3873 Mon Sep 17 00:00:00 2001
From: Gregory Ashton <gregory.ashton@ligo.org>
Date: Wed, 20 Dec 2017 11:13:23 +0100
Subject: [PATCH] Fix up SemiCoherentGlitchSearch
- Default to 1 glitch
- Only try to calculate segment 2F if length > 1800s
---
pyfstat/core.py | 12 ++++++------
1 file changed, 6 insertions(+), 6 deletions(-)
diff --git a/pyfstat/core.py b/pyfstat/core.py
index e336acc..5c14753 100755
--- a/pyfstat/core.py
+++ b/pyfstat/core.py
@@ -1033,7 +1033,7 @@ class SemiCoherentGlitchSearch(ComputeFstat):
@helper_functions.initializer
def __init__(self, label, outdir, tref, minStartTime, maxStartTime,
- nglitch=0, sftfilepattern=None, theta0_idx=0, BSGL=False,
+ nglitch=1, sftfilepattern=None, theta0_idx=0, BSGL=False,
minCoverFreq=None, maxCoverFreq=None, assumeSqrtSX=None,
detectors=None, SSBprec=None, injectSources=None):
"""
@@ -1085,11 +1085,11 @@ class SemiCoherentGlitchSearch(ComputeFstat):
twoFSum = 0
for i, theta_i_at_tref in enumerate(thetas):
ts, te = tboundaries[i], tboundaries[i+1]
-
- twoFVal = self.get_fullycoherent_twoF(
- ts, te, theta_i_at_tref[1], theta_i_at_tref[2],
- theta_i_at_tref[3], Alpha, Delta)
- twoFSum += twoFVal
+ if te - te > 1800:
+ twoFVal = self.get_fullycoherent_twoF(
+ ts, te, theta_i_at_tref[1], theta_i_at_tref[2],
+ theta_i_at_tref[3], Alpha, Delta)
+ twoFSum += twoFVal
if np.isfinite(twoFSum):
return twoFSum
--
GitLab