Update 1DC HowTo

Jürgen Knödlseder
1 parent 67fee531
Showing 4 changed files with 82 additions and 19 deletions
doc/source/users/tutorials/1dc/data_organisation.rst
doc/source/users/tutorials/1dc/howto.rst
doc/source/users/tutorials/1dc/howto_lightcurve.rst
doc/source/users/tutorials/1dc/howto_phasecurve.rst
@@ -224,21 +224,21 @@ An example of an ``EVENTS`` table is shown below.
 Event files exist for
- +-----------------+---------------+--------------------+
- | Survey          | Configuration | Prefix             |
- +=================+===============+====================+
- | Galactic Centre | Baseline      | ``gc_baseline``    |
- +-----------------+---------------+--------------------+
- | Galactic Plane  | Baseline      | ``gps_baseline``   |
- +-----------------+---------------+--------------------+
- | Extragalactic   | Baseline      | ``egal_baseline``  |
- +-----------------+---------------+--------------------+
- | Galactic Centre | Threshold     | ``gc_threshold``   |
- +-----------------+---------------+--------------------+
- | Galactic Plane  | Threshold     | ``gps_threshold``  |
- +-----------------+---------------+--------------------+
- | Extragalactic   | Threshold     | ``egal_threshold`` |
- +-----------------+---------------+--------------------+
+ +-----------------+---------------+--------------------+------------------+
+ | Survey          | Configuration | Prefix             | Energy range     |
+ +=================+===============+====================+==================+
+ | Galactic Centre | Baseline      | ``gc_baseline``    |                  |
+ +-----------------+---------------+--------------------+------------------+
+ | Galactic Plane  | Baseline      | ``gps_baseline``   | 30 GeV - 120 TeV |
+ +-----------------+---------------+--------------------+------------------+
+ | Extragalactic   | Baseline      | ``egal_baseline``  |                  |
+ +-----------------+---------------+--------------------+------------------+
+ | Galactic Centre | Threshold     | ``gc_threshold``   |                  |
+ +-----------------+---------------+--------------------+------------------+
+ | Galactic Plane  | Threshold     | ``gps_threshold``  |                  |
+ +-----------------+---------------+--------------------+------------------+
+ | Extragalactic   | Threshold     | ``egal_threshold`` |                  |
+ +-----------------+---------------+--------------------+------------------+
 Observation Definition Files
@@ -3,9 +3,8 @@
 Analysis cookbook
 =================
-The following sections provide some working examples for the analysis
-of data from the
-:ref:`first CTA Data Challenge <glossary_1dc>`.
+The following sections explain how you can address some specific analysis
+issues.
 .. toctree::
    :maxdepth: 2
@@ -14,4 +13,18 @@ of data from the
    howto_tsmap
    howto_ulimit
    howto_lightcurve
+   howto_phasecurve
    howto_exclude
+
+If you cannot find the answer to your question, have a look at the
+:ref:`first steps with ctools <quickstart>`
+and the
+:ref:`beyond the first steps <beyond>`
+sections of the User Manual. You may also consult the
+:ref:`Reference Manual <reference>`
+as well as the
+:ref:`frequently asked questions <faq>`.
+If after this you still need some help, please create a
+`support ticket <https://cta-redmine.irap.omp.eu/projects/ctools/issues/new>`_.
+Specify ``Support`` as Tracker, choose a short and meaningful subject,
+and add an exhaustive description of your request.
-.. _1dc_howto_ligutcurve:
+.. _1dc_howto_ligthcurve:
 How to generate a light curve?
 ------------------------------
+.. _1dc_howto_phasecurve:
+
+How to generate a phase curve of a pulsar or a gamma-ray binary?
+----------------------------------------------------------------
+
+So far there are no tools in ctools to generate phase curves of pulsars or
+gamma-ray binaries. A plugin for
+`tempo2 <http://www.atnf.csiro.au/research/pulsar/tempo2/>`_
+is under development by the ctools developers, and you may get in touch
+with Thierry Reposeur (reposeur@in2p3.fr) for information about this
+development.
+
+If you know the phase information of the pulsar or gamma-ray binary,
+you may specify a ``temporal`` model component of type ``PhaseCurve`` in the
+:ref:`model definition XML file <glossary_moddef>`
+for the source and fit the information using an unbinned analysis.
+An example of a
+:ref:`model definition XML file <glossary_moddef>`
+including a ``PhaseCurve`` model component is shown below.
+For more information, please read the following
+:ref:`description of the temporal model components <sec_temporal_models>`.
+
+.. code-block:: xml
+
+   <?xml version="1.0" standalone="no"?>
+   <source_library title="source library">
+     <source name="Crab" type="PointSource">
+       <spatialModel type="PointSource">
+         <parameter name="RA"  scale="1" value="83.6331" min="-360" max="360" free="0"/>
+         <parameter name="DEC" scale="1" value="22.0145" min="-90"  max="90"  free="0"/>
+       </spatialModel>
+       <spectrum type="PowerLaw">
+         <parameter name="Prefactor"   scale="1e-16" value="5.7"  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="PivotEnergy" scale="1e6"   value="0.3"  min="0.01"  max="1000.0" free="0"/>
+       </spectrum>
+       <temporal type="PhaseCurve" file="model_temporal_phasecurve.fits">
+         <parameter name="Normalization" scale="1" value="1.0"     min="0.0" max="1000.0"   free="0"/>
+         <parameter name="MJD"           scale="1" value="51544.5" min="0.0" max="100000.0" free="0"/>
+         <parameter name="Phase"         scale="1" value="0.0"     min="0.0" max="1.0"      free="0"/>
+         <parameter name="F0"            scale="1" value="1.0"     min="0.0" max="1000.0"   free="0"/>
+         <parameter name="F1"            scale="1" value="0.1"     min="0.0" max="1000.0"   free="0"/>
+         <parameter name="F2"            scale="1" value="0.01"    min="0.0" max="1000.0"   free="0"/>
+       </temporal>
+     </source>
+   </source_library>
+
+
+
+