From 0873cb8eebcf86edcae88001e5b555f059194b9b Mon Sep 17 00:00:00 2001 From: Jürgen Knödlseder Date: Thu, 4 Apr 2019 17:10:22 +0200 Subject: [PATCH] Update Fermi How To --- doc/source/users/tutorials/howto/fermi/howto_fermi_butterfly.png | Bin 41508 -> 0 bytes doc/source/users/tutorials/howto/fermi/howto_fermi_errors.rst | 154 +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++----------------------------------------------------------------------------- doc/source/users/tutorials/howto/fermi/howto_fermi_fitting.rst | 184 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++-------------------------------------------------------------------------------------------- doc/source/users/tutorials/howto/fermi/howto_fermi_spectrum.png | Bin 22496 -> 0 bytes doc/source/users/tutorials/howto/fermi/howto_fermi_tsmap.rst | 6 +++--- doc/source/users/tutorials/howto/fermi/howto_fermi_ulimit.rst | 14 +++++++------- 6 files changed, 179 insertions(+), 179 deletions(-) diff --git a/doc/source/users/tutorials/howto/fermi/howto_fermi_butterfly.png b/doc/source/users/tutorials/howto/fermi/howto_fermi_butterfly.png index 007abed..fa5416a 100755 Binary files a/doc/source/users/tutorials/howto/fermi/howto_fermi_butterfly.png and b/doc/source/users/tutorials/howto/fermi/howto_fermi_butterfly.png differ diff --git a/doc/source/users/tutorials/howto/fermi/howto_fermi_errors.rst b/doc/source/users/tutorials/howto/fermi/howto_fermi_errors.rst index f24b0d1..a00f65e 100755 --- a/doc/source/users/tutorials/howto/fermi/howto_fermi_errors.rst +++ b/doc/source/users/tutorials/howto/fermi/howto_fermi_errors.rst @@ -33,80 +33,80 @@ covariance matrix. .. code-block:: none - 2018-01-25T20:45:22: +================================================+ - 2018-01-25T20:45:22: | Compute error for source "Vela" parameter "RA" | - 2018-01-25T20:45:22: +================================================+ - 2018-01-25T20:45:22: Confidence level ..........: 68 % - 2018-01-25T20:45:22: Log-likelihood difference .: 0.494473 - 2018-01-25T20:45:22: Initial factor range ......: [128.819554251926, 128.851989886434] - 2018-01-25T20:45:41: Lower parameter factor ....: 128.834 - 2018-01-25T20:45:41: Upper parameter factor ....: 128.837 - 2018-01-25T20:45:41: Error from curvature ......: 0.00162178172537938 deg - 2018-01-25T20:45:41: Error from profile ........: 0.00169067577327553 deg - 2018-01-25T20:45:41: Negative profile error ....: 0.00165900034895117 deg - 2018-01-25T20:45:41: Positive profile error ....: 0.00172235119759989 deg - 2018-01-25T20:45:41: - 2018-01-25T20:45:41: +=================================================+ - 2018-01-25T20:45:41: | Compute error for source "Vela" parameter "DEC" | - 2018-01-25T20:45:41: +=================================================+ - 2018-01-25T20:45:41: Confidence level ..........: 68 % - 2018-01-25T20:45:41: Log-likelihood difference .: 0.494473 - 2018-01-25T20:45:41: Initial factor range ......: [-45.1949881923928, -45.1721741961801] - 2018-01-25T20:46:01: Lower parameter factor ....: -45.1847 - 2018-01-25T20:46:01: Upper parameter factor ....: -45.1824 - 2018-01-25T20:46:01: Error from curvature ......: 0.00114069981063169 deg - 2018-01-25T20:46:01: Error from profile ........: 0.00116548552428952 deg - 2018-01-25T20:46:01: Negative profile error ....: 0.00115852324517363 deg - 2018-01-25T20:46:01: Positive profile error ....: 0.00117244780340542 deg - 2018-01-25T20:46:01: - 2018-01-25T20:46:01: +=======================================================+ - 2018-01-25T20:46:01: | Compute error for source "Vela" parameter "Prefactor" | - 2018-01-25T20:46:01: +=======================================================+ - 2018-01-25T20:46:01: Confidence level ..........: 68 % - 2018-01-25T20:46:01: Log-likelihood difference .: 0.494473 - 2018-01-25T20:46:01: Initial factor range ......: [1e-07, 9.09784331905335] - 2018-01-25T20:46:45: Lower parameter factor ....: 3.93674 - 2018-01-25T20:46:45: Upper parameter factor ....: 4.89558 - 2018-01-25T20:46:45: Error from curvature ......: 4.74564602380296e-10 ph/cm2/s/MeV - 2018-01-25T20:46:45: Error from profile ........: 4.79421031311828e-10 ph/cm2/s/MeV - 2018-01-25T20:46:45: Negative profile error ....: 4.15456324058327e-10 ph/cm2/s/MeV - 2018-01-25T20:46:45: Positive profile error ....: 5.4338573856533e-10 ph/cm2/s/MeV - 2018-01-25T20:46:45: - 2018-01-25T20:46:45: +====================================================+ - 2018-01-25T20:46:45: | Compute error for source "Vela" parameter "Index1" | - 2018-01-25T20:46:45: +====================================================+ - 2018-01-25T20:46:45: Confidence level ..........: 68 % - 2018-01-25T20:46:45: Log-likelihood difference .: 0.494473 - 2018-01-25T20:46:45: Initial factor range ......: [1.0323933217421, 1.65461803910864] - 2018-01-25T20:47:31: Lower parameter factor ....: 1.31084 - 2018-01-25T20:47:31: Upper parameter factor ....: 1.37275 - 2018-01-25T20:47:31: Error from curvature ......: 0.0311112358683271 - 2018-01-25T20:47:31: Error from profile ........: 0.0309517300203693 - 2018-01-25T20:47:31: Negative profile error ....: 0.032660721248488 - 2018-01-25T20:47:31: Positive profile error ....: 0.0292427387922507 - 2018-01-25T20:47:31: - 2018-01-25T20:47:31: +==========================================================+ - 2018-01-25T20:47:31: | Compute error for source "Vela" parameter "CutoffEnergy" | - 2018-01-25T20:47:31: +==========================================================+ - 2018-01-25T20:47:31: Confidence level ..........: 68 % - 2018-01-25T20:47:31: Log-likelihood difference .: 0.494473 - 2018-01-25T20:47:31: Initial factor range ......: [0.001, 2.80853448949703] - 2018-01-25T20:48:14: Lower parameter factor ....: 0.81239 - 2018-01-25T20:48:14: Upper parameter factor ....: 1.17294 - 2018-01-25T20:48:14: Error from curvature ......: 181.899981380157 MeV - 2018-01-25T20:48:14: Error from profile ........: 180.277267463442 MeV - 2018-01-25T20:48:14: Negative profile error ....: 177.144641591912 MeV - 2018-01-25T20:48:14: Positive profile error ....: 183.409893334972 MeV - 2018-01-25T20:48:14: - 2018-01-25T20:48:14: +====================================================+ - 2018-01-25T20:48:14: | Compute error for source "Vela" parameter "Index2" | - 2018-01-25T20:48:14: +====================================================+ - 2018-01-25T20:48:14: Confidence level ..........: 68 % - 2018-01-25T20:48:14: Log-likelihood difference .: 0.494473 - 2018-01-25T20:48:14: Initial factor range ......: [0.308028525758507, 0.866644263529458] - 2018-01-25T20:49:01: Lower parameter factor ....: 0.559037 - 2018-01-25T20:49:01: Upper parameter factor ....: 0.615056 - 2018-01-25T20:49:01: Error from curvature ......: 0.0279307868885475 - 2018-01-25T20:49:01: Error from profile ........: 0.0280092058458762 - 2018-01-25T20:49:01: Negative profile error ....: 0.028299015036004 - 2018-01-25T20:49:01: Positive profile error ....: 0.0277193966557485 + 2019-04-04T15:04:51: +================================================+ + 2019-04-04T15:04:51: | Compute error for source "Vela" parameter "RA" | + 2019-04-04T15:04:51: +================================================+ + 2019-04-04T15:04:51: Confidence level ..........: 68 % + 2019-04-04T15:04:51: Log-likelihood difference .: 0.494473 + 2019-04-04T15:04:51: Initial factor range ......: [128.819554251886, 128.851989886469] + 2019-04-04T15:05:10: Lower parameter factor ....: 128.834 + 2019-04-04T15:05:10: Upper parameter factor ....: 128.837 + 2019-04-04T15:05:10: Error from curvature ......: 0.00162178172918532 deg + 2019-04-04T15:05:10: Error from profile ........: 0.00169067577726878 deg + 2019-04-04T15:05:10: Negative profile error ....: 0.00165900035287336 deg + 2019-04-04T15:05:10: Positive profile error ....: 0.0017223512016642 deg + 2019-04-04T15:05:10: + 2019-04-04T15:05:10: +=================================================+ + 2019-04-04T15:05:10: | Compute error for source "Vela" parameter "DEC" | + 2019-04-04T15:05:10: +=================================================+ + 2019-04-04T15:05:10: Confidence level ..........: 68 % + 2019-04-04T15:05:10: Log-likelihood difference .: 0.494473 + 2019-04-04T15:05:10: Initial factor range ......: [-45.1949881925207, -45.1721741960936] + 2019-04-04T15:05:31: Lower parameter factor ....: -45.1847 + 2019-04-04T15:05:31: Upper parameter factor ....: -45.1824 + 2019-04-04T15:05:31: Error from curvature ......: 0.00114069982135549 deg + 2019-04-04T15:05:31: Error from profile ........: 0.0011654855352532 deg + 2019-04-04T15:05:31: Negative profile error ....: 0.00115852325606625 deg + 2019-04-04T15:05:31: Positive profile error ....: 0.00117244781444015 deg + 2019-04-04T15:05:31: + 2019-04-04T15:05:31: +=======================================================+ + 2019-04-04T15:05:31: | Compute error for source "Vela" parameter "Prefactor" | + 2019-04-04T15:05:31: +=======================================================+ + 2019-04-04T15:05:31: Confidence level ..........: 68 % + 2019-04-04T15:05:31: Log-likelihood difference .: 0.494473 + 2019-04-04T15:05:31: Initial factor range ......: [1e-07, 9.09784331943508] + 2019-04-04T15:06:15: Lower parameter factor ....: 3.93674 + 2019-04-04T15:06:15: Upper parameter factor ....: 4.89558 + 2019-04-04T15:06:15: Error from curvature ......: 4.74564602405264e-10 ph/cm2/s/MeV + 2019-04-04T15:06:15: Error from profile ........: 4.79421031332425e-10 ph/cm2/s/MeV + 2019-04-04T15:06:15: Negative profile error ....: 4.15456324070932e-10 ph/cm2/s/MeV + 2019-04-04T15:06:15: Positive profile error ....: 5.43385738593917e-10 ph/cm2/s/MeV + 2019-04-04T15:06:15: + 2019-04-04T15:06:15: +====================================================+ + 2019-04-04T15:06:15: | Compute error for source "Vela" parameter "Index1" | + 2019-04-04T15:06:15: +====================================================+ + 2019-04-04T15:06:15: Confidence level ..........: 68 % + 2019-04-04T15:06:15: Log-likelihood difference .: 0.494473 + 2019-04-04T15:06:15: Initial factor range ......: [1.03239332172832, 1.65461803910225] + 2019-04-04T15:07:01: Lower parameter factor ....: 1.31084 + 2019-04-04T15:07:01: Upper parameter factor ....: 1.37275 + 2019-04-04T15:07:01: Error from curvature ......: 0.0311112358686962 + 2019-04-04T15:07:01: Error from profile ........: 0.0309517300207367 + 2019-04-04T15:07:01: Negative profile error ....: 0.0326607212488756 + 2019-04-04T15:07:01: Positive profile error ....: 0.0292427387925978 + 2019-04-04T15:07:01: + 2019-04-04T15:07:01: +==========================================================+ + 2019-04-04T15:07:01: | Compute error for source "Vela" parameter "CutoffEnergy" | + 2019-04-04T15:07:01: +==========================================================+ + 2019-04-04T15:07:01: Confidence level ..........: 68 % + 2019-04-04T15:07:01: Log-likelihood difference .: 0.494473 + 2019-04-04T15:07:01: Initial factor range ......: [0.001, 2.80853448937906] + 2019-04-04T15:07:45: Lower parameter factor ....: 0.81239 + 2019-04-04T15:07:45: Upper parameter factor ....: 1.17294 + 2019-04-04T15:07:45: Error from curvature ......: 181.899981373502 MeV + 2019-04-04T15:07:45: Error from profile ........: 180.277267455479 MeV + 2019-04-04T15:07:45: Negative profile error ....: 177.144641582696 MeV + 2019-04-04T15:07:45: Positive profile error ....: 183.409893328263 MeV + 2019-04-04T15:07:45: + 2019-04-04T15:07:45: +====================================================+ + 2019-04-04T15:07:45: | Compute error for source "Vela" parameter "Index2" | + 2019-04-04T15:07:45: +====================================================+ + 2019-04-04T15:07:45: Confidence level ..........: 68 % + 2019-04-04T15:07:45: Log-likelihood difference .: 0.494473 + 2019-04-04T15:07:45: Initial factor range ......: [0.3080285257564, 0.866644263519715] + 2019-04-04T15:08:32: Lower parameter factor ....: 0.559037 + 2019-04-04T15:08:32: Upper parameter factor ....: 0.615056 + 2019-04-04T15:08:32: Error from curvature ......: 0.0279307868881658 + 2019-04-04T15:08:32: Error from profile ........: 0.0280092058454933 + 2019-04-04T15:08:32: Negative profile error ....: 0.028299015035617 + 2019-04-04T15:08:32: Positive profile error ....: 0.0277193966553695 diff --git a/doc/source/users/tutorials/howto/fermi/howto_fermi_fitting.rst b/doc/source/users/tutorials/howto/fermi/howto_fermi_fitting.rst index f8abbc1..7668232 100755 --- a/doc/source/users/tutorials/howto/fermi/howto_fermi_fitting.rst +++ b/doc/source/users/tutorials/howto/fermi/howto_fermi_fitting.rst @@ -92,95 +92,95 @@ log file ``ctlike.log`` that was also created by :ref:`ctlike`: .. code-block:: none - 2018-01-25T20:07:55: +=================================+ - 2018-01-25T20:07:55: | Maximum likelihood optimisation | - 2018-01-25T20:07:55: +=================================+ - 2018-01-25T20:07:57: >Iteration 0: -logL=-177232.327, Lambda=1.0e-03 - 2018-01-25T20:07:58: >Iteration 1: -logL=-198931.136, Lambda=1.0e-03, delta=21698.809, step=1.0e+00, max(|grad|)=-22717.063466 [CutoffEnergy:4] - 2018-01-25T20:07:58: >Iteration 2: -logL=-200269.888, Lambda=1.0e-04, delta=1338.753, step=1.0e+00, max(|grad|)=-67629.683771 [CutoffEnergy:4] - 2018-01-25T20:07:59: >Iteration 3: -logL=-212351.965, Lambda=1.0e-05, delta=12082.076, step=1.0e+00, max(|grad|)=23026.802421 [Index2:6] - 2018-01-25T20:07:59: >Iteration 4: -logL=-213954.259, Lambda=1.0e-06, delta=1602.294, step=1.0e+00, max(|grad|)=-8297.380912 [CutoffEnergy:4] - 2018-01-25T20:07:59: >Iteration 5: -logL=-215347.973, Lambda=1.0e-07, delta=1393.714, step=1.0e+00, max(|grad|)=2746.754042 [Index2:6] - 2018-01-25T20:08:00: >Iteration 6: -logL=-215497.807, Lambda=1.0e-08, delta=149.834, step=1.0e+00, max(|grad|)=699.848484 [Index2:6] - 2018-01-25T20:08:00: Iteration 7: -logL=-215497.807, Lambda=1.0e-09, delta=-13.845, step=1.0e+00, max(|grad|)=-1201.982875 [Index1:3] (stalled) - 2018-01-25T20:08:01: Iteration 8: -logL=-215497.807, Lambda=1.0e-08, delta=-13.844, step=1.0e+00, max(|grad|)=-1201.959647 [Index1:3] (stalled) - 2018-01-25T20:08:01: Iteration 9: -logL=-215497.807, Lambda=1.0e-07, delta=-13.835, step=1.0e+00, max(|grad|)=-1201.727401 [Index1:3] (stalled) - 2018-01-25T20:08:02: Iteration 10: -logL=-215497.807, Lambda=1.0e-06, delta=-13.744, step=1.0e+00, max(|grad|)=-1199.408717 [Index1:3] (stalled) - 2018-01-25T20:08:02: Iteration 11: -logL=-215497.807, Lambda=1.0e-05, delta=-12.858, step=1.0e+00, max(|grad|)=-1176.594021 [Index1:3] (stalled) - 2018-01-25T20:08:03: Iteration 12: -logL=-215497.807, Lambda=1.0e-04, delta=-6.126, step=1.0e+00, max(|grad|)=-980.833963 [Index1:3] (stalled) - 2018-01-25T20:08:04: >Iteration 13: -logL=-215503.803, Lambda=1.0e-03, delta=5.995, step=1.0e+00, max(|grad|)=-285.964915 [Index1:3] - 2018-01-25T20:08:04: Iteration 14: -logL=-215503.803, Lambda=1.0e-04, delta=-8.762, step=1.0e+00, max(|grad|)=-730.461326 [Index1:3] (stalled) - 2018-01-25T20:08:05: >Iteration 15: -logL=-215504.976, Lambda=1.0e-03, delta=1.173, step=1.0e+00, max(|grad|)=-153.941336 [Index1:3] - 2018-01-25T20:08:05: Iteration 16: -logL=-215503.593, Lambda=1.0e-04, delta=-1.383, step=1.0e+00, max(|grad|)=-327.368445 [Index1:3] (stalled) - 2018-01-25T20:08:06: >Iteration 17: -logL=-215505.817, Lambda=1.0e-03, delta=2.225, step=1.0e+00, max(|grad|)=3.256216 [Index2:6] - 2018-01-25T20:08:06: >Iteration 18: -logL=-215505.824, Lambda=1.0e-04, delta=0.007, step=1.0e+00, max(|grad|)=-8.190316 [CutoffEnergy:4] - 2018-01-25T20:08:06: >Iteration 19: -logL=-215505.826, Lambda=1.0e-05, delta=0.002, step=1.0e+00, max(|grad|)=0.635263 [Index2:6] - 2018-01-25T20:08:07: - 2018-01-25T20:08:07: +=========================================+ - 2018-01-25T20:08:07: | Maximum likelihood optimisation results | - 2018-01-25T20:08:07: +=========================================+ - 2018-01-25T20:08:07: === GOptimizerLM === - 2018-01-25T20:08:07: Optimized function value ..: -215505.826 - 2018-01-25T20:08:07: Absolute precision ........: 0.005 - 2018-01-25T20:08:07: Acceptable value decrease .: 2 - 2018-01-25T20:08:07: Optimization status .......: converged - 2018-01-25T20:08:07: Number of parameters ......: 14 - 2018-01-25T20:08:07: Number of free parameters .: 7 - 2018-01-25T20:08:07: Number of iterations ......: 19 - 2018-01-25T20:08:07: Lambda ....................: 1e-06 - 2018-01-25T20:08:07: Maximum log likelihood ....: 215505.826 - 2018-01-25T20:08:07: Observed events (Nobs) ...: 202330.000 - 2018-01-25T20:08:07: Predicted events (Npred) ..: 202318.925 (Nobs - Npred = 11.0749561789562) - 2018-01-25T20:08:07: === GModels === - 2018-01-25T20:08:07: Number of models ..........: 3 - 2018-01-25T20:08:07: Number of parameters ......: 14 - 2018-01-25T20:08:07: === GModelSky === - 2018-01-25T20:08:07: Name ......................: Vela - 2018-01-25T20:08:07: Instruments ...............: all - 2018-01-25T20:08:07: Instrument scale factors ..: unity - 2018-01-25T20:08:07: Observation identifiers ...: all - 2018-01-25T20:08:07: Model type ................: PointSource - 2018-01-25T20:08:07: Model components ..........: "PointSource" * "SuperExponentialCutoffPowerLaw" * "Constant" - 2018-01-25T20:08:07: Number of parameters ......: 8 - 2018-01-25T20:08:07: Number of spatial par's ...: 2 - 2018-01-25T20:08:07: RA .......................: 128.835772118429 +/- 0.00162177813558897 [-360,360] deg (free,scale=1) - 2018-01-25T20:08:07: DEC ......................: -45.1835812138884 +/- 0.00114069728792873 [-90,90] deg (free,scale=1) - 2018-01-25T20:08:07: Number of spectral par's ..: 5 - 2018-01-25T20:08:07: Prefactor ................: 4.35032346821964e-09 +/- 4.73845090178596e-10 [1e-16,1e-06] ph/cm2/s/MeV (free,scale=1e-09,gradient) - 2018-01-25T20:08:07: Index1 ...................: -1.34367731366247 +/- 0.031088811599731 [-0,-5] (free,scale=-1,gradient) - 2018-01-25T20:08:07: CutoffEnergy .............: 990.275379109663 +/- 181.814299583364 [1,1000000] MeV (free,scale=1000,gradient) - 2018-01-25T20:08:07: PivotEnergy ..............: 1000 [10,1000000] MeV (fixed,scale=1000,gradient) - 2018-01-25T20:08:07: Index2 ...................: 0.587550117179991 +/- 0.0279192539614008 [0.1,5] (free,scale=1,gradient) - 2018-01-25T20:08:07: Number of temporal par's ..: 1 - 2018-01-25T20:08:07: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) - 2018-01-25T20:08:07: === GModelSky === - 2018-01-25T20:08:07: Name ......................: Galactic_diffuse - 2018-01-25T20:08:07: Instruments ...............: all - 2018-01-25T20:08:07: Instrument scale factors ..: unity - 2018-01-25T20:08:07: Observation identifiers ...: all - 2018-01-25T20:08:07: Model type ................: DiffuseSource - 2018-01-25T20:08:07: Model components ..........: "DiffuseMapCube" * "Constant" * "Constant" - 2018-01-25T20:08:07: Number of parameters ......: 3 - 2018-01-25T20:08:07: Number of spatial par's ...: 1 - 2018-01-25T20:08:07: Normalization ............: 1 [0.1,10] (fixed,scale=1,gradient) - 2018-01-25T20:08:07: Number of spectral par's ..: 1 - 2018-01-25T20:08:07: Normalization ............: 1.10941973066973 +/- 0.00503412002598653 [0.1,1000] (free,scale=1,gradient) - 2018-01-25T20:08:07: Number of temporal par's ..: 1 - 2018-01-25T20:08:07: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) - 2018-01-25T20:08:07: === GModelSky === - 2018-01-25T20:08:07: Name ......................: Extragalactic_diffuse - 2018-01-25T20:08:07: Instruments ...............: all - 2018-01-25T20:08:07: Instrument scale factors ..: unity - 2018-01-25T20:08:07: Observation identifiers ...: all - 2018-01-25T20:08:07: Model type ................: DiffuseSource - 2018-01-25T20:08:07: Model components ..........: "DiffuseIsotropic" * "FileFunction" * "Constant" - 2018-01-25T20:08:07: Number of parameters ......: 3 - 2018-01-25T20:08:07: Number of spatial par's ...: 1 - 2018-01-25T20:08:07: Value ....................: 1 [0,10] (fixed,scale=1,gradient) - 2018-01-25T20:08:07: Number of spectral par's ..: 1 - 2018-01-25T20:08:07: Normalization ............: 1 [0,1000] (fixed,scale=1,gradient) - 2018-01-25T20:08:07: Number of temporal par's ..: 1 - 2018-01-25T20:08:07: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) + 2019-04-04T14:34:18: +=================================+ + 2019-04-04T14:34:18: | Maximum likelihood optimisation | + 2019-04-04T14:34:18: +=================================+ + 2019-04-04T14:34:20: >Iteration 0: -logL=-177232.327, Lambda=1.0e-03 + 2019-04-04T14:34:20: >Iteration 1: -logL=-198931.136, Lambda=1.0e-03, delta=21698.809, step=1.0e+00, max(|grad|)=-22717.063466 [CutoffEnergy:4] + 2019-04-04T14:34:21: >Iteration 2: -logL=-200269.888, Lambda=1.0e-04, delta=1338.753, step=1.0e+00, max(|grad|)=-67629.683813 [CutoffEnergy:4] + 2019-04-04T14:34:21: >Iteration 3: -logL=-212351.965, Lambda=1.0e-05, delta=12082.076, step=1.0e+00, max(|grad|)=23026.802441 [Index2:6] + 2019-04-04T14:34:22: >Iteration 4: -logL=-213954.259, Lambda=1.0e-06, delta=1602.294, step=1.0e+00, max(|grad|)=-8297.380921 [CutoffEnergy:4] + 2019-04-04T14:34:22: >Iteration 5: -logL=-215347.973, Lambda=1.0e-07, delta=1393.714, step=1.0e+00, max(|grad|)=2746.754045 [Index2:6] + 2019-04-04T14:34:23: >Iteration 6: -logL=-215497.807, Lambda=1.0e-08, delta=149.834, step=1.0e+00, max(|grad|)=699.848487 [Index2:6] + 2019-04-04T14:34:23: Iteration 7: -logL=-215497.807, Lambda=1.0e-09, delta=-13.845, step=1.0e+00, max(|grad|)=-1201.982876 [Index1:3] (stalled) + 2019-04-04T14:34:24: Iteration 8: -logL=-215497.807, Lambda=1.0e-08, delta=-13.844, step=1.0e+00, max(|grad|)=-1201.959648 [Index1:3] (stalled) + 2019-04-04T14:34:24: Iteration 9: -logL=-215497.807, Lambda=1.0e-07, delta=-13.835, step=1.0e+00, max(|grad|)=-1201.727402 [Index1:3] (stalled) + 2019-04-04T14:34:24: Iteration 10: -logL=-215497.807, Lambda=1.0e-06, delta=-13.744, step=1.0e+00, max(|grad|)=-1199.408718 [Index1:3] (stalled) + 2019-04-04T14:34:25: Iteration 11: -logL=-215497.807, Lambda=1.0e-05, delta=-12.858, step=1.0e+00, max(|grad|)=-1176.594022 [Index1:3] (stalled) + 2019-04-04T14:34:26: Iteration 12: -logL=-215497.807, Lambda=1.0e-04, delta=-6.126, step=1.0e+00, max(|grad|)=-980.833964 [Index1:3] (stalled) + 2019-04-04T14:34:26: >Iteration 13: -logL=-215503.803, Lambda=1.0e-03, delta=5.995, step=1.0e+00, max(|grad|)=-285.964915 [Index1:3] + 2019-04-04T14:34:27: Iteration 14: -logL=-215503.803, Lambda=1.0e-04, delta=-8.762, step=1.0e+00, max(|grad|)=-730.461327 [Index1:3] (stalled) + 2019-04-04T14:34:27: >Iteration 15: -logL=-215504.976, Lambda=1.0e-03, delta=1.173, step=1.0e+00, max(|grad|)=-153.941337 [Index1:3] + 2019-04-04T14:34:28: Iteration 16: -logL=-215503.593, Lambda=1.0e-04, delta=-1.383, step=1.0e+00, max(|grad|)=-327.368446 [Index1:3] (stalled) + 2019-04-04T14:34:28: >Iteration 17: -logL=-215505.817, Lambda=1.0e-03, delta=2.225, step=1.0e+00, max(|grad|)=3.256215 [Index2:6] + 2019-04-04T14:34:29: >Iteration 18: -logL=-215505.824, Lambda=1.0e-04, delta=0.007, step=1.0e+00, max(|grad|)=-8.190315 [CutoffEnergy:4] + 2019-04-04T14:34:29: >Iteration 19: -logL=-215505.826, Lambda=1.0e-05, delta=0.002, step=1.0e+00, max(|grad|)=0.635264 [Index2:6] + 2019-04-04T14:34:30: + 2019-04-04T14:34:30: +=========================================+ + 2019-04-04T14:34:30: | Maximum likelihood optimisation results | + 2019-04-04T14:34:30: +=========================================+ + 2019-04-04T14:34:30: === GOptimizerLM === + 2019-04-04T14:34:30: Optimized function value ..: -215505.826 + 2019-04-04T14:34:30: Absolute precision ........: 0.005 + 2019-04-04T14:34:30: Acceptable value decrease .: 2 + 2019-04-04T14:34:30: Optimization status .......: converged + 2019-04-04T14:34:30: Number of parameters ......: 14 + 2019-04-04T14:34:30: Number of free parameters .: 7 + 2019-04-04T14:34:30: Number of iterations ......: 19 + 2019-04-04T14:34:30: Lambda ....................: 1e-06 + 2019-04-04T14:34:30: Maximum log likelihood ....: 215505.826 + 2019-04-04T14:34:30: Observed events (Nobs) ...: 202330.000 + 2019-04-04T14:34:30: Predicted events (Npred) ..: 202318.925 (Nobs - Npred = 11.0749567862367) + 2019-04-04T14:34:30: === GModels === + 2019-04-04T14:34:30: Number of models ..........: 3 + 2019-04-04T14:34:30: Number of parameters ......: 14 + 2019-04-04T14:34:30: === GModelSky === + 2019-04-04T14:34:30: Name ......................: Vela + 2019-04-04T14:34:30: Instruments ...............: all + 2019-04-04T14:34:30: Instrument scale factors ..: unity + 2019-04-04T14:34:30: Observation identifiers ...: all + 2019-04-04T14:34:30: Model type ................: PointSource + 2019-04-04T14:34:30: Model components ..........: "PointSource" * "SuperExponentialCutoffPowerLaw" * "Constant" + 2019-04-04T14:34:30: Number of parameters ......: 8 + 2019-04-04T14:34:30: Number of spatial par's ...: 2 + 2019-04-04T14:34:30: RA .......................: 128.835772118406 +/- 0.00162177814111823 [-360,360] deg (free,scale=1) + 2019-04-04T14:34:30: DEC ......................: -45.1835812138759 +/- 0.00114069728372229 [-90,90] deg (free,scale=1) + 2019-04-04T14:34:30: Number of spectral par's ..: 5 + 2019-04-04T14:34:30: Prefactor ................: 4.35032346843698e-09 +/- 4.73845090232622e-10 [1e-16,1e-06] ph/cm2/s/MeV (free,scale=1e-09,gradient) + 2019-04-04T14:34:30: Index1 ...................: -1.34367731364905 +/- 0.0310888116010815 [-0,-5] (free,scale=-1,gradient) + 2019-04-04T14:34:30: CutoffEnergy .............: 990.275379021926 +/- 181.814299577731 [1,1000000] MeV (free,scale=1000,gradient) + 2019-04-04T14:34:30: PivotEnergy ..............: 1000 [10,1000000] MeV (fixed,scale=1000,gradient) + 2019-04-04T14:34:30: Index2 ...................: 0.58755011716699 +/- 0.027919253961318 [0.1,5] (free,scale=1,gradient) + 2019-04-04T14:34:30: Number of temporal par's ..: 1 + 2019-04-04T14:34:30: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) + 2019-04-04T14:34:30: === GModelSky === + 2019-04-04T14:34:30: Name ......................: Galactic_diffuse + 2019-04-04T14:34:30: Instruments ...............: all + 2019-04-04T14:34:30: Instrument scale factors ..: unity + 2019-04-04T14:34:30: Observation identifiers ...: all + 2019-04-04T14:34:30: Model type ................: DiffuseSource + 2019-04-04T14:34:30: Model components ..........: "DiffuseMapCube" * "Constant" * "Constant" + 2019-04-04T14:34:30: Number of parameters ......: 3 + 2019-04-04T14:34:30: Number of spatial par's ...: 1 + 2019-04-04T14:34:30: Normalization ............: 1 [0.1,10] (fixed,scale=1,gradient) + 2019-04-04T14:34:30: Number of spectral par's ..: 1 + 2019-04-04T14:34:30: Normalization ............: 1.10941973066989 +/- 0.00503412002600817 [0.1,1000] (free,scale=1,gradient) + 2019-04-04T14:34:30: Number of temporal par's ..: 1 + 2019-04-04T14:34:30: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) + 2019-04-04T14:34:30: === GModelSky === + 2019-04-04T14:34:30: Name ......................: Extragalactic_diffuse + 2019-04-04T14:34:30: Instruments ...............: all + 2019-04-04T14:34:30: Instrument scale factors ..: unity + 2019-04-04T14:34:30: Observation identifiers ...: all + 2019-04-04T14:34:30: Model type ................: DiffuseSource + 2019-04-04T14:34:30: Model components ..........: "DiffuseIsotropic" * "FileFunction" * "Constant" + 2019-04-04T14:34:30: Number of parameters ......: 3 + 2019-04-04T14:34:30: Number of spatial par's ...: 1 + 2019-04-04T14:34:30: Value ....................: 1 [0,10] (fixed,scale=1,gradient) + 2019-04-04T14:34:30: Number of spectral par's ..: 1 + 2019-04-04T14:34:30: Normalization ............: 1 [0,1000] (fixed,scale=1,gradient) + 2019-04-04T14:34:30: Number of temporal par's ..: 1 + 2019-04-04T14:34:30: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) The fit converged after 19 iterations with spectral parameters that are reasonably close to those found in @@ -241,9 +241,9 @@ using the Fermi-LAT Science Tools. .. code-block:: none - 2018-01-25T20:05:57: Name ......................: Vela - 2018-01-25T20:05:57: Instruments ...............: all - 2018-01-25T20:05:57: Test Statistic ............: 326341.439917148 + 2019-04-04T14:37:50: Name ......................: Vela + 2019-04-04T14:37:50: Instruments ...............: all + 2019-04-04T14:37:50: Test Statistic ............: 326341.43991719 which is also close to the value determined by gtlike. diff --git a/doc/source/users/tutorials/howto/fermi/howto_fermi_spectrum.png b/doc/source/users/tutorials/howto/fermi/howto_fermi_spectrum.png index 529f6c8..9e1aca0 100755 Binary files a/doc/source/users/tutorials/howto/fermi/howto_fermi_spectrum.png and b/doc/source/users/tutorials/howto/fermi/howto_fermi_spectrum.png differ diff --git a/doc/source/users/tutorials/howto/fermi/howto_fermi_tsmap.rst b/doc/source/users/tutorials/howto/fermi/howto_fermi_tsmap.rst index f7a795d..853f373 100644 --- a/doc/source/users/tutorials/howto/fermi/howto_fermi_tsmap.rst +++ b/doc/source/users/tutorials/howto/fermi/howto_fermi_tsmap.rst @@ -20,14 +20,14 @@ source. You do this by using the :ref:`cttsmap` tool as follows: Input event list, counts cube or observation definition XML file [events.fits] obs.xml Test source name [Crab] Vela Input model definition XML file [$CTOOLS/share/models/crab.xml] vela_results.xml + Coordinate system (CEL - celestial, GAL - galactic) (CEL|GAL) [CEL] + Projection method (AIT|AZP|CAR|GLS|MER|MOL|SFL|SIN|STG|TAN) [CAR] First coordinate of image center in degrees (RA or galactic l) (0-360) [83.63] 128.84 Second coordinate of image center in degrees (DEC or galactic b) (-90-90) [22.01] -45.18 - Projection method (AIT|AZP|CAR|GLS|MER|MOL|SFL|SIN|STG|TAN) [CAR] - Coordinate system (CEL - celestial, GAL - galactic) (CEL|GAL) [CEL] Image scale (in degrees/pixel) [0.02] 0.3 Size of the X axis in pixels [200] 10 Size of the Y axis in pixels [200] 10 - Output Test Statistic map file [tsmap.fits] + Output Test Statistic map file [tsmap.fits] :ref:`cttsmap` created a 10 x 10 pixels map with a pixel size of 0.3 degrees around the location of the Vela pulsar and saved the result in the FITS file diff --git a/doc/source/users/tutorials/howto/fermi/howto_fermi_ulimit.rst b/doc/source/users/tutorials/howto/fermi/howto_fermi_ulimit.rst index 037f4c5..858a172 100755 --- a/doc/source/users/tutorials/howto/fermi/howto_fermi_ulimit.rst +++ b/doc/source/users/tutorials/howto/fermi/howto_fermi_ulimit.rst @@ -20,7 +20,7 @@ for the source. To do this you first have to add an additional source to the - + @@ -53,9 +53,9 @@ The results of the upper limit computation can then be extracted from the .. code-block:: bash - 2018-01-25T20:43:02: +=====================+ - 2018-01-25T20:43:02: | Upper limit results | - 2018-01-25T20:43:02: +=====================+ - 2018-01-25T20:43:02: Differential flux limit ...: 1.00100099727198e-13 ph/cm2/s/MeV at 0.001 TeV - 2018-01-25T20:43:02: Integral flux limit .......: 9.9999999627471e-10 ph/cm2/s within [0.0001-0.1] TeV - 2018-01-25T20:43:02: Energy flux limit .........: 1.10785216561628e-12 erg/cm2/s within [0.0001-0.1] TeV + 2019-04-04T15:03:22: +=====================+ + 2019-04-04T15:03:22: | Upper limit results | + 2019-04-04T15:03:22: +=====================+ + 2019-04-04T15:03:22: Differential flux limit ...: 1.77270316013682e-13 ph/cm2/s/MeV at 0.001 TeV + 2019-04-04T15:03:22: Integral flux limit .......: 1.77093045697669e-09 ph/cm2/s within [0.0001-0.1] TeV + 2019-04-04T15:03:22: Energy flux limit .........: 1.96192914922622e-12 erg/cm2/s within [0.0001-0.1] TeV -- libgit2 0.21.4