Commit 9538a7c2464df9dafc8ca514a9885b9b158ab813
1 parent
b7ab3459
Update 1DC tutorials
Showing
3 changed files
with
40 additions
and
29 deletions
doc/source/users/tutorials/1dc/first_fitting.rst
| ... | ... | @@ -30,34 +30,45 @@ converged properly, as illustrated in the example below: |
| 30 | 30 | |
| 31 | 31 | .. code-block:: bash |
| 32 | 32 | |
| 33 | - 2017-03-03T18:06:30: +=================================+ | |
| 34 | - 2017-03-03T18:06:30: | Maximum likelihood optimisation | | |
| 35 | - 2017-03-03T18:06:30: +=================================+ | |
| 36 | - 2017-03-03T18:06:56: >Iteration 0: -logL=18221386.597, Lambda=1.0e-03 | |
| 37 | - 2017-03-03T18:07:19: >Iteration 1: -logL=18206762.996, Lambda=1.0e-03, delta=14623.601, max(|grad|)=45198.417696 [Index:13] | |
| 38 | - 2017-03-03T18:07:42: >Iteration 2: -logL=18203712.440, Lambda=1.0e-04, delta=3050.557, max(|grad|)=10013.593417 [Index:3] | |
| 39 | - 2017-03-03T18:08:04: >Iteration 3: -logL=18202302.390, Lambda=1.0e-05, delta=1410.050, max(|grad|)=7461.719387 [Index:3] | |
| 40 | - 2017-03-03T18:08:27: >Iteration 4: -logL=18201825.720, Lambda=1.0e-06, delta=476.669, max(|grad|)=4076.653047 [RA:0] | |
| 41 | - 2017-03-03T18:08:50: >Iteration 5: -logL=18201757.606, Lambda=1.0e-07, delta=68.114, max(|grad|)=1547.953907 [RA:0] | |
| 42 | - 2017-03-03T18:09:13: >Iteration 6: -logL=18201754.932, Lambda=1.0e-08, delta=2.674, max(|grad|)=388.888160 [RA:0] | |
| 43 | - 2017-03-03T18:09:36: >Iteration 7: -logL=18201754.869, Lambda=1.0e-09, delta=0.063, max(|grad|)=81.488909 [RA:0] | |
| 44 | - 2017-03-03T18:09:59: >Iteration 8: -logL=18201754.866, Lambda=1.0e-10, delta=0.003, max(|grad|)=16.367278 [RA:0] | |
| 45 | - 2017-03-03T18:10:22: | |
| 46 | - 2017-03-03T18:10:22: +=========================================+ | |
| 47 | - 2017-03-03T18:10:22: | Maximum likelihood optimisation results | | |
| 48 | - 2017-03-03T18:10:22: +=========================================+ | |
| 49 | - 2017-03-03T18:10:22: === GOptimizerLM === | |
| 50 | - 2017-03-03T18:10:22: Optimized function value ..: 18201754.866 | |
| 51 | - 2017-03-03T18:10:22: Absolute precision ........: 0.005 | |
| 52 | - 2017-03-03T18:10:22: Acceptable value decrease .: 2 | |
| 53 | - 2017-03-03T18:10:22: Optimization status .......: converged | |
| 54 | - 2017-03-03T18:10:22: Number of parameters ......: 16 | |
| 55 | - 2017-03-03T18:10:22: Number of free parameters .: 10 | |
| 56 | - 2017-03-03T18:10:22: Number of iterations ......: 8 | |
| 57 | - 2017-03-03T18:10:22: Lambda ....................: 1e-11 | |
| 58 | - 2017-03-03T18:10:22: Maximum log likelihood ....: -18201754.866 | |
| 59 | - 2017-03-03T18:10:22: Observed events (Nobs) ...: 2841472.000 | |
| 60 | - 2017-03-03T18:10:22: Predicted events (Npred) ..: 2841471.996 (Nobs - Npred = 0.00435844948515296) | |
| 33 | + 2017-05-30T14:04:49: +=================================+ | |
| 34 | + 2017-05-30T14:04:49: | Maximum likelihood optimisation | | |
| 35 | + 2017-05-30T14:04:49: +=================================+ | |
| 36 | + 2017-05-30T14:05:14: >Iteration 0: -logL=22930903.429, Lambda=1.0e-03 | |
| 37 | + 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] | |
| 38 | + 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] | |
| 39 | + 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] | |
| 40 | + 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] | |
| 41 | + 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] | |
| 42 | + 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] | |
| 43 | + 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] | |
| 44 | + 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] | |
| 45 | + 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] | |
| 46 | + 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] | |
| 47 | + 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] | |
| 48 | + 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] | |
| 49 | + 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] | |
| 50 | + 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] | |
| 51 | + 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] | |
| 52 | + 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] | |
| 53 | + 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] | |
| 54 | + 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] | |
| 55 | + 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] | |
| 56 | + 2017-05-30T14:12:58: | |
| 57 | + 2017-05-30T14:12:58: +=========================================+ | |
| 58 | + 2017-05-30T14:12:58: | Maximum likelihood optimisation results | | |
| 59 | + 2017-05-30T14:12:58: +=========================================+ | |
| 60 | + 2017-05-30T14:12:58: === GOptimizerLM === | |
| 61 | + 2017-05-30T14:12:58: Optimized function value ..: 22902055.599 | |
| 62 | + 2017-05-30T14:12:58: Absolute precision ........: 0.005 | |
| 63 | + 2017-05-30T14:12:58: Acceptable value decrease .: 2 | |
| 64 | + 2017-05-30T14:12:58: Optimization status .......: converged | |
| 65 | + 2017-05-30T14:12:58: Number of parameters ......: 16 | |
| 66 | + 2017-05-30T14:12:58: Number of free parameters .: 10 | |
| 67 | + 2017-05-30T14:12:58: Number of iterations ......: 19 | |
| 68 | + 2017-05-30T14:12:58: Lambda ....................: 1e-22 | |
| 69 | + 2017-05-30T14:12:58: Maximum log likelihood ....: -22902055.599 | |
| 70 | + 2017-05-30T14:12:58: Observed events (Nobs) ...: 3348255.000 | |
| 71 | + 2017-05-30T14:12:58: Predicted events (Npred) ..: 3348253.996 (Nobs - Npred = 1.00355448853225) | |
| 61 | 72 | |
| 62 | 73 | You may also convert the fitted model positions into a `ds9 <http://ds9.si.edu>`_ |
| 63 | 74 | 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. |
| 84 | 95 | :width: 600px |
| 85 | 96 | :align: center |
| 86 | 97 | |
| 87 | - *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* | |
| 98 | + *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* | |
| 88 | 99 | |
| 89 | 100 | You can also convert the spectral parameters of the point sources into a |
| 90 | 101 | butterfly diagram for each source using the :ref:`ctbutterfly` tool. | ... | ... |
doc/source/users/tutorials/1dc/first_skymap_fitted.png
doc/source/users/tutorials/1dc/first_spectrum.png
