From d57238122ac96c95925c51a7d3538124e5bdc8ec Mon Sep 17 00:00:00 2001
From: Jürgen Knödlseder <jurgen.knodlseder@irap.omp.eu>
Date: Tue, 26 Feb 2019 11:12:55 +0100
Subject: [PATCH] Add n_off=0 special case to significance computation in ctskymap

---
 src/ctskymap/ctskymap.cpp | 5 +++++
 1 file changed, 5 insertions(+), 0 deletions(-)

diff --git a/src/ctskymap/ctskymap.cpp b/src/ctskymap/ctskymap.cpp
index 440b3e6..ab63f44 100755
--- a/src/ctskymap/ctskymap.cpp
+++ b/src/ctskymap/ctskymap.cpp
@@ -1631,6 +1631,11 @@ double ctskymap::sigma_li_ma(const double& n_on,
         sigma = -2.0 * n_off * std::log(1.0+alpha);
     }
 
+    // Handle special case of no Off-counts
+    else if (n_off == 0.0) {
+        sigma = 2.0 * n_on * std::log((1.0+alpha)/alpha);
+    }
+
     // ... otherwise handle general case
     else {
         sigma = (n_on < (alpha*n_off) ? -2.0 : 2.0) *
--
libgit2 0.21.4