Luigi Tibaldo / ctools

Commit fd51695778e80debfeb79e6b2fef36f670399ec3

Authored by Jürgen Knödlseder
1 parent a7ee33f5

Update 1DC documentation

Inline Side-by-side
Showing 3 changed files with 50 additions and 49 deletions
doc/source/users/tutorials/howto/howto_extent.rst
  View file @fd51695
... ... @@ -14,7 +14,7 @@ has to be used in the
14 14 for the source of interest.
15 15 The example below is based on the source model that was determined
16 16 during the
17   -:ref:`source detection step<1dc_select_models>`
  17 +:ref:`iterative model improvement <1dc_first_improving>`
18 18 of the
19 19 :ref:`first CTA Data Challenge <glossary_1dc>`
20 20 tutorial where two point sources were detected.
... ... @@ -52,30 +52,30 @@ will then fit the disk radius of ``Src003``. The results can be seen in the
52 52 log file and will also be written into the output
53 53 :ref:`model definition file <glossary_moddef>`.
54 54 For illustration, an excerpt of the :ref:`ctlike` log file is shown below.
55   -The disk with of ``Src003`` has been fitted to 0.334 +/- 0.003 degrees.
56   -The source is detected with a TS value of 5376.9.
57   -
58   -.. code-block:: xml
59   -
60   - 2017-07-28T05:53:53: === GModelSky ===
61   - 2017-07-28T05:53:53: Name ......................: Src003
62   - 2017-07-28T05:53:53: Instruments ...............: all
63   - 2017-07-28T05:53:53: Test Statistic ............: 5376.92376418598
64   - 2017-07-28T05:53:53: Instrument scale factors ..: unity
65   - 2017-07-28T05:53:53: Observation identifiers ...: all
66   - 2017-07-28T05:53:53: Model type ................: ExtendedSource
67   - 2017-07-28T05:53:53: Model components ..........: "RadialDisk" * "PowerLaw" * "Constant"
68   - 2017-07-28T05:53:53: Number of parameters ......: 7
69   - 2017-07-28T05:53:53: Number of spatial par's ...: 3
70   - 2017-07-28T05:53:53: RA .......................: 266.309961354122 +/- 0.00363413551819907 [-360,360] deg (free,scale=1)
71   - 2017-07-28T05:53:53: DEC ......................: -30.1920697862131 +/- 0.00315430125250908 [-90,90] deg (free,scale=1)
72   - 2017-07-28T05:53:53: Radius ...................: 0.334199097561604 +/- 0.00257183907696455 [0.001,10] deg (free,scale=1)
73   - 2017-07-28T05:53:53: Number of spectral par's ..: 3
74   - 2017-07-28T05:53:53: Prefactor ................: 6.21523600805409e-17 +/- 1.05933757692423e-18 [0,5.7e-15] ph/cm2/s/MeV (free,scale=5.7e-18,gradient)
75   - 2017-07-28T05:53:53: Index ....................: -2.73943047192846 +/- 0.0151745151814222 [9.92,-9.92] (free,scale=-2.48,gradient)
76   - 2017-07-28T05:53:53: PivotEnergy ..............: 300000 MeV (fixed,scale=300000,gradient)
77   - 2017-07-28T05:53:53: Number of temporal par's ..: 1
78   - 2017-07-28T05:53:53: Normalization ............: 1 (relative value) (fixed,scale=1,gradient)
  55 +The disk with of ``Src003`` has been fitted to 0.335 +/- 0.003 degrees.
  56 +The source is detected with a TS value of 5393.8.
  57 +
  58 +.. code-block:: bash
  59 +
  60 + 2017-08-25T19:37:47: === GModelSky ===
  61 + 2017-08-25T19:37:47: Name ......................: Src003
  62 + 2017-08-25T19:37:47: Instruments ...............: all
  63 + 2017-08-25T19:37:47: Test Statistic ............: 5393.78645148844
  64 + 2017-08-25T19:37:47: Instrument scale factors ..: unity
  65 + 2017-08-25T19:37:47: Observation identifiers ...: all
  66 + 2017-08-25T19:37:47: Model type ................: ExtendedSource
  67 + 2017-08-25T19:37:47: Model components ..........: "RadialDisk" * "PowerLaw" * "Constant"
  68 + 2017-08-25T19:37:47: Number of parameters ......: 7
  69 + 2017-08-25T19:37:47: Number of spatial par's ...: 3
  70 + 2017-08-25T19:37:47: RA .......................: 266.306975228126 +/- 0.00363322258453142 [-360,360] deg (free,scale=1)
  71 + 2017-08-25T19:37:47: DEC ......................: -30.1946569442618 +/- 0.0031545381615373 [-90,90] deg (free,scale=1)
  72 + 2017-08-25T19:37:47: Radius ...................: 0.335173109280516 +/- 0.00257226134023012 [0.001,10] deg (free,scale=1)
  73 + 2017-08-25T19:37:47: Number of spectral par's ..: 3
  74 + 2017-08-25T19:37:47: Prefactor ................: 6.23930379041412e-17 +/- 1.06204679984262e-18 [0,5.7e-15] ph/cm2/s/MeV (free,scale=5.7e-18,gradient)
  75 + 2017-08-25T19:37:47: Index ....................: -2.73810250489576 +/- 0.0151068737641297 [9.92,-9.92] (free,scale=-2.48,gradient)
  76 + 2017-08-25T19:37:47: PivotEnergy ..............: 300000 MeV (fixed,scale=300000,gradient)
  77 + 2017-08-25T19:37:47: Number of temporal par's ..: 1
  78 + 2017-08-25T19:37:47: Normalization ............: 1 (relative value) (fixed,scale=1,gradient)
79 79  
80 80 .. warning::
81 81 The parameter errors returned by :ref:`ctlike` are purely statistical. In
... ... @@ -83,7 +83,8 @@ The source is detected with a TS value of 5376.9.
83 83 as for example the limits on the knowledge of the point spread function.
84 84 These systematic uncertainties are **not** determined by ctools.
85 85  
86   -The figure below shows the residual map after subtracting the three fitted
  86 +The figure below shows the residual map, generated using the ``SIGNIFICANCE``
  87 +method, after subtracting the three fitted
87 88 sources and the diffuse emission model from the data. There is a ring-like
88 89 residual at the position of ``Src003`` which suggests that an axisymmetric
89 90 disk is not an accurate description of the data.
... ... @@ -116,36 +117,36 @@ looks as follows:
116 117  
117 118 Running :ref:`ctlike` again with that model results in a Gaussian sigma of
118 119 0.198 +/- 0.003 degrees for ``Src003``. The source is detected with a TS value
119   -of 5659.2 which is considerably larger than the TS value of 5376.9 that is found
  120 +of 5665.7 which is considerably larger than the TS value of 5393.8 that is found
120 121 above for the radial disk model.
121 122 Below an excerpt of the :ref:`ctlike` log file:
122 123  
123   -.. code-block:: xml
124   -
125   - 2017-07-28T05:37:55: === GModelSky ===
126   - 2017-07-28T05:37:55: Name ......................: Src003
127   - 2017-07-28T05:37:55: Instruments ...............: all
128   - 2017-07-28T05:37:55: Test Statistic ............: 5659.21767961094
129   - 2017-07-28T05:37:55: Instrument scale factors ..: unity
130   - 2017-07-28T05:37:55: Observation identifiers ...: all
131   - 2017-07-28T05:37:55: Model type ................: ExtendedSource
132   - 2017-07-28T05:37:55: Model components ..........: "RadialGaussian" * "PowerLaw" * "Constant"
133   - 2017-07-28T05:37:55: Number of parameters ......: 7
134   - 2017-07-28T05:37:55: Number of spatial par's ...: 3
135   - 2017-07-28T05:37:55: RA .......................: 266.302363126862 +/- 0.00476126847653064 [-360,360] deg (free,scale=1)
136   - 2017-07-28T05:37:55: DEC ......................: -30.1977095245935 +/- 0.00412840925183781 [-90,90] deg (free,scale=1)
137   - 2017-07-28T05:37:55: Sigma ....................: 0.198403793257533 +/- 0.00298520334678877 [0.001,10] deg (free,scale=1)
138   - 2017-07-28T05:37:55: Number of spectral par's ..: 3
139   - 2017-07-28T05:37:55: Prefactor ................: 7.15243704872173e-17 +/- 1.40858695554198e-18 [0,5.7e-15] ph/cm2/s/MeV (free,scale=5.7e-18,gradient)
140   - 2017-07-28T05:37:55: Index ....................: -2.71991083144987 +/- 0.0147325731477118 [9.92,-9.92] (free,scale=-2.48,gradient)
141   - 2017-07-28T05:37:55: PivotEnergy ..............: 300000 MeV (fixed,scale=300000,gradient)
142   - 2017-07-28T05:37:55: Number of temporal par's ..: 1
143   - 2017-07-28T05:37:55: Normalization ............: 1 (relative value) (fixed,scale=1,gradient)
  124 +.. code-block:: bash
  125 +
  126 + 2017-08-25T20:32:55: === GModelSky ===
  127 + 2017-08-25T20:32:55: Name ......................: Src003
  128 + 2017-08-25T20:32:55: Instruments ...............: all
  129 + 2017-08-25T20:32:55: Test Statistic ............: 5665.73173788155
  130 + 2017-08-25T20:32:55: Instrument scale factors ..: unity
  131 + 2017-08-25T20:32:55: Observation identifiers ...: all
  132 + 2017-08-25T20:32:55: Model type ................: ExtendedSource
  133 + 2017-08-25T20:32:55: Model components ..........: "RadialGaussian" * "PowerLaw" * "Constant"
  134 + 2017-08-25T20:32:55: Number of parameters ......: 7
  135 + 2017-08-25T20:32:55: Number of spatial par's ...: 3
  136 + 2017-08-25T20:32:55: RA .......................: 266.30022317901 +/- 0.00475655134090925 [-360,360] deg (free,scale=1)
  137 + 2017-08-25T20:32:55: DEC ......................: -30.1993627725406 +/- 0.0041251924859714 [-90,90] deg (free,scale=1)
  138 + 2017-08-25T20:32:55: Sigma ....................: 0.198343910095786 +/- 0.00298438230006778 [0.001,10] deg (free,scale=1)
  139 + 2017-08-25T20:32:55: Number of spectral par's ..: 3
  140 + 2017-08-25T20:32:55: Prefactor ................: 7.15427332759485e-17 +/- 1.40842792023593e-18 [0,5.7e-15] ph/cm2/s/MeV (free,scale=5.7e-18,gradient)
  141 + 2017-08-25T20:32:55: Index ....................: -2.71899753120954 +/- 0.0146676775081257 [9.92,-9.92] (free,scale=-2.48,gradient)
  142 + 2017-08-25T20:32:55: PivotEnergy ..............: 300000 MeV (fixed,scale=300000,gradient)
  143 + 2017-08-25T20:32:55: Number of temporal par's ..: 1
  144 + 2017-08-25T20:32:55: Normalization ............: 1 (relative value) (fixed,scale=1,gradient)
144 145  
145 146 .. note::
146 147 While the TS values can formally not be convert into a statistical
147 148 significance between different spatial model hypotheses, a TS improvement
148   - of 282.3 indicates a considerably better fit of the axisymmetric Gaussian
  149 + of 271.9 indicates a considerably better fit of the axisymmetric Gaussian
149 150 model with respect to the radial disk model to the data.
150 151  
151 152 The figure below shows the residual map for the fit of ``Src003`` with an
... ...
doc/source/users/tutorials/howto/howto_extent_disk.png
View file @fd51695

317 KB | W: | H:

289 KB | W: | H:

  • 2-up
  • Swipe
  • Onion skin
doc/source/users/tutorials/howto/howto_extent_gauss.png
View file @fd51695

321 KB | W: | H:

287 KB | W: | H:

  • 2-up
  • Swipe
  • Onion skin