Modify parameter limits to comply with GammaLib parameter range checking

Jürgen Knödlseder
1 parent 12d0bd2b
Showing 13 changed files with 33 additions and 33 deletions
scripts/cspull.py
scripts/csresmap.py
scripts/cssens.py
scripts/cstsdist.py
src/ctbin/ctbin.par
src/ctbkgcube/ctbkgcube.par
src/ctcubemask/ctcubemask.par
src/ctexpcube/ctexpcube.par
src/ctmodel/ctmodel.par
src/ctobssim/ctobssim.par
src/ctpsfcube/ctpsfcube.par
src/ctselect/ctselect.par
src/ctskymap/ctskymap.par
@@ -109,8 +109,8 @@ class cspull(ctools.cscript):
             pars.append(gammalib.GApplicationPar("enumbins","i","a","0","","","Number of energy bins (0=unbinned)"))
             pars.append(gammalib.GApplicationPar("npix","i","h","200","","","Number of pixels for binned"))
             pars.append(gammalib.GApplicationPar("binsz","r","h","0.05","","","Pixel size for binned (deg/pixel)"))
-            pars.append(gammalib.GApplicationPar("deadc","r","h","0.95","","","Deadtime correction factor"))
-            pars.append(gammalib.GApplicationPar("rad","r","h","5.0","","","Radius of ROI (deg)"))         
+            pars.append(gammalib.GApplicationPar("deadc","r","h","0.95","0","1","Deadtime correction factor"))
+            pars.append(gammalib.GApplicationPar("rad","r","h","5.0","0","180","Radius of ROI (deg)"))         
             pars.append(gammalib.GApplicationPar("pattern","s","h","single","","","Observation pattern (single/four)"))
             pars.append(gammalib.GApplicationPar("offset","r","h","1.5","","","Observation pattern offset (deg)"))
             pars.append(gammalib.GApplicationPar("ntrials","i","a","10","","","Number of trials"))
@@ -105,8 +105,8 @@ class csresmap(ctools.cscript):
             pars.append(gammalib.GApplicationPar("outmap","f","a","resmap.fits","","","Output residual map"))
             pars.append(gammalib.GApplicationPar("caldb","s","a","dummy","","","Calibration database"))
             pars.append(gammalib.GApplicationPar("irf","s","a","cta_dummy_irf","","","Instrument response function"))
-            pars.append(gammalib.GApplicationPar("emin","r","h","0.1","0","","Lower energy limit (TeV)"))
-            pars.append(gammalib.GApplicationPar("emax","r","h","100.0","0","","Upper energy limit (TeV)"))
+            pars.append(gammalib.GApplicationPar("emin","r","h","0.1","","","Lower energy limit (TeV)"))
+            pars.append(gammalib.GApplicationPar("emax","r","h","100.0","","","Upper energy limit (TeV)"))
             pars.append(gammalib.GApplicationPar("enumbins","i","h","20","","","Number of energy bins"))
             pars.append(gammalib.GApplicationPar("ebinalg","s","h","LOG","LIN|LOG|FILE","","Binning algorithm"))
             pars.append(gammalib.GApplicationPar("coordsys","s","a","CEL","CEL|GAL","","Coordinate System"))
@@ -109,7 +109,7 @@ class cssens(ctools.cscript):
             pars.append(gammalib.GApplicationPar("caldb","s","a","$CTOOLS/share/caldb/data/cta/dummy","","","Calibration database"))
             pars.append(gammalib.GApplicationPar("irf","s","a","cta_dummy_irf","","","Instrument response function"))
             pars.append(gammalib.GApplicationPar("type","s","a","point","","","Source model type (point/gauss/shell/disk)"))
-            pars.append(gammalib.GApplicationPar("offset","r","a","0.0","0.0","","Source offset angle (deg)"))
+            pars.append(gammalib.GApplicationPar("offset","r","a","0.0","","","Source offset angle (deg)"))
             pars.append(gammalib.GApplicationPar("bkg","s","a","$CTOOLS/share/models/bkg_dummy.txt","","","Background model file function (none=power law for E)"))
             pars.append(gammalib.GApplicationPar("duration","r","a","180000.0","","","Effective exposure time (s)"))
             pars.append(gammalib.GApplicationPar("rad","r","a","5.0","","","Radius of ROI (deg)"))
@@ -100,8 +100,8 @@ class cstsdist(ctools.cscript):
             pars.append(gammalib.GApplicationPar("psfcube","s","a","NONE","","","PSF cube file (only needed for stacked analysis)"))
             pars.append(gammalib.GApplicationPar("caldb","s","a","dummy","","","Calibration database"))
             pars.append(gammalib.GApplicationPar("irf","s","a","cta_dummy_irf","","","Instrument response function"))
-            pars.append(gammalib.GApplicationPar("ra","r","h","83.6331","0","360","RA of pointing (deg)"))
-            pars.append(gammalib.GApplicationPar("dec","r","h","22.0145","-90","90","Dec of pointing (deg)"))
+            pars.append(gammalib.GApplicationPar("ra","r","h","83.63","0","360","RA of pointing (deg)"))
+            pars.append(gammalib.GApplicationPar("dec","r","h","22.01","-90","90","Dec of pointing (deg)"))
             pars.append(gammalib.GApplicationPar("emin","r","a","0.1","","","Lower energy limit (TeV)"))
             pars.append(gammalib.GApplicationPar("emax","r","a","100.0","","","Upper energy limit (TeV)"))
             pars.append(gammalib.GApplicationPar("tmin","r","h","0.0","","","Start time (MET in s)"))
@@ -114,7 +114,7 @@ class cstsdist(ctools.cscript):
             pars.append(gammalib.GApplicationPar("ntrials","i","a","10","","","Number of trials"))
             pars.append(gammalib.GApplicationPar("type","s","a","point","","","Source model type (point/gauss/shell/disk)"))
             pars.append(gammalib.GApplicationPar("index","r","h","-2.48","","","Spectral index"))       
-            pars.append(gammalib.GApplicationPar("offset","r","a","0.0","0.0","","Source offset angle (deg)"))
+            pars.append(gammalib.GApplicationPar("offset","r","a","0.0","","","Source offset angle (deg)"))
             pars.append(gammalib.GApplicationPar("bkg","s","a","$CTOOLS/share/models/bkg_dummy.txt","","","Background model file function"))
             pars.append_standard()
             pars.append(gammalib.GApplicationPar("logfile","f","h","cstsdist.log","","","Log filename"))
@@ -23,8 +23,8 @@ outcube, f, a, &quot;cntcube.fits&quot;,,, &quot;Output counts cube&quot;
 # Energy binning parameters
 #==========================
 ebinalg,  s, a, "LOG",FILE|LIN|LOG,, "Algorithm for defining energy bins"
-emin,     r, a, 0.1,0,, "Start value for first energy bin in TeV"
-emax,     r, a, 100.0,0,, "Stop value for last energy bin in TeV"
+emin,     r, a, 0.1,,, "Start value for first energy bin in TeV"
+emax,     r, a, 100.0,,, "Stop value for last energy bin in TeV"
 enumbins, i, a, 20,,, "Number of energy bins"
 ebinfile, f, a, "NONE",,, "Name of the file containing the energy bin definition"
 #
@@ -28,8 +28,8 @@ irf,      s, a, &quot;cta_dummy_irf&quot;,,, &quot;Instrument response function&quot;
 # Energy binning parameters
 #==========================
 ebinalg,  s, h, "LOG",FILE|LIN|LOG,, "Algorithm for defining energy bins"
-emin,     r, a, 0.1,0,, "Start value for first energy bin in TeV"
-emax,     r, a, 100.0,0,, "Stop value for last energy bin in TeV"
+emin,     r, a, 0.1,,, "Lower energy limit (TeV)"
+emax,     r, a, 100.0,,, "Upper energy limit (TeV)"
 enumbins, i, a, 20,,, "Number of energy bins"
 ebinfile, f, h, "NONE",,, "Name of the file containing the energy bin definition"
 #
@@ -28,8 +28,8 @@ usepnt, b, h, no,,, &quot;Use pointing instead of RA/DEC parameters?&quot;
 ra,     r, a, 83.63,0,360, "RA for ROI centre (degrees)"
 dec,    r, a, 22.01,-90,90, "Dec for ROI centre (degrees)"
 rad,    r, a, 3.0,0,180, "Radius of ROI (degrees)"
-emin,   r, a, 0.1,0,, "Lower energy limit (TeV)"
-emax,   r, a, 100.0,0,, "Upper energy limit (TeV)"
+emin,   r, a, 0.1,,, "Lower energy limit (TeV)"
+emax,   r, a, 100.0,,, "Upper energy limit (TeV)"
 #
 # Standard parameters
 #====================
@@ -27,8 +27,8 @@ edisp,    b, h, no,,, &quot;Apply energy dispersion?&quot;
 # Energy binning parameters
 #==========================
 ebinalg,  s, h, "LOG", FILE|LIN|LOG,,"Algorithm for defining energy bins"
-emin,     r, a, 0.1,0,,"Start value for first energy bin in TeV"
-emax,     r, a, 100.0,0,, "Stop value for last energy bin in TeV"
+emin,     r, a, 0.1,,,"Lower energy limit (TeV)"
+emax,     r, a, 100.0,,, "Upper energy limit (TeV)"
 enumbins, i, a, 20,,, "Number of energy bins"
 ebinfile, f, h, "NONE",,, "Name of the file containing the energy bin definition"
 #
@@ -32,15 +32,15 @@ irf,      s, a, &quot;cta_dummy_irf&quot;,,, &quot;Instrument response function&quot;
 ra,       r, a, 83.63,0,360, "RA of pointing (degrees)"
 dec,      r, a, 22.01,-90,90, "Dec of pointing (degrees)"
 rad,      r, a, 5.0,0,180, "Radius of FOV (degrees)"
-tmin,     r, a, 0.0,0,, "Start time (MET in s)"
-tmax,     r, a, 1800.0,0,, "End time (MET in s)"
+tmin,     r, a, 0.0,,, "Start time (MET in s)"
+tmax,     r, a, 1800.0,,, "End time (MET in s)"
 deadc,    r, h, 0.95,0,1, "Average deadtime correction factor"
 #
 # Energy binning parameters
 #==========================
 ebinalg,  s, h, "LOG",FILE|LIN|LOG,, "Algorithm for defining energy bins"
-emin,     r, a, 0.1,0,, "Start value for first energy bin in TeV"
-emax,     r, a, 100.0,0,, "Stop value for last energy bin in TeV"
+emin,     r, a, 0.1,,, "Start value for first energy bin in TeV"
+emax,     r, a, 100.0,,, "Stop value for last energy bin in TeV"
 enumbins, i, a, 20,,, "Number of energy bins"
 ebinfile, f, h, "NONE",,, "Name of the file containing the energy bin definition"
 #
@@ -33,10 +33,10 @@ seed,   i, h, 1,,, &quot;Random number generator seed&quot;
 ra,     r, a, 83.63,0,360, "RA of pointing (degrees)"
 dec,    r, a, 22.01,-90,90, "Dec of pointing (degrees)"
 rad,    r, a, 5.0,0,180, "Radius of FOV (degrees)"
-tmin,   r, a, 0.0,0,, "Start time (MET in s)"
-tmax,   r, a, 1800.0,0,, "End time (MET in s)"
-emin,   r, a, 0.1,0,, "Lower energy limit (TeV)"
-emax,   r, a, 100.0,0,, "Upper energy limit (TeV)"
+tmin,   r, a, 0.0,,, "Start time (MET in s)"
+tmax,   r, a, 1800.0,,, "End time (MET in s)"
+emin,   r, a, 0.1,,, "Lower energy limit (TeV)"
+emax,   r, a, 100.0,,, "Upper energy limit (TeV)"
 edisp,  b, h, no,,, "Apply energy dispersion?"
 deadc,  r, h, 0.95,0,1, "Average deadtime correction factor"
 #
@@ -27,8 +27,8 @@ edisp,    b, h, no,,, &quot;Apply energy dispersion?&quot;
 # Energy binning parameters
 #==========================
 ebinalg,  s, h, "LOG", FILE|LIN|LOG,, "Algorithm for defining energy bins"
-emin,     r, a, 0.1,,, "Start value for first energy bin in TeV"
-emax,     r, a, 100.0,,, "Stop value for last energy bin in TeV"
+emin,     r, a, 0.1,,, "Lower energy limit (TeV)"
+emax,     r, a, 100.0,,, "Upper energy limit (TeV)"
 enumbins, i, a, 20,,, "Number of energy bins"
 ebinfile, f, h, "NONE",,, "Name of the file containing the energy bin definition"
 #
@@ -45,7 +45,7 @@ proj,     s, a, &quot;CAR&quot;, AIT|AZP|CAR|MER|MOL|STG|TAN,, &quot;Projection method&quot;
 #
 # Angular separation binning parameters
 #======================================
-amax,     r, h, 0.3,0,, "Upper bound of angular separation between true and measued photon direction (degrees)"
+amax,     r, h, 0.3,,, "Upper bound of angular separation between true and measued photon direction (degrees)"
 anumbins, i, h, 200,,, "Number of angular separation bins"
 #
 # Standard parameters
@@ -27,10 +27,10 @@ usepnt,  b, h, no,,, &quot;Use pointing instead of RA/DEC parameters?&quot;
 ra,      r, a, 83.63,0,360, "RA for ROI centre (degrees)"
 dec,     r, a, 22.01,-90,90, "Dec for ROI centre (degrees)"
 rad,     r, a, 3.0,0,180, "Radius of ROI (degrees)"
-tmin,    r, a, 0.0,0,, "Start time (CTA MET in seconds)"
-tmax,    r, a, 0.0,0,, "End time (CTA MET in seconds)"
-emin,    r, a, 0.1,0,, "Lower energy limit (TeV)"
-emax,    r, a, 100.0,0,, "Upper energy limit (TeV)"
+tmin,    r, a, 0.0,,, "Start time (CTA MET in seconds)"
+tmax,    r, a, 0.0,,, "End time (CTA MET in seconds)"
+emin,    r, a, 0.1,,, "Lower energy limit (TeV)"
+emax,    r, a, 100.0,,, "Upper energy limit (TeV)"
 expr,    s, h, "",,, "Additional selection expression"
 usethres,s, h, "NONE",NONE|DEFAULT|USER,, "Energy threshold type"
 #
@@ -22,8 +22,8 @@ outmap, f, a, &quot;skymap.fits&quot;,,, &quot;Output file name&quot;
 #
 # Event selection parameters
 #===========================
-emin, r, a, 0.1,0,, "Minimum energy in TeV"
-emax, r, a, 100.0,0,, "Maximum energy in TeV"
+emin, r, a, 0.1,,, "Lower energy limit (TeV)"
+emax, r, a, 100.0,,, "Upper energy limit (TeV)"
 #
 # Spatial binning parameters
 #===========================