Commit e948b1d582c4b8e785c80fac671c70ef8040ab71
1 parent
39a15056
Update first Data Challenge tutorial
Showing
15 changed files
with
112 additions
and
188 deletions
doc/source/users/tutorials/1dc/first_detect.rst
| ... | ... | @@ -76,8 +76,8 @@ for a model fitting with :ref:`ctlike`. |
| 76 | 76 | <parameter name="PivotEnergy" value="1" scale="300000" free="0" /> |
| 77 | 77 | </spectrum> |
| 78 | 78 | <spatialModel type="PointSource"> |
| 79 | - <parameter name="RA" value="266.4120906928" error="0" scale="1" free="1" /> | |
| 80 | - <parameter name="DEC" value="-29.0219729468991" error="0" scale="1" free="1" /> | |
| 79 | + <parameter name="RA" value="266.424004498437" error="0" scale="1" free="1" /> | |
| 80 | + <parameter name="DEC" value="-29.0049010253548" error="0" scale="1" free="1" /> | |
| 81 | 81 | </spatialModel> |
| 82 | 82 | </source> |
| 83 | 83 | <source name="Src002" type="PointSource"> |
| ... | ... | @@ -87,9 +87,9 @@ for a model fitting with :ref:`ctlike`. |
| 87 | 87 | <parameter name="PivotEnergy" value="1" scale="300000" free="0" /> |
| 88 | 88 | </spectrum> |
| 89 | 89 | <spatialModel type="PointSource"> |
| 90 | - <parameter name="RA" value="266.793148201606" error="0" scale="1" free="1" /> | |
| 91 | - <parameter name="DEC" value="-28.1253038992376" error="0" scale="1" free="1" /> | |
| 92 | - </spatialModel> | |
| 90 | + <parameter name="RA" value="266.831945177213" error="0" scale="1" free="1" /> | |
| 91 | + <parameter name="DEC" value="-28.1460284439951" error="0" scale="1" free="1" /> | |
| 92 | + </spatialModel> | |
| 93 | 93 | </source> |
| 94 | 94 | <source name="Background" type="CTAIrfBackground"> |
| 95 | 95 | <spectrum type="PowerLaw"> | ... | ... |
doc/source/users/tutorials/1dc/first_fitting.rst
| ... | ... | @@ -48,49 +48,34 @@ converged properly, as illustrated in the example below: |
| 48 | 48 | |
| 49 | 49 | .. code-block:: bash |
| 50 | 50 | |
| 51 | - 2017-06-08T08:20:32: +=================================+ | |
| 52 | - 2017-06-08T08:20:32: | Maximum likelihood optimisation | | |
| 53 | - 2017-06-08T08:20:32: +=================================+ | |
| 54 | - 2017-06-08T08:20:46: >Iteration 0: -logL=242767.783, Lambda=1.0e-03 | |
| 55 | - 2017-06-08T08:21:00: >Iteration 1: -logL=222518.221, Lambda=1.0e-03, delta=20249.562, step=1.0e+00, max(|grad|)=50533.167265 [Index:13] | |
| 56 | - 2017-06-08T08:21:14: >Iteration 2: -logL=218997.783, Lambda=1.0e-04, delta=3520.439, step=1.0e+00, max(|grad|)=-13548.755375 [DEC:1] | |
| 57 | - 2017-06-08T08:21:29: >Iteration 3: -logL=217965.881, Lambda=1.0e-05, delta=1031.901, step=1.0e+00, max(|grad|)=-11344.447813 [DEC:1] | |
| 58 | - 2017-06-08T08:21:43: >Iteration 4: -logL=217775.661, Lambda=1.0e-06, delta=190.220, step=1.0e+00, max(|grad|)=-7778.127593 [DEC:1] | |
| 59 | - 2017-06-08T08:21:57: >Iteration 5: -logL=217725.868, Lambda=1.0e-07, delta=49.793, step=1.0e+00, max(|grad|)=-4796.822100 [DEC:1] | |
| 60 | - 2017-06-08T08:22:11: >Iteration 6: -logL=217701.511, Lambda=1.0e-08, delta=24.357, step=1.0e+00, max(|grad|)=-3267.160620 [RA:6] | |
| 61 | - 2017-06-08T08:22:25: >Iteration 7: -logL=217688.626, Lambda=1.0e-09, delta=12.885, step=1.0e+00, max(|grad|)=-2438.542445 [RA:6] | |
| 62 | - 2017-06-08T08:22:40: >Iteration 8: -logL=217681.375, Lambda=1.0e-10, delta=7.250, step=1.0e+00, max(|grad|)=-1905.466691 [RA:6] | |
| 63 | - 2017-06-08T08:22:54: >Iteration 9: -logL=217677.367, Lambda=1.0e-11, delta=4.008, step=1.0e+00, max(|grad|)=-1338.131602 [RA:6] | |
| 64 | - 2017-06-08T08:23:08: >Iteration 10: -logL=217675.391, Lambda=1.0e-12, delta=1.977, step=1.0e+00, max(|grad|)=-1117.908664 [RA:6] | |
| 65 | - 2017-06-08T08:23:22: >Iteration 11: -logL=217674.202, Lambda=1.0e-13, delta=1.188, step=1.0e+00, max(|grad|)=-683.041271 [RA:6] | |
| 66 | - 2017-06-08T08:23:36: >Iteration 12: -logL=217673.712, Lambda=1.0e-14, delta=0.490, step=1.0e+00, max(|grad|)=-481.729534 [RA:6] | |
| 67 | - 2017-06-08T08:23:50: >Iteration 13: -logL=217673.495, Lambda=1.0e-15, delta=0.217, step=1.0e+00, max(|grad|)=-385.504388 [RA:6] | |
| 68 | - 2017-06-08T08:24:04: >Iteration 14: -logL=217673.327, Lambda=1.0e-16, delta=0.168, step=1.0e+00, max(|grad|)=-314.151335 [RA:6] | |
| 69 | - 2017-06-08T08:24:18: >Iteration 15: -logL=217673.231, Lambda=1.0e-17, delta=0.096, step=1.0e+00, max(|grad|)=-227.509872 [RA:6] | |
| 70 | - 2017-06-08T08:24:32: >Iteration 16: -logL=217673.170, Lambda=1.0e-18, delta=0.060, step=1.0e+00, max(|grad|)=-176.781601 [RA:6] | |
| 71 | - 2017-06-08T08:24:47: >Iteration 17: -logL=217673.138, Lambda=1.0e-19, delta=0.032, step=1.0e+00, max(|grad|)=-124.109305 [RA:6] | |
| 72 | - 2017-06-08T08:25:01: >Iteration 18: -logL=217673.124, Lambda=1.0e-20, delta=0.014, step=1.0e+00, max(|grad|)=-88.895131 [RA:6] | |
| 73 | - 2017-06-08T08:25:15: >Iteration 19: -logL=217673.117, Lambda=1.0e-21, delta=0.007, step=1.0e+00, max(|grad|)=-63.724112 [RA:6] | |
| 74 | - 2017-06-08T08:25:29: >Iteration 20: -logL=217673.112, Lambda=1.0e-22, delta=0.006, step=1.0e+00, max(|grad|)=61.601206 [DEC:7] | |
| 75 | - 2017-06-08T08:25:43: >Iteration 21: -logL=217673.107, Lambda=1.0e-23, delta=0.005, step=1.0e+00, max(|grad|)=58.318115 [DEC:7] | |
| 76 | - 2017-06-08T08:25:57: >Iteration 22: -logL=217673.101, Lambda=1.0e-24, delta=0.006, step=1.0e+00, max(|grad|)=50.936104 [DEC:7] | |
| 77 | - 2017-06-08T08:26:12: >Iteration 23: -logL=217673.097, Lambda=1.0e-25, delta=0.004, step=1.0e+00, max(|grad|)=46.423925 [DEC:7] | |
| 78 | - 2017-06-08T08:26:26: | |
| 79 | - 2017-06-08T08:26:26: +=========================================+ | |
| 80 | - 2017-06-08T08:26:26: | Maximum likelihood optimisation results | | |
| 81 | - 2017-06-08T08:26:26: +=========================================+ | |
| 82 | - 2017-06-08T08:26:26: === GOptimizerLM === | |
| 83 | - 2017-06-08T08:26:26: Optimized function value ..: 217673.097 | |
| 84 | - 2017-06-08T08:26:26: Absolute precision ........: 0.005 | |
| 85 | - 2017-06-08T08:26:26: Acceptable value decrease .: 2 | |
| 86 | - 2017-06-08T08:26:26: Optimization status .......: converged | |
| 87 | - 2017-06-08T08:26:26: Number of parameters ......: 16 | |
| 88 | - 2017-06-08T08:26:26: Number of free parameters .: 10 | |
| 89 | - 2017-06-08T08:26:26: Number of iterations ......: 23 | |
| 90 | - 2017-06-08T08:26:26: Lambda ....................: 1e-26 | |
| 91 | - 2017-06-08T08:26:26: Maximum log likelihood ....: -217673.097 | |
| 92 | - 2017-06-08T08:26:26: Observed events (Nobs) ...: 2205594.000 | |
| 93 | - 2017-06-08T08:26:26: Predicted events (Npred) ..: 2205593.994 (Nobs - Npred = 0.00562562886625528) | |
| 51 | + 2017-07-27T14:05:58: +=================================+ | |
| 52 | + 2017-07-27T14:05:58: | Maximum likelihood optimisation | | |
| 53 | + 2017-07-27T14:05:58: +=================================+ | |
| 54 | + 2017-07-27T14:06:07: >Iteration 0: -logL=291743.621, Lambda=1.0e-03 | |
| 55 | + 2017-07-27T14:06:17: >Iteration 1: -logL=266150.955, Lambda=1.0e-03, delta=25592.666, step=1.0e+00, max(|grad|)=65223.588855 [Index:13] | |
| 56 | + 2017-07-27T14:06:26: >Iteration 2: -logL=260665.430, Lambda=1.0e-04, delta=5485.525, step=1.0e+00, max(|grad|)=10907.804073 [Index:3] | |
| 57 | + 2017-07-27T14:06:36: >Iteration 3: -logL=258904.835, Lambda=1.0e-05, delta=1760.595, step=1.0e+00, max(|grad|)=6694.957421 [Index:3] | |
| 58 | + 2017-07-27T14:06:45: >Iteration 4: -logL=258545.054, Lambda=1.0e-06, delta=359.781, step=1.0e+00, max(|grad|)=2242.517259 [RA:0] | |
| 59 | + 2017-07-27T14:06:55: >Iteration 5: -logL=258531.804, Lambda=1.0e-07, delta=13.250, step=1.0e+00, max(|grad|)=631.774929 [RA:0] | |
| 60 | + 2017-07-27T14:07:05: >Iteration 6: -logL=258531.685, Lambda=1.0e-08, delta=0.120, step=1.0e+00, max(|grad|)=139.793278 [RA:0] | |
| 61 | + 2017-07-27T14:07:14: >Iteration 7: -logL=258531.673, Lambda=1.0e-09, delta=0.012, step=1.0e+00, max(|grad|)=62.479910 [DEC:1] | |
| 62 | + 2017-07-27T14:07:24: >Iteration 8: -logL=258531.670, Lambda=1.0e-10, delta=0.002, step=1.0e+00, max(|grad|)=29.031154 [DEC:1] | |
| 63 | + 2017-07-27T14:07:33: | |
| 64 | + 2017-07-27T14:07:33: +=========================================+ | |
| 65 | + 2017-07-27T14:07:33: | Maximum likelihood optimisation results | | |
| 66 | + 2017-07-27T14:07:33: +=========================================+ | |
| 67 | + 2017-07-27T14:07:33: === GOptimizerLM === | |
| 68 | + 2017-07-27T14:07:33: Optimized function value ..: 258531.670 | |
| 69 | + 2017-07-27T14:07:33: Absolute precision ........: 0.005 | |
| 70 | + 2017-07-27T14:07:33: Acceptable value decrease .: 2 | |
| 71 | + 2017-07-27T14:07:33: Optimization status .......: converged | |
| 72 | + 2017-07-27T14:07:33: Number of parameters ......: 16 | |
| 73 | + 2017-07-27T14:07:33: Number of free parameters .: 10 | |
| 74 | + 2017-07-27T14:07:33: Number of iterations ......: 8 | |
| 75 | + 2017-07-27T14:07:33: Lambda ....................: 1e-11 | |
| 76 | + 2017-07-27T14:07:33: Maximum log likelihood ....: -258531.670 | |
| 77 | + 2017-07-27T14:07:33: Observed events (Nobs) ...: 2069142.000 | |
| 78 | + 2017-07-27T14:07:33: Predicted events (Npred) ..: 2069141.997 (Nobs - Npred = 0.00287596974521875) | |
| 94 | 79 | |
| 95 | 80 | You may also convert the fitted model positions into a `ds9 <http://ds9.si.edu>`_ |
| 96 | 81 | region file using the :ref:`csmodelinfo` script so that you can overlay the | ... | ... |
doc/source/users/tutorials/1dc/first_improving.rst
| ... | ... | @@ -24,7 +24,7 @@ You do this with the :ref:`csresmap` script by providing the output |
| 24 | 24 | Input background cube file (only needed for stacked analysis) [NONE] bkgcube.fits |
| 25 | 25 | Input model definition XML file [$CTOOLS/share/models/crab.xml] stacked_results_cutoff.xml |
| 26 | 26 | Residual map computation algorithm (SUB|SUBDIV|SUBDIVSQRT|SIGNIFICANCE) [SUBDIV] SUB |
| 27 | - Output residual map file [resmap.fits] resmap.fits | |
| 27 | + Output residual map file [resmap.fits] | |
| 28 | 28 | |
| 29 | 29 | This produces the file ``resmap.fits`` that contains a residual map that |
| 30 | 30 | you can display for example with `ds9 <http://ds9.si.edu>`_. |
| ... | ... | @@ -52,12 +52,12 @@ as illustrated below: |
| 52 | 52 | <spectrum type="ExponentialCutoffPowerLaw"> |
| 53 | 53 | <parameter name="Prefactor" scale="1e-18" value="5.7" min="1e-07" max="1000.0" free="1"/> |
| 54 | 54 | <parameter name="Index" scale="-1" value="2.48" min="0.0" max="+5.0" free="1"/> |
| 55 | - <parameter name="CutoffEnergy" scale="1e6" value="10.0" min="0.01" max="1000.0" free="1"/> | |
| 55 | + <parameter name="CutoffEnergy" scale="1e7" value="1.0" min="0.01" max="1000.0" free="1"/> | |
| 56 | 56 | <parameter name="PivotEnergy" scale="1e6" value="0.3" min="0.01" max="1000.0" free="0"/> |
| 57 | 57 | </spectrum> |
| 58 | 58 | <spatialModel type="PointSource"> |
| 59 | - <parameter name="RA" value="266.4120906928" error="0" scale="1" free="1" /> | |
| 60 | - <parameter name="DEC" value="-29.0219729468991" error="0" scale="1" free="1" /> | |
| 59 | + <parameter name="RA" value="266.424004498437" error="0" scale="1" free="1" /> | |
| 60 | + <parameter name="DEC" value="-29.0049010253548" error="0" scale="1" free="1" /> | |
| 61 | 61 | </spatialModel> |
| 62 | 62 | </source> |
| 63 | 63 | <source name="Src002" type="PointSource"> |
| ... | ... | @@ -67,8 +67,8 @@ as illustrated below: |
| 67 | 67 | <parameter name="PivotEnergy" value="1" scale="300000" free="0" /> |
| 68 | 68 | </spectrum> |
| 69 | 69 | <spatialModel type="PointSource"> |
| 70 | - <parameter name="RA" value="266.793148201606" error="0" scale="1" free="1" /> | |
| 71 | - <parameter name="DEC" value="-28.1253038992376" error="0" scale="1" free="1" /> | |
| 70 | + <parameter name="RA" value="266.831945177213" error="0" scale="1" free="1" /> | |
| 71 | + <parameter name="DEC" value="-28.1460284439951" error="0" scale="1" free="1" /> | |
| 72 | 72 | </spatialModel> |
| 73 | 73 | </source> |
| 74 | 74 | <source name="IEM" type="DiffuseSource"> |
| ... | ... | @@ -79,10 +79,10 @@ as illustrated below: |
| 79 | 79 | <parameter name="Normalization" value="1" scale="1" min="0.001" max="1000" free="0" /> |
| 80 | 80 | </spatialModel> |
| 81 | 81 | </source> |
| 82 | - <source name="Background" type="CTAIrfBackground"> | |
| 82 | + <source name="BackgroundModel" type="CTACubeBackground"> | |
| 83 | 83 | <spectrum type="PowerLaw"> |
| 84 | - <parameter name="Prefactor" value="1" scale="1" min="0.1" max="10" free="1" /> | |
| 85 | - <parameter name="Index" value="0" scale="1" min="-10" max="10" free="1" /> | |
| 84 | + <parameter name="Prefactor" value="1" error="0" scale="1" min="0.01" max="100" free="1" /> | |
| 85 | + <parameter name="Index" value="0" error="0" scale="1" min="-5" max="5" free="1" /> | |
| 86 | 86 | <parameter name="PivotEnergy" value="1" scale="1000000" free="0" /> |
| 87 | 87 | </spectrum> |
| 88 | 88 | </source> |
| ... | ... | @@ -90,77 +90,16 @@ as illustrated below: |
| 90 | 90 | |
| 91 | 91 | Repeating the fit with this model and producing a corresponding residual map |
| 92 | 92 | produces the map shown below. Now, the residuals near the two point sources |
| 93 | -are flatter and the diffuse emission has disappeared. However, there | |
| 94 | -is still a clear emission hole near ``Src001`` which suggests that this source | |
| 95 | -may be extended. | |
| 93 | +are flatter and the diffuse emission has largely disappeared. The diffuse | |
| 94 | +emission is in fact a bit oversubtracted, which is probably due to the | |
| 95 | +additional source components visible in the residual map. It is left as an | |
| 96 | +exercise to the User to include model components for these additional | |
| 97 | +components in the | |
| 98 | +:ref:`model definition file <glossary_moddef>` | |
| 99 | +in order to reduce the fit residuals. | |
| 96 | 100 | |
| 97 | 101 | .. figure:: first_skymap_residual_iem.png |
| 98 | 102 | :width: 400px |
| 99 | 103 | :align: center |
| 100 | 104 | |
| 101 | 105 | *Residual sky map after subtraction of the fitted model including a diffuse emission component* |
| 102 | - | |
| 103 | -In our attempt to iteratively refine the source model, we therefore use in our next | |
| 104 | -iteration a radial disk source instead of a point source for ``Src001`` by | |
| 105 | -specifying the following | |
| 106 | -:ref:`model definition file <glossary_moddef>` | |
| 107 | -on input to the :ref:`ctlike` tool: | |
| 108 | - | |
| 109 | -.. code-block:: xml | |
| 110 | - | |
| 111 | - <?xml version="1.0" encoding="UTF-8" standalone="no"?> | |
| 112 | - <source_library title="source library"> | |
| 113 | - <source name="Src001" type="PointSource"> | |
| 114 | - <spectrum type="ExponentialCutoffPowerLaw"> | |
| 115 | - <parameter name="Prefactor" scale="1e-18" value="5.7" min="1e-07" max="1000.0" free="1"/> | |
| 116 | - <parameter name="Index" scale="-1" value="2.48" min="0.0" max="+5.0" free="1"/> | |
| 117 | - <parameter name="CutoffEnergy" scale="1e6" value="10.0" min="0.01" max="1000.0" free="1"/> | |
| 118 | - <parameter name="PivotEnergy" scale="1e6" value="0.3" min="0.01" max="1000.0" free="0"/> | |
| 119 | - </spectrum> | |
| 120 | - <spatialModel type="RadialDisk"> | |
| 121 | - <parameter name="RA" value="266.4121" scale="1.0" min="-360" max="360" free="1"/> | |
| 122 | - <parameter name="DEC" value="-29.0220" scale="1.0" min="-90" max="90" free="1"/> | |
| 123 | - <parameter name="Radius" value="0.01" scale="1.0" min="0.001" max="10" free="1"/> | |
| 124 | - </spatialModel> | |
| 125 | - </source> | |
| 126 | - <source name="Src002" type="PointSource"> | |
| 127 | - <spectrum type="PowerLaw"> | |
| 128 | - <parameter name="Prefactor" value="1" error="0" scale="5.7e-18" min="0" free="1" /> | |
| 129 | - <parameter name="Index" value="1" error="-0" scale="-2.48" min="-4.03225806451613" max="4.03225806451613" free="1" /> | |
| 130 | - <parameter name="PivotEnergy" value="1" scale="300000" free="0" /> | |
| 131 | - </spectrum> | |
| 132 | - <spatialModel type="PointSource"> | |
| 133 | - <parameter name="RA" value="266.793148201606" error="0" scale="1" free="1" /> | |
| 134 | - <parameter name="DEC" value="-28.1253038992376" error="0" scale="1" free="1" /> | |
| 135 | - </spatialModel> | |
| 136 | - </source> | |
| 137 | - <source name="IEM" type="DiffuseSource"> | |
| 138 | - <spectrum type="ConstantValue"> | |
| 139 | - <parameter name="Value" value="1" scale="1" min="1e-05" max="100000" free="1" /> | |
| 140 | - </spectrum> | |
| 141 | - <spatialModel type="MapCubeFunction" file="$CTADATA/models/cube_iem.fits.gz"> | |
| 142 | - <parameter name="Normalization" value="1" scale="1" min="0.001" max="1000" free="0" /> | |
| 143 | - </spatialModel> | |
| 144 | - </source> | |
| 145 | - <source name="Background" type="CTAIrfBackground"> | |
| 146 | - <spectrum type="PowerLaw"> | |
| 147 | - <parameter name="Prefactor" value="1" scale="1" min="0.1" max="10" free="1" /> | |
| 148 | - <parameter name="Index" value="0" scale="1" min="-10" max="10" free="1" /> | |
| 149 | - <parameter name="PivotEnergy" value="1" scale="1000000" free="0" /> | |
| 150 | - </spectrum> | |
| 151 | - </source> | |
| 152 | - </source_library> | |
| 153 | - | |
| 154 | -As shown in the residual map below, the emission hole at the position of | |
| 155 | -``Src001`` has now disappeared. The disk model has in fact been fitted with | |
| 156 | -a radius of 5.3 +/- 0.1 arcmin, indicating that ``Src001`` is significantly | |
| 157 | -extended. The remaining residuals suggest that also ``Src002`` may be | |
| 158 | -slightly extended, but it is left as an exercise to the User to test this | |
| 159 | -hypothesis. | |
| 160 | - | |
| 161 | -.. figure:: first_skymap_residual_iem_disk.png | |
| 162 | - :width: 400px | |
| 163 | - :align: center | |
| 164 | - | |
| 165 | - *Residual sky map after subtraction of the fitted model including a diffuse emission component and a disk model for Src001* | |
| 166 | - | ... | ... |
doc/source/users/tutorials/1dc/first_select_events.rst
| ... | ... | @@ -46,33 +46,33 @@ displays a summary of the content of this new |
| 46 | 46 | |
| 47 | 47 | .. code-block:: bash |
| 48 | 48 | |
| 49 | - 2017-06-08T07:42:01: +=========+ | |
| 50 | - 2017-06-08T07:42:01: | Summary | | |
| 51 | - 2017-06-08T07:42:01: +=========+ | |
| 52 | - 2017-06-08T07:42:01: === Observations === | |
| 53 | - 2017-06-08T07:42:01: Unbinned observations .....: 57 | |
| 54 | - 2017-06-08T07:42:01: Binned observations .......: 0 | |
| 55 | - 2017-06-08T07:42:01: === Events === | |
| 56 | - 2017-06-08T07:42:01: Number of events ..........: 3349183 | |
| 57 | - 2017-06-08T07:42:01: Number of bins ............: 0 | |
| 58 | - 2017-06-08T07:42:01: === Pointings === | |
| 59 | - 2017-06-08T07:42:01: Mean offset angle .........: Unknown | |
| 60 | - 2017-06-08T07:42:01: Mean zenith angle .........: 0.00 deg | |
| 61 | - 2017-06-08T07:42:01: Mean azimuth angle ........: 0.00 deg | |
| 62 | - 2017-06-08T07:42:01: === Energy range === | |
| 63 | - 2017-06-08T07:42:01: Minimum energy ............: 100 GeV | |
| 64 | - 2017-06-08T07:42:01: Maximum energy ............: 100 TeV | |
| 65 | - 2017-06-08T07:42:01: === Time range === | |
| 66 | - 2017-06-08T07:42:01: MJD (days) ................: 59235.500 - 59276.921 | |
| 67 | - 2017-06-08T07:42:01: UTC .......................: 2021-01-21T11:58:51 - 2021-03-03T22:04:51 | |
| 68 | - 2017-06-08T07:42:01: MET (seconds) .............: 664502400.000 - 668081160.000 | |
| 69 | - 2017-06-08T07:42:01: Total ontime ..............: 102600.00 s = 1710.00 min = 28.50 h | |
| 70 | - 2017-06-08T07:42:01: Total livetime ............: 100548.00 s = 1675.80 min = 27.93 h | |
| 49 | + 2017-07-27T10:02:26: +=========+ | |
| 50 | + 2017-07-27T10:02:26: | Summary | | |
| 51 | + 2017-07-27T10:02:26: +=========+ | |
| 52 | + 2017-07-27T10:02:26: === Observations === | |
| 53 | + 2017-07-27T10:02:26: Unbinned observations .....: 57 | |
| 54 | + 2017-07-27T10:02:26: Binned observations .......: 0 | |
| 55 | + 2017-07-27T10:02:26: === Events === | |
| 56 | + 2017-07-27T10:02:26: Number of events ..........: 3084639 | |
| 57 | + 2017-07-27T10:02:26: Number of bins ............: 0 | |
| 58 | + 2017-07-27T10:02:26: === Pointings === | |
| 59 | + 2017-07-27T10:02:26: Mean offset angle .........: Unknown | |
| 60 | + 2017-07-27T10:02:26: Mean zenith angle .........: 0.00 deg | |
| 61 | + 2017-07-27T10:02:26: Mean azimuth angle ........: 0.00 deg | |
| 62 | + 2017-07-27T10:02:26: === Energy range === | |
| 63 | + 2017-07-27T10:02:26: Minimum energy ............: 100 GeV | |
| 64 | + 2017-07-27T10:02:26: Maximum energy ............: 100 TeV | |
| 65 | + 2017-07-27T10:02:26: === Time range === | |
| 66 | + 2017-07-27T10:02:26: MJD (days) ................: 59235.500 - 59276.921 | |
| 67 | + 2017-07-27T10:02:26: UTC .......................: 2021-01-21T11:58:51 - 2021-03-03T22:04:51 | |
| 68 | + 2017-07-27T10:02:26: MET (seconds) .............: 664502400.000 - 668081160.000 | |
| 69 | + 2017-07-27T10:02:26: Total ontime ..............: 102600.00 s = 1710.00 min = 28.50 h | |
| 70 | + 2017-07-27T10:02:26: Total livetime ............: 100548.00 s = 1675.80 min = 27.93 h | |
| 71 | 71 | |
| 72 | 72 | There are still 57 observations in the |
| 73 | 73 | :ref:`observation definition file <glossary_obsdef>` |
| 74 | 74 | but now the events are restricted to the energy interval 100 GeV - 100 TeV. In |
| 75 | -total there are 3 349 183 events within this interval in the selected observations. | |
| 75 | +total there are 3 084 639 events within this interval in the selected observations. | |
| 76 | 76 | |
| 77 | 77 | The content of ``obs_selected.xml`` will look similar to the content of |
| 78 | 78 | ``obs.xml`` with the original event list file names replaced by the names of |
| ... | ... | @@ -82,17 +82,17 @@ the new event files: |
| 82 | 82 | |
| 83 | 83 | <?xml version="1.0" encoding="UTF-8" standalone="no"?> |
| 84 | 84 | <observation_list title="observation list"> |
| 85 | - <observation name="GPS" id="120380" instrument="CTA"> | |
| 86 | - <parameter name="EventList" file="selected_gps_baseline_120380.fits" /> | |
| 85 | + <observation name="GPS" id="110380" instrument="CTA"> | |
| 86 | + <parameter name="EventList" file="selected_gps_baseline_110380.fits" /> | |
| 87 | 87 | <parameter name="Calibration" database="1dc" response="South_z20_50h" /> |
| 88 | 88 | </observation> |
| 89 | - <observation name="GPS" id="120381" instrument="CTA"> | |
| 90 | - <parameter name="EventList" file="selected_gps_baseline_120381.fits" /> | |
| 89 | + <observation name="GPS" id="110381" instrument="CTA"> | |
| 90 | + <parameter name="EventList" file="selected_gps_baseline_110381.fits" /> | |
| 91 | 91 | <parameter name="Calibration" database="1dc" response="South_z20_50h" /> |
| 92 | 92 | </observation> |
| 93 | 93 | ... |
| 94 | - <observation name="GPS" id="121177" instrument="CTA"> | |
| 95 | - <parameter name="EventList" file="selected_gps_baseline_121177.fits" /> | |
| 94 | + <observation name="GPS" id="111177" instrument="CTA"> | |
| 95 | + <parameter name="EventList" file="selected_gps_baseline_111177.fits" /> | |
| 96 | 96 | <parameter name="Calibration" database="1dc" response="South_z20_50h" /> |
| 97 | 97 | </observation> |
| 98 | 98 | </observation_list> | ... | ... |
doc/source/users/tutorials/1dc/first_select_obs.rst
| ... | ... | @@ -67,32 +67,32 @@ observations into the console: |
| 67 | 67 | |
| 68 | 68 | The output of :ref:`csobsinfo` is shown below. :ref:`csobsselect` selected |
| 69 | 69 | 57 observations from the Galactic Plane Survey data containing a total of |
| 70 | -6 209 281 events. | |
| 70 | +6 046 100 events. | |
| 71 | 71 | |
| 72 | 72 | .. code-block:: bash |
| 73 | 73 | |
| 74 | - 2017-06-08T07:23:37: +=========+ | |
| 75 | - 2017-06-08T07:23:37: | Summary | | |
| 76 | - 2017-06-08T07:23:37: +=========+ | |
| 77 | - 2017-06-08T07:23:37: === Observations === | |
| 78 | - 2017-06-08T07:23:37: Unbinned observations .....: 57 | |
| 79 | - 2017-06-08T07:23:37: Binned observations .......: 0 | |
| 80 | - 2017-06-08T07:23:37: === Events === | |
| 81 | - 2017-06-08T07:23:37: Number of events ..........: 6209281 | |
| 82 | - 2017-06-08T07:23:37: Number of bins ............: 0 | |
| 83 | - 2017-06-08T07:23:37: === Pointings === | |
| 84 | - 2017-06-08T07:23:37: Mean offset angle .........: Unknown | |
| 85 | - 2017-06-08T07:23:37: Mean zenith angle .........: 0.00 deg | |
| 86 | - 2017-06-08T07:23:37: Mean azimuth angle ........: 0.00 deg | |
| 87 | - 2017-06-08T07:23:37: === Energy range === | |
| 88 | - 2017-06-08T07:23:37: Minimum energy ............: 30 GeV | |
| 89 | - 2017-06-08T07:23:37: Maximum energy ............: 160 TeV | |
| 90 | - 2017-06-08T07:23:37: === Time range === | |
| 91 | - 2017-06-08T07:23:37: MJD (days) ................: 59235.500 - 59276.921 | |
| 92 | - 2017-06-08T07:23:37: UTC .......................: 2021-01-21T11:58:51 - 2021-03-03T22:04:51 | |
| 93 | - 2017-06-08T07:23:37: MET (seconds) .............: 664502400.000 - 668081160.000 | |
| 94 | - 2017-06-08T07:23:37: Total ontime ..............: 102600.00 s = 1710.00 min = 28.50 h | |
| 95 | - 2017-06-08T07:23:37: Total livetime ............: 100548.00 s = 1675.80 min = 27.93 h | |
| 74 | + 2017-07-27T09:58:12: +=========+ | |
| 75 | + 2017-07-27T09:58:12: | Summary | | |
| 76 | + 2017-07-27T09:58:12: +=========+ | |
| 77 | + 2017-07-27T09:58:12: === Observations === | |
| 78 | + 2017-07-27T09:58:12: Unbinned observations .....: 57 | |
| 79 | + 2017-07-27T09:58:12: Binned observations .......: 0 | |
| 80 | + 2017-07-27T09:58:12: === Events === | |
| 81 | + 2017-07-27T09:58:12: Number of events ..........: 6046100 | |
| 82 | + 2017-07-27T09:58:12: Number of bins ............: 0 | |
| 83 | + 2017-07-27T09:58:12: === Pointings === | |
| 84 | + 2017-07-27T09:58:12: Mean offset angle .........: Unknown | |
| 85 | + 2017-07-27T09:58:12: Mean zenith angle .........: 0.00 deg | |
| 86 | + 2017-07-27T09:58:12: Mean azimuth angle ........: 0.00 deg | |
| 87 | + 2017-07-27T09:58:12: === Energy range === | |
| 88 | + 2017-07-27T09:58:12: Minimum energy ............: 30 GeV | |
| 89 | + 2017-07-27T09:58:12: Maximum energy ............: 160 TeV | |
| 90 | + 2017-07-27T09:58:12: === Time range === | |
| 91 | + 2017-07-27T09:58:12: MJD (days) ................: 59235.500 - 59276.921 | |
| 92 | + 2017-07-27T09:58:12: UTC .......................: 2021-01-21T11:58:51 - 2021-03-03T22:04:51 | |
| 93 | + 2017-07-27T09:58:12: MET (seconds) .............: 664502400.000 - 668081160.000 | |
| 94 | + 2017-07-27T09:58:12: Total ontime ..............: 102600.00 s = 1710.00 min = 28.50 h | |
| 95 | + 2017-07-27T09:58:12: Total livetime ............: 100548.00 s = 1675.80 min = 27.93 h | |
| 96 | 96 | |
| 97 | 97 | The resulting |
| 98 | 98 | :ref:`observation definition file <glossary_obsdef>` |
| ... | ... | @@ -102,17 +102,17 @@ will look as follows: |
| 102 | 102 | |
| 103 | 103 | <?xml version="1.0" encoding="UTF-8" standalone="no"?> |
| 104 | 104 | <observation_list title="observation list"> |
| 105 | - <observation name="GPS" id="120380" instrument="CTA"> | |
| 106 | - <parameter name="EventList" file="/Users/jurgen/analysis/cta/dc/1dc/1dc.final/validation/1dc.south/data/baseline/gps/gps_baseline_120380.fits" /> | |
| 105 | + <observation name="GPS" id="110380" instrument="CTA"> | |
| 106 | + <parameter name="EventList" file="/Users/jurgen/analysis/cta/dc/1dc/1dc/tutorial/1dc/data/baseline/gps/gps_baseline_110380.fits" /> | |
| 107 | 107 | <parameter name="Calibration" database="1dc" response="South_z20_50h" /> |
| 108 | 108 | </observation> |
| 109 | - <observation name="GPS" id="120381" instrument="CTA"> | |
| 110 | - <parameter name="EventList" file="/Users/jurgen/analysis/cta/dc/1dc/1dc.final/validation/1dc.south/data/baseline/gps/gps_baseline_120381.fits" /> | |
| 109 | + <observation name="GPS" id="110381" instrument="CTA"> | |
| 110 | + <parameter name="EventList" file="/Users/jurgen/analysis/cta/dc/1dc/1dc/tutorial/1dc/data/baseline/gps/gps_baseline_110381.fits" /> | |
| 111 | 111 | <parameter name="Calibration" database="1dc" response="South_z20_50h" /> |
| 112 | 112 | </observation> |
| 113 | 113 | ... |
| 114 | - <observation name="GPS" id="121177" instrument="CTA"> | |
| 115 | - <parameter name="EventList" file="/Users/jurgen/analysis/cta/dc/1dc/1dc.final/validation/1dc.south/data/baseline/gps/gps_baseline_121177.fits" /> | |
| 114 | + <observation name="GPS" id="111177" instrument="CTA"> | |
| 115 | + <parameter name="EventList" file="/Users/jurgen/analysis/cta/dc/1dc/1dc/tutorial/1dc/data/baseline/gps/gps_baseline_111177.fits" /> | |
| 116 | 116 | <parameter name="Calibration" database="1dc" response="South_z20_50h" /> |
| 117 | 117 | </observation> |
| 118 | 118 | </observation_list> | ... | ... |
doc/source/users/tutorials/1dc/first_skymap.png
doc/source/users/tutorials/1dc/first_skymap.rst
| ... | ... | @@ -44,7 +44,7 @@ All events between 100 GeV and 100 TeV are collected in the sky map. |
| 44 | 44 | The sky map is written into the FITS file ``skymap.fits`` that is created in |
| 45 | 45 | the working directory. |
| 46 | 46 | The sky map, displayed using |
| 47 | -`ds9 <http://ds9.si.edu>`_, | |
| 47 | +`ds9 <http://ds9.si.edu>`_ with a small Gaussian kernel smoothing applied, | |
| 48 | 48 | is shown below: |
| 49 | 49 | |
| 50 | 50 | .. figure:: first_skymap.png | ... | ... |
doc/source/users/tutorials/1dc/first_skymap_bkgsubtract.png
doc/source/users/tutorials/1dc/first_skymap_fitted.png
doc/source/users/tutorials/1dc/first_skymap_residual.png
doc/source/users/tutorials/1dc/first_skymap_residual_iem.png
doc/source/users/tutorials/1dc/first_skymap_sources.png
doc/source/users/tutorials/1dc/first_spectrum_cutoff_stacked.png
doc/source/users/tutorials/1dc/first_spectrum_stacked.png
doc/source/users/tutorials/1dc/first_stacked.rst
| ... | ... | @@ -150,8 +150,8 @@ that was generated by :ref:`ctbkgcube`. |
| 150 | 150 | <parameter name="PivotEnergy" value="1" scale="300000" free="0" /> |
| 151 | 151 | </spectrum> |
| 152 | 152 | <spatialModel type="PointSource"> |
| 153 | - <parameter name="RA" value="266.4120906928" error="0" scale="1" free="1" /> | |
| 154 | - <parameter name="DEC" value="-29.0219729468991" error="0" scale="1" free="1" /> | |
| 153 | + <parameter name="RA" value="266.424004498437" error="0" scale="1" free="1" /> | |
| 154 | + <parameter name="DEC" value="-29.0049010253548" error="0" scale="1" free="1" /> | |
| 155 | 155 | </spatialModel> |
| 156 | 156 | </source> |
| 157 | 157 | <source name="Src002" type="PointSource"> |
| ... | ... | @@ -161,8 +161,8 @@ that was generated by :ref:`ctbkgcube`. |
| 161 | 161 | <parameter name="PivotEnergy" value="1" scale="300000" free="0" /> |
| 162 | 162 | </spectrum> |
| 163 | 163 | <spatialModel type="PointSource"> |
| 164 | - <parameter name="RA" value="266.793148201606" error="0" scale="1" free="1" /> | |
| 165 | - <parameter name="DEC" value="-28.1253038992376" error="0" scale="1" free="1" /> | |
| 164 | + <parameter name="RA" value="266.831945177213" error="0" scale="1" free="1" /> | |
| 165 | + <parameter name="DEC" value="-28.1460284439951" error="0" scale="1" free="1" /> | |
| 166 | 166 | </spatialModel> |
| 167 | 167 | </source> |
| 168 | 168 | <source name="BackgroundModel" type="CTACubeBackground"> | ... | ... |