developing test

Luigi Tibaldo
1 parent 474ea1c6
Showing 3 changed files with 34 additions and 16 deletions
cscripts/obsutils.py
test/data/crab_onoff.xml
test/test_cslightcrv.py
@@ -815,4 +815,6 @@ def get_onoff_obs(cls, obs):
  
     onoff_obs = phagen.obs().copy()
  
+    #need to reset statistic for fit???
+
     return onoff_obs
+<?xml version="1.0" standalone="no"?>
+<source_library title="source library">
+  <source name="Crab" type="PointSource">
+    <spectrum type="PowerLaw">
+       <parameter name="Prefactor"   scale="1e-16" value="5.7"  min="1e-07" max="1000.0" free="1"/>
+       <parameter name="Index"       scale="-1"    value="2.48" min="0.0"   max="+5.0"   free="1"/>
+       <parameter name="PivotEnergy" scale="1e6"   value="0.3"  min="0.01"  max="1000.0" free="0"/>
+    </spectrum>
+    <spatialModel type="PointSource">
+      <parameter name="RA"  scale="1.0" value="83.6331" min="-360" max="360" free="0"/>
+      <parameter name="DEC" scale="1.0" value="22.0145" min="-90"  max="90"  free="0"/>
+    </spatialModel>
+  </source>
+  <source name="Background" type="RadialAcceptance" instrument="CTAOnOff">
+    <spectrum type="PowerLaw">
+       <parameter name="Prefactor"   scale="1e-6" value="61.8" min="0.0"  max="1000.0" free="1"/>
+       <parameter name="Index"       scale="-1"   value="1.85" min="0.0"  max="+5.0"   free="1"/>
+       <parameter name="PivotEnergy" scale="1e6"  value="1.0"  min="0.01" max="1000.0" free="0"/>
+    </spectrum>
+    <radialModel type="Gaussian">
+       <parameter name="Sigma" scale="1.0" value="3.0" min="0.01" max="10.0" free="1"/>
+    </radialModel>
+  </source>
+</source_library>
@@ -43,7 +43,8 @@ class Test(test):
         test.__init__(self)
  
         # Add off-axis events for classical analysis
-        self._myevents1 = self._datadir + '/crab_offaxis1.fits'
+        self._offaxis_events = self._datadir + '/crab_offaxis1.fits'
+        self._model_onoff    = self._datadir + '/crab_onoff.xml'
  
         # Return
         return
@@ -240,24 +241,15 @@ class Test(test):
         self._check_light_curve('cslightcrv_py4.fits', 2)
  
         # cslightcrv with classical analysis
-
-        # first we build an observation container to pass the model
-        # to csphagen that will change the background to OnOff
-        obs = gammalib.GObservations()
-        for s, events in enumerate([self._myevents1]):
-            run = gammalib.GCTAObservation(events)
-            run.id(str(s + 1))
-            run.response(self._irf, gammalib.GCaldb('cta', self._caldb))
-            run.model(gammalib.GModels(self._model))
-            obs.append(run)
-
-        # set up cslightcrv
         lcrv = cscripts.cslightcrv()
-        lcrv['inobs']     = obs
+        lcrv['inobs']     = self._offaxis_events
+        lcrv['inmodel']   = self._model_onoff
         lcrv['srcname']   = 'Crab'
+        lcrv['caldb']    = self._caldb
+        lcrv['irf']      = self._irf
         lcrv['tbinalg']   = 'LIN'
-        lcrv['tmin']      = 6.3110887e8
-        lcrv['tmax']      = 6.3110977e8
+        lcrv['tmin']      = 'MJD 58849.0008'
+        lcrv['tmax']      = 'MJD 58849.0112'
         lcrv['tbins']     = 2
         lcrv['method']    = 'ONOFF'
         lcrv['emin']      = 0.1