More and more docs

Michael Mayer
1 parent d92c32be
Showing 3 changed files with 85 additions and 0 deletions
doc/source/user_manual/iact_analysis/analysis.rst
doc/source/user_manual/iact_analysis/high_level.rst
doc/source/user_manual/iact_analysis/tipps.rst → doc/source/user_manual/iact_analysis/tips.rst
@@ -385,6 +385,37 @@ sky cube. In the output model the background models per obseervation will be rem
 for the newly created background cube is included. Sky models present in the input model XML file will also be included in the 
 new XML file, which subsequently can be used for binned :ref:`ctlike`.
+Example for stacked model XML file
+^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
+The output model of :ref:`ctbkgcube` looks the following:
+
+.. code-block:: xml
+
+	<?xml version="1.0" encoding="UTF-8" standalone="no"?>
+	<source_library title="source library">
+	  <source name="BackgroundModel" type="CTACubeBackground" instrument="CTA,HESS,MAGIC,VERITAS">
+	    <spectrum type="PowerLaw">
+	      <parameter name="Prefactor" value="1" error="0" scale="1" min="0.01" max="100" free="1" />
+	      <parameter name="Index" value="0" error="0" scale="1" min="-5" max="5" free="1" />
+	      <parameter name="Scale" value="1" scale="1e+06" free="0" />
+	    </spectrum>
+	  </source>
+	  <source name="Crab" type="PointSource">
+	    <spectrum type="PowerLaw">
+	       <parameter name="Prefactor" scale="1e-17" value="3.0"  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="Scale"     scale="1e6"   value="1."  min="0.01"  max="1000.0" free="0"/>
+	    </spectrum>
+	    <spatialModel type="SkyDirFunction">
+	      <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_library>
+	
+The background model with ``type=CTACubeBackground`` is used to scale the background cube stored in the FITS file
+created by :ref:`ctbkgcube`.
+
 Run ctlike
 ^^^^^^^^^^
 Having all the intermediate data products ready, a binned analysis can be conducted using :ref:`ctlike`.
@@ -111,13 +111,64 @@ source of interest. For this purpose, the tool :ref:`cterror` can be used:
 	Input event list, counts cube or observation definition XML file [selected_obs.xml] 
 	Input model XML file [crab_models.xml] 
 	Source of interest [Crab] 
+	Output model XML file [cterror_results.xml] 
+The output model does not contain asymmetric errors yet. The positive and negative uncertainties can be read from the
+logfile (or from screen) if ``debug=yes`` was specified.
 Compute spectral points
 -----------------------
+A very common task in astronomy is to compute spectral data points. To determine a spectral data point, a small energy
+range is considered and the model prefactor and its uncertainty is evaluated. The tool :ref:`csspec` works on both, binned
+and unbinned data. 
+
+Unbinned
+^^^^^^^^
+
+.. code-block:: bash
+
+	$ csspec debug=yes
+	Input event list, counts cube, or observation definition XML file [events.fits] selected_obs.xml 
+	Input model XML file [$CTOOLS/share/models/crab.xml] crab_models.xml 
+	Source name [Crab] 
+	Algorithm for defining energy bins (FILE|LIN|LOG) [LOG] 
+	Lower energy limit for spectral points (TeV) [0.1] 0.5
+	Upper energy limit for spectral points (TeV) [100.0] 50.0
+	Number of spectral points (1-10000) [20] 10
+	Output spectrum file [spectrum.fits] 
+	
+Binned
+^^^^^^
+
+.. code-block:: bash
+
+	$ csspec debug=yes
+	Input event list, counts cube, or observation definition XML file [cntcube.fits] 
+	Input exposure cube file (only needed for stacked analysis) [expcube.fits] 
+	Input PSF cube file (only needed for stacked analysis) [psfcube.fits] 
+	Input background cube file (only needed for stacked analysis) [bkgcube.fits] 
+	Input model XML file [binned_models.xml] 
+	Source name [Crab] 
+	Number of spectral points (1-10000) [10] 
+	Lower energy limit for spectral points (TeV) [0.5] 
+	Upper energy limit for spectral points (TeV) [50.0] 
+	Output spectrum file [spectrum.fits] 
+
+
+Plot spectral points
+^^^^^^^^^^^^^^^^^^^^
+Instead of supporting a plotting library, simple example scripts are available to visualise data products.
+Have a look at ``$CTOOLS/examples/show_spectrum.py`` how a spectrum can be plotted using ``matplotlib``.
+To have first glance at the above computed spectrum one can use this script in the following way:
+
+.. code-block:: bash
+  
+  $ python $CTOOLS/examples/show_spectrum.py spectrum.fits
+
 Compute light curves
 --------------------
+Euqally common in astronomy are light curves from a time-series analysis.
 Compute a residual map
 ----------------------
@@ -9,6 +9,9 @@ Visualise observations
 Visualise models
 ----------------
+Increase analysis speed
+-----------------------
+
 Create count map with one energy bin
 ------------------------------------