Update default cssens.py parameters.

Jürgen Knödlseder
1 parent 7057753b
Showing 2 changed files with 7 additions and 7 deletions
cscripts/cssens.par
cscripts/cssens.py
@@ -31,9 +31,9 @@ outfile, f, h, sensitivity.dat,,, &quot;Output sensitivity ASCII file&quot;
 offset,   r, h, 0.0,,, "Source offset angle (deg)"
 duration, r, a, 180000.0,,, "Effective exposure time (s)"
 rad,      r, a, 5.0,,, "Radius of ROI (deg)"
-emin,     r, a, 0.030,,, "Lower energy limit (TeV)"
-emax,     r, a, 180.0,,, "Upper energy limit (TeV)"
-bins,     i, a, 20,,, "Number of energy bins for differential sensitivity computation"
+emin,     r, a, 0.03163,,, "Lower energy limit (TeV)"
+emax,     r, a, 125.8925,,, "Upper energy limit (TeV)"
+bins,     i, a, 18,,, "Number of energy bins for differential sensitivity computation"
 enumbins, i, h, 0,,, "Number of energy bins for binned analysis (0=unbinned)"
 npix,     i, h, 200,,, "Number of pixels for binned analysis"
 binsz,    r, h, 0.05,,, "Pixel size for binned analysis (deg/pixel)"
@@ -289,10 +289,10 @@ class cssens(ctools.cscript):
         # Set energy boundaries
         self._set_obs_ebounds(emin, emax)
  
-        # Determine energy boundaries from first observation in the container
-        loge      = math.log10(math.sqrt(emin.TeV()*emax.TeV()))
-        e_mean    = math.pow(10.0, loge)
-        erg_mean  = e_mean * tev2erg
+        # Determine mean energy for energy boundary
+        e_mean   = math.sqrt(emin.TeV()*emax.TeV())
+        loge     = math.log10(e_mean)
+        erg_mean = e_mean * tev2erg
  
         # Compute Crab unit. This is the factor with which the Prefactor needs
         # to be multiplied to get 1 Crab.