Andreas Specovius / ctools

Commit 7bad2878e14856288c80b6db8788b4a8672af1b1

Authored by Jürgen Knödlseder
1 parent ba0f65c1

Update first Data Challenge tutorial

Inline Side-by-side
Showing 21 changed files with 216 additions and 117 deletions
doc/source/users/tutorials/1dc/first_improving.rst
  View file @7bad287
... ... @@ -94,9 +94,11 @@ are flatter and the diffuse emission has largely disappeared. The diffuse
94 94 emission is in fact a bit oversubtracted, which is probably due to the
95 95 additional source components visible in the residual map. It is left as an
96 96 exercise to the User to include model components for these additional
97   -components in the
  97 +sources in the
98 98 :ref:`model definition file <glossary_moddef>`
99 99 in order to reduce the fit residuals.
  100 +An example for such an exercise is shown
  101 +:ref:`here <1dc_howto_extent>`.
100 102  
101 103 .. figure:: first_skymap_residual_iem.png
102 104 :width: 400px
... ...
doc/source/users/tutorials/1dc/first_spectrum_cutoff_unbinned.png
View file @7bad287

55.4 KB | W: | H:

55.2 KB | W: | H:

  • 2-up
  • Swipe
  • Onion skin
doc/source/users/tutorials/1dc/first_unbinned.rst
  View file @7bad287
... ... @@ -11,7 +11,7 @@ Fitting the model components directly to the events
11 11 Use an unbinned maximum likelihood analysis if you have doubts about the
12 12 impact of the selected binning on your analysis and if the number of
13 13 observations is not too large (typically a few tens of 30 minute
14   - observations). Or if you simply to not worry about computation time.
  14 + observations). Or if you simply do not worry about computation time.
15 15  
16 16 Instead of binning the events into a counts cube and computing the effective
17 17 :ref:`instrument response functions <glossary_irf>`
... ... @@ -47,47 +47,35 @@ converged properly, as illustrated in the example below:
47 47  
48 48 .. code-block:: bash
49 49  
50   - 2017-06-08T09:56:42: +=================================+
51   - 2017-06-08T09:56:42: | Maximum likelihood optimisation |
52   - 2017-06-08T09:56:42: +=================================+
53   - 2017-06-08T09:57:14: >Iteration 0: -logL=22938810.009, Lambda=1.0e-03
54   - 2017-06-08T09:57:43: >Iteration 1: -logL=22914874.739, Lambda=1.0e-03, delta=23935.270, step=1.0e+00, max(|grad|)=52907.253765 [Index:14]
55   - 2017-06-08T09:58:11: >Iteration 2: -logL=22911216.020, Lambda=1.0e-04, delta=3658.719, step=1.0e+00, max(|grad|)=-13199.549393 [DEC:1]
56   - 2017-06-08T09:58:41: >Iteration 3: -logL=22909769.414, Lambda=1.0e-05, delta=1446.606, step=1.0e+00, max(|grad|)=-11987.289647 [DEC:1]
57   - 2017-06-08T09:59:10: >Iteration 4: -logL=22909405.587, Lambda=1.0e-06, delta=363.826, step=1.0e+00, max(|grad|)=-8985.645939 [DEC:1]
58   - 2017-06-08T09:59:39: >Iteration 5: -logL=22909337.513, Lambda=1.0e-07, delta=68.074, step=1.0e+00, max(|grad|)=-5891.649861 [DEC:1]
59   - 2017-06-08T10:00:07: >Iteration 6: -logL=22909309.777, Lambda=1.0e-08, delta=27.736, step=1.0e+00, max(|grad|)=-3762.571320 [DEC:1]
60   - 2017-06-08T10:00:36: >Iteration 7: -logL=22909295.909, Lambda=1.0e-09, delta=13.868, step=1.0e+00, max(|grad|)=-2458.462512 [RA:7]
61   - 2017-06-08T10:01:05: >Iteration 8: -logL=22909288.709, Lambda=1.0e-10, delta=7.200, step=1.0e+00, max(|grad|)=-1835.957201 [RA:7]
62   - 2017-06-08T10:01:34: >Iteration 9: -logL=22909284.904, Lambda=1.0e-11, delta=3.805, step=1.0e+00, max(|grad|)=-1368.259068 [RA:7]
63   - 2017-06-08T10:02:04: >Iteration 10: -logL=22909282.863, Lambda=1.0e-12, delta=2.041, step=1.0e+00, max(|grad|)=-1022.653174 [RA:7]
64   - 2017-06-08T10:02:35: >Iteration 11: -logL=22909281.751, Lambda=1.0e-13, delta=1.112, step=1.0e+00, max(|grad|)=-765.437910 [RA:7]
65   - 2017-06-08T10:03:05: >Iteration 12: -logL=22909281.140, Lambda=1.0e-14, delta=0.611, step=1.0e+00, max(|grad|)=-572.313029 [RA:7]
66   - 2017-06-08T10:03:33: >Iteration 13: -logL=22909280.802, Lambda=1.0e-15, delta=0.338, step=1.0e+00, max(|grad|)=-430.018412 [RA:7]
67   - 2017-06-08T10:04:02: >Iteration 14: -logL=22909280.613, Lambda=1.0e-16, delta=0.189, step=1.0e+00, max(|grad|)=-321.908506 [RA:7]
68   - 2017-06-08T10:04:30: >Iteration 15: -logL=22909280.507, Lambda=1.0e-17, delta=0.106, step=1.0e+00, max(|grad|)=-240.266994 [RA:7]
69   - 2017-06-08T10:04:58: >Iteration 16: -logL=22909280.448, Lambda=1.0e-18, delta=0.059, step=1.0e+00, max(|grad|)=-180.004153 [RA:7]
70   - 2017-06-08T10:05:27: >Iteration 17: -logL=22909280.414, Lambda=1.0e-19, delta=0.033, step=1.0e+00, max(|grad|)=-135.633545 [RA:7]
71   - 2017-06-08T10:05:55: >Iteration 18: -logL=22909280.395, Lambda=1.0e-20, delta=0.019, step=1.0e+00, max(|grad|)=-102.119554 [RA:7]
72   - 2017-06-08T10:06:24: >Iteration 19: -logL=22909280.385, Lambda=1.0e-21, delta=0.011, step=1.0e+00, max(|grad|)=-76.657216 [RA:7]
73   - 2017-06-08T10:06:52: >Iteration 20: -logL=22909280.379, Lambda=1.0e-22, delta=0.006, step=1.0e+00, max(|grad|)=-57.439799 [RA:7]
74   - 2017-06-08T10:07:20: >Iteration 21: -logL=22909280.375, Lambda=1.0e-23, delta=0.003, step=1.0e+00, max(|grad|)=-43.064181 [RA:7]
75   - 2017-06-08T10:07:49:
76   - 2017-06-08T10:07:49: +=========================================+
77   - 2017-06-08T10:07:49: | Maximum likelihood optimisation results |
78   - 2017-06-08T10:07:49: +=========================================+
79   - 2017-06-08T10:07:49: === GOptimizerLM ===
80   - 2017-06-08T10:07:49: Optimized function value ..: 22909280.375
81   - 2017-06-08T10:07:49: Absolute precision ........: 0.005
82   - 2017-06-08T10:07:49: Acceptable value decrease .: 2
83   - 2017-06-08T10:07:49: Optimization status .......: converged
84   - 2017-06-08T10:07:49: Number of parameters ......: 17
85   - 2017-06-08T10:07:49: Number of free parameters .: 11
86   - 2017-06-08T10:07:49: Number of iterations ......: 21
87   - 2017-06-08T10:07:49: Lambda ....................: 1e-24
88   - 2017-06-08T10:07:49: Maximum log likelihood ....: -22909280.375
89   - 2017-06-08T10:07:49: Observed events (Nobs) ...: 3349183.000
90   - 2017-06-08T10:07:49: Predicted events (Npred) ..: 3349177.996 (Nobs - Npred = 5.00393125554547)
  50 + 2017-07-27T16:28:26: +=================================+
  51 + 2017-07-27T16:28:26: | Maximum likelihood optimisation |
  52 + 2017-07-27T16:28:26: +=================================+
  53 + 2017-07-27T16:28:54: >Iteration 0: -logL=20949416.217, Lambda=1.0e-03
  54 + 2017-07-27T16:29:21: >Iteration 1: -logL=20918674.227, Lambda=1.0e-03, delta=30741.990, step=1.0e+00, max(|grad|)=74249.049334 [Index:14]
  55 + 2017-07-27T16:29:49: >Iteration 2: -logL=20913114.563, Lambda=1.0e-04, delta=5559.664, step=1.0e+00, max(|grad|)=6728.499004 [RA:0]
  56 + 2017-07-27T16:30:16: >Iteration 3: -logL=20910994.152, Lambda=1.0e-05, delta=2120.411, step=1.0e+00, max(|grad|)=4949.846735 [RA:0]
  57 + 2017-07-27T16:30:44: >Iteration 4: -logL=20910115.048, Lambda=1.0e-06, delta=879.104, step=1.0e+00, max(|grad|)=2889.871293 [RA:0]
  58 + 2017-07-27T16:31:12: >Iteration 5: -logL=20909966.036, Lambda=1.0e-07, delta=149.012, step=1.0e+00, max(|grad|)=1190.543059 [RA:0]
  59 + 2017-07-27T16:31:39: >Iteration 6: -logL=20909948.012, Lambda=1.0e-08, delta=18.024, step=1.0e+00, max(|grad|)=550.745751 [DEC:1]
  60 + 2017-07-27T16:32:06: >Iteration 7: -logL=20909946.600, Lambda=1.0e-09, delta=1.412, step=1.0e+00, max(|grad|)=208.667347 [DEC:1]
  61 + 2017-07-27T16:32:33: >Iteration 8: -logL=20909946.575, Lambda=1.0e-10, delta=0.025, step=1.0e+00, max(|grad|)=51.069713 [DEC:1]
  62 + 2017-07-27T16:33:00: >Iteration 9: -logL=20909946.573, Lambda=1.0e-11, delta=0.002, step=1.0e+00, max(|grad|)=16.914402 [DEC:1]
  63 + 2017-07-27T16:33:26:
  64 + 2017-07-27T16:33:26: +=========================================+
  65 + 2017-07-27T16:33:26: | Maximum likelihood optimisation results |
  66 + 2017-07-27T16:33:26: +=========================================+
  67 + 2017-07-27T16:33:26: === GOptimizerLM ===
  68 + 2017-07-27T16:33:26: Optimized function value ..: 20909946.573
  69 + 2017-07-27T16:33:26: Absolute precision ........: 0.005
  70 + 2017-07-27T16:33:26: Acceptable value decrease .: 2
  71 + 2017-07-27T16:33:26: Optimization status .......: converged
  72 + 2017-07-27T16:33:26: Number of parameters ......: 17
  73 + 2017-07-27T16:33:26: Number of free parameters .: 11
  74 + 2017-07-27T16:33:26: Number of iterations ......: 9
  75 + 2017-07-27T16:33:26: Lambda ....................: 1e-12
  76 + 2017-07-27T16:33:26: Maximum log likelihood ....: -20909946.573
  77 + 2017-07-27T16:33:26: Observed events (Nobs) ...: 3084639.000
  78 + 2017-07-27T16:33:26: Predicted events (Npred) ..: 3084635.972 (Nobs - Npred = 3.0282306545414)
91 79  
92 80 Also the :ref:`ctbutterfly` tool and the :ref:`csspec` and :ref:`csresmap`
93 81 scripts can be run into unbinned mode by providing the
... ...
doc/source/users/tutorials/howto/howto_display_butterfly.png
View file @7bad287

37.8 KB | W: | H:

33.2 KB | W: | H:

  • 2-up
  • Swipe
  • Onion skin
doc/source/users/tutorials/howto/howto_display_spectrum.png
View file @7bad287

22.4 KB | W: | H:

22.5 KB | W: | H:

  • 2-up
  • Swipe
  • Onion skin
doc/source/users/tutorials/howto/howto_extent.rst
  View file @7bad287
... ... @@ -11,63 +11,172 @@ How to determine the extension of a source?
11 11 To determine the spatial extension of a source, an extended source model
12 12 has to be used in the
13 13 :ref:`model definition file <glossary_moddef>`
14   -of the source of interest.
15   -In the example below, a ``RadialGaussian`` model is used for the source
16   -``Src001`` to model the spatial morphology of the source by an axisymmetric
17   -Gaussian profile with the Gaussian width ``sigma`` as free parameter.
18   -Note that at the same time the source type was set to
19   -``ExtendedSource``.
  14 +for the source of interest.
  15 +The example below is based on the source model that was determined
  16 +during the
  17 +:ref:`source detection step<1dc_select_models>`
  18 +of the
  19 +:ref:`first CTA Data Challenge <glossary_1dc>`
  20 +tutorial where two point sources were detected.
  21 +As in the tutorial, ``Src001`` is modelled with an exponentially cut-off power
  22 +law spectrum, ``Src002`` with a simple power law.
  23 +In addition, a diffuse emission component was added to the model.
  24 +
  25 +As noted during the
  26 +:ref:`iterative improvement of the source model <1dc_first_improving>`
  27 +there is a bright extended source south-west of the Galactic Centre, and
  28 +we add now an additional ``Src003`` component with a ``RadialDisk`` model to
  29 +describe the spatial morphology of the source.
  30 +The
  31 +:ref:`model definition <glossary_moddef>`
  32 +for ``Src003`` is shown below:
20 33  
21 34 .. code-block:: xml
22 35  
23   - <source name="Src001" type="PointSource">
24   - <spectrum type="ExponentialCutoffPowerLaw">
25   - <parameter name="Prefactor" scale="1e-18" value="5.7" min="1e-07" max="1000.0" free="1"/>
26   - <parameter name="Index" scale="-1" value="2.48" min="0.0" max="+5.0" free="1"/>
27   - <parameter name="CutoffEnergy" scale="1e6" value="10.0" min="0.01" max="1000.0" free="1"/>
28   - <parameter name="PivotEnergy" scale="1e6" value="0.3" min="0.01" max="1000.0" free="0"/>
  36 + <source name="Src003" type="ExtendedSource" tscalc="1">
  37 + <spectrum type="PowerLaw">
  38 + <parameter name="Prefactor" scale="5.7e-18" value="1.0" min="0" max="1000.0" free="1"/>
  39 + <parameter name="Index" scale="-2.48" value="1.0" min="-4.0" max="4.0" free="1"/>
  40 + <parameter name="PivotEnergy" scale="300000" value="1.0" free="0" />
29 41 </spectrum>
30   - <spatialModel type="RadialGaussian">
31   - <parameter name="RA" scale="1.0" value="266.4121" min="-360" max="360" free="1"/>
32   - <parameter name="DEC" scale="1.0" value="-29.0220" min="-90" max="90" free="1"/>
33   - <parameter name="Sigma" scale="1.0" value="0.01" min="0.001" max="10" free="1"/>
  42 + <spatialModel type="RadialDisk">
  43 + <parameter name="RA" scale="1.0" value="266.3070" min="-360" max="360" free="1"/>
  44 + <parameter name="DEC" scale="1.0" value="-30.1876" min="-90" max="90" free="1"/>
  45 + <parameter name="Radius" scale="1.0" value="0.3" min="0.001" max="10" free="1"/>
34 46 </spatialModel>
35 47 </source>
36 48  
37   -Running :ref:`ctlike` using the modified
  49 +Running :ref:`ctlike` using the
38 50 :ref:`model definition file <glossary_moddef>`
39   -will then fit the Gaussian width of ``Src001``. The results can be seen in the
  51 +will then fit the disk radius of ``Src003``. The results can be seen in the
40 52 log file and will also be written into the output
41 53 :ref:`model definition file <glossary_moddef>`.
42   -
43 54 For illustration, an excerpt of the :ref:`ctlike` log file is shown below.
44   -The Gaussian width of ``Src001`` has been fitted to 0.048 +/- 0.001 degrees.
  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.
45 57  
46 58 .. code-block:: xml
47 59  
48   - 2017-06-08T15:07:46: === GModelSky ===
49   - 2017-06-08T15:07:46: Name ......................: Src001
50   - 2017-06-08T15:07:46: Instruments ...............: all
51   - 2017-06-08T15:07:46: Instrument scale factors ..: unity
52   - 2017-06-08T15:07:46: Observation identifiers ...: all
53   - 2017-06-08T15:07:46: Model type ................: ExtendedSource
54   - 2017-06-08T15:07:46: Model components ..........: "RadialGaussian" * "ExponentialCutoffPowerLaw" * "Constant"
55   - 2017-06-08T15:07:46: Number of parameters ......: 8
56   - 2017-06-08T15:07:46: Number of spatial par's ...: 3
57   - 2017-06-08T15:07:46: RA .......................: 266.417513130306 +/- 0.00132553002208653 [-360,360] deg (free,scale=1)
58   - 2017-06-08T15:07:46: DEC ......................: -29.0076435535415 +/- 0.00117490755596776 [-90,90] deg (free,scale=1)
59   - 2017-06-08T15:07:46: Sigma ....................: 0.0484844490479065 +/- 0.000997962295674141 [0.001,10] deg (free,scale=1)
60   - 2017-06-08T15:07:46: Number of spectral par's ..: 4
61   - 2017-06-08T15:07:46: Prefactor ................: 3.26956513671326e-17 +/- 6.26524553994195e-19 [1e-25,1e-15] ph/cm2/s/MeV (free,scale=1e-18,gradient)
62   - 2017-06-08T15:07:46: Index ....................: -2.32005423294123 +/- 0.0173504097487491 [-0,-5] (free,scale=-1,gradient)
63   - 2017-06-08T15:07:46: CutoffEnergy .............: 162314916.632706 +/- 182782066.654979 [10000,1000000000] MeV (free,scale=1000000,gradient)
64   - 2017-06-08T15:07:46: PivotEnergy ..............: 300000 [10000,1000000000] MeV (fixed,scale=1000000,gradient)
65   - 2017-06-08T15:07:46: Number of temporal par's ..: 1
66   - 2017-06-08T15:07:46: Normalization ............: 1 (relative value) (fixed,scale=1,gradient)
  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)
67 79  
68 80 .. warning::
69 81 The parameter errors returned by :ref:`ctlike` are purely statistical. In
70 82 addition to the statistical errors there are systematic uncertainties, such
71 83 as for example the limits on the knowledge of the point spread function.
  84 + These systematic uncertainties are **not** determined by ctools.
  85 +
  86 +The figure below shows the residual map after subtracting the three fitted
  87 +sources and the diffuse emission model from the data. There is a ring-like
  88 +residual at the position of ``Src003`` which suggests that an axisymmetric
  89 +disk is not an accurate description of the data.
  90 +
  91 +.. figure:: howto_extent_disk.png
  92 + :width: 400px
  93 + :align: center
  94 +
  95 + *Residual sky map for a radial disk spatial shape for Src003*
  96 +
  97 +There are other spatial models in ctools, and we try in a second step an
  98 +axisymmetric Gaussian spatial shape. The corresponding
  99 +:ref:`model definition file <glossary_moddef>`
  100 +looks as follows:
  101 +
  102 +.. code-block:: xml
  103 +
  104 + <source name="Src003" type="ExtendedSource" tscalc="1">
  105 + <spectrum type="PowerLaw">
  106 + <parameter name="Prefactor" scale="5.7e-18" value="1.0" min="0" max="1000.0" free="1"/>
  107 + <parameter name="Index" scale="-2.48" value="1.0" min="-4.0" max="4.0" free="1"/>
  108 + <parameter name="PivotEnergy" scale="300000" value="1.0" free="0" />
  109 + </spectrum>
  110 + <spatialModel type="RadialGaussian">
  111 + <parameter name="RA" scale="1.0" value="266.3070" min="-360" max="360" free="1"/>
  112 + <parameter name="DEC" scale="1.0" value="-30.1876" min="-90" max="90" free="1"/>
  113 + <parameter name="Sigma" scale="1.0" value="0.2" min="0.001" max="10" free="1"/>
  114 + </spatialModel>
  115 + </source>
  116 +
  117 +Running :ref:`ctlike` again with that model results in a Gaussian sigma of
  118 +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 +above for the radial disk model.
  121 +Below an excerpt of the :ref:`ctlike` log file:
  122 +
  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)
  144 +
  145 +.. note::
  146 + While the TS values can formally not be convert into a statistical
  147 + significance between different spatial model hypotheses, a TS improvement
  148 + of 282.3 indicates a considerably better fit of the axisymmetric Gaussian
  149 + model with respect to the radial disk model to the data.
  150 +
  151 +The figure below shows the residual map for the fit of ``Src003`` with an
  152 +axisymmetric Gaussian model.
  153 +The map now looks pretty flat around ``Src003``, suggesting that an axisymmetric
  154 +Gaussian model is an appropriate description for the morphology of the gamma-ray
  155 +source.
  156 +
  157 +.. figure:: howto_extent_gauss.png
  158 + :width: 400px
  159 + :align: center
  160 +
  161 + *Residual sky map for an axisymmetric Gaussian spatial shape for Src003*
  162 +
  163 +.. tip::
  164 + The region overlays for the residual sky maps were generated using the
  165 + :ref:`csmodelinfo` script.
  166 +
  167 +.. warning::
  168 + The fitting of extended spatial models takes more computing time
  169 + than the fitting of point sources. The computing time is related to the
  170 + spatial extent of the source and to the spatial shape, with a Gaussian
  171 + disk model taking considerably more computing time than a radial disk
  172 + model due to the tails of the Gaussian function. It is therefore **recommended
  173 + to use by default radial disk models for the extension fitting**, and only
  174 + switch to a Gaussian disk models when really needed, or for the determination
  175 + of final values for a publication.
72 176  
  177 + For reference, here the computing times on Mac OS X for the example using
  178 + different spatial morphology hypotheses for ``Src003``:
73 179  
  180 + * Point source: 12 min
  181 + * Radial disk source: 15 min
  182 + * Axisymmetric Gaussian source: 40 min
... ...
doc/source/users/tutorials/howto/howto_extent_disk.png 0 → 100644
View file @7bad287

317 KB
doc/source/users/tutorials/howto/howto_extent_gauss.png 0 → 100644
View file @7bad287

321 KB
doc/source/users/tutorials/howto/howto_filtered.png
View file @7bad287

59.8 KB | W: | H:

60 KB | W: | H:

  • 2-up
  • Swipe
  • Onion skin
doc/source/users/tutorials/howto/howto_filtered_exclude.png
View file @7bad287

56.5 KB | W: | H:

53.1 KB | W: | H:

  • 2-up
  • Swipe
  • Onion skin
doc/source/users/tutorials/howto/howto_lightcurve.png
View file @7bad287

23.9 KB | W: | H:

24.2 KB | W: | H:

  • 2-up
  • Swipe
  • Onion skin
doc/source/users/tutorials/howto/howto_phasecurve.rst
  View file @7bad287
... ... @@ -195,8 +195,9 @@ probability threshold to &gt;20% by adding the attribute
195 195 Upper energy limit (TeV) [160.0]
196 196 Output event list or observation definition XML file [obs_vela_phased_selected.xml] obs_vela_phased_prob_selected.xml
197 197  
198   -The selected events as a function of phase are displayed below. Unfortunately
199   -there is still no clear signal for a phase dependence of the event rate.
  198 +The selected events as a function of phase are displayed below. There is now a
  199 +hint for a phase dependence of the event rate, although this hint is not very
  200 +significant.
200 201  
201 202 .. figure:: howto_phasecurve_probphases.png
202 203 :width: 600px
... ...
doc/source/users/tutorials/howto/howto_phasecurve_events.png
View file @7bad287

129 KB | W: | H:

135 KB | W: | H:

  • 2-up
  • Swipe
  • Onion skin
doc/source/users/tutorials/howto/howto_phasecurve_phases.png
View file @7bad287

16 KB | W: | H:

16 KB | W: | H:

  • 2-up
  • Swipe
  • Onion skin
doc/source/users/tutorials/howto/howto_phasecurve_prob.png
View file @7bad287

153 KB | W: | H:

158 KB | W: | H:

  • 2-up
  • Swipe
  • Onion skin
doc/source/users/tutorials/howto/howto_phasecurve_probphases.png
View file @7bad287

13.3 KB | W: | H:

12.8 KB | W: | H:

  • 2-up
  • Swipe
  • Onion skin
doc/source/users/tutorials/howto/howto_ts.rst
  View file @7bad287
... ... @@ -23,13 +23,12 @@ The structure of the XML file with significance computation requested for
23 23 <spectrum type="ExponentialCutoffPowerLaw">
24 24 <parameter name="Prefactor" scale="1e-18" value="5.7" min="1e-07" max="1000.0" free="1"/>
25 25 <parameter name="Index" scale="-1" value="2.48" min="0.0" max="+5.0" free="1"/>
26   - <parameter name="CutoffEnergy" scale="1e6" value="10.0" min="0.01" max="1000.0" free="1"/>
  26 + <parameter name="CutoffEnergy" scale="1e7" value="1.0" min="0.01" max="1000.0" free="1"/>
27 27 <parameter name="PivotEnergy" scale="1e6" value="0.3" min="0.01" max="1000.0" free="0"/>
28 28 </spectrum>
29   - <spatialModel type="RadialDisk">
30   - <parameter name="RA" scale="1.0" value="266.4121" min="-360" max="360" free="1"/>
31   - <parameter name="DEC" scale="1.0" value="-29.0220" min="-90" max="90" free="1"/>
32   - <parameter name="Radius" scale="1.0" value="0.01" min="0.001" max="10" free="1"/>
  29 + <spatialModel type="PointSource">
  30 + <parameter name="RA" value="266.424004498437" error="0" scale="1" free="1" />
  31 + <parameter name="DEC" value="-29.0049010253548" error="0" scale="1" free="1" />
33 32 </spatialModel>
34 33 </source>
35 34 ...
... ... @@ -49,15 +48,15 @@ done using a stacked analysis.
49 48  
50 49 .. code-block:: bash
51 50  
52   - 2017-06-02T10:54:10: === GModelSky ===
53   - 2017-06-02T10:54:10: Name ......................: Src001
54   - 2017-06-02T10:54:10: Instruments ...............: all
55   - 2017-06-02T10:54:10: Test Statistic ............: 8794.57365422026
  51 + 2017-07-28T01:26:56: === GModelSky ===
  52 + 2017-07-28T01:26:56: Name ......................: Src001
  53 + 2017-07-28T01:26:56: Instruments ...............: all
  54 + 2017-07-28T01:26:56: Test Statistic ............: 15404.1562257314
56 55 ...
57   - 2017-06-02T10:54:10: === GModelSky ===
58   - 2017-06-02T10:54:10: Name ......................: Src002
59   - 2017-06-02T10:54:10: Instruments ...............: all
60   - 2017-06-02T10:54:10: Test Statistic ............: 1439.46021600731
  56 + 2017-07-28T01:26:56: === GModelSky ===
  57 + 2017-07-28T01:26:56: Name ......................: Src002
  58 + 2017-07-28T01:26:56: Instruments ...............: all
  59 + 2017-07-28T01:26:56: Test Statistic ............: 3419.3020299729
61 60  
62 61 The Test Statistic values are also written into the
63 62 :ref:`model definition file <glossary_moddef>`
... ... @@ -67,10 +66,10 @@ created by :ref:`ctlike`:
67 66  
68 67 <?xml version="1.0" encoding="UTF-8" standalone="no"?>
69 68 <source_library title="source library">
70   - <source name="Src001" type="ExtendedSource" ts="8794.574" tscalc="1">
  69 + <source name="Src001" type="PointSource" ts="15404.156" tscalc="1">
71 70 ...
72 71 </source>
73   - <source name="Src002" type="PointSource" ts="1439.460" tscalc="1">
  72 + <source name="Src002" type="PointSource" ts="3419.302" tscalc="1">
74 73 ...
75 74 </source>
76 75 ...
... ...
doc/source/users/tutorials/howto/howto_tsmap.png
View file @7bad287

24.5 KB | W: | H:

24.8 KB | W: | H:

  • 2-up
  • Swipe
  • Onion skin
doc/source/users/tutorials/howto/howto_tsmap.rst
  View file @7bad287
... ... @@ -90,16 +90,16 @@ The figure below shows the resulting Test Statistic map.
90 90  
91 91 .. code-block:: bash
92 92  
93   - cttsmap inobs=cntcube.fits inmodel=stacked_models.xml srcname=Src001 expcube=expcube.fits psfcube=psfcube.fits bkgcube=bkgcube.fits edisp=no errors=no usepnt=no nxpix=10 nypix=10 binsz=0.05 coordsys=GAL proj=CAR xref=0.0 yref=0.0 publish=no chatter=2 clobber=yes debug=no mode=ql logL0=-227864.58395644117 binmin=0 binmax=10 outmap=tsmap_split_0.fits logfile=tsmap_split_0.log &
94   - cttsmap inobs=cntcube.fits inmodel=stacked_models.xml srcname=Src001 expcube=expcube.fits psfcube=psfcube.fits bkgcube=bkgcube.fits edisp=no errors=no usepnt=no nxpix=10 nypix=10 binsz=0.05 coordsys=GAL proj=CAR xref=0.0 yref=0.0 publish=no chatter=2 clobber=yes debug=no mode=ql logL0=-227864.58395644117 binmin=10 binmax=20 outmap=tsmap_split_1.fits logfile=tsmap_split_1.log &
95   - cttsmap inobs=cntcube.fits inmodel=stacked_models.xml srcname=Src001 expcube=expcube.fits psfcube=psfcube.fits bkgcube=bkgcube.fits edisp=no errors=no usepnt=no nxpix=10 nypix=10 binsz=0.05 coordsys=GAL proj=CAR xref=0.0 yref=0.0 publish=no chatter=2 clobber=yes debug=no mode=ql logL0=-227864.58395644117 binmin=20 binmax=30 outmap=tsmap_split_2.fits logfile=tsmap_split_2.log &
96   - cttsmap inobs=cntcube.fits inmodel=stacked_models.xml srcname=Src001 expcube=expcube.fits psfcube=psfcube.fits bkgcube=bkgcube.fits edisp=no errors=no usepnt=no nxpix=10 nypix=10 binsz=0.05 coordsys=GAL proj=CAR xref=0.0 yref=0.0 publish=no chatter=2 clobber=yes debug=no mode=ql logL0=-227864.58395644117 binmin=30 binmax=40 outmap=tsmap_split_3.fits logfile=tsmap_split_3.log &
97   - cttsmap inobs=cntcube.fits inmodel=stacked_models.xml srcname=Src001 expcube=expcube.fits psfcube=psfcube.fits bkgcube=bkgcube.fits edisp=no errors=no usepnt=no nxpix=10 nypix=10 binsz=0.05 coordsys=GAL proj=CAR xref=0.0 yref=0.0 publish=no chatter=2 clobber=yes debug=no mode=ql logL0=-227864.58395644117 binmin=40 binmax=50 outmap=tsmap_split_4.fits logfile=tsmap_split_4.log &
98   - cttsmap inobs=cntcube.fits inmodel=stacked_models.xml srcname=Src001 expcube=expcube.fits psfcube=psfcube.fits bkgcube=bkgcube.fits edisp=no errors=no usepnt=no nxpix=10 nypix=10 binsz=0.05 coordsys=GAL proj=CAR xref=0.0 yref=0.0 publish=no chatter=2 clobber=yes debug=no mode=ql logL0=-227864.58395644117 binmin=50 binmax=60 outmap=tsmap_split_5.fits logfile=tsmap_split_5.log &
99   - cttsmap inobs=cntcube.fits inmodel=stacked_models.xml srcname=Src001 expcube=expcube.fits psfcube=psfcube.fits bkgcube=bkgcube.fits edisp=no errors=no usepnt=no nxpix=10 nypix=10 binsz=0.05 coordsys=GAL proj=CAR xref=0.0 yref=0.0 publish=no chatter=2 clobber=yes debug=no mode=ql logL0=-227864.58395644117 binmin=60 binmax=70 outmap=tsmap_split_6.fits logfile=tsmap_split_6.log &
100   - cttsmap inobs=cntcube.fits inmodel=stacked_models.xml srcname=Src001 expcube=expcube.fits psfcube=psfcube.fits bkgcube=bkgcube.fits edisp=no errors=no usepnt=no nxpix=10 nypix=10 binsz=0.05 coordsys=GAL proj=CAR xref=0.0 yref=0.0 publish=no chatter=2 clobber=yes debug=no mode=ql logL0=-227864.58395644117 binmin=70 binmax=80 outmap=tsmap_split_7.fits logfile=tsmap_split_7.log &
101   - cttsmap inobs=cntcube.fits inmodel=stacked_models.xml srcname=Src001 expcube=expcube.fits psfcube=psfcube.fits bkgcube=bkgcube.fits edisp=no errors=no usepnt=no nxpix=10 nypix=10 binsz=0.05 coordsys=GAL proj=CAR xref=0.0 yref=0.0 publish=no chatter=2 clobber=yes debug=no mode=ql logL0=-227864.58395644117 binmin=80 binmax=90 outmap=tsmap_split_8.fits logfile=tsmap_split_8.log &
102   - cttsmap inobs=cntcube.fits inmodel=stacked_models.xml srcname=Src001 expcube=expcube.fits psfcube=psfcube.fits bkgcube=bkgcube.fits edisp=no errors=no usepnt=no nxpix=10 nypix=10 binsz=0.05 coordsys=GAL proj=CAR xref=0.0 yref=0.0 publish=no chatter=2 clobber=yes debug=no mode=ql logL0=-227864.58395644117 binmin=90 binmax=100 outmap=tsmap_split_9.fits logfile=tsmap_split_9.log &
  93 + cttsmap inobs=cntcube.fits inmodel=stacked_models.xml srcname=Src001 expcube=expcube.fits psfcube=psfcube.fits bkgcube=bkgcube.fits edisp=no errors=no usepnt=no nxpix=10 nypix=10 binsz=0.05 coordsys=GAL proj=CAR xref=0.0 yref=0.0 publish=no chatter=2 clobber=yes debug=no mode=ql logL0=-274462.5864105515 binmin=0 binmax=10 outmap=tsmap_split_0.fits logfile=tsmap_split_0.log &
  94 + cttsmap inobs=cntcube.fits inmodel=stacked_models.xml srcname=Src001 expcube=expcube.fits psfcube=psfcube.fits bkgcube=bkgcube.fits edisp=no errors=no usepnt=no nxpix=10 nypix=10 binsz=0.05 coordsys=GAL proj=CAR xref=0.0 yref=0.0 publish=no chatter=2 clobber=yes debug=no mode=ql logL0=-274462.5864105515 binmin=10 binmax=20 outmap=tsmap_split_1.fits logfile=tsmap_split_1.log &
  95 + cttsmap inobs=cntcube.fits inmodel=stacked_models.xml srcname=Src001 expcube=expcube.fits psfcube=psfcube.fits bkgcube=bkgcube.fits edisp=no errors=no usepnt=no nxpix=10 nypix=10 binsz=0.05 coordsys=GAL proj=CAR xref=0.0 yref=0.0 publish=no chatter=2 clobber=yes debug=no mode=ql logL0=-274462.5864105515 binmin=20 binmax=30 outmap=tsmap_split_2.fits logfile=tsmap_split_2.log &
  96 + cttsmap inobs=cntcube.fits inmodel=stacked_models.xml srcname=Src001 expcube=expcube.fits psfcube=psfcube.fits bkgcube=bkgcube.fits edisp=no errors=no usepnt=no nxpix=10 nypix=10 binsz=0.05 coordsys=GAL proj=CAR xref=0.0 yref=0.0 publish=no chatter=2 clobber=yes debug=no mode=ql logL0=-274462.5864105515 binmin=30 binmax=40 outmap=tsmap_split_3.fits logfile=tsmap_split_3.log &
  97 + cttsmap inobs=cntcube.fits inmodel=stacked_models.xml srcname=Src001 expcube=expcube.fits psfcube=psfcube.fits bkgcube=bkgcube.fits edisp=no errors=no usepnt=no nxpix=10 nypix=10 binsz=0.05 coordsys=GAL proj=CAR xref=0.0 yref=0.0 publish=no chatter=2 clobber=yes debug=no mode=ql logL0=-274462.5864105515 binmin=40 binmax=50 outmap=tsmap_split_4.fits logfile=tsmap_split_4.log &
  98 + cttsmap inobs=cntcube.fits inmodel=stacked_models.xml srcname=Src001 expcube=expcube.fits psfcube=psfcube.fits bkgcube=bkgcube.fits edisp=no errors=no usepnt=no nxpix=10 nypix=10 binsz=0.05 coordsys=GAL proj=CAR xref=0.0 yref=0.0 publish=no chatter=2 clobber=yes debug=no mode=ql logL0=-274462.5864105515 binmin=50 binmax=60 outmap=tsmap_split_5.fits logfile=tsmap_split_5.log &
  99 + cttsmap inobs=cntcube.fits inmodel=stacked_models.xml srcname=Src001 expcube=expcube.fits psfcube=psfcube.fits bkgcube=bkgcube.fits edisp=no errors=no usepnt=no nxpix=10 nypix=10 binsz=0.05 coordsys=GAL proj=CAR xref=0.0 yref=0.0 publish=no chatter=2 clobber=yes debug=no mode=ql logL0=-274462.5864105515 binmin=60 binmax=70 outmap=tsmap_split_6.fits logfile=tsmap_split_6.log &
  100 + cttsmap inobs=cntcube.fits inmodel=stacked_models.xml srcname=Src001 expcube=expcube.fits psfcube=psfcube.fits bkgcube=bkgcube.fits edisp=no errors=no usepnt=no nxpix=10 nypix=10 binsz=0.05 coordsys=GAL proj=CAR xref=0.0 yref=0.0 publish=no chatter=2 clobber=yes debug=no mode=ql logL0=-274462.5864105515 binmin=70 binmax=80 outmap=tsmap_split_7.fits logfile=tsmap_split_7.log &
  101 + cttsmap inobs=cntcube.fits inmodel=stacked_models.xml srcname=Src001 expcube=expcube.fits psfcube=psfcube.fits bkgcube=bkgcube.fits edisp=no errors=no usepnt=no nxpix=10 nypix=10 binsz=0.05 coordsys=GAL proj=CAR xref=0.0 yref=0.0 publish=no chatter=2 clobber=yes debug=no mode=ql logL0=-274462.5864105515 binmin=80 binmax=90 outmap=tsmap_split_8.fits logfile=tsmap_split_8.log &
  102 + cttsmap inobs=cntcube.fits inmodel=stacked_models.xml srcname=Src001 expcube=expcube.fits psfcube=psfcube.fits bkgcube=bkgcube.fits edisp=no errors=no usepnt=no nxpix=10 nypix=10 binsz=0.05 coordsys=GAL proj=CAR xref=0.0 yref=0.0 publish=no chatter=2 clobber=yes debug=no mode=ql logL0=-274462.5864105515 binmin=90 binmax=100 outmap=tsmap_split_9.fits logfile=tsmap_split_9.log &
103 103  
104 104 Now run the ASCII file by typing
105 105  
... ...
doc/source/users/tutorials/howto/howto_ulimit.rst
  View file @7bad287
... ... @@ -51,12 +51,12 @@ extracted from the :ref:`ctulimit` log file:
51 51  
52 52 .. code-block:: bash
53 53  
54   - 2017-06-08T15:39:54: +=====================+
55   - 2017-06-08T15:39:54: | Upper limit results |
56   - 2017-06-08T15:39:54: +=====================+
57   - 2017-06-08T15:39:54: Differential flux limit ...: 2.4438738640373e-20 ph/cm2/s/MeV at 1 TeV
58   - 2017-06-08T15:39:54: Integral flux limit .......: 1.64945554704815e-14 ph/cm2/s within [1-100] TeV
59   - 2017-06-08T15:39:54: Energy flux limit .........: 7.26289446286815e-14 erg/cm2/s within [1-100] TeV
  54 + 2017-07-28T07:03:03: +=====================+
  55 + 2017-07-28T07:03:03: | Upper limit results |
  56 + 2017-07-28T07:03:03: +=====================+
  57 + 2017-07-28T07:03:03: Differential flux limit ...: 6.39388221102118e-21 ph/cm2/s/MeV at 1 TeV
  58 + 2017-07-28T07:03:03: Integral flux limit .......: 4.31545368823519e-15 ph/cm2/s within [1-100] TeV
  59 + 2017-07-28T07:03:03: Energy flux limit .........: 1.90018365473006e-14 erg/cm2/s within [1-100] TeV
60 60  
61 61 .. warning::
62 62 It is important to fix the spectral index of the test source in the
... ...
doc/source/users/tutorials/index.rst
  View file @7bad287
... ... @@ -8,6 +8,6 @@ Tutorials
8 8  
9 9 quickstart/index
10 10 beyond/index
11   - howto/index
12 11 1dc/index
13   - iact_analysis/index
14 12 \ No newline at end of file
  13 + howto/index
  14 + iact_analysis/index
... ...