diff --git a/pyfstat.py b/pyfstat.py
index 4b044ea23b8f88f1b9a2d9f08c0908f627a0df87..e1a8c6b9465d294014d62108675c4451f01ba0de 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):