From 6a937bc455857f6dfd9c0e32f1a22edb5859bb2c Mon Sep 17 00:00:00 2001
From: Jürgen Knödlseder <jurgen.knodlseder@irap.omp.eu>
Date: Thu, 11 Aug 2016 01:58:00 +0200
Subject: [PATCH] Use at least 30 energy bins per decade in binned cspull

---
 cscripts/cspull.py | 16 +++++++++++++---
 1 file changed, 13 insertions(+), 3 deletions(-)

diff --git a/cscripts/cspull.py b/cscripts/cspull.py
index 4f8fa79..11369e8 100755
--- a/cscripts/cspull.py
+++ b/cscripts/cspull.py
@@ -182,17 +182,27 @@ class cspull(ctools.cscript):
         # If we have a binned obeservation then specify the lower and
         # upper energy limit
         if self._enumbins > 0:
+
+            # Get energy range in TeV
             emin = self['emin'].real()
             emax = self['emax'].real()
+
+            # Set number of energy bins to at least 30 per energy decade
+            enumbins = int((math.log10(self['emax'].real()) -
+                            math.log10(self['emin'].real())) * 30.0)
+            if self._enumbins > enumbins:
+                enumbins = self._enumbins
+
         else:
-            emin = None
-            emax = None
+            emin     = None
+            emax     = None
+            enumbins = 0
 
         # Simulate events
         obs = obsutils.sim(self._obs,
                            emin=emin,
                            emax=emax,
-                           nbins=self._enumbins,
+                           nbins=enumbins,
                            seed=seed,
                            binsz=self._binsz,
                            npix=self._npix,
--
libgit2 0.21.4