Commit 809bc8f6c86ebd2516315ef11b1a24e1cacae138

Authored by Jürgen Knödlseder
1 parent 16621160

Add bkgcube parameter

doc/source/reference_manual/cspull.rst
@@ -26,9 +26,15 @@ General parameters @@ -26,9 +26,15 @@ General parameters
26 ``outfile [file]`` 26 ``outfile [file]``
27 ASCII file containing the individual pull values. 27 ASCII file containing the individual pull values.
28 28
29 -``ntrials [integer]``  
30 - Number of samples for generating the pull distribution.  
31 - 29 +``expcube [file]``
  30 + Exposure cube file (only needed for stacked analysis).
  31 +
  32 +``psfcube [file]``
  33 + PSF cube file (only needed for stacked analysis).
  34 +
  35 +``bkgcube [file]``
  36 + Background cube file (only needed for stacked analysis).
  37 +
32 ``caldb [string]`` 38 ``caldb [string]``
33 Calibration database. 39 Calibration database.
34 40
@@ -38,6 +44,9 @@ General parameters @@ -38,6 +44,9 @@ General parameters
38 ``(edisp = no) [boolean]`` 44 ``(edisp = no) [boolean]``
39 Apply energy dispersion to response computation. 45 Apply energy dispersion to response computation.
40 46
  47 +``ntrials [integer]``
  48 + Number of samples for generating the pull distribution.
  49 +
41 ``ra [real]`` 50 ``ra [real]``
42 Right Ascension of CTA pointing (J2000, in degrees). 51 Right Ascension of CTA pointing (J2000, in degrees).
43 52
scripts/cspull.py
@@ -98,8 +98,11 @@ class cspull(ctools.cscript): @@ -98,8 +98,11 @@ class cspull(ctools.cscript):
98 pars.append(gammalib.GApplicationPar("outfile","f","a","pull.dat","","","Output file name")) 98 pars.append(gammalib.GApplicationPar("outfile","f","a","pull.dat","","","Output file name"))
99 pars.append(gammalib.GApplicationPar("expcube","s","a","NONE","","","Exposure cube file (only needed for stacked analysis)")) 99 pars.append(gammalib.GApplicationPar("expcube","s","a","NONE","","","Exposure cube file (only needed for stacked analysis)"))
100 pars.append(gammalib.GApplicationPar("psfcube","s","a","NONE","","","PSF cube file (only needed for stacked analysis)")) 100 pars.append(gammalib.GApplicationPar("psfcube","s","a","NONE","","","PSF cube file (only needed for stacked analysis)"))
  101 + pars.append(gammalib.GApplicationPar("bkgcube","s","a","NONE","","","Background cube file (only needed for stacked analysis)"))
101 pars.append(gammalib.GApplicationPar("caldb","s","a","dummy","","","Calibration database")) 102 pars.append(gammalib.GApplicationPar("caldb","s","a","dummy","","","Calibration database"))
102 pars.append(gammalib.GApplicationPar("irf","s","a","cta_dummy_irf","","","Instrument response function")) 103 pars.append(gammalib.GApplicationPar("irf","s","a","cta_dummy_irf","","","Instrument response function"))
  104 + pars.append(gammalib.GApplicationPar("edisp","b","h","no","","","Apply energy dispersion?"))
  105 + pars.append(gammalib.GApplicationPar("ntrials","i","a","10","","","Number of trials"))
103 pars.append(gammalib.GApplicationPar("ra","r","a","83.6331","0","360","RA of pointing (deg)")) 106 pars.append(gammalib.GApplicationPar("ra","r","a","83.6331","0","360","RA of pointing (deg)"))
104 pars.append(gammalib.GApplicationPar("dec","r","a","22.0145","-90","90","Dec of pointing (deg)")) 107 pars.append(gammalib.GApplicationPar("dec","r","a","22.0145","-90","90","Dec of pointing (deg)"))
105 pars.append(gammalib.GApplicationPar("emin","r","a","0.1","","","Lower energy limit (TeV)")) 108 pars.append(gammalib.GApplicationPar("emin","r","a","0.1","","","Lower energy limit (TeV)"))
@@ -113,8 +116,6 @@ class cspull(ctools.cscript): @@ -113,8 +116,6 @@ class cspull(ctools.cscript):
113 pars.append(gammalib.GApplicationPar("rad","r","h","5.0","0","180","Radius of ROI (deg)")) 116 pars.append(gammalib.GApplicationPar("rad","r","h","5.0","0","180","Radius of ROI (deg)"))
114 pars.append(gammalib.GApplicationPar("pattern","s","h","single","","","Observation pattern (single/four)")) 117 pars.append(gammalib.GApplicationPar("pattern","s","h","single","","","Observation pattern (single/four)"))
115 pars.append(gammalib.GApplicationPar("offset","r","h","1.5","","","Observation pattern offset (deg)")) 118 pars.append(gammalib.GApplicationPar("offset","r","h","1.5","","","Observation pattern offset (deg)"))
116 - pars.append(gammalib.GApplicationPar("ntrials","i","a","10","","","Number of trials"))  
117 - pars.append(gammalib.GApplicationPar("edisp","b","h","no","","","Apply energy dispersion?"))  
118 pars.append_standard() 119 pars.append_standard()
119 pars.append(gammalib.GApplicationPar("logfile","f","h","cspull.log","","","Log filename")) 120 pars.append(gammalib.GApplicationPar("logfile","f","h","cspull.log","","","Log filename"))
120 pars.save(parfile) 121 pars.save(parfile)