From 8e3068f32ba5b4efe52779d7e7bbdabc8fae93cf Mon Sep 17 00:00:00 2001
From: Lucie Gerard <lgerard@wgs18.zeuthen.desy.de>
Date: Fri, 18 Sep 2015 14:00:34 +0200
Subject: [PATCH] new branck 1526-sensitivity_tev2erg_conversion, modified cssens.py. Sensitivity definition, one of the erg_mean is replaced by e_mean

---
 cscripts/cssens.py | 12 ++++++++++--
 1 file changed, 10 insertions(+), 2 deletions(-)

diff --git a/cscripts/cssens.py b/cscripts/cssens.py
index c916953..62be69a 100755
--- a/cscripts/cssens.py
+++ b/cscripts/cssens.py
@@ -52,6 +52,7 @@ class cssens(ctools.cscript):
         self.m_dec  = None
 
         # Make sure that parfile exists
+        print "making sure the parfile exists"
         file = self.parfile()
 
         # Initialise application
@@ -88,6 +89,8 @@ class cssens(ctools.cscript):
         """
         # Set parfile name
         parfile = self.name+".par"
+        
+        print "parfile ",parfile
 
         try:
             pars = gammalib.GApplicationPars(parfile)
@@ -233,9 +236,11 @@ class cssens(ctools.cscript):
 
             # Set energies
             if self.m_type == "Differential":
+                print "Differential "
                 emin  = self.m_ebounds.emin(ieng)
                 emax  = self.m_ebounds.emax(ieng)
             elif self.m_type == "Integral":
+                print "Integral "
                 emin  = self.m_ebounds.emin(ieng)
                 emax  = self.m_ebounds.emax()
             else:
@@ -469,6 +474,9 @@ class cssens(ctools.cscript):
             sim = obsutils.sim(obs, nbins=self.m_enumbins, seed=iter,
                                binsz=self.m_binsz, npix=self.m_npix,
                                log=self.m_log, debug=self.m_debug, edisp=self.m_edisp)
+            sim = obsutils.sim(obs, nbins=self.m_enumbins, seed=iter,
+                               binsz=self.m_binsz, npix=self.m_npix,
+                               log=self.m_log, debug=self.m_debug, edisp=self.m_edisp)
 
             # Determine number of events in simulation
             nevents = 0.0
@@ -576,7 +584,7 @@ class cssens(ctools.cscript):
             energy      = gammalib.GEnergy(e_mean, "TeV")
             time        = gammalib.GTime()
             sensitivity = result_all[self.m_srcname].spectral().eval(energy, time) * \
-                          erg_mean*erg_mean * 1.0e6
+                          e_mean*erg_mean * 1.0e6
 
             # Compute flux correction factor based on average TS
             correct = 1.0
@@ -727,7 +735,7 @@ class cssens(ctools.cscript):
                     self.log(str(par)+"\n")
 
         # Store result
-        result = {'loge': loge, 'emin': emin.TeV(), 'emax': emax.TeV(), \
+        result = {'loge': loge, 'emin': math.log10(emin.TeV()), 'emax': math.log10(emax.TeV()), \
                   'crab_flux': crab_flux, 'photon_flux': photon_flux, \
                   'energy_flux': energy_flux, \
                   'sensitivity': sensitivity}
--
libgit2 0.21.4