From bbd56b7e5afdb996e04a7e38017e87d1a8857409 Mon Sep 17 00:00:00 2001
From: "gregory.ashton" <gregory.ashton@ligo.org>
Date: Wed, 23 Nov 2016 17:08:29 +0100
Subject: [PATCH] Fixes error in DeltaOmega (solid angle) caluclatoin

---
 pyfstat.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/pyfstat.py b/pyfstat.py
index 111f160..391f49c 100755
--- a/pyfstat.py
+++ b/pyfstat.py
@@ -2173,7 +2173,7 @@ class MCMCFollowUpSearch(MCMCSemiCoherentSearch):
             DeltaAlpha = self.get_width_from_prior(self.theta_prior, 'Alpha')
             DeltaDelta = self.get_width_from_prior(self.theta_prior, 'Delta')
             DeltaMid = self.get_mid_from_prior(self.theta_prior, 'Delta')
-            DeltaOmega = np.sin(DeltaMid) * DeltaDelta * DeltaAlpha
+            DeltaOmega = np.sin(np.pi/2 - DeltaMid) * DeltaDelta * DeltaAlpha
             logging.info('Search over Alpha and Delta')
         else:
             logging.info('No sky search requested')
-- 
GitLab