Use new response functions in example scripts

Jürgen Knödlseder
1 parent a5f21c87
Showing 7 changed files with 14 additions and 12 deletions
examples/check_models.py
examples/pipeline_binned_disk.py
examples/pipeline_binned_mem.py
examples/pipeline_stacked_disk.py
examples/pipeline_stacked_mem.py
examples/pipeline_unbinned_disk.py
examples/pipeline_unbinned_mem.py
@@ -41,8 +41,8 @@ def pipeline(model_name):
 	FITS files on disk. All data is only hold in memory.
 	"""
 	# Set script parameters
-	caldb       = "dummy"
-	irf         = "cta_dummy_irf"
+	caldb       = "prod2"
+	irf         = "South_50h"
 	ra          =   83.63
 	dec         =   22.01
 	rad_sim     =   10.0
@@ -36,7 +36,7 @@ from ctools import obsutils
 def setup_observations(pattern="four", ra=83.63, dec=22.01, offset=1.5, \
                        emin=0.1, emax=100.0, rad=5.0, duration=1800.0, \
                        deadc=0.95, \
-                       caldb="dummy", irf="cta_dummy_irf"):
+                       caldb="prod2", irf="South_50h"):
     """
     Returns an observation container.
@@ -96,7 +96,7 @@ def run_pipeline(obs, emin=0.1, emax=100.0, \
                  enumbins=20, nxpix=200, nypix=200, binsz=0.02, \
                  coordsys="CEL", proj="CAR", \
                  model="${CTOOLS}/share/models/crab.xml", \
-                 caldb="dummy", irf="cta_dummy_irf", \
+                 caldb="prod2", irf="South_50h", \
                  debug=False):
     """
     Simulation and binned analysis pipeline.
@@ -165,6 +165,7 @@ def run_pipeline(obs, emin=0.1, emax=100.0, \
     like["outmodel"].filename("fit_results.xml")
     like["expcube"].filename("NONE")
     like["psfcube"].filename("NONE")
+    like["bkgcube"].filename("NONE")
     like["caldb"].string(caldb)
     like["irf"].string(irf)
     like["debug"].boolean(True) # Switch this always on for results in console
@@ -36,7 +36,7 @@ from ctools import obsutils
 def setup_observations(pattern="four", ra=83.63, dec=22.01, offset=1.5, \
                        emin=0.1, emax=100.0, rad=5.0, duration=1800.0, \
                        deadc=0.95, \
-                       caldb="dummy", irf="cta_dummy_irf"):
+                       caldb="prod2", irf="South_50h"):
     """
     Returns an observation container.
@@ -36,7 +36,7 @@ from ctools import obsutils
 def setup_observations(pattern="four", ra=83.63, dec=22.01, offset=1.5, \
                        emin=0.1, emax=100.0, rad=5.0, duration=1800.0, \
                        deadc=0.95, \
-                       caldb="dummy", irf="cta_dummy_irf"):
+                       caldb="prod2", irf="South_50h"):
     """
     Returns an observation container.
@@ -96,7 +96,7 @@ def run_pipeline(obs, ra=83.63, dec=22.01, emin=0.1, emax=100.0, \
                  enumbins=20, nxpix=200, nypix=200, binsz=0.02, \
                  coordsys="CEL", proj="CAR", \
                  model="${CTOOLS}/share/models/crab.xml", \
-                 caldb="dummy", irf="cta_dummy_irf", \
+                 caldb="prod2", irf="South_50h", \
                  debug=False):
     """
     Simulation and stacked analysis pipeline.
@@ -210,6 +210,7 @@ def run_pipeline(obs, ra=83.63, dec=22.01, emin=0.1, emax=100.0, \
     like["outmodel"].filename("fit_results.xml")
     like["expcube"].filename("expcube.fits")
     like["psfcube"].filename("psfcube.fits")
+    like["bkgcube"].filename("bkgcube.fits")
     like["caldb"].string(caldb)
     like["irf"].string(irf)
     like["debug"].boolean(True) # Switch this always on for results in console
@@ -36,7 +36,7 @@ from ctools import obsutils
 def setup_observations(pattern="four", ra=83.63, dec=22.01, offset=1.5, \
                        emin=0.1, emax=100.0, rad=5.0, duration=1800.0, \
                        deadc=0.95, \
-                       caldb="dummy", irf="cta_dummy_irf"):
+                       caldb="prod2", irf="South_50h"):
     """
     Returns an observation container.
@@ -188,7 +188,7 @@ def run_pipeline(obs, ra=83.63, dec=22.01, emin=0.1, emax=100.0, \
     # Set Exposure and Psf cube for first CTA observation
     # (ctbin will create an observation with a single container)
-    bin.obs()[0].response(expcube.expcube(), psfcube.psfcube())
+    bin.obs()[0].response(expcube.expcube(), psfcube.psfcube(), bkgcube.bkgcube())
     # Perform maximum likelihood fitting
     like = ctools.ctlike(bin.obs())
@@ -36,7 +36,7 @@ from ctools import obsutils
 def setup_observations(pattern="four", ra=83.63, dec=22.01, offset=1.5, \
                        emin=0.1, emax=100.0, rad=5.0, duration=1800.0, \
                        deadc=0.95, \
-                       caldb="dummy", irf="cta_dummy_irf"):
+                       caldb="prod2", irf="South_50h"):
     """
     Returns an observation container.
@@ -96,7 +96,7 @@ def run_pipeline(obs, ra=83.63, dec=22.01, rad=3.0, \
                  emin=0.1, emax=100.0, \
                  tmin=0.0, tmax=0.0, \
                  model="${CTOOLS}/share/models/crab.xml", \
-                 caldb="dummy", irf="cta_dummy_irf", \
+                 caldb="prod2", irf="South_50h", \
                  debug=False):
     """
     Simulation and unbinned analysis pipeline.
@@ -36,7 +36,7 @@ from ctools import obsutils
 def setup_observations(pattern="four", ra=83.63, dec=22.01, offset=1.5, \
                        emin=0.1, emax=100.0, rad=5.0, duration=1800.0, \
                        deadc=0.95, \
-                       caldb="dummy", irf="cta_dummy_irf"):
+                       caldb="prod2", irf="South_50h"):
     """
     Returns an observation container.