Added phasecut parameters and expscape>0 check

fermiPy · Sep 29, 2016 · 7cb86c9 · 7cb86c9
1 parent 7f8c3cf
commit 7cb86c9
Show file tree

Hide file tree

Showing 2 changed files with 7 additions and 1 deletion.
diff --git a/fermipy/defaults.py b/fermipy/defaults.py
@@ -37,6 +37,8 @@ def make_default_dict(d):
     'evclass': (None, 'Event class selection.', int),
     'evtype': (None, 'Event type selection.', int),
     'convtype': (None, 'Conversion type selection.', int),
+    'phasemin': (None, 'Minimum pulsar phase', float),
+    'phasemax': (None, 'Maximum pulsar phase', float),
     'target': (None, 'Choose an object on which to center the ROI.  '
                      'This option takes precendence over ra/dec or glon/glat.', str),
     'ra': (None, '', float),

diff --git a/fermipy/gtanalysis.py b/fermipy/gtanalysis.py
@@ -4822,6 +4822,8 @@ def _select_data(self, overwrite=False):
                            dec=self.roi.skydir.dec.deg,
                            rad=self.config['selection']['radius'],
                            convtype=self.config['selection']['convtype'],
+                           phasemin=self.config['selection']['phasemin'],
+                           phasemax=self.config['selection']['phasemax'],
                            evtype=self.config['selection']['evtype'],
                            evclass=self.config['selection']['evclass'],
                            tmin=self.config['selection']['tmin'],
@@ -4909,11 +4911,13 @@ def _scale_srcmap(self, scale_map):
                 continue
 
             scale = scale_map[hdu.name]
+            if scale<1e-20:
+                self.logger.warning("The expscale parameter was zero, setting it to 1e-8")
+                scale = 1e-8
             if 'SRCMAP_SCALE' in hdu.header:
                 old_scale = hdu.header['SRCMAP_SCALE']
             else:
                 old_scale = 1.0
-
             hdu.data *= scale/old_scale
             hdu.header['SRCMAP_SCALE'] = scale