Nathan Kelley-Hoskins / ctools

Commit 0873cb8eebcf86edcae88001e5b555f059194b9b

Authored by Jürgen Knödlseder
1 parent 19ce1581

Update Fermi How To

Inline Side-by-side
Showing 6 changed files with 179 additions and 179 deletions
doc/source/users/tutorials/howto/fermi/howto_fermi_butterfly.png
View file @0873cb8

40.5 KB | W: | H:

77.1 KB | W: | H:

  • 2-up
  • Swipe
  • Onion skin
doc/source/users/tutorials/howto/fermi/howto_fermi_errors.rst
  View file @0873cb8
@@ -33,80 +33,80 @@ covariance matrix. @@ -33,80 +33,80 @@ covariance matrix.
33 33
34 .. code-block:: none 34 .. code-block:: none
35 35
36 - 2018-01-25T20:45:22: +================================================+  
37 - 2018-01-25T20:45:22: | Compute error for source "Vela" parameter "RA" |  
38 - 2018-01-25T20:45:22: +================================================+  
39 - 2018-01-25T20:45:22: Confidence level ..........: 68 %  
40 - 2018-01-25T20:45:22: Log-likelihood difference .: 0.494473  
41 - 2018-01-25T20:45:22: Initial factor range ......: [128.819554251926, 128.851989886434]  
42 - 2018-01-25T20:45:41: Lower parameter factor ....: 128.834  
43 - 2018-01-25T20:45:41: Upper parameter factor ....: 128.837  
44 - 2018-01-25T20:45:41: Error from curvature ......: 0.00162178172537938 deg  
45 - 2018-01-25T20:45:41: Error from profile ........: 0.00169067577327553 deg  
46 - 2018-01-25T20:45:41: Negative profile error ....: 0.00165900034895117 deg  
47 - 2018-01-25T20:45:41: Positive profile error ....: 0.00172235119759989 deg  
48 - 2018-01-25T20:45:41:  
49 - 2018-01-25T20:45:41: +=================================================+  
50 - 2018-01-25T20:45:41: | Compute error for source "Vela" parameter "DEC" |  
51 - 2018-01-25T20:45:41: +=================================================+  
52 - 2018-01-25T20:45:41: Confidence level ..........: 68 %  
53 - 2018-01-25T20:45:41: Log-likelihood difference .: 0.494473  
54 - 2018-01-25T20:45:41: Initial factor range ......: [-45.1949881923928, -45.1721741961801]  
55 - 2018-01-25T20:46:01: Lower parameter factor ....: -45.1847  
56 - 2018-01-25T20:46:01: Upper parameter factor ....: -45.1824  
57 - 2018-01-25T20:46:01: Error from curvature ......: 0.00114069981063169 deg  
58 - 2018-01-25T20:46:01: Error from profile ........: 0.00116548552428952 deg  
59 - 2018-01-25T20:46:01: Negative profile error ....: 0.00115852324517363 deg  
60 - 2018-01-25T20:46:01: Positive profile error ....: 0.00117244780340542 deg  
61 - 2018-01-25T20:46:01:  
62 - 2018-01-25T20:46:01: +=======================================================+  
63 - 2018-01-25T20:46:01: | Compute error for source "Vela" parameter "Prefactor" |  
64 - 2018-01-25T20:46:01: +=======================================================+  
65 - 2018-01-25T20:46:01: Confidence level ..........: 68 %  
66 - 2018-01-25T20:46:01: Log-likelihood difference .: 0.494473  
67 - 2018-01-25T20:46:01: Initial factor range ......: [1e-07, 9.09784331905335]  
68 - 2018-01-25T20:46:45: Lower parameter factor ....: 3.93674  
69 - 2018-01-25T20:46:45: Upper parameter factor ....: 4.89558  
70 - 2018-01-25T20:46:45: Error from curvature ......: 4.74564602380296e-10 ph/cm2/s/MeV  
71 - 2018-01-25T20:46:45: Error from profile ........: 4.79421031311828e-10 ph/cm2/s/MeV  
72 - 2018-01-25T20:46:45: Negative profile error ....: 4.15456324058327e-10 ph/cm2/s/MeV  
73 - 2018-01-25T20:46:45: Positive profile error ....: 5.4338573856533e-10 ph/cm2/s/MeV  
74 - 2018-01-25T20:46:45:  
75 - 2018-01-25T20:46:45: +====================================================+  
76 - 2018-01-25T20:46:45: | Compute error for source "Vela" parameter "Index1" |  
77 - 2018-01-25T20:46:45: +====================================================+  
78 - 2018-01-25T20:46:45: Confidence level ..........: 68 %  
79 - 2018-01-25T20:46:45: Log-likelihood difference .: 0.494473  
80 - 2018-01-25T20:46:45: Initial factor range ......: [1.0323933217421, 1.65461803910864]  
81 - 2018-01-25T20:47:31: Lower parameter factor ....: 1.31084  
82 - 2018-01-25T20:47:31: Upper parameter factor ....: 1.37275  
83 - 2018-01-25T20:47:31: Error from curvature ......: 0.0311112358683271  
84 - 2018-01-25T20:47:31: Error from profile ........: 0.0309517300203693  
85 - 2018-01-25T20:47:31: Negative profile error ....: 0.032660721248488  
86 - 2018-01-25T20:47:31: Positive profile error ....: 0.0292427387922507  
87 - 2018-01-25T20:47:31:  
88 - 2018-01-25T20:47:31: +==========================================================+  
89 - 2018-01-25T20:47:31: | Compute error for source "Vela" parameter "CutoffEnergy" |  
90 - 2018-01-25T20:47:31: +==========================================================+  
91 - 2018-01-25T20:47:31: Confidence level ..........: 68 %  
92 - 2018-01-25T20:47:31: Log-likelihood difference .: 0.494473  
93 - 2018-01-25T20:47:31: Initial factor range ......: [0.001, 2.80853448949703]  
94 - 2018-01-25T20:48:14: Lower parameter factor ....: 0.81239  
95 - 2018-01-25T20:48:14: Upper parameter factor ....: 1.17294  
96 - 2018-01-25T20:48:14: Error from curvature ......: 181.899981380157 MeV  
97 - 2018-01-25T20:48:14: Error from profile ........: 180.277267463442 MeV  
98 - 2018-01-25T20:48:14: Negative profile error ....: 177.144641591912 MeV  
99 - 2018-01-25T20:48:14: Positive profile error ....: 183.409893334972 MeV  
100 - 2018-01-25T20:48:14:  
101 - 2018-01-25T20:48:14: +====================================================+  
102 - 2018-01-25T20:48:14: | Compute error for source "Vela" parameter "Index2" |  
103 - 2018-01-25T20:48:14: +====================================================+  
104 - 2018-01-25T20:48:14: Confidence level ..........: 68 %  
105 - 2018-01-25T20:48:14: Log-likelihood difference .: 0.494473  
106 - 2018-01-25T20:48:14: Initial factor range ......: [0.308028525758507, 0.866644263529458]  
107 - 2018-01-25T20:49:01: Lower parameter factor ....: 0.559037  
108 - 2018-01-25T20:49:01: Upper parameter factor ....: 0.615056  
109 - 2018-01-25T20:49:01: Error from curvature ......: 0.0279307868885475  
110 - 2018-01-25T20:49:01: Error from profile ........: 0.0280092058458762  
111 - 2018-01-25T20:49:01: Negative profile error ....: 0.028299015036004  
112 - 2018-01-25T20:49:01: Positive profile error ....: 0.0277193966557485 36 + 2019-04-04T15:04:51: +================================================+
  37 + 2019-04-04T15:04:51: | Compute error for source "Vela" parameter "RA" |
  38 + 2019-04-04T15:04:51: +================================================+
  39 + 2019-04-04T15:04:51: Confidence level ..........: 68 %
  40 + 2019-04-04T15:04:51: Log-likelihood difference .: 0.494473
  41 + 2019-04-04T15:04:51: Initial factor range ......: [128.819554251886, 128.851989886469]
  42 + 2019-04-04T15:05:10: Lower parameter factor ....: 128.834
  43 + 2019-04-04T15:05:10: Upper parameter factor ....: 128.837
  44 + 2019-04-04T15:05:10: Error from curvature ......: 0.00162178172918532 deg
  45 + 2019-04-04T15:05:10: Error from profile ........: 0.00169067577726878 deg
  46 + 2019-04-04T15:05:10: Negative profile error ....: 0.00165900035287336 deg
  47 + 2019-04-04T15:05:10: Positive profile error ....: 0.0017223512016642 deg
  48 + 2019-04-04T15:05:10:
  49 + 2019-04-04T15:05:10: +=================================================+
  50 + 2019-04-04T15:05:10: | Compute error for source "Vela" parameter "DEC" |
  51 + 2019-04-04T15:05:10: +=================================================+
  52 + 2019-04-04T15:05:10: Confidence level ..........: 68 %
  53 + 2019-04-04T15:05:10: Log-likelihood difference .: 0.494473
  54 + 2019-04-04T15:05:10: Initial factor range ......: [-45.1949881925207, -45.1721741960936]
  55 + 2019-04-04T15:05:31: Lower parameter factor ....: -45.1847
  56 + 2019-04-04T15:05:31: Upper parameter factor ....: -45.1824
  57 + 2019-04-04T15:05:31: Error from curvature ......: 0.00114069982135549 deg
  58 + 2019-04-04T15:05:31: Error from profile ........: 0.0011654855352532 deg
  59 + 2019-04-04T15:05:31: Negative profile error ....: 0.00115852325606625 deg
  60 + 2019-04-04T15:05:31: Positive profile error ....: 0.00117244781444015 deg
  61 + 2019-04-04T15:05:31:
  62 + 2019-04-04T15:05:31: +=======================================================+
  63 + 2019-04-04T15:05:31: | Compute error for source "Vela" parameter "Prefactor" |
  64 + 2019-04-04T15:05:31: +=======================================================+
  65 + 2019-04-04T15:05:31: Confidence level ..........: 68 %
  66 + 2019-04-04T15:05:31: Log-likelihood difference .: 0.494473
  67 + 2019-04-04T15:05:31: Initial factor range ......: [1e-07, 9.09784331943508]
  68 + 2019-04-04T15:06:15: Lower parameter factor ....: 3.93674
  69 + 2019-04-04T15:06:15: Upper parameter factor ....: 4.89558
  70 + 2019-04-04T15:06:15: Error from curvature ......: 4.74564602405264e-10 ph/cm2/s/MeV
  71 + 2019-04-04T15:06:15: Error from profile ........: 4.79421031332425e-10 ph/cm2/s/MeV
  72 + 2019-04-04T15:06:15: Negative profile error ....: 4.15456324070932e-10 ph/cm2/s/MeV
  73 + 2019-04-04T15:06:15: Positive profile error ....: 5.43385738593917e-10 ph/cm2/s/MeV
  74 + 2019-04-04T15:06:15:
  75 + 2019-04-04T15:06:15: +====================================================+
  76 + 2019-04-04T15:06:15: | Compute error for source "Vela" parameter "Index1" |
  77 + 2019-04-04T15:06:15: +====================================================+
  78 + 2019-04-04T15:06:15: Confidence level ..........: 68 %
  79 + 2019-04-04T15:06:15: Log-likelihood difference .: 0.494473
  80 + 2019-04-04T15:06:15: Initial factor range ......: [1.03239332172832, 1.65461803910225]
  81 + 2019-04-04T15:07:01: Lower parameter factor ....: 1.31084
  82 + 2019-04-04T15:07:01: Upper parameter factor ....: 1.37275
  83 + 2019-04-04T15:07:01: Error from curvature ......: 0.0311112358686962
  84 + 2019-04-04T15:07:01: Error from profile ........: 0.0309517300207367
  85 + 2019-04-04T15:07:01: Negative profile error ....: 0.0326607212488756
  86 + 2019-04-04T15:07:01: Positive profile error ....: 0.0292427387925978
  87 + 2019-04-04T15:07:01:
  88 + 2019-04-04T15:07:01: +==========================================================+
  89 + 2019-04-04T15:07:01: | Compute error for source "Vela" parameter "CutoffEnergy" |
  90 + 2019-04-04T15:07:01: +==========================================================+
  91 + 2019-04-04T15:07:01: Confidence level ..........: 68 %
  92 + 2019-04-04T15:07:01: Log-likelihood difference .: 0.494473
  93 + 2019-04-04T15:07:01: Initial factor range ......: [0.001, 2.80853448937906]
  94 + 2019-04-04T15:07:45: Lower parameter factor ....: 0.81239
  95 + 2019-04-04T15:07:45: Upper parameter factor ....: 1.17294
  96 + 2019-04-04T15:07:45: Error from curvature ......: 181.899981373502 MeV
  97 + 2019-04-04T15:07:45: Error from profile ........: 180.277267455479 MeV
  98 + 2019-04-04T15:07:45: Negative profile error ....: 177.144641582696 MeV
  99 + 2019-04-04T15:07:45: Positive profile error ....: 183.409893328263 MeV
  100 + 2019-04-04T15:07:45:
  101 + 2019-04-04T15:07:45: +====================================================+
  102 + 2019-04-04T15:07:45: | Compute error for source "Vela" parameter "Index2" |
  103 + 2019-04-04T15:07:45: +====================================================+
  104 + 2019-04-04T15:07:45: Confidence level ..........: 68 %
  105 + 2019-04-04T15:07:45: Log-likelihood difference .: 0.494473
  106 + 2019-04-04T15:07:45: Initial factor range ......: [0.3080285257564, 0.866644263519715]
  107 + 2019-04-04T15:08:32: Lower parameter factor ....: 0.559037
  108 + 2019-04-04T15:08:32: Upper parameter factor ....: 0.615056
  109 + 2019-04-04T15:08:32: Error from curvature ......: 0.0279307868881658
  110 + 2019-04-04T15:08:32: Error from profile ........: 0.0280092058454933
  111 + 2019-04-04T15:08:32: Negative profile error ....: 0.028299015035617
  112 + 2019-04-04T15:08:32: Positive profile error ....: 0.0277193966553695
doc/source/users/tutorials/howto/fermi/howto_fermi_fitting.rst
  View file @0873cb8
@@ -92,95 +92,95 @@ log file ``ctlike.log`` that was also created by :ref:`ctlike`: @@ -92,95 +92,95 @@ log file ``ctlike.log`` that was also created by :ref:`ctlike`:
92 92
93 .. code-block:: none 93 .. code-block:: none
94 94
95 - 2018-01-25T20:07:55: +=================================+  
96 - 2018-01-25T20:07:55: | Maximum likelihood optimisation |  
97 - 2018-01-25T20:07:55: +=================================+  
98 - 2018-01-25T20:07:57: >Iteration 0: -logL=-177232.327, Lambda=1.0e-03  
99 - 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]  
100 - 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]  
101 - 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]  
102 - 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]  
103 - 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]  
104 - 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]  
105 - 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)  
106 - 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)  
107 - 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)  
108 - 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)  
109 - 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)  
110 - 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)  
111 - 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]  
112 - 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)  
113 - 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]  
114 - 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)  
115 - 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]  
116 - 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]  
117 - 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]  
118 - 2018-01-25T20:08:07:  
119 - 2018-01-25T20:08:07: +=========================================+  
120 - 2018-01-25T20:08:07: | Maximum likelihood optimisation results |  
121 - 2018-01-25T20:08:07: +=========================================+  
122 - 2018-01-25T20:08:07: === GOptimizerLM ===  
123 - 2018-01-25T20:08:07: Optimized function value ..: -215505.826  
124 - 2018-01-25T20:08:07: Absolute precision ........: 0.005  
125 - 2018-01-25T20:08:07: Acceptable value decrease .: 2  
126 - 2018-01-25T20:08:07: Optimization status .......: converged  
127 - 2018-01-25T20:08:07: Number of parameters ......: 14  
128 - 2018-01-25T20:08:07: Number of free parameters .: 7  
129 - 2018-01-25T20:08:07: Number of iterations ......: 19  
130 - 2018-01-25T20:08:07: Lambda ....................: 1e-06  
131 - 2018-01-25T20:08:07: Maximum log likelihood ....: 215505.826  
132 - 2018-01-25T20:08:07: Observed events (Nobs) ...: 202330.000  
133 - 2018-01-25T20:08:07: Predicted events (Npred) ..: 202318.925 (Nobs - Npred = 11.0749561789562)  
134 - 2018-01-25T20:08:07: === GModels ===  
135 - 2018-01-25T20:08:07: Number of models ..........: 3  
136 - 2018-01-25T20:08:07: Number of parameters ......: 14  
137 - 2018-01-25T20:08:07: === GModelSky ===  
138 - 2018-01-25T20:08:07: Name ......................: Vela  
139 - 2018-01-25T20:08:07: Instruments ...............: all  
140 - 2018-01-25T20:08:07: Instrument scale factors ..: unity  
141 - 2018-01-25T20:08:07: Observation identifiers ...: all  
142 - 2018-01-25T20:08:07: Model type ................: PointSource  
143 - 2018-01-25T20:08:07: Model components ..........: "PointSource" * "SuperExponentialCutoffPowerLaw" * "Constant"  
144 - 2018-01-25T20:08:07: Number of parameters ......: 8  
145 - 2018-01-25T20:08:07: Number of spatial par's ...: 2  
146 - 2018-01-25T20:08:07: RA .......................: 128.835772118429 +/- 0.00162177813558897 [-360,360] deg (free,scale=1)  
147 - 2018-01-25T20:08:07: DEC ......................: -45.1835812138884 +/- 0.00114069728792873 [-90,90] deg (free,scale=1)  
148 - 2018-01-25T20:08:07: Number of spectral par's ..: 5  
149 - 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)  
150 - 2018-01-25T20:08:07: Index1 ...................: -1.34367731366247 +/- 0.031088811599731 [-0,-5] (free,scale=-1,gradient)  
151 - 2018-01-25T20:08:07: CutoffEnergy .............: 990.275379109663 +/- 181.814299583364 [1,1000000] MeV (free,scale=1000,gradient)  
152 - 2018-01-25T20:08:07: PivotEnergy ..............: 1000 [10,1000000] MeV (fixed,scale=1000,gradient)  
153 - 2018-01-25T20:08:07: Index2 ...................: 0.587550117179991 +/- 0.0279192539614008 [0.1,5] (free,scale=1,gradient)  
154 - 2018-01-25T20:08:07: Number of temporal par's ..: 1  
155 - 2018-01-25T20:08:07: Normalization ............: 1 (relative value) (fixed,scale=1,gradient)  
156 - 2018-01-25T20:08:07: === GModelSky ===  
157 - 2018-01-25T20:08:07: Name ......................: Galactic_diffuse  
158 - 2018-01-25T20:08:07: Instruments ...............: all  
159 - 2018-01-25T20:08:07: Instrument scale factors ..: unity  
160 - 2018-01-25T20:08:07: Observation identifiers ...: all  
161 - 2018-01-25T20:08:07: Model type ................: DiffuseSource  
162 - 2018-01-25T20:08:07: Model components ..........: "DiffuseMapCube" * "Constant" * "Constant"  
163 - 2018-01-25T20:08:07: Number of parameters ......: 3  
164 - 2018-01-25T20:08:07: Number of spatial par's ...: 1  
165 - 2018-01-25T20:08:07: Normalization ............: 1 [0.1,10] (fixed,scale=1,gradient)  
166 - 2018-01-25T20:08:07: Number of spectral par's ..: 1  
167 - 2018-01-25T20:08:07: Normalization ............: 1.10941973066973 +/- 0.00503412002598653 [0.1,1000] (free,scale=1,gradient)  
168 - 2018-01-25T20:08:07: Number of temporal par's ..: 1  
169 - 2018-01-25T20:08:07: Normalization ............: 1 (relative value) (fixed,scale=1,gradient)  
170 - 2018-01-25T20:08:07: === GModelSky ===  
171 - 2018-01-25T20:08:07: Name ......................: Extragalactic_diffuse  
172 - 2018-01-25T20:08:07: Instruments ...............: all  
173 - 2018-01-25T20:08:07: Instrument scale factors ..: unity  
174 - 2018-01-25T20:08:07: Observation identifiers ...: all  
175 - 2018-01-25T20:08:07: Model type ................: DiffuseSource  
176 - 2018-01-25T20:08:07: Model components ..........: "DiffuseIsotropic" * "FileFunction" * "Constant"  
177 - 2018-01-25T20:08:07: Number of parameters ......: 3  
178 - 2018-01-25T20:08:07: Number of spatial par's ...: 1  
179 - 2018-01-25T20:08:07: Value ....................: 1 [0,10] (fixed,scale=1,gradient)  
180 - 2018-01-25T20:08:07: Number of spectral par's ..: 1  
181 - 2018-01-25T20:08:07: Normalization ............: 1 [0,1000] (fixed,scale=1,gradient)  
182 - 2018-01-25T20:08:07: Number of temporal par's ..: 1  
183 - 2018-01-25T20:08:07: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) 95 + 2019-04-04T14:34:18: +=================================+
  96 + 2019-04-04T14:34:18: | Maximum likelihood optimisation |
  97 + 2019-04-04T14:34:18: +=================================+
  98 + 2019-04-04T14:34:20: >Iteration 0: -logL=-177232.327, Lambda=1.0e-03
  99 + 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]
  100 + 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]
  101 + 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]
  102 + 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]
  103 + 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]
  104 + 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]
  105 + 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)
  106 + 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)
  107 + 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)
  108 + 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)
  109 + 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)
  110 + 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)
  111 + 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]
  112 + 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)
  113 + 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]
  114 + 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)
  115 + 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]
  116 + 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]
  117 + 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]
  118 + 2019-04-04T14:34:30:
  119 + 2019-04-04T14:34:30: +=========================================+
  120 + 2019-04-04T14:34:30: | Maximum likelihood optimisation results |
  121 + 2019-04-04T14:34:30: +=========================================+
  122 + 2019-04-04T14:34:30: === GOptimizerLM ===
  123 + 2019-04-04T14:34:30: Optimized function value ..: -215505.826
  124 + 2019-04-04T14:34:30: Absolute precision ........: 0.005
  125 + 2019-04-04T14:34:30: Acceptable value decrease .: 2
  126 + 2019-04-04T14:34:30: Optimization status .......: converged
  127 + 2019-04-04T14:34:30: Number of parameters ......: 14
  128 + 2019-04-04T14:34:30: Number of free parameters .: 7
  129 + 2019-04-04T14:34:30: Number of iterations ......: 19
  130 + 2019-04-04T14:34:30: Lambda ....................: 1e-06
  131 + 2019-04-04T14:34:30: Maximum log likelihood ....: 215505.826
  132 + 2019-04-04T14:34:30: Observed events (Nobs) ...: 202330.000
  133 + 2019-04-04T14:34:30: Predicted events (Npred) ..: 202318.925 (Nobs - Npred = 11.0749567862367)
  134 + 2019-04-04T14:34:30: === GModels ===
  135 + 2019-04-04T14:34:30: Number of models ..........: 3
  136 + 2019-04-04T14:34:30: Number of parameters ......: 14
  137 + 2019-04-04T14:34:30: === GModelSky ===
  138 + 2019-04-04T14:34:30: Name ......................: Vela
  139 + 2019-04-04T14:34:30: Instruments ...............: all
  140 + 2019-04-04T14:34:30: Instrument scale factors ..: unity
  141 + 2019-04-04T14:34:30: Observation identifiers ...: all
  142 + 2019-04-04T14:34:30: Model type ................: PointSource
  143 + 2019-04-04T14:34:30: Model components ..........: "PointSource" * "SuperExponentialCutoffPowerLaw" * "Constant"
  144 + 2019-04-04T14:34:30: Number of parameters ......: 8
  145 + 2019-04-04T14:34:30: Number of spatial par's ...: 2
  146 + 2019-04-04T14:34:30: RA .......................: 128.835772118406 +/- 0.00162177814111823 [-360,360] deg (free,scale=1)
  147 + 2019-04-04T14:34:30: DEC ......................: -45.1835812138759 +/- 0.00114069728372229 [-90,90] deg (free,scale=1)
  148 + 2019-04-04T14:34:30: Number of spectral par's ..: 5
  149 + 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)
  150 + 2019-04-04T14:34:30: Index1 ...................: -1.34367731364905 +/- 0.0310888116010815 [-0,-5] (free,scale=-1,gradient)
  151 + 2019-04-04T14:34:30: CutoffEnergy .............: 990.275379021926 +/- 181.814299577731 [1,1000000] MeV (free,scale=1000,gradient)
  152 + 2019-04-04T14:34:30: PivotEnergy ..............: 1000 [10,1000000] MeV (fixed,scale=1000,gradient)
  153 + 2019-04-04T14:34:30: Index2 ...................: 0.58755011716699 +/- 0.027919253961318 [0.1,5] (free,scale=1,gradient)
  154 + 2019-04-04T14:34:30: Number of temporal par's ..: 1
  155 + 2019-04-04T14:34:30: Normalization ............: 1 (relative value) (fixed,scale=1,gradient)
  156 + 2019-04-04T14:34:30: === GModelSky ===
  157 + 2019-04-04T14:34:30: Name ......................: Galactic_diffuse
  158 + 2019-04-04T14:34:30: Instruments ...............: all
  159 + 2019-04-04T14:34:30: Instrument scale factors ..: unity
  160 + 2019-04-04T14:34:30: Observation identifiers ...: all
  161 + 2019-04-04T14:34:30: Model type ................: DiffuseSource
  162 + 2019-04-04T14:34:30: Model components ..........: "DiffuseMapCube" * "Constant" * "Constant"
  163 + 2019-04-04T14:34:30: Number of parameters ......: 3
  164 + 2019-04-04T14:34:30: Number of spatial par's ...: 1
  165 + 2019-04-04T14:34:30: Normalization ............: 1 [0.1,10] (fixed,scale=1,gradient)
  166 + 2019-04-04T14:34:30: Number of spectral par's ..: 1
  167 + 2019-04-04T14:34:30: Normalization ............: 1.10941973066989 +/- 0.00503412002600817 [0.1,1000] (free,scale=1,gradient)
  168 + 2019-04-04T14:34:30: Number of temporal par's ..: 1
  169 + 2019-04-04T14:34:30: Normalization ............: 1 (relative value) (fixed,scale=1,gradient)
  170 + 2019-04-04T14:34:30: === GModelSky ===
  171 + 2019-04-04T14:34:30: Name ......................: Extragalactic_diffuse
  172 + 2019-04-04T14:34:30: Instruments ...............: all
  173 + 2019-04-04T14:34:30: Instrument scale factors ..: unity
  174 + 2019-04-04T14:34:30: Observation identifiers ...: all
  175 + 2019-04-04T14:34:30: Model type ................: DiffuseSource
  176 + 2019-04-04T14:34:30: Model components ..........: "DiffuseIsotropic" * "FileFunction" * "Constant"
  177 + 2019-04-04T14:34:30: Number of parameters ......: 3
  178 + 2019-04-04T14:34:30: Number of spatial par's ...: 1
  179 + 2019-04-04T14:34:30: Value ....................: 1 [0,10] (fixed,scale=1,gradient)
  180 + 2019-04-04T14:34:30: Number of spectral par's ..: 1
  181 + 2019-04-04T14:34:30: Normalization ............: 1 [0,1000] (fixed,scale=1,gradient)
  182 + 2019-04-04T14:34:30: Number of temporal par's ..: 1
  183 + 2019-04-04T14:34:30: Normalization ............: 1 (relative value) (fixed,scale=1,gradient)
184 184
185 The fit converged after 19 iterations with spectral parameters that are 185 The fit converged after 19 iterations with spectral parameters that are
186 reasonably close to those found in 186 reasonably close to those found in
@@ -241,9 +241,9 @@ using the Fermi-LAT Science Tools. @@ -241,9 +241,9 @@ using the Fermi-LAT Science Tools.
241 241
242 .. code-block:: none 242 .. code-block:: none
243 243
244 - 2018-01-25T20:05:57: Name ......................: Vela  
245 - 2018-01-25T20:05:57: Instruments ...............: all  
246 - 2018-01-25T20:05:57: Test Statistic ............: 326341.439917148 244 + 2019-04-04T14:37:50: Name ......................: Vela
  245 + 2019-04-04T14:37:50: Instruments ...............: all
  246 + 2019-04-04T14:37:50: Test Statistic ............: 326341.43991719
247 247
248 which is also close to the value determined by gtlike. 248 which is also close to the value determined by gtlike.
249 249
doc/source/users/tutorials/howto/fermi/howto_fermi_spectrum.png
View file @0873cb8

22 KB | W: | H:

44.2 KB | W: | H:

  • 2-up
  • Swipe
  • Onion skin
doc/source/users/tutorials/howto/fermi/howto_fermi_tsmap.rst
  View file @0873cb8
@@ -20,14 +20,14 @@ source. You do this by using the :ref:`cttsmap` tool as follows: @@ -20,14 +20,14 @@ source. You do this by using the :ref:`cttsmap` tool as follows:
20 Input event list, counts cube or observation definition XML file [events.fits] obs.xml 20 Input event list, counts cube or observation definition XML file [events.fits] obs.xml
21 Test source name [Crab] Vela 21 Test source name [Crab] Vela
22 Input model definition XML file [$CTOOLS/share/models/crab.xml] vela_results.xml 22 Input model definition XML file [$CTOOLS/share/models/crab.xml] vela_results.xml
  23 + Coordinate system (CEL - celestial, GAL - galactic) (CEL|GAL) [CEL]
  24 + Projection method (AIT|AZP|CAR|GLS|MER|MOL|SFL|SIN|STG|TAN) [CAR]
23 First coordinate of image center in degrees (RA or galactic l) (0-360) [83.63] 128.84 25 First coordinate of image center in degrees (RA or galactic l) (0-360) [83.63] 128.84
24 Second coordinate of image center in degrees (DEC or galactic b) (-90-90) [22.01] -45.18 26 Second coordinate of image center in degrees (DEC or galactic b) (-90-90) [22.01] -45.18
25 - Projection method (AIT|AZP|CAR|GLS|MER|MOL|SFL|SIN|STG|TAN) [CAR]  
26 - Coordinate system (CEL - celestial, GAL - galactic) (CEL|GAL) [CEL]  
27 Image scale (in degrees/pixel) [0.02] 0.3 27 Image scale (in degrees/pixel) [0.02] 0.3
28 Size of the X axis in pixels [200] 10 28 Size of the X axis in pixels [200] 10
29 Size of the Y axis in pixels [200] 10 29 Size of the Y axis in pixels [200] 10
30 - Output Test Statistic map file [tsmap.fits] 30 + Output Test Statistic map file [tsmap.fits]
31 31
32 :ref:`cttsmap` created a 10 x 10 pixels map with a pixel size of 0.3 degrees 32 :ref:`cttsmap` created a 10 x 10 pixels map with a pixel size of 0.3 degrees
33 around the location of the Vela pulsar and saved the result in the FITS file 33 around the location of the Vela pulsar and saved the result in the FITS file
doc/source/users/tutorials/howto/fermi/howto_fermi_ulimit.rst
  View file @0873cb8
@@ -20,7 +20,7 @@ for the source. To do this you first have to add an additional source to the @@ -20,7 +20,7 @@ for the source. To do this you first have to add an additional source to the
20 <source_library title="source library"> 20 <source_library title="source library">
21 <source type="PointSource" name="Test"> 21 <source type="PointSource" name="Test">
22 <spectrum type="PowerLaw"> 22 <spectrum type="PowerLaw">
23 - <parameter name="PhotonFlux" scale="1e-12" value="1.0" min="1e-07" max="1000.0" free="1"/> 23 + <parameter name="PhotonFlux" scale="1e-12" value="1.0" min="1e-07" max="1e+07" free="1"/>
24 <parameter name="Index" scale="1.0" value="-2.0" min="-5.0" max="+5.0" free="0"/> 24 <parameter name="Index" scale="1.0" value="-2.0" min="-5.0" max="+5.0" free="0"/>
25 <parameter name="LowerLimit" scale="1.0" value="100.0" min="10.0" max="1000000.0" free="0"/> 25 <parameter name="LowerLimit" scale="1.0" value="100.0" min="10.0" max="1000000.0" free="0"/>
26 <parameter name="UpperLimit" scale="1.0" value="100000.0" min="10.0" max="1000000.0" free="0"/> 26 <parameter name="UpperLimit" scale="1.0" value="100000.0" min="10.0" max="1000000.0" free="0"/>
@@ -53,9 +53,9 @@ The results of the upper limit computation can then be extracted from the @@ -53,9 +53,9 @@ The results of the upper limit computation can then be extracted from the
53 53
54 .. code-block:: bash 54 .. code-block:: bash
55 55
56 - 2018-01-25T20:43:02: +=====================+  
57 - 2018-01-25T20:43:02: | Upper limit results |  
58 - 2018-01-25T20:43:02: +=====================+  
59 - 2018-01-25T20:43:02: Differential flux limit ...: 1.00100099727198e-13 ph/cm2/s/MeV at 0.001 TeV  
60 - 2018-01-25T20:43:02: Integral flux limit .......: 9.9999999627471e-10 ph/cm2/s within [0.0001-0.1] TeV  
61 - 2018-01-25T20:43:02: Energy flux limit .........: 1.10785216561628e-12 erg/cm2/s within [0.0001-0.1] TeV 56 + 2019-04-04T15:03:22: +=====================+
  57 + 2019-04-04T15:03:22: | Upper limit results |
  58 + 2019-04-04T15:03:22: +=====================+
  59 + 2019-04-04T15:03:22: Differential flux limit ...: 1.77270316013682e-13 ph/cm2/s/MeV at 0.001 TeV
  60 + 2019-04-04T15:03:22: Integral flux limit .......: 1.77093045697669e-09 ph/cm2/s within [0.0001-0.1] TeV
  61 + 2019-04-04T15:03:22: Energy flux limit .........: 1.96192914922622e-12 erg/cm2/s within [0.0001-0.1] TeV