From c2f959532557391d8490604ea7f53614c853fda4 Mon Sep 17 00:00:00 2001
From: Gregory Ashton <gregory.ashton@ligo.org>
Date: Tue, 11 Oct 2016 21:58:47 +0200
Subject: [PATCH] Adds hack factor to magnify BSGL: needs to be discussed and
understood
---
pyfstat.py | 8 +++++---
1 file changed, 5 insertions(+), 3 deletions(-)
diff --git a/pyfstat.py b/pyfstat.py
index 4b044ea..e1a8c6b 100755
--- a/pyfstat.py
+++ b/pyfstat.py
@@ -283,7 +283,7 @@ class ComputeFstat(object):
self.BSGLSetup = lalpulsar.CreateBSGLSetup(numDetectors,
Fstar0sc,
oLGX,
- False,
+ True,
1)
self.twoFX = np.zeros(10)
self.whatToCompute = (self.whatToCompute +
@@ -304,6 +304,8 @@ class ComputeFstat(object):
argp=None):
""" Returns the twoF fully-coherently at a single point """
+ BSGL_PREFACTOR = 10 * 1 / np.log10(np.exp(1))
+
self.PulsarDopplerParams.fkdot = np.array([F0, F1, F2, 0, 0, 0, 0])
self.PulsarDopplerParams.Alpha = Alpha
self.PulsarDopplerParams.Delta = Delta
@@ -330,7 +332,7 @@ class ComputeFstat(object):
self.twoFX[1] = self.FstatResults.twoFPerDet(1)
BSGL = lalpulsar.ComputeBSGL(twoF, self.twoFX,
self.BSGLSetup)
- return BSGL
+ return BSGL_PREFACTOR * BSGL
self.windowRange.t0 = int(tstart) # TYPE UINT4
self.windowRange.tau = int(tend - tstart) # TYPE UINT4
@@ -355,7 +357,7 @@ class ComputeFstat(object):
BSGL = lalpulsar.ComputeBSGL(2*FS.F_mn.data[0][0], self.twoFX,
self.BSGLSetup)
- return BSGL
+ return BSGL_PREFACTOR * BSGL
class SemiCoherentGlitchSearch(BaseSearchClass, ComputeFstat):
--
GitLab