diff --git a/doc/source/users/tutorials/1dc/first_fitting.rst b/doc/source/users/tutorials/1dc/first_fitting.rst index 00ce287..7f6a340 100644 --- a/doc/source/users/tutorials/1dc/first_fitting.rst +++ b/doc/source/users/tutorials/1dc/first_fitting.rst @@ -30,34 +30,45 @@ converged properly, as illustrated in the example below: .. code-block:: bash - 2017-03-03T18:06:30: +=================================+ - 2017-03-03T18:06:30: | Maximum likelihood optimisation | - 2017-03-03T18:06:30: +=================================+ - 2017-03-03T18:06:56: >Iteration 0: -logL=18221386.597, Lambda=1.0e-03 - 2017-03-03T18:07:19: >Iteration 1: -logL=18206762.996, Lambda=1.0e-03, delta=14623.601, max(|grad|)=45198.417696 [Index:13] - 2017-03-03T18:07:42: >Iteration 2: -logL=18203712.440, Lambda=1.0e-04, delta=3050.557, max(|grad|)=10013.593417 [Index:3] - 2017-03-03T18:08:04: >Iteration 3: -logL=18202302.390, Lambda=1.0e-05, delta=1410.050, max(|grad|)=7461.719387 [Index:3] - 2017-03-03T18:08:27: >Iteration 4: -logL=18201825.720, Lambda=1.0e-06, delta=476.669, max(|grad|)=4076.653047 [RA:0] - 2017-03-03T18:08:50: >Iteration 5: -logL=18201757.606, Lambda=1.0e-07, delta=68.114, max(|grad|)=1547.953907 [RA:0] - 2017-03-03T18:09:13: >Iteration 6: -logL=18201754.932, Lambda=1.0e-08, delta=2.674, max(|grad|)=388.888160 [RA:0] - 2017-03-03T18:09:36: >Iteration 7: -logL=18201754.869, Lambda=1.0e-09, delta=0.063, max(|grad|)=81.488909 [RA:0] - 2017-03-03T18:09:59: >Iteration 8: -logL=18201754.866, Lambda=1.0e-10, delta=0.003, max(|grad|)=16.367278 [RA:0] - 2017-03-03T18:10:22: - 2017-03-03T18:10:22: +=========================================+ - 2017-03-03T18:10:22: | Maximum likelihood optimisation results | - 2017-03-03T18:10:22: +=========================================+ - 2017-03-03T18:10:22: === GOptimizerLM === - 2017-03-03T18:10:22: Optimized function value ..: 18201754.866 - 2017-03-03T18:10:22: Absolute precision ........: 0.005 - 2017-03-03T18:10:22: Acceptable value decrease .: 2 - 2017-03-03T18:10:22: Optimization status .......: converged - 2017-03-03T18:10:22: Number of parameters ......: 16 - 2017-03-03T18:10:22: Number of free parameters .: 10 - 2017-03-03T18:10:22: Number of iterations ......: 8 - 2017-03-03T18:10:22: Lambda ....................: 1e-11 - 2017-03-03T18:10:22: Maximum log likelihood ....: -18201754.866 - 2017-03-03T18:10:22: Observed events (Nobs) ...: 2841472.000 - 2017-03-03T18:10:22: Predicted events (Npred) ..: 2841471.996 (Nobs - Npred = 0.00435844948515296) + 2017-05-30T14:04:49: +=================================+ + 2017-05-30T14:04:49: | Maximum likelihood optimisation | + 2017-05-30T14:04:49: +=================================+ + 2017-05-30T14:05:14: >Iteration 0: -logL=22930903.429, Lambda=1.0e-03 + 2017-05-30T14:05:38: >Iteration 1: -logL=22906772.079, Lambda=1.0e-03, delta=24131.350, step=1.0e+00, max(|grad|)=53111.535756 [Index:13] + 2017-05-30T14:06:00: >Iteration 2: -logL=22903418.347, Lambda=1.0e-04, delta=3353.732, step=1.0e+00, max(|grad|)=-11638.756860 [RA:0] + 2017-05-30T14:06:23: >Iteration 3: -logL=22902335.786, Lambda=1.0e-05, delta=1082.561, step=1.0e+00, max(|grad|)=-9993.636957 [RA:0] + 2017-05-30T14:06:46: >Iteration 4: -logL=22902144.463, Lambda=1.0e-06, delta=191.323, step=1.0e+00, max(|grad|)=-7053.128713 [RA:0] + 2017-05-30T14:07:09: >Iteration 5: -logL=22902098.908, Lambda=1.0e-07, delta=45.555, step=1.0e+00, max(|grad|)=4767.092710 [RA:6] + 2017-05-30T14:07:33: >Iteration 6: -logL=22902077.092, Lambda=1.0e-08, delta=21.816, step=1.0e+00, max(|grad|)=3512.238312 [RA:6] + 2017-05-30T14:07:56: >Iteration 7: -logL=22902066.319, Lambda=1.0e-09, delta=10.774, step=1.0e+00, max(|grad|)=2564.648439 [RA:6] + 2017-05-30T14:08:21: >Iteration 8: -logL=22902060.969, Lambda=1.0e-10, delta=5.350, step=1.0e+00, max(|grad|)=1861.818754 [RA:6] + 2017-05-30T14:08:45: >Iteration 9: -logL=22902058.299, Lambda=1.0e-11, delta=2.670, step=1.0e+00, max(|grad|)=1346.424825 [RA:6] + 2017-05-30T14:09:08: >Iteration 10: -logL=22902056.961, Lambda=1.0e-12, delta=1.338, step=1.0e+00, max(|grad|)=971.176097 [RA:6] + 2017-05-30T14:09:31: >Iteration 11: -logL=22902056.287, Lambda=1.0e-13, delta=0.674, step=1.0e+00, max(|grad|)=699.274062 [RA:6] + 2017-05-30T14:09:54: >Iteration 12: -logL=22902055.947, Lambda=1.0e-14, delta=0.340, step=1.0e+00, max(|grad|)=502.893171 [RA:6] + 2017-05-30T14:10:18: >Iteration 13: -logL=22902055.774, Lambda=1.0e-15, delta=0.172, step=1.0e+00, max(|grad|)=361.372646 [RA:6] + 2017-05-30T14:10:41: >Iteration 14: -logL=22902055.687, Lambda=1.0e-16, delta=0.088, step=1.0e+00, max(|grad|)=259.519694 [RA:6] + 2017-05-30T14:11:04: >Iteration 15: -logL=22902055.642, Lambda=1.0e-17, delta=0.045, step=1.0e+00, max(|grad|)=186.309625 [RA:6] + 2017-05-30T14:11:27: >Iteration 16: -logL=22902055.620, Lambda=1.0e-18, delta=0.023, step=1.0e+00, max(|grad|)=133.719821 [RA:6] + 2017-05-30T14:11:50: >Iteration 17: -logL=22902055.608, Lambda=1.0e-19, delta=0.012, step=1.0e+00, max(|grad|)=95.962554 [RA:6] + 2017-05-30T14:12:13: >Iteration 18: -logL=22902055.602, Lambda=1.0e-20, delta=0.006, step=1.0e+00, max(|grad|)=68.863372 [RA:6] + 2017-05-30T14:12:35: >Iteration 19: -logL=22902055.599, Lambda=1.0e-21, delta=0.003, step=1.0e+00, max(|grad|)=49.416455 [RA:6] + 2017-05-30T14:12:58: + 2017-05-30T14:12:58: +=========================================+ + 2017-05-30T14:12:58: | Maximum likelihood optimisation results | + 2017-05-30T14:12:58: +=========================================+ + 2017-05-30T14:12:58: === GOptimizerLM === + 2017-05-30T14:12:58: Optimized function value ..: 22902055.599 + 2017-05-30T14:12:58: Absolute precision ........: 0.005 + 2017-05-30T14:12:58: Acceptable value decrease .: 2 + 2017-05-30T14:12:58: Optimization status .......: converged + 2017-05-30T14:12:58: Number of parameters ......: 16 + 2017-05-30T14:12:58: Number of free parameters .: 10 + 2017-05-30T14:12:58: Number of iterations ......: 19 + 2017-05-30T14:12:58: Lambda ....................: 1e-22 + 2017-05-30T14:12:58: Maximum log likelihood ....: -22902055.599 + 2017-05-30T14:12:58: Observed events (Nobs) ...: 3348255.000 + 2017-05-30T14:12:58: Predicted events (Npred) ..: 3348253.996 (Nobs - Npred = 1.00355448853225) You may also convert the fitted model positions into a `ds9 `_ region file using the :ref:`csmodelinfo` script so that you can overlay the @@ -84,7 +95,7 @@ which is required to assure the proper convergence of the fit. :width: 600px :align: center - *Background subtracted sky map of the events recorded around the Galactic Centre during the Galactic Centre Survey with the fitted positions of the sources shown as black circles* + *Background subtracted sky map of the events recorded around the Galactic Centre during the Galactic Plane Survey with the fitted positions of the sources shown as black circles* You can also convert the spectral parameters of the point sources into a butterfly diagram for each source using the :ref:`ctbutterfly` tool. diff --git a/doc/source/users/tutorials/1dc/first_skymap_fitted.png b/doc/source/users/tutorials/1dc/first_skymap_fitted.png index 25a381c..e9dc8c4 100644 --- a/doc/source/users/tutorials/1dc/first_skymap_fitted.png +++ b/doc/source/users/tutorials/1dc/first_skymap_fitted.png diff --git a/doc/source/users/tutorials/1dc/first_spectrum.png b/doc/source/users/tutorials/1dc/first_spectrum.png index c9268eb..204ad0f 100644 --- a/doc/source/users/tutorials/1dc/first_spectrum.png +++ b/doc/source/users/tutorials/1dc/first_spectrum.png