first crack at csphasecrv

Luigi Tibaldo
1 parent 5e827cf5
Showing 3 changed files with 25 additions and 6 deletions
cscripts/cslightcrv.py
cscripts/csphasecrv.par
cscripts/csphasecrv.py
@@ -74,8 +74,8 @@ class cslightcrv(ctools.csobservation):
         # Initialise some members
         self._srcname = ''
         self._tbins   = gammalib.GGti()
+        self._onoff   = False
         self._stacked = False
-        self._onoff = False
         self._fits    = gammalib.GFits()
         # Return
@@ -21,7 +21,8 @@ inmodel,   f, a, $CTOOLS/share/models/crab.xml,,, &quot;Input model definition XML fi
 srcname,   s, a, Crab,,, "Source name"
 caldb,     s, a, prod2,,, "Calibration database"
 irf,       s, a, South_0.5h,,, "Instrument response function"
-edisp,     b, h, no,,, "Apply energy dispersion?"
+inexclusion,  f, h, NONE,,, "Sky map of exclusion regions for On/Off analysis"
+edisp,     b, h, no,,, "Apply energy dispersion in Cube analysis?"
 outfile,   f, a, phasecrv.fits,,, "Output fits file name."
 #
@@ -30,9 +31,10 @@ outfile,   f, a, phasecrv.fits,,, &quot;Output fits file name.&quot;
 phbinalg,  s, a, LIN,LIN|FILE,, "Algorithm for defining phase bins"
 phbins,    i, a, 5,1,10000, "Number of phase bins"
 phbinfile, f, a, NONE,,, "File defining the phase binning"
+method,    s, a, CUBE, CUBE|ONOFF,, "Analysis method"
 emin,      r, a, 0.1,,, "Lower energy limit of events (TeV)"
 emax,      r, a, 100.0,,, "Upper energy limit of events (TeV)"
-enumbins,  i, a, 0,0,100, "Number of energy bins for binned (0=unbinned)"
+enumbins,  i, a, 0,0,100, "Number of energy bins for binned (0=unbinned for Cube analysis only))"
 coordsys,  s, a, CEL,CEL|GAL,, "Coordinate System"
 proj,      s, a, CAR,AIT|AZP|CAR|GLS|MER|MOL|SFL|SIN|STG|TAN,, "Projection method"
 xref,      r, a, 83.63,0,360, "First coordinate of image center in degrees (RA or galactic l)"
@@ -40,6 +42,14 @@ yref,      r, a, 22.01,-90,90, &quot;Second coordinate of image center in degrees (DE
 nxpix,     i, a, 200,,, "Size of the X axis in pixels"
 nypix,     i, a, 200,,, "Size of the Y axis in pixels"
 binsz,     r, a, 0.02,,, "Pixel size (deg/pixel)"
+srcshape,  s, h, CIRCLE,CIRCLE,, "Shape of the source region"
+rad,       r, a, 0.2,0,180, "Radius of source region circle (deg)"
+bkgmethod, s, h, REFLECTED,REFLECTED,, "Method for background estimation"
+bkgregmin, i, h, 2,1,1000, "Minimum number of background regions"
+maxoffset, r, h, 4.0,,, "Maximum offset of source from camera center to accept observation"
+etruemin,  r, h, 0.01,,, "Minimum true energy (TeV)"
+etruemax,  r, h, 300.0,,, "Maximum true energy (TeV)"
+etruebins, i, h, 30,,, "Number of bins per decade for true energy bins"
 statistic, s, h, DEFAULT,DEFAULT|CSTAT|WSTAT|CHI2,, "Optimization statistic"
 #
@@ -63,6 +63,7 @@ class csphasecrv(ctools.csobservation):
         # [[ph1min,ph1max], [ph2min,ph2max],..]
         self._srcname   = ''
         self._phbins    = [[0.0,1.0]]
+        self._onoff     = False
         self._stacked   = False
         self._fits      = gammalib.GFits()
         self._fitmodels = {}
@@ -92,8 +93,13 @@ class csphasecrv(ctools.csobservation):
         # Get phase boundaries
         self._phbins = self._create_tbounds()
+        # Set On/Off analysis flag and query relevant user parameters
+        self._onoff = self._is_onoff()
+
+        # If cube analysis is selected
         # Set stacked analysis flag and query relevant user parameters
-        self._stacked = self._is_stacked()
+        if not self._onoff:
+            self._stacked = self._is_stacked()
         # Query the hidden parameters, just in case
         self['edisp'].boolean()
@@ -355,10 +361,13 @@ class csphasecrv(ctools.csobservation):
                 select.obs()[i].id(oldid+'_'+str(phbin[0])+'-'+str(phbin[1]))
             obs = select.obs()
-            # If a stacked analysis is requested then bin the events
+            # If an On/Off analysis is requested generate the On/Off observations
+            if self._onoff:
+                obs = obsutils.get_onoff_obs(self, select.obs())
+            # Otherwise, if stacked analysis is requested then bin the events
             # and compute the stacked response functions and setup
             # an observation container with a single stacked observation.
-            if self._stacked:
+            elif self._stacked:
                 obs = obsutils.get_stacked_obs(self, select.obs())
             # Header