Domenico Tiziani / ctools

Commit b107f9bdb85de512dfbdfad7a9abce3b05fd5b15

Authored by Jürgen Knödlseder
1 parent 2edf36eb

Update with latest default values

Inline Side-by-side
Showing 1 changed file with 230 additions and 221 deletions
doc/source/user_manual/getting_started/quickstart.rst
  View file @b107f9b
... ... @@ -34,17 +34,17 @@ will query for a number of parameters:
34 34 .. code-block:: bash
35 35  
36 36 $ ctobssim
37   - Model [test/data/crab.xml] $GAMMALIB/share/models/crab.xml
38   - Calibration database [test/irf] $GAMMALIB/share/caldb/cta
39   - Instrument response function [cta_dummy_irf]
  37 + Model [$CTOOLS/share/models/crab.xml]
40 38 RA of pointing (degrees) (0-360) [83.63]
41 39 Dec of pointing (degrees) (-90-90) [22.01]
42   - Radius of FOV (degrees) (0-180) [10.0]
  40 + Calibration database [dummy]
  41 + Instrument response function [cta_dummy_irf]
  42 + Radius of FOV (degrees) (0-180) [5.0]
43 43 Start time (MET in s) (0) [0.0]
44 44 End time (MET in s) (0) [1800.0]
45 45 Lower energy limit (TeV) (0) [0.1]
46 46 Upper energy limit (TeV) (0) [100.0]
47   - Output event data file or observation definition file [events.fits]
  47 + Output event data file or observation definition file [events.fits]
48 48  
49 49 Each line represents a query for one parameter. The line starts with a short
50 50 description of the parameter, followed by the default parameter value proposed
... ... @@ -53,13 +53,13 @@ is the case for the majority of parameters shown here), the default parameter
53 53 will be used. Otherwise, the specified parameter will overwrite the
54 54 default parameter.
55 55  
56   -You may have recognised that the environment variable ``$GAMMALIB`` has
  56 +You may have recognised that the environment variable ``$CTOOLS`` has
57 57 been used in the path names of the first two parameters. ctools will
58 58 automatically expand the environment variables.
59 59  
60 60 The CTA instrument properties (effective area, PSF width) are taken for
61   -the moment from Konrad Bernloehrs sensitivity files which are stored in
62   -``$GAMMALIB/share/caldb/cta``.
  61 +the moment from a dummy performance table that is located in
  62 +``$CTOOLS/share/caldb/cta``.
63 63  
64 64 Events are simulated based on the instrument properties and based on a
65 65 source and background model. Only events that fall within the specified
... ... @@ -69,13 +69,13 @@ event data file. The duration of the simulation is taken here to 30 minutes
69 69 (or 1800 seconds). Events are simulated for energies between 0.1 and 100 TeV.
70 70  
71 71 The source and background model is defined by the XML file
72   -``$GAMMALIB/share/models/crab.xml``:
  72 +``$CTOOLS/share/models/crab.xml``:
73 73  
74 74 .. code-block:: xml
75 75  
76 76 <?xml version="1.0" standalone="no"?>
77 77 <source_library title="source library">
78   - <source name="MAGIC Crab (Albert et al. 2008, ApJ, 674, 1037)" type="PointSource">
  78 + <source name="Crab" type="PointSource">
79 79 <spectrum type="PowerLaw">
80 80 <parameter name="Prefactor" scale="1e-16" value="5.7" min="1e-07" max="1000.0" free="1"/>
81 81 <parameter name="Index" scale="-1" value="2.48" min="0.0" max="+5.0" free="1"/>
... ... @@ -87,7 +87,7 @@ The source and background model is defined by the XML file
87 87 </spatialModel>
88 88 </source>
89 89 <source name="Background" type="RadialAcceptance" instrument="CTA">
90   - <spectrum type="FileFunction" file="$GAMMALIB/share/models/bkg_dummy.txt">
  90 + <spectrum type="FileFunction" file="$CTOOLS/share/models/bkg_dummy.txt">
91 91 <parameter scale="1.0" name="Normalization" min="0.0" max="1000.0" value="1.0" free="1"/>
92 92 </spectrum>
93 93 <radialModel type="Gaussian">
... ... @@ -161,24 +161,26 @@ from this file are shown here:
161 161  
162 162 .. code-block:: xml
163 163  
164   - 2014-01-09T20:51:36: +======================+
165   - 2014-01-09T20:51:36: | Simulate observation |
166   - 2014-01-09T20:51:36: +======================+
167   - 2014-01-09T20:51:36: === Observation ===
168   - 2014-01-09T20:51:36: Simulation area ...........: 1.9635e+11 cm2
169   - 2014-01-09T20:51:36: Simulation cone ...........: RA=83.63 deg, Dec=22.01 deg, r=10.5 deg
170   - 2014-01-09T20:51:36: Time interval .............: 0 - 1800 s
171   - 2014-01-09T20:51:36: Energy range ..............: 100 GeV - 100 TeV
172   - 2014-01-09T20:51:36: MC source photons .........: 207547 [Crab]
173   - 2014-01-09T20:51:36: MC source events ..........: 995 [Crab]
174   - 2014-01-09T20:51:36: MC source events ..........: 995 (all source models)
175   - 2014-01-09T20:51:36: MC background events ......: 5146
176   - 2014-01-09T20:51:36:
177   - 2014-01-09T20:51:36: +==================+
178   - 2014-01-09T20:51:36: | Save observation |
179   - 2014-01-09T20:51:36: +==================+
180   - 2014-01-09T20:51:37:
181   - 2014-01-09T20:51:37: Application "ctobssim" terminated after 24 wall clock seconds, consuming 0.307116 seconds of CPU time.
  164 + 2014-10-30T22:35:06: +======================+
  165 + 2014-10-30T22:35:06: | Simulate observation |
  166 + 2014-10-30T22:35:06: +======================+
  167 + 2014-10-30T22:35:06: === Observation ===
  168 + 2014-10-30T22:35:06: Simulation area ...........: 1.9635e+11 cm2
  169 + 2014-10-30T22:35:06: Simulation cone ...........: RA=83.63 deg, Dec=22.01 deg, r=5.5 deg
  170 + 2014-10-30T22:35:06: Time interval .............: 0 - 1800 s
  171 + 2014-10-30T22:35:06: Photon energy range .......: 100 GeV - 100 TeV
  172 + 2014-10-30T22:35:06: Event energy range ........: 100 GeV - 100 TeV
  173 + 2014-10-30T22:35:06: MC source photons .........: 207547 [Crab]
  174 + 2014-10-30T22:35:06: MC source events ..........: 995 [Crab]
  175 + 2014-10-30T22:35:06: MC source events ..........: 995 (all source models)
  176 + 2014-10-30T22:35:06: MC background events ......: 5146
  177 + 2014-10-30T22:35:06: MC events .................: 6141 (all models)
  178 + 2014-10-30T22:35:06:
  179 + 2014-10-30T22:35:06: +==================+
  180 + 2014-10-30T22:35:06: | Save observation |
  181 + 2014-10-30T22:35:06: +==================+
  182 + 2014-10-30T22:35:06:
  183 + 2014-10-30T22:35:06: Application "ctobssim" terminated after 10 wall clock seconds, consuming 0.3604 seconds of CPU time.
182 184  
183 185 Each line starts with the UTC time at which the line has been written. In
184 186 this run, 207547 Crab photons have been thrown over an area of 19.6 square
... ... @@ -224,7 +226,7 @@ ctbin is executed by typing:
224 226 .. code-block:: bash
225 227  
226 228 $ ctbin
227   - Input event list or observation definition file [test/data/crab_events.fits.gz] events.fits
  229 + Input event list or observation definition file [events.fits]
228 230 First coordinate of image center in degrees (RA or galactic l) [83.63]
229 231 Second coordinate of image center in degrees (DEC or galactic b) [22.01]
230 232 Projection method e.g. AIT|AZP|CAR|MER|STG|TAN (AIT|AZP|CAR|MER|STG|TAN) [CAR]
... ... @@ -277,40 +279,37 @@ reproduced below:
277 279  
278 280 .. code-block:: xml
279 281  
280   - 2014-01-09T21:16:18: +====================+
281   - 2014-01-09T21:16:18: | Binned observation |
282   - 2014-01-09T21:16:18: +====================+
283   - 2014-01-09T21:16:18: === GObservations ===
284   - 2014-01-09T21:16:18: Number of observations ....: 1
285   - 2014-01-09T21:16:18: Number of predicted events : 0
286   - 2014-01-09T21:16:18: === GCTAObservation ===
287   - 2014-01-09T21:16:18: Name ......................:
288   - 2014-01-09T21:16:18: Identifier ................:
289   - 2014-01-09T21:16:18: Instrument ................: CTA
290   - 2014-01-09T21:16:18: Statistics ................: Poisson
291   - 2014-01-09T21:16:18: Ontime ....................: 1800 s
292   - 2014-01-09T21:16:18: Livetime ..................: 1710 s
293   - 2014-01-09T21:16:18: Deadtime correction .......: 0.95
294   - 2014-01-09T21:16:18: === GCTAPointing ===
295   - 2014-01-09T21:16:18: Pointing direction ........:
296   - 2014-01-09T21:16:18: === GCTAResponse ===
297   - 2014-01-09T21:16:18: Calibration database ......:
298   - 2014-01-09T21:16:18: Response name .............:
299   - 2014-01-09T21:16:18: RMF file name .............:
300   - 2014-01-09T21:16:18: === GCTAEventCube ===
301   - 2014-01-09T21:16:18: Number of events ..........: 5542
302   - 2014-01-09T21:16:18: Number of elements ........: 800000
303   - 2014-01-09T21:16:18: Number of pixels ..........: 40000
304   - 2014-01-09T21:16:18: Number of energy bins .....: 20
305   - 2014-01-09T21:16:18: Time interval .............: -3.15576e+08 - -3.15574e+08 sec
306   - 2014-01-09T21:16:18:
307   - 2014-01-09T21:16:18:
308   - 2014-01-09T21:16:18:
309   - 2014-01-09T21:16:18: +==================+
310   - 2014-01-09T21:16:18: | Save observation |
311   - 2014-01-09T21:16:18: +==================+
312   - 2014-01-09T21:16:18:
313   - 2014-01-09T21:16:18: Application "ctbin" terminated after 37 wall clock seconds, consuming 0.242495 seconds of CPU time.
  282 + 2014-10-30T22:39:35: +====================+
  283 + 2014-10-30T22:39:35: | Binned observation |
  284 + 2014-10-30T22:39:35: +====================+
  285 + 2014-10-30T22:39:35: === GObservations ===
  286 + 2014-10-30T22:39:35: Number of observations ....: 1
  287 + 2014-10-30T22:39:35: Number of predicted events : 0
  288 + 2014-10-30T22:39:35: === GCTAObservation ===
  289 + 2014-10-30T22:39:35: Name ......................:
  290 + 2014-10-30T22:39:35: Identifier ................:
  291 + 2014-10-30T22:39:35: Instrument ................: CTA
  292 + 2014-10-30T22:39:35: Event file ................: events.fits
  293 + 2014-10-30T22:39:35: Event type ................: CountsMap
  294 + 2014-10-30T22:39:35: Statistics ................: Poisson
  295 + 2014-10-30T22:39:35: Ontime ....................: 1800 s
  296 + 2014-10-30T22:39:35: Livetime ..................: 1710 s
  297 + 2014-10-30T22:39:35: Deadtime correction .......: 0.95
  298 + 2014-10-30T22:39:35: === GCTAPointing ===
  299 + 2014-10-30T22:39:35: Pointing direction ........: (RA,Dec)=(83.63,22.01)
  300 + 2014-10-30T22:39:35: Response function .........: undefined
  301 + 2014-10-30T22:39:35: === GCTAEventCube ===
  302 + 2014-10-30T22:39:35: Number of events ..........: 5542
  303 + 2014-10-30T22:39:35: Number of elements ........: 800000
  304 + 2014-10-30T22:39:35: Number of pixels ..........: 40000
  305 + 2014-10-30T22:39:35: Number of energy bins .....: 20
  306 + 2014-10-30T22:39:35: Time interval .............: 0 - 1800 sec
  307 + 2014-10-30T22:39:35:
  308 + 2014-10-30T22:39:35: +==================+
  309 + 2014-10-30T22:39:35: | Save observation |
  310 + 2014-10-30T22:39:35: +==================+
  311 + 2014-10-30T22:39:36:
  312 + 2014-10-30T22:39:36: Application "ctbin" terminated after 18 wall clock seconds, consuming 0.308684 seconds of CPU time.
314 313  
315 314 From the 6141 events that have been simulated and stored in the
316 315 ``events.fits`` file, 5542 lie within the cube boundaries and are thus put
... ... @@ -331,10 +330,10 @@ do the fit by typing:
331 330 .. code-block:: bash
332 331  
333 332 $ ctlike
334   - Event list, counts cube or observation definition file [test/data/crab_events.fits.gz] cntmap.fits
335   - Calibration database [test/irf] $GAMMALIB/share/caldb/cta
  333 + Event list, counts cube or observation definition file [events.fits] cntmap.fits
  334 + Calibration database [dummy]
336 335 Instrument response function [cta_dummy_irf]
337   - Source model [test/data/crab.xml] $GAMMALIB/share/models/crab.xml
  336 + Source model [$CTOOLS/share/models/crab.xml]
338 337 Source model output file [crab_results.xml]
339 338  
340 339 Fitting of the data is done in *binned* mode, which means that the events
... ... @@ -357,8 +356,8 @@ result file that has been produced by the run:
357 356 <source_library title="source library">
358 357 <source name="Crab" type="PointSource">
359 358 <spectrum type="PowerLaw">
360   - <parameter name="Prefactor" value="6.07927" error="0.204581" scale="1e-16" min="1e-07" max="1000" free="1" />
361   - <parameter name="Index" value="2.5009" error="0.0252056" scale="-1" min="0" max="5" free="1" />
  359 + <parameter name="Prefactor" value="6.07928" error="0.204582" scale="1e-16" min="1e-07" max="1000" free="1" />
  360 + <parameter name="Index" value="2.5009" error="0.0252057" scale="-1" min="0" max="5" free="1" />
362 361 <parameter name="Scale" value="0.3" scale="1e+06" min="0.01" max="1000" free="0" />
363 362 </spectrum>
364 363 <spatialModel type="SkyDirFunction">
... ... @@ -367,11 +366,11 @@ result file that has been produced by the run:
367 366 </spatialModel>
368 367 </source>
369 368 <source name="Background" type="RadialAcceptance" instrument="CTA">
370   - <spectrum type="FileFunction" file="$GAMMALIB/share/models/bkg_dummy.txt">
371   - <parameter name="Normalization" value="0.99069" error="0.0192887" scale="1" min="0" max="1000" free="1" />
  369 + <spectrum type="FileFunction" file="$CTOOLS/share/models/bkg_dummy.txt">
  370 + <parameter name="Normalization" value="0.990708" error="0.0192887" scale="1" min="0" max="1000" free="1" />
372 371 </spectrum>
373 372 <radialModel type="Gaussian">
374   - <parameter name="Sigma" value="3.02527" error="0.0594142" scale="1" min="0.01" max="10" free="1" />
  373 + <parameter name="Sigma" value="3.02523" error="0.0594117" scale="1" min="0.01" max="10" free="1" />
375 374 </radialModel>
376 375 </source>
377 376 </source_library>
... ... @@ -387,74 +386,81 @@ this run:
387 386  
388 387 .. code-block:: xml
389 388  
390   - 2014-01-09T21:26:56: +=================================+
391   - 2014-01-09T21:26:56: | Maximum likelihood optimisation |
392   - 2014-01-09T21:26:56: +=================================+
393   - 2014-01-09T21:26:57: Initial iteration: func=22752.5, Lambda=0.001
394   - 2014-01-09T21:26:59: Iteration 1: func=22750.6, Lambda=0.0001, delta=1.91297, max(grad)=0 [-1]
395   - 2014-01-09T21:27:00: Iteration 2: func=22750.6, Lambda=1e-05, delta=0.00782419, max(grad)=0.000463786 [6]
396   - 2014-01-09T21:27:01: Iteration 3: func=22750.6, Lambda=1e-06, delta=2.49478e-06, max(grad)=0 [-1]
397   - 2014-01-09T21:27:03: Iteration 4: func=22750.6, Lambda=1e-05, delta=0, max(grad)=7.93613e-08 [6] (stalled)
398   - 2014-01-09T21:27:04: Iteration 5: func=22750.6, Lambda=0.0001, delta=0, max(grad)=7.91125e-08 [6] (stalled)
399   - 2014-01-09T21:27:05: Iteration 6: func=22750.6, Lambda=0.001, delta=0, max(grad)=7.66376e-08 [6] (stalled)
400   - 2014-01-09T21:27:07: Iteration 7: func=22750.6, Lambda=0.01, delta=0, max(grad)=1.28271e-07 [7] (stalled)
401   - 2014-01-09T21:27:08: Iteration 8: func=22750.6, Lambda=0.001, delta=2.98678e-09, max(grad)=1.23881e-06 [7]
402   - 2014-01-09T21:27:09:
403   - 2014-01-09T21:27:09: +=========================================+
404   - 2014-01-09T21:27:09: | Maximum likelihood optimization results |
405   - 2014-01-09T21:27:09: +=========================================+
406   - 2014-01-09T21:27:09: === GOptimizerLM ===
407   - 2014-01-09T21:27:09: Optimized function value ..: 22750.6
408   - 2014-01-09T21:27:09: Absolute precision ........: 1e-06
409   - 2014-01-09T21:27:09: Optimization status .......: converged
410   - 2014-01-09T21:27:09: Number of parameters ......: 9
411   - 2014-01-09T21:27:09: Number of free parameters .: 4
412   - 2014-01-09T21:27:09: Number of iterations ......: 8
413   - 2014-01-09T21:27:09: Lambda ....................: 0.001
414   - 2014-01-09T21:27:09: Maximum log likelihood ....: -22750.6
415   - 2014-01-09T21:27:09: Observed events (Nobs) ...: 5542
416   - 2014-01-09T21:27:09: Predicted events (Npred) ..: 5542 (Nobs - Npred = 5.52808e-06)
417   - 2014-01-09T21:27:09: === GModels ===
418   - 2014-01-09T21:27:09: Number of models ..........: 2
419   - 2014-01-09T21:27:09: Number of parameters ......: 9
420   - 2014-01-09T21:27:09: === GModelSky ===
421   - 2014-01-09T21:27:09: Name ......................: Crab
422   - 2014-01-09T21:27:09: Instruments ...............: all
423   - 2014-01-09T21:27:09: Instrument scale factors ..: unity
424   - 2014-01-09T21:27:09: Observation identifiers ...: all
425   - 2014-01-09T21:27:09: Model type ................: PointSource
426   - 2014-01-09T21:27:09: Model components ..........: "SkyDirFunction" * "PowerLaw" * "Constant"
427   - 2014-01-09T21:27:09: Number of parameters ......: 6
428   - 2014-01-09T21:27:09: Number of spatial par's ...: 2
429   - 2014-01-09T21:27:09: RA .......................: 83.6331 [-360,360] deg (fixed,scale=1)
430   - 2014-01-09T21:27:09: DEC ......................: 22.0145 [-90,90] deg (fixed,scale=1)
431   - 2014-01-09T21:27:09: Number of spectral par's ..: 3
432   - 2014-01-09T21:27:09: Prefactor ................: 6.07927e-16 +/- 2.04581e-17 [1e-23,1e-13] ph/cm2/s/MeV (free,scale=1e-16,gradient)
433   - 2014-01-09T21:27:09: Index ....................: -2.5009 +/- 0.0252056 [-0,-5] (free,scale=-1,gradient)
434   - 2014-01-09T21:27:09: PivotEnergy ..............: 300000 [10000,1e+09] MeV (fixed,scale=1e+06,gradient)
435   - 2014-01-09T21:27:09: Number of temporal par's ..: 1
436   - 2014-01-09T21:27:09: Constant .................: 1 (relative value) (fixed,scale=1,gradient)
437   - 2014-01-09T21:27:09: === GCTAModelRadialAcceptance ===
438   - 2014-01-09T21:27:09: Name ......................: Background
439   - 2014-01-09T21:27:09: Instruments ...............: CTA
440   - 2014-01-09T21:27:09: Instrument scale factors ..: unity
441   - 2014-01-09T21:27:09: Observation identifiers ...: all
442   - 2014-01-09T21:27:09: Model type ................: "Gaussian" * "FileFunction" * "Constant"
443   - 2014-01-09T21:27:09: Number of parameters ......: 3
444   - 2014-01-09T21:27:09: Number of radial par's ....: 1
445   - 2014-01-09T21:27:09: Sigma ....................: 3.02527 +/- 0.0594142 [0.01,10] deg2 (free,scale=1,gradient)
446   - 2014-01-09T21:27:09: Number of spectral par's ..: 1
447   - 2014-01-09T21:27:09: Normalization ............: 0.99069 +/- 0.0192887 [0,1000] (free,scale=1,gradient)
448   - 2014-01-09T21:27:09: Number of temporal par's ..: 1
449   - 2014-01-09T21:27:09: Constant .................: 1 (relative value) (fixed,scale=1,gradient)
450   - 2014-01-09T21:27:09:
451   - 2014-01-09T21:27:09: +==============+
452   - 2014-01-09T21:27:09: | Save results |
453   - 2014-01-09T21:27:09: +==============+
454   - 2014-01-09T21:27:09:
455   - 2014-01-09T21:27:09: Application "ctlike" terminated after 55 wall clock seconds, consuming 13.4985 seconds of CPU time.
456   -
457   -The maximum likelihood optimizer required 8 iterations to converge. This
  389 + 2014-10-30T22:41:08: +=================================+
  390 + 2014-10-30T22:41:08: | Maximum likelihood optimisation |
  391 + 2014-10-30T22:41:08: +=================================+
  392 + 2014-10-30T22:41:09: >Iteration 0: -logL=22752.553, Lambda=1.0e-03
  393 + 2014-10-30T22:41:11: >Iteration 1: -logL=22750.640, Lambda=1.0e-03, delta=1.913, max(|grad|)=-0.941594 [Index:3]
  394 + 2014-10-30T22:41:12: >Iteration 2: -logL=22750.632, Lambda=1.0e-04, delta=0.008, max(|grad|)=-0.079270 [Index:3]
  395 + 2014-10-30T22:41:14: >Iteration 3: -logL=22750.632, Lambda=1.0e-05, delta=0.000, max(|grad|)=-0.003073 [Index:3]
  396 + 2014-10-30T22:41:16:
  397 + 2014-10-30T22:41:16: +=========================================+
  398 + 2014-10-30T22:41:16: | Maximum likelihood optimization results |
  399 + 2014-10-30T22:41:16: +=========================================+
  400 + 2014-10-30T22:41:16: === GOptimizerLM ===
  401 + 2014-10-30T22:41:16: Optimized function value ..: 22750.632
  402 + 2014-10-30T22:41:16: Absolute precision ........: 0.005
  403 + 2014-10-30T22:41:16: Optimization status .......: converged
  404 + 2014-10-30T22:41:16: Number of parameters ......: 9
  405 + 2014-10-30T22:41:16: Number of free parameters .: 4
  406 + 2014-10-30T22:41:16: Number of iterations ......: 3
  407 + 2014-10-30T22:41:16: Lambda ....................: 1e-06
  408 + 2014-10-30T22:41:16:
  409 + 2014-10-30T22:41:16: +=========================================+
  410 + 2014-10-30T22:41:16: | Maximum likelihood optimization results |
  411 + 2014-10-30T22:41:16: +=========================================+
  412 + 2014-10-30T22:41:16: === GOptimizerLM ===
  413 + 2014-10-30T22:41:16: Optimized function value ..: 22750.632
  414 + 2014-10-30T22:41:16: Absolute precision ........: 0.005
  415 + 2014-10-30T22:41:16: Optimization status .......: converged
  416 + 2014-10-30T22:41:16: Number of parameters ......: 9
  417 + 2014-10-30T22:41:16: Number of free parameters .: 4
  418 + 2014-10-30T22:41:16: Number of iterations ......: 3
  419 + 2014-10-30T22:41:16: Lambda ....................: 1e-06
  420 + 2014-10-30T22:41:16: Maximum log likelihood ....: -22750.632
  421 + 2014-10-30T22:41:16: Observed events (Nobs) ...: 5542.000
  422 + 2014-10-30T22:41:16: Predicted events (Npred) ..: 5542.000 (Nobs - Npred = 3.3498e-06)
  423 + 2014-10-30T22:41:16: === GModels ===
  424 + 2014-10-30T22:41:16: Number of models ..........: 2
  425 + 2014-10-30T22:41:16: Number of parameters ......: 9
  426 + 2014-10-30T22:41:16: === GModelSky ===
  427 + 2014-10-30T22:41:16: Name ......................: Crab
  428 + 2014-10-30T22:41:16: Instruments ...............: all
  429 + 2014-10-30T22:41:16: Instrument scale factors ..: unity
  430 + 2014-10-30T22:41:16: Observation identifiers ...: all
  431 + 2014-10-30T22:41:16: Model type ................: PointSource
  432 + 2014-10-30T22:41:16: Model components ..........: "SkyDirFunction" * "PowerLaw" * "Constant"
  433 + 2014-10-30T22:41:16: Number of parameters ......: 6
  434 + 2014-10-30T22:41:16: Number of spatial par's ...: 2
  435 + 2014-10-30T22:41:16: RA .......................: 83.6331 [-360,360] deg (fixed,scale=1)
  436 + 2014-10-30T22:41:16: DEC ......................: 22.0145 [-90,90] deg (fixed,scale=1)
  437 + 2014-10-30T22:41:16: Number of spectral par's ..: 3
  438 + 2014-10-30T22:41:16: Prefactor ................: 6.07928e-16 +/- 2.04582e-17 [1e-23,1e-13] ph/cm2/s/MeV (free,scale=1e-16,gradient)
  439 + 2014-10-30T22:41:16: Index ....................: -2.5009 +/- 0.0252057 [-0,-5] (free,scale=-1,gradient)
  440 + 2014-10-30T22:41:16: PivotEnergy ..............: 300000 [10000,1e+09] MeV (fixed,scale=1e+06,gradient)
  441 + 2014-10-30T22:41:16: Number of temporal par's ..: 1
  442 + 2014-10-30T22:41:16: Constant .................: 1 (relative value) (fixed,scale=1,gradient)
  443 + 2014-10-30T22:41:16: === GCTAModelRadialAcceptance ===
  444 + 2014-10-30T22:41:16: Name ......................: Background
  445 + 2014-10-30T22:41:16: Instruments ...............: CTA
  446 + 2014-10-30T22:41:16: Instrument scale factors ..: unity
  447 + 2014-10-30T22:41:16: Observation identifiers ...: all
  448 + 2014-10-30T22:41:16: Model type ................: "Gaussian" * "FileFunction" * "Constant"
  449 + 2014-10-30T22:41:16: Number of parameters ......: 3
  450 + 2014-10-30T22:41:16: Number of radial par's ....: 1
  451 + 2014-10-30T22:41:16: Sigma ....................: 3.02523 +/- 0.0594117 [0.01,10] deg2 (free,scale=1,gradient)
  452 + 2014-10-30T22:41:16: Number of spectral par's ..: 1
  453 + 2014-10-30T22:41:16: Normalization ............: 0.990708 +/- 0.0192887 [0,1000] (free,scale=1,gradient)
  454 + 2014-10-30T22:41:16: Number of temporal par's ..: 1
  455 + 2014-10-30T22:41:16: Constant .................: 1 (relative value) (fixed,scale=1,gradient)
  456 + 2014-10-30T22:41:16:
  457 + 2014-10-30T22:41:16: +==============+
  458 + 2014-10-30T22:41:16: | Save results |
  459 + 2014-10-30T22:41:16: +==============+
  460 + 2014-10-30T22:41:16:
  461 + 2014-10-30T22:41:16: Application "ctlike" terminated after 22 wall clock seconds, consuming 8.54182 seconds of CPU time.
  462 +
  463 +The maximum likelihood optimizer required 3 iterations to converge. This
458 464 is pretty fast, but recall that we used the same model file for the simulation
459 465 and for fitting, hence the initial parameter values were already very close
460 466 to the best fitting values. To see the impact of the initial parameters on
... ... @@ -498,7 +504,7 @@ following example shows how to run ctselect:
498 504 .. code-block:: bash
499 505  
500 506 $ ctselect
501   - Input event list or observation definition file [test/data/crab_events.fits.gz] events.fits
  507 + Input event list or observation definition file [events.fits]
502 508 RA for ROI centre (degrees) (0-360) [83.63]
503 509 Dec for ROI centre (degrees) (-90-90) [22.01]
504 510 Radius of ROI (degrees) (0-180) [3.0]
... ... @@ -532,15 +538,18 @@ part:
532 538  
533 539 .. code-block:: xml
534 540  
535   - 2014-01-09T21:36:20: +=================+
536   - 2014-01-09T21:36:20: | Event selection |
537   - 2014-01-09T21:36:20: +=================+
538   - 2014-01-09T21:36:20: === Observation ===
539   - 2014-01-09T21:36:21: Energy range ..............: 0.1 - 100 TeV
540   - 2014-01-09T21:36:21: Acceptance cone centre ....: RA=83.63, DEC=22.01 deg
541   - 2014-01-09T21:36:21: Acceptance cone radius ....: 3 deg
542   - 2014-01-09T21:36:21: cfitsio selection .........: ENERGY >= 0.10000000 && ENERGY <= 100.00000000 && ANGSEP(83.630000,22.010000,RA,DEC) <= 3.000000
543   - 2014-01-09T21:36:21: FITS filename .............: ctselectKKHEiF[EVENTS][ENERGY >= 0.10000000 && ENERGY <= 100.00000000 && ANGSEP(83.630000,22.010000,RA,DEC) <= 3.000000]
  541 + 2014-10-30T22:44:01: +=================+
  542 + 2014-10-30T22:44:01: | Event selection |
  543 + 2014-10-30T22:44:01: +=================+
  544 + 2014-10-30T22:44:01: === Observation ===
  545 + 2014-10-30T22:44:02: Requested energy range ....: 0.1 - 100 TeV
  546 + 2014-10-30T22:44:02: Energy range of data ......: 0.1 - 100 TeV
  547 + 2014-10-30T22:44:02: Selected energy range .....: 0.1 - 100 TeV
  548 + 2014-10-30T22:44:02: Requested ROI .............: Centre(RA,DEC)=(83.63, 22.01) deg, Radius=3 deg
  549 + 2014-10-30T22:44:02: ROI of data ...............: Centre(RA,DEC)=(83.63, 22.01) deg, Radius=5 deg
  550 + 2014-10-30T22:44:02: Selected ROI ..............: Centre(RA,DEC)=(83.63, 22.01) deg, Radius=3 deg
  551 + 2014-10-30T22:44:02: cfitsio selection .........: ENERGY >= 0.10000000 && ENERGY <= 100.00000000 && ANGSEP(83.630000,22.010000,RA,DEC) <= 3.000000
  552 + 2014-10-30T22:44:02: FITS filename .............: /var/tmp/tmp.0.VSaNXN[EVENTS][ENERGY >= 0.10000000 && ENERGY <= 100.00000000 && ANGSEP(83.630000,22.010000,RA,DEC) <= 3.000000]
544 553  
545 554 Now that you have selected the events of interest, you can run ctlike in
546 555 unbinned mode. To do this you have to specify the selected event list
... ... @@ -550,9 +559,9 @@ instead of the counts cube:
550 559  
551 560 $ ctlike
552 561 Event list, counts cube or observation definition file [cntmap.fits] selected_events.fits
553   - Calibration database [$GAMMALIB/share/caldb/cta]
  562 + Calibration database [dummy]
554 563 Instrument response function [cta_dummy_irf]
555   - Source model [$GAMMALIB/share/models/crab.xml]
  564 + Source model [$CTOOLS/share/models/crab.xml]
556 565 Source model output file [crab_results.xml]
557 566  
558 567 You will recognise that ctlike runs much faster in unbinned mode compared
... ... @@ -578,68 +587,67 @@ the same events for both analyses.
578 587  
579 588 .. code-block:: xml
580 589  
581   - 2014-01-09T21:47:54: +=================================+
582   - 2014-01-09T21:47:54: | Maximum likelihood optimisation |
583   - 2014-01-09T21:47:54: +=================================+
584   - 2014-01-09T21:47:54: Initial iteration: func=44393.9, Lambda=0.001
585   - 2014-01-09T21:47:54: Iteration 1: func=44392.2, Lambda=0.001, delta=1.7034, max(grad)=9.60142 [6]
586   - 2014-01-09T21:47:54: Iteration 2: func=44392.1, Lambda=0.0001, delta=0.099482, max(grad)=0 [-1]
587   - 2014-01-09T21:47:54: Iteration 3: func=44392.1, Lambda=1e-05, delta=6.63429e-05, max(grad)=0.00559314 [6]
588   - 2014-01-09T21:47:54: Iteration 4: func=44392.1, Lambda=1e-06, delta=4.7803e-08, max(grad)=0 [-1]
589   - 2014-01-09T21:47:54:
590   - 2014-01-09T21:47:54: +=========================================+
591   - 2014-01-09T21:47:54: | Maximum likelihood optimization results |
592   - 2014-01-09T21:47:54: +=========================================+
593   - 2014-01-09T21:47:54: === GOptimizerLM ===
594   - 2014-01-09T21:47:54: Optimized function value ..: 44392.1
595   - 2014-01-09T21:47:54: Absolute precision ........: 1e-06
596   - 2014-01-09T21:47:54: Optimization status .......: converged
597   - 2014-01-09T21:47:54: Number of parameters ......: 9
598   - 2014-01-09T21:47:54: Number of free parameters .: 4
599   - 2014-01-09T21:47:54: Number of iterations ......: 4
600   - 2014-01-09T21:47:54: Lambda ....................: 1e-07
601   - 2014-01-09T21:47:54: Maximum log likelihood ....: -44392.1
602   - 2014-01-09T21:47:54: Observed events (Nobs) ...: 6127
603   - 2014-01-09T21:47:54: Predicted events (Npred) ..: 6127 (Nobs - Npred = 1.01613e-07)
604   - 2014-01-09T21:47:54: === GModels ===
605   - 2014-01-09T21:47:54: Number of models ..........: 2
606   - 2014-01-09T21:47:54: Number of parameters ......: 9
607   - 2014-01-09T21:47:54: === GModelSky ===
608   - 2014-01-09T21:47:54: Name ......................: Crab
609   - 2014-01-09T21:47:54: Instruments ...............: all
610   - 2014-01-09T21:47:54: Instrument scale factors ..: unity
611   - 2014-01-09T21:47:54: Observation identifiers ...: all
612   - 2014-01-09T21:47:54: Model type ................: PointSource
613   - 2014-01-09T21:47:54: Model components ..........: "SkyDirFunction" * "PowerLaw" * "Constant"
614   - 2014-01-09T21:47:54: Number of parameters ......: 6
615   - 2014-01-09T21:47:54: Number of spatial par's ...: 2
616   - 2014-01-09T21:47:54: RA .......................: 83.6331 [-360,360] deg (fixed,scale=1)
617   - 2014-01-09T21:47:54: DEC ......................: 22.0145 [-90,90] deg (fixed,scale=1)
618   - 2014-01-09T21:47:54: Number of spectral par's ..: 3
619   - 2014-01-09T21:47:54: Prefactor ................: 6.13416e-16 +/- 2.05772e-17 [1e-23,1e-13] ph/cm2/s/MeV (free,scale=1e-16,gradient)
620   - 2014-01-09T21:47:54: Index ....................: -2.50585 +/- 0.0250821 [-0,-5] (free,scale=-1,gradient)
621   - 2014-01-09T21:47:54: PivotEnergy ..............: 300000 [10000,1e+09] MeV (fixed,scale=1e+06,gradient)
622   - 2014-01-09T21:47:54: Number of temporal par's ..: 1
623   - 2014-01-09T21:47:54: Constant .................: 1 (relative value) (fixed,scale=1,gradient)
624   - 2014-01-09T21:47:54: === GCTAModelRadialAcceptance ===
625   - 2014-01-09T21:47:54: Name ......................: Background
626   - 2014-01-09T21:47:54: Instruments ...............: CTA
627   - 2014-01-09T21:47:54: Instrument scale factors ..: unity
628   - 2014-01-09T21:47:54: Observation identifiers ...: all
629   - 2014-01-09T21:47:54: Model type ................: "Gaussian" * "FileFunction" * "Constant"
630   - 2014-01-09T21:47:54: Number of parameters ......: 3
631   - 2014-01-09T21:47:54: Number of radial par's ....: 1
632   - 2014-01-09T21:47:54: Sigma ....................: 3.04429 +/- 0.0330229 [0.01,10] deg2 (free,scale=1,gradient)
633   - 2014-01-09T21:47:54: Number of spectral par's ..: 1
634   - 2014-01-09T21:47:54: Normalization ............: 0.997516 +/- 0.0175361 [0,1000] (free,scale=1,gradient)
635   - 2014-01-09T21:47:54: Number of temporal par's ..: 1
636   - 2014-01-09T21:47:54: Constant .................: 1 (relative value) (fixed,scale=1,gradient)
637   - 2014-01-09T21:47:54:
638   - 2014-01-09T21:47:54: +==============+
639   - 2014-01-09T21:47:54: | Save results |
640   - 2014-01-09T21:47:54: +==============+
641   - 2014-01-09T21:47:54:
642   - 2014-01-09T21:47:54: Application "ctlike" terminated after 6 wall clock seconds, consuming 0.163384 seconds of CPU time.
  590 + 2014-10-30T22:55:43: +=================================+
  591 + 2014-10-30T22:55:43: | Maximum likelihood optimisation |
  592 + 2014-10-30T22:55:43: +=================================+
  593 + 2014-10-30T22:55:43: >Iteration 0: -logL=44400.561, Lambda=1.0e-03
  594 + 2014-10-30T22:55:43: >Iteration 1: -logL=44396.253, Lambda=1.0e-03, delta=4.308, max(|grad|)=9.645475 [Sigma:6]
  595 + 2014-10-30T22:55:43: >Iteration 2: -logL=44396.162, Lambda=1.0e-04, delta=0.091, max(|grad|)=-0.336139 [Sigma:6]
  596 + 2014-10-30T22:55:43: >Iteration 3: -logL=44396.162, Lambda=1.0e-05, delta=0.000, max(|grad|)=0.005603 [Sigma:6]
  597 + 2014-10-30T22:55:43:
  598 + 2014-10-30T22:55:43: +=========================================+
  599 + 2014-10-30T22:55:43: | Maximum likelihood optimization results |
  600 + 2014-10-30T22:55:43: +=========================================+
  601 + 2014-10-30T22:55:43: === GOptimizerLM ===
  602 + 2014-10-30T22:55:43: Optimized function value ..: 44396.162
  603 + 2014-10-30T22:55:43: Absolute precision ........: 0.005
  604 + 2014-10-30T22:55:43: Optimization status .......: converged
  605 + 2014-10-30T22:55:43: Number of parameters ......: 9
  606 + 2014-10-30T22:55:43: Number of free parameters .: 4
  607 + 2014-10-30T22:55:43: Number of iterations ......: 3
  608 + 2014-10-30T22:55:43: Lambda ....................: 1e-06
  609 + 2014-10-30T22:55:43: Maximum log likelihood ....: -44396.162
  610 + 2014-10-30T22:55:43: Observed events (Nobs) ...: 6127.000
  611 + 2014-10-30T22:55:43: Predicted events (Npred) ..: 6127.000 (Nobs - Npred = 0.000141861)
  612 + 2014-10-30T22:55:43: === GModels ===
  613 + 2014-10-30T22:55:43: Number of models ..........: 2
  614 + 2014-10-30T22:55:43: Number of parameters ......: 9
  615 + 2014-10-30T22:55:43: === GModelSky ===
  616 + 2014-10-30T22:55:43: Name ......................: Crab
  617 + 2014-10-30T22:55:43: Instruments ...............: all
  618 + 2014-10-30T22:55:43: Instrument scale factors ..: unity
  619 + 2014-10-30T22:55:43: Observation identifiers ...: all
  620 + 2014-10-30T22:55:43: Model type ................: PointSource
  621 + 2014-10-30T22:55:43: Model components ..........: "SkyDirFunction" * "PowerLaw" * "Constant"
  622 + 2014-10-30T22:55:43: Number of parameters ......: 6
  623 + 2014-10-30T22:55:43: Number of spatial par's ...: 2
  624 + 2014-10-30T22:55:43: RA .......................: 83.6331 [-360,360] deg (fixed,scale=1)
  625 + 2014-10-30T22:55:43: DEC ......................: 22.0145 [-90,90] deg (fixed,scale=1)
  626 + 2014-10-30T22:55:43: Number of spectral par's ..: 3
  627 + 2014-10-30T22:55:43: Prefactor ................: 6.13326e-16 +/- 2.05738e-17 [1e-23,1e-13] ph/cm2/s/MeV (free,scale=1e-16,gradient)
  628 + 2014-10-30T22:55:43: Index ....................: -2.5057 +/- 0.0250818 [-0,-5] (free,scale=-1,gradient)
  629 + 2014-10-30T22:55:43: PivotEnergy ..............: 300000 [10000,1e+09] MeV (fixed,scale=1e+06,gradient)
  630 + 2014-10-30T22:55:43: Number of temporal par's ..: 1
  631 + 2014-10-30T22:55:43: Constant .................: 1 (relative value) (fixed,scale=1,gradient)
  632 + 2014-10-30T22:55:43: === GCTAModelRadialAcceptance ===
  633 + 2014-10-30T22:55:43: Name ......................: Background
  634 + 2014-10-30T22:55:43: Instruments ...............: CTA
  635 + 2014-10-30T22:55:43: Instrument scale factors ..: unity
  636 + 2014-10-30T22:55:43: Observation identifiers ...: all
  637 + 2014-10-30T22:55:43: Model type ................: "Gaussian" * "FileFunction" * "Constant"
  638 + 2014-10-30T22:55:43: Number of parameters ......: 3
  639 + 2014-10-30T22:55:43: Number of radial par's ....: 1
  640 + 2014-10-30T22:55:43: Sigma ....................: 3.04429 +/- 0.0330229 [0.01,10] deg2 (free,scale=1,gradient)
  641 + 2014-10-30T22:55:43: Number of spectral par's ..: 1
  642 + 2014-10-30T22:55:43: Normalization ............: 0.996765 +/- 0.0175227 [0,1000] (free,scale=1,gradient)
  643 + 2014-10-30T22:55:43: Number of temporal par's ..: 1
  644 + 2014-10-30T22:55:43: Constant .................: 1 (relative value) (fixed,scale=1,gradient)
  645 + 2014-10-30T22:55:43:
  646 + 2014-10-30T22:55:43: +==============+
  647 + 2014-10-30T22:55:43: | Save results |
  648 + 2014-10-30T22:55:43: +==============+
  649 + 2014-10-30T22:55:43:
  650 + 2014-10-30T22:55:43: Application "ctlike" terminated after 5 wall clock seconds, consuming 0.246747 seconds of CPU time.
643 651  
644 652  
645 653 .. _sec_butterfly_cta:
... ... @@ -658,14 +666,14 @@ line.
658 666 .. code-block:: bash
659 667  
660 668 $ ctbutterfly
661   - Input event list, cube or observation definition file [$CTOOLS/test/data/crab_events.fits.gz]
662   - Calibration database [share/caldb/data/cta/dummy//cta_dummy_irf] share/caldb/data/cta/dummy/bcf/
  669 + Input event list, cube or observation definition file [events.fits]
  670 + Calibration database [dummy]
663 671 Instrument response function [cta_dummy_irf]
664   - Source model [$CTOOLS/test/data/crab.xml]
665   - Source of interest [CrabNebula] Crab
  672 + Source model [$CTOOLS/share/models/crab.xml] crab_results.xml
  673 + Source of interest [Crab]
666 674 Start value for first energy bin in TeV [0.1]
667   - Stop value for last energy bin in TeV [100.]
668   - Output ascii file [butterfly.txt]
  675 + Stop value for last energy bin in TeV [100.0]
  676 + Output ascii file [butterfly.txt]
669 677  
670 678 Below some lines of the ``ctbutterfly.log``:
671 679  
... ... @@ -692,7 +700,8 @@ in the ctools ``example`` folder. You will need matplotlib on your system
692 700 to make this work. To launch the script, type:
693 701  
694 702 .. code-block:: bash
695   - python $CTOOLS/share/examples/python/show_butterfly.py butterfly.txt
  703 +
  704 + python $CTOOLS/share/examples/python/show_butterfly.py butterfly.txt
696 705  
697 706 This will result in a canvas which should look like the following:
698 707  
... ...