diff --git a/doc/source/users/tutorials/quickstart/binning.rst b/doc/source/users/tutorials/quickstart/binning.rst index c1a6bd1..fe55dfa 100755 --- a/doc/source/users/tutorials/quickstart/binning.rst +++ b/doc/source/users/tutorials/quickstart/binning.rst @@ -20,10 +20,10 @@ using the :ref:`ctbin` tool as follows: $ ctbin Input event list or observation definition XML file [events.fits] selected_events.fits + 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] Second coordinate of image center in degrees (DEC or galactic b) (-90-90) [22.51] - 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] Size of the X axis in pixels [200] Size of the Y axis in pixels [200] @@ -31,7 +31,7 @@ using the :ref:`ctbin` tool as follows: Lower energy limit (TeV) [0.1] Upper energy limit (TeV) [100.0] Number of energy bins (1-200) [20] - Output counts cube file [cntcube.fits] + Output counts cube file or observation definition XML file [cntcube.fits] In this example the events from file ``selected_events.fits`` will be binned into a counts cube stored into the file ``cntcube.fits``. The counts cube @@ -59,7 +59,7 @@ above. The ``EBOUNDS`` table has 20 rows, one for each energy bin, while the ``GTI`` table has just a single row, indicating the start and stop time of the simulated data. -.. figure:: cntmap-fits.jpg +.. figure:: cntmap_fits.png :width: 600px :align: center @@ -69,7 +69,7 @@ simulated data. An image of the first bin, covering the energy range 100 - 141 GeV, is shown below: -.. figure:: cntmap-map.jpg +.. figure:: cntmap_map.png :height: 400px :align: center @@ -80,17 +80,17 @@ reproduced below: .. code-block:: none - 2018-01-24T14:23:03: +=================+ - 2018-01-24T14:23:03: | Bin observation | - 2018-01-24T14:23:03: +=================+ - 2018-01-24T14:23:03: === CTA observation (id=000001) === - 2018-01-24T14:23:03: Events in list ............: 21974 - 2018-01-24T14:23:03: Events in cube ............: 18799 - 2018-01-24T14:23:03: Events outside RoI ........: 0 - 2018-01-24T14:23:03: Events with invalid WCS ...: 0 - 2018-01-24T14:23:03: Events outside cube area ..: 3175 - 2018-01-24T14:23:03: Events outside energy bins : 0 - -From the 21974 events that were simulated and stored in the -``selected_events.fits`` file, 18799 lie within the cube boundaries and are thus + 2019-04-02T13:46:34: +=================+ + 2019-04-02T13:46:34: | Bin observation | + 2019-04-02T13:46:34: +=================+ + 2019-04-02T13:46:34: === CTA observation (id=000001) === + 2019-04-02T13:46:34: Events in list ............: 22708 + 2019-04-02T13:46:34: Events in cube ............: 19452 + 2019-04-02T13:46:34: Events outside RoI ........: 0 + 2019-04-02T13:46:34: Events with invalid WCS ...: 0 + 2019-04-02T13:46:34: Events outside cube area ..: 3256 + 2019-04-02T13:46:34: Events outside energy bins : 0 + +From the 22708 events that were simulated and stored in the +``selected_events.fits`` file, 19452 lie within the cube boundaries and are thus put into the resulting counts cube. diff --git a/doc/source/users/tutorials/quickstart/butterfly.jpg b/doc/source/users/tutorials/quickstart/butterfly.jpg deleted file mode 100755 index 1201504..0000000 --- a/doc/source/users/tutorials/quickstart/butterfly.jpg +++ /dev/null diff --git a/doc/source/users/tutorials/quickstart/butterfly.png b/doc/source/users/tutorials/quickstart/butterfly.png new file mode 100644 index 0000000..dc16170 --- /dev/null +++ b/doc/source/users/tutorials/quickstart/butterfly.png diff --git a/doc/source/users/tutorials/quickstart/butterfly.rst b/doc/source/users/tutorials/quickstart/butterfly.rst index 928b2b0..2186783 100755 --- a/doc/source/users/tutorials/quickstart/butterfly.rst +++ b/doc/source/users/tutorials/quickstart/butterfly.rst @@ -38,15 +38,15 @@ default is set to 100. Below an excerpt of the first lines of the .. code-block:: none - 103514.216667934 8.08473784238882e-15 7.8714657958793e-15 8.29800988889835e-15 - 110917.48152624 6.81070529541223e-15 6.63639780809342e-15 6.98501278273104e-15 - 118850.222743702 5.73744103089563e-15 5.59501033791654e-15 5.87987172387472e-15 - 127350.308101666 4.83330700055084e-15 4.71693954114459e-15 4.94967445995709e-15 + 103514.216667934 8.17999288138647e-15 7.96480369463198e-15 8.39518206814095e-15 + 110917.48152624 6.88846691601701e-15 6.71266078582656e-15 7.06427304620745e-15 + 118850.222743702 5.80085791529665e-15 5.65725791756428e-15 5.94445791302902e-15 + 127350.308101666 4.88496975650954e-15 4.76769161306379e-15 5.00224789995529e-15 ... A graphical display of the results is shown below: -.. figure:: butterfly.jpg +.. figure:: butterfly.png :width: 600px :align: center diff --git a/doc/source/users/tutorials/quickstart/cntmap-fits.jpg b/doc/source/users/tutorials/quickstart/cntmap-fits.jpg deleted file mode 100755 index 265cf59..0000000 --- a/doc/source/users/tutorials/quickstart/cntmap-fits.jpg +++ /dev/null diff --git a/doc/source/users/tutorials/quickstart/cntmap-map.jpg b/doc/source/users/tutorials/quickstart/cntmap-map.jpg deleted file mode 100755 index e755c06..0000000 --- a/doc/source/users/tutorials/quickstart/cntmap-map.jpg +++ /dev/null diff --git a/doc/source/users/tutorials/quickstart/cntmap_fits.png b/doc/source/users/tutorials/quickstart/cntmap_fits.png new file mode 100644 index 0000000..985d377 --- /dev/null +++ b/doc/source/users/tutorials/quickstart/cntmap_fits.png diff --git a/doc/source/users/tutorials/quickstart/cntmap_map.png b/doc/source/users/tutorials/quickstart/cntmap_map.png new file mode 100644 index 0000000..ebe5ae5 --- /dev/null +++ b/doc/source/users/tutorials/quickstart/cntmap_map.png diff --git a/doc/source/users/tutorials/quickstart/ctools-from-python.ipynb b/doc/source/users/tutorials/quickstart/ctools-from-python.ipynb index 455cd71..5225a4b 100644 --- a/doc/source/users/tutorials/quickstart/ctools-from-python.ipynb +++ b/doc/source/users/tutorials/quickstart/ctools-from-python.ipynb @@ -175,7 +175,10 @@ "\r\n", "tmax [time]\r\n", " Stop time (UTC string, JD, MJD or MET in seconds).\r\n", - "\r", + "\r\n", + "(mjdref = 51544.5) [real]\r\n", + " Reference Modified Julian Day (MJD) for simulated events. The times in\r\n", + " seconds for each event are counted from this reference time on.\r\n", "\r\n", "emin [real]\r\n", " Lower energy limit of simulated events (in TeV).\r\n", @@ -197,6 +200,9 @@ "Standard parameters\r\n", "-------------------\r\n", "\r\n", + "(nthreads = 0) [integer]\r\n", + " Number of parallel processes (0=use all available CPUs).\r\n", + "\r\n", "(publish = no) [boolean]\r\n", " Specifies whether the event list(s) should be published on VO Hub.\r\n", "\r\n", @@ -275,7 +281,7 @@ "=== GObservations ===\n", " Number of observations ....: 1\n", " Number of models ..........: 2\n", - " Number of observed events .: 200741\n", + " Number of observed events .: 204568\n", " Number of predicted events : 0\n" ] } @@ -318,9 +324,9 @@ " Caldb instrument ..........: prod2\n", " Response name .............: South_0.5h\n", " Energy dispersion .........: Not used\n", - " Save energy range .........: undefined\n", + " Safe energy range .........: undefined\n", "=== GCTAEventList ===\n", - " Number of events ..........: 200741 (disposed in \"events.fits\")\n", + " Number of events ..........: 204568 (disposed in \"events.fits\")\n", " Time interval .............: 58849.0008007407 - 58849.0424674074 days\n", " Energy interval ...........: 0.03 - 150 TeV\n", " Region of interest ........: RA=83.5, DEC=22.8 [0,0] Radius=5 deg\n", @@ -373,9 +379,9 @@ "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAERCAYAAAB4jRxOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGQZJREFUeJzt3X2QZGd13/Hvb1EkLDBgsJGMxK4gslBZRcKLI16cFAPC\nloCExZgUWpkYRDCUcMnkhUJUCWd3TBUhdiABBNhJ1mtMSiyvAslYhcBoIKpFoCCkxWJXWvMiaSO0\nmFgEQchKiJM/+o621fTs9Ex33+6e/n6qprb7dvc9p+/OnH76ufeem6pCkrTxbZp0ApKkdljwJWlO\nWPAlaU5Y8CVpTljwJWlOWPAlaU5Y8CVpTljwJWlOHDOOlSYJ8GbgYcB1VfX+ccSRJA1uXCP8rcBJ\nwD3AwTHFkCStwUAFP8nOJIeS7O1Zfk6S/UluSXJR10NPAPZU1euB144wX0nSOg06wt8FnN29IMkm\n4JJm+RnAtiSnNw8fBO5qbv94BHlKkoY0UMGvqms4UsCXnQkcqKpbq+peYDedqRyAjwHnJHkH8PlR\nJStJWr9hdtqeBNzedf8gnQ8BqupHwKuO9uIktumUpHWoqqzndcPstO0XcE1FfPv27Vx99dVU1dT8\nbN++feI5mJM5zWNe5nT0n6uvvprt27cPUbKHG+EfBDZ33T8ZuGMtK9ixY8cQ4SVpfiwsLLCwsMDi\n4uK617GWEX544Kj+OuDUJFuSHAucC1y+luA7duxgaWlpLS+RpLm0tLQ09CB50MMyLwX2AKcluS3J\n+VV1H3AhcBVwE7C7qvatJfiOHTtYWFhYY8rjNW35gDkNypwGN415mdPRLSwsDF3wUzWZfadJavv2\n7fd/TZEkrWxpaYmlpSUWFxepde60nWjBn1RsSZpVSdZd8G2eJklzYqIF3522kjSYUey0dUpHkmaI\nUzqSpFU5pSNJM8ApHUmaM07pSJJW5ZSOJM0Ap3Qkac44pSNJWpUFX5LmhAVfkuaEO20laQa401aS\n5ow7bSVJq7LgS9KcsOBL0pyw4EvSnLDgS9Kc8LBMSZoBHpYpSXPGwzIlSauy4EvSnLDgS9KcOGbS\nCUyDe+65h71797YWb/PmzTz60Y9uLZ4kwZgKfpJnAW8GbgI+UFWfH0ecUXnb297O4uI7efCDHzP2\nWD/+8Q859dSf54Yb/sfYY0lSt3GN8Au4GzgOODimGCPzgx/8kMOHL+Dw4d9vIdr13H33q1qII0kP\nNNAcfpKdSQ4l2duz/Jwk+5PckuSi5eVV9fmqegHwRuAPRpuyJGk9Bt1puws4u3tBkk3AJc3yM4Bt\nSU7ved33gGOHTVKSNLyBpnSq6pokW3oWnwkcqKpbAZLsBrYC+5P8Bp0PgofT+VCQJE3YMHP4JwG3\nd90/SOdDgKq6DLhstRV0nya8sLDAwsLCEOlI0saztLQ0shY0wxT8fqf2rqlXwrB9ISRpo+sdDC8u\nLq57XcOceHUQ2Nx1/2TgjrWswOZpkjSYUTRPW0vBDw8c1V8HnJpkS5JjgXOBy4fKRpI0NoMelnkp\nsAc4LcltSc6vqvuAC4Gr6Jxgtbuq9q0l+I4dO5y3l6QBLCwsDD3CH/QonfNWWH4lcOVQGUiSWuEF\nUCRpBngBlBG5+OLf5y1vORZop7XC4x//Kr7+9etbiCVpo5nZC6A4wpekwTjCHxFH+JJmxcyO8CVJ\n7XFKR5JmgFM6I+KUjqRZ4ZSOJGlVFnxJmhMTvYj5cmuFeWuvcOutB0jW9Y1sXU44YQt33vmt1uJJ\nGr1RtEl2Dp/25/Dhqayxk/SQwrRsa0nDcQ5fkrQqC74kzQmPw5ekGeBx+CPiHL6kWeEcviRpVRZ8\nSZoTFnxJmhMWfEmaExZ8SZoTHpYpSTPAwzJHxMMyJc0KD8uUJK3Kgi9Jc8KCPxeOI0krPyeeeMqk\n36ykFUy0H77acpi29hkcOtRen39JazO2EX6S45P8zyTPH1cMSdLgxjmlcxHwwTGuX5K0BgMV/CQ7\nkxxKsrdn+TlJ9ie5JclFXcvPAr4GfAfwO74kTYFBR/i7gLO7FyTZBFzSLD8D2Jbk9ObhZwNPA84D\nXjWaVCVJwxhop21VXZNkS8/iM4EDVXUrQJLdwFZgf1W9qVn228B3R5ivJGmdhjlK5yTg9q77B+l8\nCNyvqv78aCvoPk14YWGBhYWFIdKRpI1naWlpZC1ohin4/ebm13Ts37B9ISRpo+sdDC8uLq57XcMc\npXMQ2Nx1/2TgjrWswOZpkjSYVpunJTkFuKKqntjcfxBwM3AW8G3gS8C2qto34PpsntaatBjPRm3S\nOI29eVqSS4E9wGlJbktyflXdB1wIXAXcBOwetNgvc4S/EdnGQRoH2yOPiCP82Y01Lb9DUltmtj2y\nI3xJGowj/BFxhD+7sabld0hqiyN8SdrgHOGPiCP82Y01Lb9DUltmdoQvSWqPUzqSNAOc0hkRp3Rm\nN9a0/A5JbXFKR5K0Kgu+JM0J5/AlaQY4hz8izuHPbqxp+R2S2uIcviRpVcNcAEWasE5nzraccMIW\n7rzzW63Fk0ZtogV/x44dXtpQQzhMm1Njhw619+Ei9RrFpQ6dw8c5fGMNHm9afmc1v5zDlyStyoIv\nSXPCgi9Jc8KCL0lzwoIvSXPC1gqSNANsrTAiHpZprEHjTcvvrOaXh2VKreic2dvGz4knnjLpN6sN\nyNYK0sDaO7PXs3o1Do7wJWlOWPAlaU6MZUonyenA64BHAZ+tqj8eRxxJ0uDGMsKvqv1VdQHwUuCZ\n44ghSVqbgQp+kp1JDiXZ27P8nCT7k9yS5KKex/4Z8BfAX44uXUnSeg06wt8FnN29IMkm4JJm+RnA\ntmYqB4CquqKqXgC8bES5SpKGMNAcflVdk2RLz+IzgQNVdStAkt3AVmB/kmcBLwaOAz45wnwlSes0\nzE7bk4Dbu+4fpPMhQFV9DvjcaivoPk3YK19J0k8bxZWulg1T8PudGbKms1KG7QshSRtd72B4cXFx\n3esa5iidg8DmrvsnA3esZQU2T5NWYhsHPVCrzdOSnAJcUVVPbO4/CLgZOAv4NvAlYFtV7RtwfTZP\na81GbWjWfvO0jfnebAo3S8bePC3JpcAe4LQktyU5v6ruAy4ErgJuAnYPWuyXOcKXpMHYHnlEHOEb\na/riOcJXfzPbHtkRviQNxhH+iDjCN9b0xXOEr/4c4UvSBucIf0Qc4Rtr+uI5wld/MzvClyS1xykd\nSZoBTumMiFM6xpq+eE7pqD+ndCRJq7LgS9KccA5fkmaAc/gj4hy+saYvnnP46s85fEnSqiz4kjQn\nLPiSNCfcaStJM8CdtiPiTltjTV88d9qqP3faSpJWZcGXpDlhwZekOWHBl6Q54VE6kjQDPEpnRDxK\nx1jTF6/NWA8GDrcUC044YQt33vmt1uJtNMMcpXPMqJORNGsO0+YH56FD66pVGgHn8CVpTljwJWlO\nWPAlaU6MbQ4/yVbgBcDPAn9aVZ8eVyxJ0urGVvCr6hPAJ5I8AvgjwIIvSRM08JROkp1JDiXZ27P8\nnCT7k9yS5KI+L30T8O5hE5UkDWctc/i7gLO7FyTZBFzSLD8D2Jbk9K7H3wr8ZVXdMIJcJUlDGLjg\nV9U1wF09i88EDlTVrVV1L7Ab2AqQ5ELgLOAlSV49onwlSes07Bz+ScDtXfcP0vkQoKreBbzraC/u\nPk14YWGBhYWFIdORpI1laWlpZC1ohi34/U6ZG/iUvWH7QkjSRtc7GF5cXFz3uoY9Dv8gsLnr/snA\nHYO+2OZpkjSY1punJTkFuKKqntjcfxBwM525+m8DXwK2VdW+AdZl87TWbNSmXzZPm71YnXjT8rc/\ni1q5xGGSS4E9wGlJbktyflXdB1wIXAXcBOwepNgvc4QvSYOxPfKIOMI31vTF26ixOvGm5W9/Fs3s\nRcwd4UvSYBzhj4gjfGNNX7yNGqsTb1r+9mfRzI7wJUntcUpHkmaAUzoj4pSOsaYv3kaN1Yk3LX/7\ns8gpHUnSqpzSkaQZ4JTOiDilY6zpi7dRY3XiTcvf/ixySkeStCqndCRpBjilMyJO6Rhr+uJt1Fid\neNPytz+LnNKRJK3Kgi9Jc8KCL6llx5GklZ8TTzxl0m92qgx7icOh7Nixw2vZSnPnMG3tMzh0aF1T\n3VNpFNe2dact7rQ11jTG26ix2o638XYQu9NWkrQqC74kzQkLviTNCQu+JM0JWytI0gywtcKIeJSO\nsaYv3kaN1XY8j9Lp5pSOJM0JC74kzQkLviTNCQu+JM0JC74kzYmxFPwkj0vy35J8aBzrlySt3VgK\nflV9s6peNY51S5LWZ6CCn2RnkkNJ9vYsPyfJ/iS3JLloPClKkkZh0BH+LuDs7gVJNgGXNMvPALYl\nOb3ndRunGbUkzbiBCn5VXQPc1bP4TOBAVd1aVfcCu4GtAEkemeS9wJMc+UvSdBjmilcnAbd33T9I\n50OAqvo74ILVVtDdF8IrX0nSTxvFla6WDVPw+03XrLlphYVekla2XCNHUfiHKfgHgc1d908G7ljL\nCobt/CZJ82K58C8uLq57HWs5LDM8cFR/HXBqki1JjgXOBS5fS3DbI0vSYFprj5zkUmABeBRwCNhe\nVbuSPA/4z3Q+OHZW1VsHDmx75BZiLduorXY3dlvfjRmr7Xi2R+420JROVZ23wvIrgSvXExg6I3zn\n8CVpdaOYw/cCKDjCN9Y0xtuosdqO5wi/m83TJGlOeE1bSZoBXtN2RJzSMdb0xduosdqO55RON6d0\nJGlOOKUjSTPAKZ0RcUrHWNMXb6PGajueUzrdnNKRpDnhlI6kDew4krTyc+KJp4z1nTilMyJO6Rhr\n+uJt1Fhtx9t400dO6UiSVmXBl6Q5YcGXpDnhTltJmgHutB0Rd9oaa/ribdRYbcdzp203p3QkaU5Y\n8CVpTljwJWlOWPAlaU5Y8CVpTnhYpiTNAA/LHBEPyzTW9MXbqLHajudhmd2c0pGkOWHBl6Q5YcGX\npDlhwZekOXHMOFaa5HjgPcBh4HNVdek44kiSBjeuEf6LgQ9X1WuAF44pxpgsTTqBPpYmnUAfS5NO\noI+lSSfQx9KkE1jB0qQT6GNp0gn0sTTpBEZqoIKfZGeSQ0n29iw/J8n+JLckuajroZOB25vb940o\n15YsTTqBPpYmnUAfS5NOoI+lSSfQx9KkE1jB0qQT6GNp0gn0sTTpBEZq0BH+LuDs7gVJNgGXNMvP\nALYlOb15+HY6RR86B8JKkiZsoIJfVdcAd/UsPhM4UFW3VtW9wG5ga/PYZcBLkrwbuGJUyUqS1m/g\nM22TbAGuqKp/0Nz/TeDsqnp1c/9lwJlV9XsDrm86TrOVpBmz3jNthzlKp1/AgYv4ehOWJK3PMEfp\nHAQ2d90/GbhjuHQkSeOyloIfHjiqvw44NcmWJMcC5wKXjzI5SdLoDHpY5qXAHuC0JLclOb+q7gMu\nBK4CbgJ2V9W+8aUqSRrGoEfpnFdVj6mq46pqc1XtapZfWVVPqKpfqqq3Hm0dSf4wyb4kNyT5aJKH\nrfC8lY7tH4skL0ny10nuS/KUozzvW0luTPKVJF+akpxa21ZJfi7JVUluTvKpJA9f4Xn3Jbm+2U4f\nH1MuR33fSY5NsjvJgSRfSLK533pazunlSb7TbJvrk7yyhZz6nj/T85x3NtvphiRPmnROSZ6V5Htd\n2+lNLeR0cpLPJvlakq8m6XvgSZvbapCc1rWtqqqVH+C5wKbm9luBf9/nOZuAvwG2AH8PuAE4fcx5\nPQH4JeCzwFOO8rxvAD/X0rZaNae2txXwH4A3NLcvAt66wvO+P+Zts+r7Bi4A3tPcfimdb5+Tzunl\nwDvb+P3pivmPgScBe1d4/HnAJ5vbTwOunYKcngVc3vJ2OhF4UnP7ocDNff7/Wt1WA+a05m3VWvO0\nqvpMVf2kuXstR07M6na0Y/vHldfNVXWA1U8QCy01mxswp7a31Vbgfc3t9wEvWuF54z76apD33Z3r\nR4CzpiAnaPkkxOp//ky3rcCfN8/9IvDwJCdMOCdofzvdWVU3NLd/AOwDTup5WqvbasCcYI3balLd\nMl8JXNln+UkcackAnSOB+r3JSSjgU0muS/I7k06G9rfVo6vqEHR+GYFfWOF5xyX5UpI9ScbxATTI\n+77/OdXZ1/S9JI8cQy5ryQngxc10wIeS9BvwtK037//FdPy9Pb2ZEvxkkl9uM3CSU+h8A/liz0MT\n21ZHyQnWuK1G2i0zyaeB7k+95euLXVxVVzTPuRi4t/p30Bzq2P5h8hrAM6vqziS/AHw6yb5mtDKp\nnEa+rY6S01rmUTc32+lxwGeT7K2qbw6TV2+afZb1vu/e54z7OneD5HQ5cGlV3ZvkNXS+gYz7m8dq\nxvL3NqQvA1uq6v8meR7wceC0NgIneSidb4Sva0bVD3i4z0vGvq1WyWnN22qkBb+qfu1ojyd5OfB8\n4DkrPGUsx/avlteA67iz+fdvk1xG52v8ugv+CHIa+bY6Wk7NjrYTqupQkhOB76ywjuXt9M0kS8CT\ngVEW/EHe9+3AY4E7kjwIeFhVrTaNMNaceuL/Vzr7RCbtIJ3ttGzi59J0F7WqujLJe5I8sqr+bpxx\nkxxDp7C+v6o+0ecprW+r1XJaz7ZqbUonyTnAG4AXVtXhFZ426WP7+86HJTm++aQlyUOAXwf+epI5\n0f62uhx4RXP75cBP/QImeUSTC0l+Hngm8LUR5zHI+76iyRHgn9PZ+T1Oq+bUfEgu28rot8tKes+f\n6XY58NsASZ4OfG952m5SOXXPiyc5k077l7EW+8afAl+rqnes8PgkttVRc1rXthrnnuaePcoHgFuB\n65uf5aMofhH4i67nnUNnj/QB4I0t5PUiOiPCHwHfBq7szQt4HJ0jL74CfHXceQ2SU9vbCngk8Jkm\n3qeBRzTLnwr8l+b2M4C9zXa6EXjFmHL5qfcNLAL/tLl9HPCh5vFrgVNa+D1aLae30BkkfAX4K+C0\nFnK6lM4o9DBwG3A+8Brg1V3PuYTOEUY3cpSj1NrKCfjdru20B3haCzn9Kp027st/49c3/58T21aD\n5LSebTVw8zRJ0mzzmraSNCcs+JI0Jyz4kjQnLPiSNCcs+JI0Jyz4kjQnLPgamSR3j2g9Vya5K0nv\nyUunJLk2nRbNH2jORFx+7HXpXFeZJLuSfCNH2jSv+4zoUUqyPcnBJDuSvKLJ7StJDqfTevv6JG9Z\n4bUPTfLdJMf3LL8iyYuSnNe07v1Ys/z+7SEt8zh8jUyS71dV3+scrHE9zwaOB15TVS/sWv5B4CNV\n9eEk7wVuqKo/adonXA88uap+kmQXnbaxlw2bS09em+pIx9f1vH47cHdVvb1n+TeAp9Yq7R+a9//x\nqvpAc/8RdE72emxV3ZPkLOB3q+rFzRnhn6+qp643X208jvA1Vkk2J/lM0yXy08tdIpM8Pp2Lk9yY\n5M3d3w6q6mqgt1EUdHowfbS53d2i+TnAl3uK8U/9bjcj7J1Jrk7yN0ku7Hrst5J8sRllvzdJmuV3\nJ/mPSb5CpzPh89O5kM91Sd6R5PJ03JLkUc1r0oy2B+3Q+YBWA0ke0nxLuTbJl5O8oHloN7Ct63W/\nSadH+z29K6yqHwIH08JFTTQ7LPgat0uAP6uqJ9E5rf5dzfJ3AP+pqv4hncZUR/2q2RTTu7qKencL\n4l+l0zmw2x91Tem8v2v5E4Bfo3MRi+1JHpTkdDoXSnlmVT0F+AnwW83zHwJ8oaqe3MT4Y+DsqvpH\nNC2iq/M1+f3A8hTKc+l8+1hvD5h/R6edxtPpdNR8e9Oj55PA03LkamPnAh84ynq+DPyTdeagDciC\nr3F7BkeK0vvpFOfl5R9pbvdrld3raO1pfxH4257HXl9VT6mqJ1fVv+ha/smq+nFV/W/gEJ120GcB\nTwGua0byz6HTPwk6/Uw+1tw+Hfh6Vd3W3O8utruA5TivbO6v168DFze5XA0cS6ft9D10iv6Lkzwa\n+GU6fXlW8h3gMUPkoQ1mpO2RpT56R+79RvKrXrWnqr7bdONcnkfvbk/7I+DBA+bT3an1Pjp/AwHe\nV1UX93n+j+rIjq4VuzxW1cF0Wkg/m07r7PMGzGclL6r+1xHYDbyezj6Oj62yT+HBdLaNBDjC12j1\nK4Z7ODLv/DKOXEPgC8BLmtvnrrCu3vVdTafdMTywRfM+4NQBclkp378CXpLOxW2WL9j+2J7nAOwH\nHpcjF0V/ac/6dgL/Hfhg14fEenwKeN39ST5wHv4zwBl0uiYebToHOhfDaKuNt2aABV+j9DNJbkty\ne/PvvwJ+Dzg/yQ105sWXC9m/Bv5Ns/zvA/9neSVJPg98EHhOs57lC7O8sXnNLXTaNe9sll9J54LO\n3f6waw7/+u5DOLsUQFXto3NVr6uS3AhcRWea6P7nNM/7f8BraS51CXy/O286PdMfAvzZqluqTx5d\n/gA4PsneJF8Ftnfl8BPgMuBnq2rPKut9Bkef8tGc8bBMTUSSn6mqHzW3XwqcW1W/McT6Pgq8oaq+\nPqocV4jzkOYIGJK8G7ilmgtUJPkV4G1V1fvhs/za7cAPquptY8rtucBrm8MyfwW4oKr+5ThiaTY5\nwtekPLU5VPNG4ALg3w65vjdyZFQ+Tr/TfGu4CXgY8CcASS4CPtzksZIfNK/fMeqkkpxH58in5SOD\nHknXNwMJHOFL0txwhC9Jc8KCL0lzwoIvSXPCgi9Jc8KCL0lz4v8DiQ2Pubkbe1IAAAAASUVORK5C\nYII=\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAERCAYAAAB4jRxOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGP1JREFUeJzt3XuwpHV95/H3Z2TBoFGjm0ACMuAapEK56yULXnbLVoyg\n7ooxbgnEjeIaLU0R92KJVZidmaSKdXPbNaImZsnEuDWOVxRiKNFIa1F4YUXAIAPjjWEWZ4wbXHXX\nRTJ+949+DtO0feb0OX0/z/tVNUX3093P93sezvn2r3/Pr79PqgpJ0ua3Zd4JSJJmw4IvSS1hwZek\nlrDgS1JLWPAlqSUs+JLUEhZ8SWoJC74ktcRR09hpkgC/AzwMuKGq3j2NOJKk0U1rhH8ucALwQ2D/\nlGJIktZhpIKf5PIkB5PcMrD9nCR7ktyR5OK+hx4HXF9VrwdeO8F8JUkbNOoIfydwdv+GJFuAy5rt\npwPnJzmteXg/cE9z++8nkKckaUwjFfyquo7DBXzFGcDeqrqzqu4DdtObygH4EHBOkrcAn55UspKk\njRvnpO0JwF199/fTexOgqn4AvPJIL05im05J2oCqykZeN85J22EB11XEt23bxrXXXktVLcy/bdu2\nzT0HczKnNuZlTkf+d+2117Jt27YxSvZ4I/z9wEl9908E7l7PDrZv3z5GeElqj06nQ6fTYceOHRve\nx3pG+OGBo/obgMcm2ZrkaOA84Mr1BN++fTvdbnc9L5GkVup2u2MPkkddlrkLuB44Ncm+JBdW1SHg\nIuAa4FZgd1Xdtp7g27dvp9PprDPl6Vq0fMCcRmVOo1vEvMzpyDqdztgFP1XzOXeapLZt23b/xxRJ\n0uq63S7dbpcdO3ZQGzxpO9eCP6/YkrSskmy44Ns8TZJaYq4F35O2kjSaSZy0dUpHkpaIUzqSpDU5\npSNJS8ApHUlqGad0JElrckpHkpaAUzqS1DJO6UiS1mTBl6SWsOBLUkt40laSloAnbSWpZTxpK0la\nkwVfklrCgi9JLWHBl6SWsOBLUku4LFOSloDLMiWpZVyWKUlakwVfklrCgi9JLXHUvBNYBPv27WPX\nrl0zi3fmmWfyzGc+c2bxJAmmVPCTPAP4HeBW4D1V9elpxJmUd77zT7n00i5btjx96rGqDnDKKe/j\nK1+5ceqxJKnftEb4BXwPOAbYP6UYE1MFVc/h0KHfmkG0G6l65QziSNIDjTSHn+TyJAeT3DKw/Zwk\ne5LckeTile1V9emqej7wRuC3J5uyJGkjRj1puxM4u39Dki3AZc3204Hzk5w28LrvAEePm6QkaXwj\nTelU1XVJtg5sPgPYW1V3AiTZDZwL7Enyy/TeCB5O701BkjRn48zhnwDc1Xd/P703AarqCuCKtXbQ\n/zXhTqdDp9MZIx1J2ny63e7EWtCMU/CHfbV3Xb0Sxu0LIUmb3eBgeMeOHRve1zhfvNoPnNR3/0Tg\n7vXswOZpkjSaSTRPW0/BDw8c1d8APDbJ1iRHA+cBV46VjSRpakZdlrkLuB44Ncm+JBdW1SHgIuAa\nel+w2l1Vt60n+Pbt2523l6QRdDqdsUf4o67SuWCV7VcDV4+VgSRpJrwAiiQtAS+AMiGXXPJbXHrp\n0cBsWis85jGv5KtftZeOpPVb2gugOMKXpNE4wp8QR/iSlsXSjvAlSbPjlI4kLQGndCbEKR1Jy8Ip\nHUnSmrym7Rzceedekg29QW/Iccdt5cCBb8wsnqTFNNeCv9JaoW3tFQ4d+j7rbCw6loMHZ/fmImk6\nJtEm2Tl8Zj+HD09mlgUfwqIca0njcQ5fkrQmC74ktYTr8CVpCbgOf0Kcw5e0LJzDlyStyYIvSS1h\nwZeklrDgS1JLWPAlqSVclilJS8BlmRPiskxJy8JlmZKkNVnwJaklLPitcAxJZvLv+ONPnvcPK2kV\nXgClFe5lVucM7L0vLa6pjfCTHJvkfyR53rRiSJJGN80pnYuB905x/5KkdRip4Ce5PMnBJLcMbD8n\nyZ4kdyS5uG/7WcCXgW8BfsaXpAUw6gh/J3B2/4YkW4DLmu2nA+cnOa15+JnAmcAFwCsnk6okaRwj\nnbStquuSbB3YfAawt6ruBEiyGzgX2FNVb2q2/Rrw7QnmK0naoHFW6ZwA3NV3fz+9N4H7VdVfHGkH\n/V8T7nQ6dDqdMdKRpM2n2+1OrAXNOAV/2Nz8utb+jdsXQpI2u8HB8I4dOza8r3FW6ewHTuq7fyJw\n93p2YPM0SRrNTJunJTkZuKqqHt/cfxBwO3AW8E3g88D5VXXbiPuzedrMZIbxbNQmTdPUm6cl2QVc\nD5yaZF+SC6vqEHARcA1wK7B71GK/whH+ZmQbB2kabI88IY7wlzfWovwOSbOytO2RHeFL0mgc4U+I\nI/zljbUov0PSrDjCl6RNzhH+hDjCX95Yi/I7JM3K0o7wJUmz45SOJC0Bp3QmxCmd5Y21KL9D0qw4\npSNJWpMFX5Jawjl8SVoCzuFPiHP4yxtrUX6HpFkZZw5/nH740pz1GrXNynHHbeXAgW/MLJ40aRZ8\nLbF7meUnpYMHZ/fmIk2Dc/iStAScw58Q5/CNNWq8RfmdVXu5Dl+StCYLviS1hAVfklrCgi9JLWHB\nl6SWcFmmJC0Bl2VOiMsyjTWaB9P7stf0+a1ercbWCtJMzO6bvX6rV9PgHL4ktYQFX5JawoIvSS0x\nlTn8JKcBrwMeBXyyqv54GnEkSaObygi/qvZU1WuAlwBPm0YMSdL6jFTwk1ye5GCSWwa2n5NkT5I7\nklw88Ni/BP4S+KvJpStJ2qhRR/g7gbP7NyTZAlzWbD8dOL+ZygGgqq6qqucDL51QrpKkMYw0h19V\n1yXZOrD5DGBvVd0JkGQ3cC6wJ8kzgBcBxwAfnWC+kqQNGuek7QnAXX3399N7E6CqPgV8aq0d9H9N\nuNPp0Ol0xkhHkjafbrc7sRY04xT8YV8FXNfXEMftCyFJm93gYHjHjh0b3tc4q3T2Ayf13T8RuHs9\nO7B5mrSaY0gyk3/HH3/yvH9YjWCmzdOSnAxcVVWPb+4/CLgdOAv4JvB54Pyqum3E/dk8bWY2a0Oz\n2V/TdnP+bF6rd5lM/Zq2SXYB1wOnJtmX5MKqOgRcBFwD3ArsHrXYr3CEL0mjsT3yhDjCN9bixXOE\nr+GmPsKfFkf4kjQaR/gT4gjfWIsXzxG+hnOEL0mbnCP8CXGEb6zFi+cIX8Mt7QhfkjQ7TulI0hJw\nSmdCnNIx1uLFc0pHwzmlI0lakwVfklrCOXxJWgLO4U+Ic/jGWrx4zuFrOOfwJUlrsuBLUktY8CWp\nJTxpK0lLwJO2E+JJW2MtXjxP2mo4T9pKktZkwZeklrDgS1JLWPAlqSWOmmfw7du30+l06HQ680xD\narljSDZ0DnBDjjtuKwcOfGNm8TaLbrc79qpGV+ngKh1jLWK8zRqrF29R/vaXkat0JElrsuBLUktY\n8CWpJSz4ktQSU1ulk+Rc4PnATwJ/VlUfn1YsSdLaplbwq+ojwEeSPAL4PcCCL0lzNPKUTpLLkxxM\ncsvA9nOS7ElyR5KLh7z0TcDbxk1UkjSe9czh7wTO7t+QZAtwWbP9dOD8JKf1Pf5m4K+q6qYJ5CpJ\nGsPIBb+qrgPuGdh8BrC3qu6sqvuA3cC5AEkuAs4CXpzkVRPKV5K0QePO4Z8A3NV3fz+9NwGq6q3A\nW4/04v5m/rZYkKQfN4mWCivGLfjDvt478nemx716iyRtdoOD4R07dmx4X+Ouw98PnNR3/0Tg7lFf\n7CUOJWk0M7/EYZKTgauq6vHN/QcBt9Obq/8m8Hng/Kq6bYR92TxtZjZrIy6bpy1frF68RfnbX0Yz\naZ6WZBdwPXBqkn1JLqyqQ8BFwDXArcDuUYr9Ckf4kjQaL2I+IY7wjbV48TZrrF68RfnbX0ZL2x7Z\nEb4kjcYR/oQ4wjfW4sXbrLF68Rblb38ZLe0IX5I0O07pSNIScEpnQpzSMdbixdussXrxFuVvfxk5\npSNJWpNTOpK0BJzSmRCndIy1ePE2a6xevEX5219GTulIktbklI4kLQGndCbEKR1jLV68zRqrF29R\n/vaXkVM6kqQ1WfAlqSUs+JLUEp60lTRjx5BkJv+OP/7kef+wE+NJ2wnxpK2xFi/eZo0163ib7wSx\nJ20lSWuy4EtSS1jwJaklLPiS1BKu0pGkJeAqnQlxlY6xFi/eZo0163iu0unnlI4ktYQFX5JawoIv\nSS1hwZeklrDgS1JLTKXgJzklyX9L8r5p7F+StH5TKfhV9fWqeuU09i1J2piRCn6Sy5McTHLLwPZz\nkuxJckeSi6eToiRpEkYd4e8Ezu7fkGQLcFmz/XTg/CSnDbxuQ18OkCRN3kgFv6quA+4Z2HwGsLeq\n7qyq+4DdwLkASR6Z5B3AExz5S9JiOGqM154A3NV3fz+9NwGq6u+A16y1g/6+EJ1Oh06nM0Y6krT5\ndLvdifUcG6fgD5uuWXfTCgu9JK1upUZOovCPU/D3Ayf13T8RuHs9Oxi385sktcVK4d+xY8eG97Ge\nZZnhgaP6G4DHJtma5GjgPODK9QS3PbIkjWZm7ZGT7AI6wKOAg8C2qtqZ5LnAf6X3xnF5Vb155MC2\nR55BrBWbtdXu5m7ruzljzTqe7ZH7jTSlU1UXrLL9auDqjQSG3gjfOXxJ03MMyWxWhx933FYOHPjG\n1PY/iTl8L4CCI3xjLWK8zRpr1vE236cJL4AiSVqT17SVpCXgNW0nxCkdYy1evM0aa9bxnNLp55SO\nJLWEUzqStASc0pkQp3SMtXjxNmusWcdzSqefUzqS1BJO6UjSEnBKZ0Kc0jHW4sXbrLFmHc8pnX5O\n6UhSS1jwJaklLPiS1BKetJWkJeBJ2wnxpK2xFi/eZo0163ietO3nlI4ktYQFX5JawoIvSS1hwZek\nlrDgS1JLuCxTkpaAyzInxGWZxlq8eJs11qzjuSyzn1M6ktQSFnxJagkLviS1hAVfklriqGnsNMmx\nwNuBe4FPVdWuacSRJI1uWiP8FwHvr6pXAy+YUowp6c47gSG6805giO68ExiiO+8EhujOO4FVdOed\nwBDdeScwRHfeCUzUSAU/yeVJDia5ZWD7OUn2JLkjycV9D50I3NXcPjShXGekO+8EhujOO4EhuvNO\nYIjuvBMYojvvBFbRnXcCQ3TnncAQ3XknMFGjjvB3Amf3b0iyBbis2X46cH6S05qH76JX9KG3EFaS\nNGcjFfyqug64Z2DzGcDeqrqzqu4DdgPnNo9dAbw4yduAqyaVrCRp40b+pm2SrcBVVfWPm/u/Apxd\nVa9q7r8UOKOqfnPE/S3G12wlacls9Ju246zSGRZw5CK+0YQlSRszziqd/cBJffdPBO4eLx1J0rSs\np+CHB47qbwAem2RrkqOB84ArJ5mcJGlyRl2WuQu4Hjg1yb4kF1bVIeAi4BrgVmB3Vd02vVQlSeMY\ndZXOBVX1c1V1TFWdVFU7m+1XV9Xjqurnq+rNR9pHkt9NcluSm5J8MMnDVnneamv7pyLJi5P8TZJD\nSZ50hOd9I8nNSb6Y5PMLktPMjlWSn0pyTZLbk3wsycNXed6hJDc2x+nDU8rliD93kqOT7E6yN8ln\nkpw0bD8zzullSb7VHJsbk7xiBjkN/f7MwHP+qDlONyV5wrxzSvKMJN/pO05vmkFOJyb5ZJIvJ/lS\nkqELT2Z5rEbJaUPHqqpm8g94NrCluf1m4D8Nec4W4CvAVuAfADcBp005r8cBPw98EnjSEZ73NeCn\nZnSs1sxp1scK+M/AG5rbFwNvXuV5353ysVnz5wZeA7y9uf0Sep8+553Ty4A/msXvT1/MfwY8Abhl\nlcefC3y0uX0m8NkFyOkZwJUzPk7HA09obj8UuH3I/7+ZHqsRc1r3sZpZ87Sq+kRV/ai5+1kOfzGr\n35HW9k8rr9urai9rf0EszKjZ3Ig5zfpYnQu8q7n9LuCFqzxv2quvRvm5+3P9AHDWAuQEM/4SYg3/\n/ky/c4G/aJ77OeDhSY6bc04w++N0oKpuam5/H7gNOGHgaTM9ViPmBOs8VvPqlvkK4Ooh20/gcEsG\n6K0EGvZDzkMBH0tyQ5Jfn3cyzP5Y/UxVHYTeLyPw06s875gkn09yfZJpvAGN8nPf/5zqnWv6TpJH\nTiGX9eQE8KJmOuB9SYYNeGZtMO//yWL8vT2lmRL8aJJfmGXgJCfT+wTyuYGH5nasjpATrPNYTbRb\nZpKPA/3veivXF7ukqq5qnnMJcF8N76A51tr+cfIawdOq6kCSnwY+nuS2ZrQyr5wmfqyOkNN65lFP\nao7TKcAnk9xSVV8fJ6/BNIdsG/y5B58z7evcjZLTlcCuqrovyavpfQKZ9iePtUzl721MXwC2VtX/\nTfJc4MPAqbMInOSh9D4Rvq4ZVT/g4SEvmfqxWiOndR+riRb8qvqlIz2e5GXA84BnrfKUqaztXyuv\nEfdxoPnv3ya5gt7H+A0X/AnkNPFjdaScmhNtx1XVwSTHA99aZR8rx+nrSbrAE4FJFvxRfu67gEcD\ndyd5EPCwqlprGmGqOQ3E/1N650TmbT+947Ri7t+l6S9qVXV1krcneWRV/d004yY5il5hfXdVfWTI\nU2Z+rNbKaSPHamZTOknOAd4AvKCq7l3lafNe2z90PizJsc07LUkeAjwH+Jt55sTsj9WVwMub2y8D\nfuwXMMkjmlxI8g+BpwFfnnAeo/zcVzU5Avwreie/p2nNnJo3yRXnMvnjsprB78/0uxL4NYAkTwG+\nszJtN6+c+ufFk5xBr/3LVIt948+AL1fVW1Z5fB7H6og5behYTfNM88AZ5b3AncCNzb+VVRQ/C/xl\n3/POoXdGei/wxhnk9UJ6I8IfAN8Erh7MCziF3sqLLwJfmnZeo+Q062MFPBL4RBPv48Ajmu1PBt7Z\n3H4qcEtznG4GXj6lXH7s5wZ2AP+iuX0M8L7m8c8CJ8/g92itnC6lN0j4IvDXwKkzyGkXvVHovcA+\n4ELg1cCr+p5zGb0VRjdzhFVqs8oJ+I2+43Q9cOYMcno6vTbuK3/jNzb/P+d2rEbJaSPHauTmaZKk\n5eY1bSWpJSz4ktQSFnxJagkLviS1hAVfklrCgi9JLWHB18Qk+d6E9nN1knuSDH556eQkn02vRfN7\nmm8irjz2uvSuq0ySnUm+lsNtmjf8jehJSrItyf4k25O8vMnti0nuTa/19o1JLl3ltQ9N8u0kxw5s\nvyrJC5Nc0LTu/VCz/f7jIa1wHb4mJsl3q2rodQ7WuZ9nAscCr66qF/Rtfy/wgap6f5J3ADdV1Z80\n7RNuBJ5YVT9KspNe29grxs1lIK8tdbjj60Zevw34XlX94cD2rwFPrjXaPzQ//4er6j3N/UfQ+7LX\no6vqh0nOAn6jql7UfCP801X15I3mq83HEb6mKslJST7RdIn8+EqXyCSPSe/iJDcn+Z3+TwdVdS0w\n2CgKej2YPtjc7m/R/CzgCwPF+Md+t5sR9uVJrk3ylSQX9T32q0k+14yy35EkzfbvJfn9JF+k15nw\neeldyOeGJG9JcmV67kjyqOY1aUbbo3bofECrgSQPaT6lfDbJF5I8v3loN3B+3+t+hV6P9h8O7rCq\n/g+wPzO4qImWhwVf03YZ8OdV9QR6X6t/a7P9LcB/qap/Qq8x1RE/ajbF9J6+ot7fgvjp9DoH9vu9\nvimdd/dtfxzwS/QuYrEtyYOSnEbvQilPq6onAT8CfrV5/kOAz1TVE5sYfwycXVX/lKZFdPU+Jr8b\nWJlCeTa9Tx8b7QHzH+m103gKvY6af9j06PkocGYOX23sPOA9R9jPF4B/vsEctAlZ8DVtT+VwUXo3\nveK8sv0Dze1hrbIHHak97c8Cfzvw2Our6klV9cSq+td92z9aVX9fVf8LOEivHfRZwJOAG5qR/LPo\n9U+CXj+TDzW3TwO+WlX7mvv9xXYnsBLnFc39jXoOcEmTy7XA0fTaTv+QXtF/UZKfAX6BXl+e1XwL\n+Lkx8tAmM9H2yNIQgyP3YSP5Na/aU1Xfbrpxrsyj97en/QHw4BHz6e/Ueoje30CAd1XVJUOe/4M6\nfKJr1S6PVbU/vRbSz6TXOvuCEfNZzQtr+HUEdgOvp3eO40NrnFN4ML1jIwGO8DVZw4rh9Ryed34p\nh68h8Bngxc3t81bZ1+D+rqXX7hge2KL5NuCxI+SyWr5/Dbw4vYvbrFyw/dEDzwHYA5ySwxdFf8nA\n/i4H/jvw3r43iY34GPC6+5N84Dz8J4DT6XVNPNJ0DvQuhjGrNt5aAhZ8TdJPJNmX5K7mv/8W+E3g\nwiQ30ZsXXylk/w749832fwT875WdJPk08F7gWc1+Vi7M8sbmNXfQa9d8ebP9anoXdO73u31z+Df2\nL+HsUwBVdRu9q3pdk+Rm4Bp600T3P6d53v8DXktzqUvgu/150+uZ/hDgz9c8UkPy6PPbwLFJbkny\nJWBbXw4/Aq4AfrKqrl9jv0/lyFM+ahmXZWoukvxEVf2guf0S4Lyq+uUx9vdB4A1V9dVJ5bhKnIc0\nK2BI8jbgjmouUJHkF4E/qKrBN5+V124Dvl9VfzCl3J4NvLZZlvmLwGuq6t9MI5aWkyN8zcuTm6Wa\nNwOvAf7DmPt7I4dH5dP0682nhluBhwF/ApDkYuD9TR6r+X7z+u2TTirJBfRWPq2sDHokfZ8MJHCE\nL0mt4QhfklrCgi9JLWHBl6SWsOBLUktY8CWpJf4/YSiW8RnJDJkAAAAASUVORK5CYII=\n", "text/plain": [ - "" + "" ] }, "metadata": {}, @@ -440,7 +446,7 @@ "output_type": "stream", "text": [ "=== GOptimizerLM ===\n", - " Optimized function value ..: 755051.482\n", + " Optimized function value ..: 766517.033\n", " Absolute precision ........: 0.005\n", " Acceptable value decrease .: 2\n", " Optimization status .......: converged\n", @@ -486,8 +492,8 @@ " RA .......................: 83.6331 [-360,360] deg (fixed,scale=1)\n", " DEC ......................: 22.0145 [-90,90] deg (fixed,scale=1)\n", " Number of spectral par's ..: 3\n", - " Prefactor ................: 5.69396634052842e-16 +/- 6.14427401423559e-18 [1e-23,1e-13] ph/cm2/s/MeV (free,scale=1e-16,gradient)\n", - " Index ....................: -2.4703556883579 +/- 0.00792538918685176 [-0,-5] (free,scale=-1,gradient)\n", + " Prefactor ................: 5.71641985530777e-16 +/- 6.16044729438793e-18 [1e-23,1e-13] ph/cm2/s/MeV (free,scale=1e-16,gradient)\n", + " Index ....................: -2.47207574641491 +/- 0.00792255648868446 [-0,-5] (free,scale=-1,gradient)\n", " PivotEnergy ..............: 300000 [10000,1000000000] MeV (fixed,scale=1000000,gradient)\n", " Number of temporal par's ..: 1\n", " Normalization ............: 1 (relative value) (fixed,scale=1,gradient)\n", @@ -499,8 +505,8 @@ " Model type ................: \"PowerLaw\" * \"Constant\"\n", " Number of parameters ......: 4\n", " Number of spectral par's ..: 3\n", - " Prefactor ................: 0.9951294839032 +/- 0.00725734189801964 [0.001,1000] ph/cm2/s/MeV (free,scale=1,gradient)\n", - " Index ....................: -0.00197954530750781 +/- 0.00257407441298647 [-5,5] (free,scale=1,gradient)\n", + " Prefactor ................: 0.999546835049722 +/- 0.00722086200687532 [0.001,1000] ph/cm2/s/MeV (free,scale=1,gradient)\n", + " Index ....................: -0.000149670672999089 +/- 0.00255125105704583 [-5,5] (free,scale=1,gradient)\n", " PivotEnergy ..............: 1000000 [10000,1000000000] MeV (fixed,scale=1000000,gradient)\n", " Number of temporal par's ..: 1\n", " Normalization ............: 1 (relative value) (fixed,scale=1,gradient)\n" @@ -527,7 +533,7 @@ "name": "stdout", "output_type": "stream", "text": [ - "755051.481962\n" + "766517.033411\n" ] } ], @@ -564,11 +570,11 @@ " Model components ..........: \"PointSource\" * \"PowerLaw\" * \"Constant\"\n", " Number of parameters ......: 6\n", " Number of spatial par's ...: 2\n", - " RA .......................: 83.6324977668637 +/- 0.000645834980167342 [-360,360] deg (free,scale=1)\n", - " DEC ......................: 22.0156097018479 +/- 0.000603833638486857 [-90,90] deg (free,scale=1)\n", + " RA .......................: 83.6321512928169 +/- 0.000642709898034108 [-360,360] deg (free,scale=1)\n", + " DEC ......................: 22.0153650383811 +/- 0.000605599763935102 [-90,90] deg (free,scale=1)\n", " Number of spectral par's ..: 3\n", - " Prefactor ................: 5.69123259664292e-16 +/- 6.14234382363775e-18 [1e-23,1e-13] ph/cm2/s/MeV (free,scale=1e-16,gradient)\n", - " Index ....................: -2.47037429659339 +/- 0.00792634862711166 [-0,-5] (free,scale=-1,gradient)\n", + " Prefactor ................: 5.71414829573802e-16 +/- 6.15917758147358e-18 [1e-23,1e-13] ph/cm2/s/MeV (free,scale=1e-16,gradient)\n", + " Index ....................: -2.47211135009306 +/- 0.0079230681070739 [-0,-5] (free,scale=-1,gradient)\n", " PivotEnergy ..............: 300000 [10000,1000000000] MeV (fixed,scale=1000000,gradient)\n", " Number of temporal par's ..: 1\n", " Normalization ............: 1 (relative value) (fixed,scale=1,gradient)\n", @@ -580,8 +586,8 @@ " Model type ................: \"PowerLaw\" * \"Constant\"\n", " Number of parameters ......: 4\n", " Number of spectral par's ..: 3\n", - " Prefactor ................: 0.995120900180708 +/- 0.00725727612144389 [0.001,1000] ph/cm2/s/MeV (free,scale=1,gradient)\n", - " Index ....................: -0.00198038808966019 +/- 0.00257407550475833 [-5,5] (free,scale=1,gradient)\n", + " Prefactor ................: 0.999542788067004 +/- 0.00722079542969056 [0.001,1000] ph/cm2/s/MeV (free,scale=1,gradient)\n", + " Index ....................: -0.00014887248036908 +/- 0.00255123898184119 [-5,5] (free,scale=1,gradient)\n", " PivotEnergy ..............: 1000000 [10000,1000000000] MeV (fixed,scale=1000000,gradient)\n", " Number of temporal par's ..: 1\n", " Normalization ............: 1 (relative value) (fixed,scale=1,gradient)\n" @@ -611,7 +617,7 @@ "name": "stdout", "output_type": "stream", "text": [ - "4.27571920515\n" + "4.2173620963\n" ] } ], diff --git a/doc/source/users/tutorials/quickstart/energy_dispersion.rst b/doc/source/users/tutorials/quickstart/energy_dispersion.rst index 5a13599..edfcf00 100755 --- a/doc/source/users/tutorials/quickstart/energy_dispersion.rst +++ b/doc/source/users/tutorials/quickstart/energy_dispersion.rst @@ -48,9 +48,7 @@ You then select events from the simulated data as before using the $ ctselect Input event list or observation definition XML file [events.fits] events_edisp.fits - RA for ROI centre (degrees) (0-360) [83.63] - Dec for ROI centre (degrees) (-90-90) [22.51] - Radius of ROI (degrees) (0-180) [3.0] + Radius of ROI around pointing or specified RA/DEC (degrees) (0-180) [3.0] Start time (UTC string, JD, MJD or MET in seconds) [2020-01-01T00:10:00] Stop time (UTC string, JD, MJD or MET in seconds) [2020-01-01T00:40:00] Lower energy limit (TeV) [0.1] @@ -64,10 +62,10 @@ Then you bin as before the selected events into a counts cube using the $ ctbin Input event list or observation definition XML file [selected_events.fits] selected_events_edisp.fits + 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] Second coordinate of image center in degrees (DEC or galactic b) (-90-90) [22.51] - 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] Size of the X axis in pixels [200] Size of the Y axis in pixels [200] @@ -75,7 +73,7 @@ Then you bin as before the selected events into a counts cube using the Lower energy limit (TeV) [0.1] Upper energy limit (TeV) [100.0] Number of energy bins (1-200) [20] - Output counts cube file [cntcube.fits] cntcube_edisp.fits + Output counts cube file or observation definition XML file [cntcube.fits] cntcube_edisp.fits As next step you need to compute the :ref:`energy dispersion cube ` @@ -84,22 +82,22 @@ using the :ref:`ctexpcube` tool. You run the tool as follows: .. code-block:: bash $ ctedispcube - Input event list or observation definition XML file [NONE] selected_events_edisp.fits - Calibration database [prod2] - Instrument response function [South_0.5h] - Input counts cube file to extract energy dispersion cube definition [NONE] - First coordinate of image center in degrees (RA or galactic l) (0-360) [83.63] - Second coordinate of image center in degrees (DEC or galactic b) (-90-90) [22.51] - 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) [1.0] - Size of the X axis in pixels [10] - Size of the Y axis in pixels [10] - Algorithm for defining energy bins (FILE|LIN|LOG) [LOG] - Lower energy limit (TeV) [0.1] - Upper energy limit (TeV) [100.0] - Number of energy bins (1-200) [20] - Output energy dispersion cube file [edispcube.fits] + Input event list or observation definition XML file [NONE] selected_events_edisp.fits + Calibration database [prod2] + Instrument response function [South_0.5h] + Input counts cube file to extract energy dispersion cube definition [NONE] + 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] + Second coordinate of image center in degrees (DEC or galactic b) (-90-90) [22.51] + Image scale (in degrees/pixel) [1.0] + Size of the X axis in pixels [10] + Size of the Y axis in pixels [10] + Algorithm for defining energy bins (FILE|LIN|LOG) [LOG] + Lower energy limit (TeV) [0.1] + Upper energy limit (TeV) [100.0] + Number of energy bins (1-1000) [20] + Output energy dispersion cube file [edispcube.fits] Now you are ready to perform a binned maximum likelihood analysis taking the energy dispersion into account. You do this by running the :ref:`ctlike` tool @@ -112,11 +110,71 @@ energy dispersion cube: Input event list, counts cube or observation definition XML file [selected_events.fits] cntcube_edisp.fits Input exposure cube file [expcube.fits] Input PSF cube file [psfcube.fits] - Input energy dispersion cube file [NONE] edispcube.fits Input background cube file [bkgcube.fits] + Input energy dispersion cube file [NONE] edispcube.fits Input model definition XML file [$CTOOLS/share/models/crab.xml] models.xml Output model definition XML file [crab_results.xml] crab_results_edisp.xml +And here is the output in the log file: + +.. code-block:: none + + 2019-04-02T14:20:38: +=================================+ + 2019-04-02T14:20:38: | Maximum likelihood optimisation | + 2019-04-02T14:20:38: +=================================+ + 2019-04-02T14:21:49: >Iteration 0: -logL=57889.195, Lambda=1.0e-03 + 2019-04-02T14:22:51: >Iteration 1: -logL=57887.419, Lambda=1.0e-03, delta=1.776, step=1.0e+00, max(|grad|)=1.111374 [Index:3] + 2019-04-02T14:23:54: >Iteration 2: -logL=57887.417, Lambda=1.0e-04, delta=0.002, step=1.0e+00, max(|grad|)=-0.005923 [Index:7] + 2019-04-02T14:24:57: + 2019-04-02T14:24:57: +=========================================+ + 2019-04-02T14:24:57: | Maximum likelihood optimisation results | + 2019-04-02T14:24:57: +=========================================+ + 2019-04-02T14:24:57: === GOptimizerLM === + 2019-04-02T14:24:57: Optimized function value ..: 57887.417 + 2019-04-02T14:24:57: Absolute precision ........: 0.005 + 2019-04-02T14:24:57: Acceptable value decrease .: 2 + 2019-04-02T14:24:57: Optimization status .......: converged + 2019-04-02T14:24:57: Number of parameters ......: 10 + 2019-04-02T14:24:57: Number of free parameters .: 4 + 2019-04-02T14:24:57: Number of iterations ......: 2 + 2019-04-02T14:24:57: Lambda ....................: 1e-05 + 2019-04-02T14:24:57: Maximum log likelihood ....: -57887.417 + 2019-04-02T14:24:57: Observed events (Nobs) ...: 19137.000 + 2019-04-02T14:24:57: Predicted events (Npred) ..: 19136.996 (Nobs - Npred = 0.00354148293627077) + 2019-04-02T14:24:57: === GModels === + 2019-04-02T14:24:57: Number of models ..........: 2 + 2019-04-02T14:24:57: Number of parameters ......: 10 + 2019-04-02T14:24:57: === GModelSky === + 2019-04-02T14:24:57: Name ......................: Crab + 2019-04-02T14:24:57: Instruments ...............: all + 2019-04-02T14:24:57: Instrument scale factors ..: unity + 2019-04-02T14:24:57: Observation identifiers ...: all + 2019-04-02T14:24:57: Model type ................: PointSource + 2019-04-02T14:24:57: Model components ..........: "PointSource" * "PowerLaw" * "Constant" + 2019-04-02T14:24:57: Number of parameters ......: 6 + 2019-04-02T14:24:57: Number of spatial par's ...: 2 + 2019-04-02T14:24:57: RA .......................: 83.6331 [-360,360] deg (fixed,scale=1) + 2019-04-02T14:24:57: DEC ......................: 22.0145 [-90,90] deg (fixed,scale=1) + 2019-04-02T14:24:57: Number of spectral par's ..: 3 + 2019-04-02T14:24:57: Prefactor ................: 5.52559284054621e-16 +/- 9.88229994960437e-18 [1e-23,1e-13] ph/cm2/s/MeV (free,scale=1e-16,gradient) + 2019-04-02T14:24:57: Index ....................: -2.48163444213634 +/- 0.015305403980771 [-0,-5] (free,scale=-1,gradient) + 2019-04-02T14:24:57: PivotEnergy ..............: 300000 [10000,1000000000] MeV (fixed,scale=1000000,gradient) + 2019-04-02T14:24:57: Number of temporal par's ..: 1 + 2019-04-02T14:24:57: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) + 2019-04-02T14:24:57: === GCTAModelCubeBackground === + 2019-04-02T14:24:57: Name ......................: BackgroundModel + 2019-04-02T14:24:57: Instruments ...............: CTA, HESS, MAGIC, VERITAS + 2019-04-02T14:24:57: Instrument scale factors ..: unity + 2019-04-02T14:24:57: Observation identifiers ...: all + 2019-04-02T14:24:57: Model type ................: "PowerLaw" * "Constant" + 2019-04-02T14:24:57: Number of parameters ......: 4 + 2019-04-02T14:24:57: Number of spectral par's ..: 3 + 2019-04-02T14:24:57: Prefactor ................: 1.00540991217377 +/- 0.0157241034891596 [0.01,100] ph/cm2/s/MeV (free,scale=1,gradient) + 2019-04-02T14:24:57: Index ....................: 0.00380886384530723 +/- 0.00942814666809632 [-5,5] (free,scale=1,gradient) + 2019-04-02T14:24:57: PivotEnergy ..............: 1000000 MeV (fixed,scale=1000000,gradient) + 2019-04-02T14:24:57: Number of temporal par's ..: 1 + 2019-04-02T14:24:57: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) + You can also perform an unbinned maximum likelihood analysis taking the energy dispersion into account. In that case the energy dispersion information will be directly determined from the @@ -131,3 +189,65 @@ and no energy dispersion cube is required: Instrument response function [South_0.5h] Input model definition XML file [models.xml] $CTOOLS/share/models/crab.xml Output model definition XML file [crab_results_edisp.xml] + +Here the output in the log file: + +.. code-block:: none + + 2019-04-02T14:26:58: +=================================+ + 2019-04-02T14:26:58: | Maximum likelihood optimisation | + 2019-04-02T14:26:58: +=================================+ + 2019-04-02T14:27:00: >Iteration 0: -logL=143060.165, Lambda=1.0e-03 + 2019-04-02T14:27:02: >Iteration 1: -logL=143059.529, Lambda=1.0e-03, delta=0.636, step=1.0e+00, max(|grad|)=-0.938341 [Prefactor:6] + 2019-04-02T14:27:04: >Iteration 2: -logL=143059.529, Lambda=1.0e-04, delta=0.000, step=1.0e+00, max(|grad|)=-0.002461 [Index:7] + 2019-04-02T14:27:05: + 2019-04-02T14:27:05: +=========================================+ + 2019-04-02T14:27:05: | Maximum likelihood optimisation results | + 2019-04-02T14:27:05: +=========================================+ + 2019-04-02T14:27:05: === GOptimizerLM === + 2019-04-02T14:27:05: Optimized function value ..: 143059.529 + 2019-04-02T14:27:05: Absolute precision ........: 0.005 + 2019-04-02T14:27:05: Acceptable value decrease .: 2 + 2019-04-02T14:27:05: Optimization status .......: converged + 2019-04-02T14:27:05: Number of parameters ......: 10 + 2019-04-02T14:27:05: Number of free parameters .: 4 + 2019-04-02T14:27:05: Number of iterations ......: 2 + 2019-04-02T14:27:05: Lambda ....................: 1e-05 + 2019-04-02T14:27:05: Maximum log likelihood ....: -143059.529 + 2019-04-02T14:27:05: Observed events (Nobs) ...: 22407.000 + 2019-04-02T14:27:05: Predicted events (Npred) ..: 22406.999 (Nobs - Npred = 0.000615512442891486) + 2019-04-02T14:27:05: === GModels === + 2019-04-02T14:27:05: Number of models ..........: 2 + 2019-04-02T14:27:05: Number of parameters ......: 10 + 2019-04-02T14:27:05: === GModelSky === + 2019-04-02T14:27:05: Name ......................: Crab + 2019-04-02T14:27:05: Instruments ...............: all + 2019-04-02T14:27:05: Instrument scale factors ..: unity + 2019-04-02T14:27:05: Observation identifiers ...: all + 2019-04-02T14:27:05: Model type ................: PointSource + 2019-04-02T14:27:05: Model components ..........: "PointSource" * "PowerLaw" * "Constant" + 2019-04-02T14:27:05: Number of parameters ......: 6 + 2019-04-02T14:27:05: Number of spatial par's ...: 2 + 2019-04-02T14:27:05: RA .......................: 83.6331 [-360,360] deg (fixed,scale=1) + 2019-04-02T14:27:05: DEC ......................: 22.0145 [-90,90] deg (fixed,scale=1) + 2019-04-02T14:27:05: Number of spectral par's ..: 3 + 2019-04-02T14:27:05: Prefactor ................: 5.61701723486666e-16 +/- 1.00237442767986e-17 [1e-23,1e-13] ph/cm2/s/MeV (free,scale=1e-16,gradient) + 2019-04-02T14:27:05: Index ....................: -2.48356027289941 +/- 0.0152626158555975 [-0,-5] (free,scale=-1,gradient) + 2019-04-02T14:27:05: PivotEnergy ..............: 300000 [10000,1000000000] MeV (fixed,scale=1000000,gradient) + 2019-04-02T14:27:05: Number of temporal par's ..: 1 + 2019-04-02T14:27:05: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) + 2019-04-02T14:27:05: === GCTAModelIrfBackground === + 2019-04-02T14:27:05: Name ......................: CTABackgroundModel + 2019-04-02T14:27:05: Instruments ...............: CTA + 2019-04-02T14:27:05: Instrument scale factors ..: unity + 2019-04-02T14:27:05: Observation identifiers ...: all + 2019-04-02T14:27:05: Model type ................: "PowerLaw" * "Constant" + 2019-04-02T14:27:05: Number of parameters ......: 4 + 2019-04-02T14:27:05: Number of spectral par's ..: 3 + 2019-04-02T14:27:05: Prefactor ................: 1.0079566950951 +/- 0.0133706835965654 [0.001,1000] ph/cm2/s/MeV (free,scale=1,gradient) + 2019-04-02T14:27:05: Index ....................: 0.00303753293882809 +/- 0.00807470190154737 [-5,5] (free,scale=1,gradient) + 2019-04-02T14:27:05: PivotEnergy ..............: 1000000 [10000,1000000000] MeV (fixed,scale=1000000,gradient) + 2019-04-02T14:27:05: Number of temporal par's ..: 1 + 2019-04-02T14:27:05: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) + + diff --git a/doc/source/users/tutorials/quickstart/fitting.rst b/doc/source/users/tutorials/quickstart/fitting.rst index 227d56a..557803b 100755 --- a/doc/source/users/tutorials/quickstart/fitting.rst +++ b/doc/source/users/tutorials/quickstart/fitting.rst @@ -53,8 +53,8 @@ that was produced by the run: - - + + @@ -64,8 +64,8 @@ that was produced by the run: - - + + @@ -100,61 +100,61 @@ Below the last lines of the log file that was produced by this run: .. code-block:: none - 2018-01-24T14:29:51: +=================================+ - 2018-01-24T14:29:51: | Maximum likelihood optimisation | - 2018-01-24T14:29:51: +=================================+ - 2018-01-24T14:29:51: >Iteration 0: -logL=56469.945, Lambda=1.0e-03 - 2018-01-24T14:29:52: >Iteration 1: -logL=56468.674, Lambda=1.0e-03, delta=1.271, step=1.0e+00, max(|grad|)=-1.725465 [Prefactor:6] - 2018-01-24T14:29:53: >Iteration 2: -logL=56468.674, Lambda=1.0e-04, delta=0.000, step=1.0e+00, max(|grad|)=-0.020159 [Index:7] - 2018-01-24T14:29:53: - 2018-01-24T14:29:53: +=========================================+ - 2018-01-24T14:29:53: | Maximum likelihood optimisation results | - 2018-01-24T14:29:53: +=========================================+ - 2018-01-24T14:29:53: === GOptimizerLM === - 2018-01-24T14:29:53: Optimized function value ..: 56468.674 - 2018-01-24T14:29:53: Absolute precision ........: 0.005 - 2018-01-24T14:29:53: Acceptable value decrease .: 2 - 2018-01-24T14:29:53: Optimization status .......: converged - 2018-01-24T14:29:53: Number of parameters ......: 10 - 2018-01-24T14:29:53: Number of free parameters .: 4 - 2018-01-24T14:29:53: Number of iterations ......: 2 - 2018-01-24T14:29:53: Lambda ....................: 1e-05 - 2018-01-24T14:29:53: Maximum log likelihood ....: -56468.674 - 2018-01-24T14:29:53: Observed events (Nobs) ...: 18799.000 - 2018-01-24T14:29:53: Predicted events (Npred) ..: 18798.998 (Nobs - Npred = 0.00197278879204532) - 2018-01-24T14:29:53: === GModels === - 2018-01-24T14:29:53: Number of models ..........: 2 - 2018-01-24T14:29:53: Number of parameters ......: 10 - 2018-01-24T14:29:53: === GModelSky === - 2018-01-24T14:29:53: Name ......................: Crab - 2018-01-24T14:29:53: Instruments ...............: all - 2018-01-24T14:29:53: Instrument scale factors ..: unity - 2018-01-24T14:29:53: Observation identifiers ...: all - 2018-01-24T14:29:53: Model type ................: PointSource - 2018-01-24T14:29:53: Model components ..........: "PointSource" * "PowerLaw" * "Constant" - 2018-01-24T14:29:53: Number of parameters ......: 6 - 2018-01-24T14:29:53: Number of spatial par's ...: 2 - 2018-01-24T14:29:53: RA .......................: 83.6331 [-360,360] deg (fixed,scale=1) - 2018-01-24T14:29:53: DEC ......................: 22.0145 [-90,90] deg (fixed,scale=1) - 2018-01-24T14:29:53: Number of spectral par's ..: 3 - 2018-01-24T14:29:53: Prefactor ................: 5.76059824355157e-16 +/- 1.00804418203193e-17 [1e-23,1e-13] ph/cm2/s/MeV (free,scale=1e-16,gradient) - 2018-01-24T14:29:53: Index ....................: -2.48246178370535 +/- 0.0150306767905499 [-0,-5] (free,scale=-1,gradient) - 2018-01-24T14:29:53: PivotEnergy ..............: 300000 [10000,1000000000] MeV (fixed,scale=1000000,gradient) - 2018-01-24T14:29:53: Number of temporal par's ..: 1 - 2018-01-24T14:29:53: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) - 2018-01-24T14:29:53: === GCTAModelCubeBackground === - 2018-01-24T14:29:53: Name ......................: BackgroundModel - 2018-01-24T14:29:53: Instruments ...............: CTA, HESS, MAGIC, VERITAS - 2018-01-24T14:29:53: Instrument scale factors ..: unity - 2018-01-24T14:29:53: Observation identifiers ...: all - 2018-01-24T14:29:53: Model type ................: "PowerLaw" * "Constant" - 2018-01-24T14:29:53: Number of parameters ......: 4 - 2018-01-24T14:29:53: Number of spectral par's ..: 3 - 2018-01-24T14:29:53: Prefactor ................: 1.00863986588954 +/- 0.0159030370642366 [0.01,100] ph/cm2/s/MeV (free,scale=1,gradient) - 2018-01-24T14:29:53: Index ....................: 0.0113292722406164 +/- 0.00951929088334377 [-5,5] (free,scale=1,gradient) - 2018-01-24T14:29:53: PivotEnergy ..............: 1000000 MeV (fixed,scale=1000000,gradient) - 2018-01-24T14:29:53: Number of temporal par's ..: 1 - 2018-01-24T14:29:53: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) + 2019-04-02T13:55:43: +=================================+ + 2019-04-02T13:55:43: | Maximum likelihood optimisation | + 2019-04-02T13:55:43: +=================================+ + 2019-04-02T13:55:44: >Iteration 0: -logL=57895.557, Lambda=1.0e-03 + 2019-04-02T13:55:46: >Iteration 1: -logL=57893.741, Lambda=1.0e-03, delta=1.815, step=1.0e+00, max(|grad|)=5.091294 [Index:7] + 2019-04-02T13:55:47: >Iteration 2: -logL=57893.740, Lambda=1.0e-04, delta=0.001, step=1.0e+00, max(|grad|)=0.076392 [Index:7] + 2019-04-02T13:55:48: + 2019-04-02T13:55:48: +=========================================+ + 2019-04-02T13:55:48: | Maximum likelihood optimisation results | + 2019-04-02T13:55:48: +=========================================+ + 2019-04-02T13:55:48: === GOptimizerLM === + 2019-04-02T13:55:48: Optimized function value ..: 57893.740 + 2019-04-02T13:55:48: Absolute precision ........: 0.005 + 2019-04-02T13:55:48: Acceptable value decrease .: 2 + 2019-04-02T13:55:48: Optimization status .......: converged + 2019-04-02T13:55:48: Number of parameters ......: 10 + 2019-04-02T13:55:48: Number of free parameters .: 4 + 2019-04-02T13:55:48: Number of iterations ......: 2 + 2019-04-02T13:55:48: Lambda ....................: 1e-05 + 2019-04-02T13:55:48: Maximum log likelihood ....: -57893.740 + 2019-04-02T13:55:48: Observed events (Nobs) ...: 19452.000 + 2019-04-02T13:55:48: Predicted events (Npred) ..: 19451.999 (Nobs - Npred = 0.00144846074545057) + 2019-04-02T13:55:48: === GModels === + 2019-04-02T13:55:48: Number of models ..........: 2 + 2019-04-02T13:55:48: Number of parameters ......: 10 + 2019-04-02T13:55:48: === GModelSky === + 2019-04-02T13:55:48: Name ......................: Crab + 2019-04-02T13:55:48: Instruments ...............: all + 2019-04-02T13:55:48: Instrument scale factors ..: unity + 2019-04-02T13:55:48: Observation identifiers ...: all + 2019-04-02T13:55:48: Model type ................: PointSource + 2019-04-02T13:55:48: Model components ..........: "PointSource" * "PowerLaw" * "Constant" + 2019-04-02T13:55:48: Number of parameters ......: 6 + 2019-04-02T13:55:48: Number of spatial par's ...: 2 + 2019-04-02T13:55:48: RA .......................: 83.6331 [-360,360] deg (fixed,scale=1) + 2019-04-02T13:55:48: DEC ......................: 22.0145 [-90,90] deg (fixed,scale=1) + 2019-04-02T13:55:48: Number of spectral par's ..: 3 + 2019-04-02T13:55:48: Prefactor ................: 5.79620732342312e-16 +/- 1.01207859506902e-17 [1e-23,1e-13] ph/cm2/s/MeV (free,scale=1e-16,gradient) + 2019-04-02T13:55:48: Index ....................: -2.48767829376908 +/- 0.0150183262185825 [-0,-5] (free,scale=-1,gradient) + 2019-04-02T13:55:48: PivotEnergy ..............: 300000 [10000,1000000000] MeV (fixed,scale=1000000,gradient) + 2019-04-02T13:55:48: Number of temporal par's ..: 1 + 2019-04-02T13:55:48: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) + 2019-04-02T13:55:48: === GCTAModelCubeBackground === + 2019-04-02T13:55:48: Name ......................: BackgroundModel + 2019-04-02T13:55:48: Instruments ...............: CTA, HESS, MAGIC, VERITAS + 2019-04-02T13:55:48: Instrument scale factors ..: unity + 2019-04-02T13:55:48: Observation identifiers ...: all + 2019-04-02T13:55:48: Model type ................: "PowerLaw" * "Constant" + 2019-04-02T13:55:48: Number of parameters ......: 4 + 2019-04-02T13:55:48: Number of spectral par's ..: 3 + 2019-04-02T13:55:48: Prefactor ................: 0.997287677101311 +/- 0.0156090791023617 [0.01,100] ph/cm2/s/MeV (free,scale=1,gradient) + 2019-04-02T13:55:48: Index ....................: -0.00905917021700502 +/- 0.00939297053758322 [-5,5] (free,scale=1,gradient) + 2019-04-02T13:55:48: PivotEnergy ..............: 1000000 MeV (fixed,scale=1000000,gradient) + 2019-04-02T13:55:48: Number of temporal par's ..: 1 + 2019-04-02T13:55:48: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) The maximum likelihood optimizer required 2 iterations to converge. This is pretty fast, but recall that you used the same model file for the simulation @@ -200,13 +200,13 @@ are not too far from the best fitting values). .. code-block:: none - 2018-01-24T14:32:49: Name ......................: Crab - 2018-01-24T14:32:49: Instruments ...............: all - 2018-01-24T14:32:49: Test Statistic ............: 21626.0036468824 + 2019-04-02T13:58:58: Name ......................: Crab + 2019-04-02T13:58:58: Instruments ...............: all + 2019-04-02T13:58:58: Test Statistic ............: 21529.9862965408 The Test Statistic value will also be added as new attribute ``ts`` to the XML result file: .. code-block:: xml - + diff --git a/doc/source/users/tutorials/quickstart/onoff.rst b/doc/source/users/tutorials/quickstart/onoff.rst index 8e11279..fdfbaf2 100644 --- a/doc/source/users/tutorials/quickstart/onoff.rst +++ b/doc/source/users/tutorials/quickstart/onoff.rst @@ -49,12 +49,15 @@ To derive On/Off observations from the ``events_edisp.fits`` event list, type: Input event list or observation definition XML file [obs.xml] events_edisp.fits Calibration database [prod2] Instrument response function [South_0.5h] + Input model definition XML file (if NONE, use point source) [NONE] $CTOOLS/share/models/crab.xml + Source name [Crab] Algorithm for defining energy bins (FILE|LIN|LOG) [LOG] Start value for first energy bin in TeV [0.1] Stop value for last energy bin in TeV [100.0] Number of energy bins [120] 30 Stack multiple observations into single PHA, ARF and RMF files? [no] Output observation definition XML file [onoff_obs.xml] + Output model definition XML file [onoff_model.xml] Method for background estimation (REFLECTED|CUSTOM) [REFLECTED] Coordinate system (CEL - celestial, GAL - galactic) (CEL|GAL) [CEL] Right Ascension of source region centre (deg) (0-360) [83.63] @@ -121,6 +124,34 @@ into the spectral response computation. location) can be set using the hidden ``prefix`` parameter of :ref:`csphagen`. +:ref:`csphagen` also produced the +:ref:`model definition XML file ` +``onoff_model.xml`` that can be directly used for model fitting: + +.. code-block:: xml + + + + + + + + + + + + + + + + + + + + + + + There are also come ancillary `ds9 `_ region files, that show the On region and the Off regions, ``onoff_on.reg`` and ``onoff_off.reg``, respectively. Below there is @@ -135,43 +166,8 @@ the On and Off regions. the green circles the Off regions, and the white circle the On region.* -Now you are ready to fit a model to the On/Off data. As usual, you need to -create for this purpose a -:ref:`model definition file `. -So far, only fitting of point sources is supported, hence you need a single -point source as celestial component of the -:ref:`model definition file `. -As background component you need to add a ``CTAIrfBackground`` model for -the ``CTAOnOff`` instrument, which is the instrument code that defines an -On/Off analysis. Your -:ref:`model definition XML file ` -should look as follows: - -.. code-block:: bash - - - - - - - - - - - - - - - - - - - - - - - -You can now fit this model using an On/Off analysis by specifying the +You can now fit the model ``onoff_model.xml`` using an On/Off analysis by +specifying the :ref:`output observation definition file ` and the :ref:`model definition file ` @@ -180,9 +176,9 @@ to :ref:`ctlike`: .. code-block:: bash $ ctlike - Input event list, counts cube or observation definition XML file [selected_events.fits] onoff_obs.xml - Input model definition XML file [$CTOOLS/share/models/crab.xml] crab_onoff.xml - Output model definition XML file [crab_results.xml] + Input event list, counts cube or observation definition XML file [selected_events_edisp.fits] onoff_obs.xml + Input model definition XML file [$CTOOLS/share/models/crab.xml] onoff_model.xml + Output model definition XML file [crab_results_edisp.xml] crab_results.xml Below you see the corresponding output from the ``ctlike.log`` file. The fitted parameters are still the same within statistical uncertainties as the ones @@ -191,61 +187,62 @@ the background is not well known a priori. .. code-block:: none - 2018-01-24T15:03:20: +=================================+ - 2018-01-24T15:03:20: | Maximum likelihood optimisation | - 2018-01-24T15:03:20: +=================================+ - 2018-01-24T15:03:20: >Iteration 0: -logL=-44981.922, Lambda=1.0e-03 - 2018-01-24T15:03:20: >Iteration 1: -logL=-44983.914, Lambda=1.0e-03, delta=1.992, step=1.0e+00, max(|grad|)=2.056476 [Index:7] - 2018-01-24T15:03:20: >Iteration 2: -logL=-44983.915, Lambda=1.0e-04, delta=0.001, step=1.0e+00, max(|grad|)=0.001765 [Index:3] - 2018-01-24T15:03:20: - 2018-01-24T15:03:20: +=========================================+ - 2018-01-24T15:03:20: | Maximum likelihood optimisation results | - 2018-01-24T15:03:20: +=========================================+ - 2018-01-24T15:03:20: === GOptimizerLM === - 2018-01-24T15:03:20: Optimized function value ..: -44983.915 - 2018-01-24T15:03:20: Absolute precision ........: 0.005 - 2018-01-24T15:03:20: Acceptable value decrease .: 2 - 2018-01-24T15:03:20: Optimization status .......: converged - 2018-01-24T15:03:20: Number of parameters ......: 10 - 2018-01-24T15:03:20: Number of free parameters .: 4 - 2018-01-24T15:03:20: Number of iterations ......: 2 - 2018-01-24T15:03:20: Lambda ....................: 1e-05 - 2018-01-24T15:03:20: Maximum log likelihood ....: 44983.915 - 2018-01-24T15:03:20: Observed events (Nobs) ...: 7309.000 - 2018-01-24T15:03:20: Predicted events (Npred) ..: 7306.626 (Nobs - Npred = 2.37371839911157) - 2018-01-24T15:03:20: === GModels === - 2018-01-24T15:03:20: Number of models ..........: 2 - 2018-01-24T15:03:20: Number of parameters ......: 10 - 2018-01-24T15:03:20: === GModelSky === - 2018-01-24T15:03:20: Name ......................: Crab - 2018-01-24T15:03:20: Instruments ...............: all - 2018-01-24T15:03:20: Instrument scale factors ..: unity - 2018-01-24T15:03:20: Observation identifiers ...: all - 2018-01-24T15:03:20: Model type ................: PointSource - 2018-01-24T15:03:20: Model components ..........: "PointSource" * "PowerLaw" * "Constant" - 2018-01-24T15:03:20: Number of parameters ......: 6 - 2018-01-24T15:03:20: Number of spatial par's ...: 2 - 2018-01-24T15:03:20: RA .......................: 83.6331 [-360,360] deg (fixed,scale=1) - 2018-01-24T15:03:20: DEC ......................: 22.0145 [-90,90] deg (fixed,scale=1) - 2018-01-24T15:03:20: Number of spectral par's ..: 3 - 2018-01-24T15:03:20: Prefactor ................: 5.59354639816441e-16 +/- 7.3028712529757e-18 [1e-23,1e-13] ph/cm2/s/MeV (free,scale=1e-16,gradient) - 2018-01-24T15:03:20: Index ....................: -2.47773254234859 +/- 0.0110401748178683 [-0,-5] (free,scale=-1,gradient) - 2018-01-24T15:03:20: PivotEnergy ..............: 300000 [10000,1000000000] MeV (fixed,scale=1000000,gradient) - 2018-01-24T15:03:20: Number of temporal par's ..: 1 - 2018-01-24T15:03:20: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) - 2018-01-24T15:03:20: === GCTAModelIrfBackground === - 2018-01-24T15:03:20: Name ......................: BackgroundModel - 2018-01-24T15:03:20: Instruments ...............: CTAOnOff - 2018-01-24T15:03:20: Instrument scale factors ..: unity - 2018-01-24T15:03:20: Observation identifiers ...: all - 2018-01-24T15:03:20: Model type ................: "PowerLaw" * "Constant" - 2018-01-24T15:03:20: Number of parameters ......: 4 - 2018-01-24T15:03:20: Number of spectral par's ..: 3 - 2018-01-24T15:03:20: Prefactor ................: 0.977042261939201 +/- 0.0498175569832557 [0.01,100] ph/cm2/s/MeV (free,scale=1,gradient) - 2018-01-24T15:03:20: Index ....................: 0.00289056644049351 +/- 0.0299910588265984 [-5,5] (free,scale=1,gradient) - 2018-01-24T15:03:20: PivotEnergy ..............: 1000000 MeV (fixed,scale=1000000,gradient) - 2018-01-24T15:03:20: Number of temporal par's ..: 1 - 2018-01-24T15:03:20: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) + 2019-04-02T14:55:29: +=================================+ + 2019-04-02T14:55:29: | Maximum likelihood optimisation | + 2019-04-02T14:55:29: +=================================+ + 2019-04-02T14:55:29: >Iteration 0: -logL=-47436.484, Lambda=1.0e-03 + 2019-04-02T14:55:29: >Iteration 1: -logL=-47439.247, Lambda=1.0e-03, delta=2.762, step=1.0e+00, max(|grad|)=14.136296 [Index:7] + 2019-04-02T14:55:29: >Iteration 2: -logL=-47439.266, Lambda=1.0e-04, delta=0.020, step=1.0e+00, max(|grad|)=0.089563 [Index:7] + 2019-04-02T14:55:29: >Iteration 3: -logL=-47439.266, Lambda=1.0e-05, delta=0.000, step=1.0e+00, max(|grad|)=0.001727 [Index:7] + 2019-04-02T14:55:29: + 2019-04-02T14:55:29: +=========================================+ + 2019-04-02T14:55:29: | Maximum likelihood optimisation results | + 2019-04-02T14:55:29: +=========================================+ + 2019-04-02T14:55:29: === GOptimizerLM === + 2019-04-02T14:55:29: Optimized function value ..: -47439.266 + 2019-04-02T14:55:29: Absolute precision ........: 0.005 + 2019-04-02T14:55:29: Acceptable value decrease .: 2 + 2019-04-02T14:55:29: Optimization status .......: converged + 2019-04-02T14:55:29: Number of parameters ......: 10 + 2019-04-02T14:55:29: Number of free parameters .: 4 + 2019-04-02T14:55:29: Number of iterations ......: 3 + 2019-04-02T14:55:29: Lambda ....................: 1e-06 + 2019-04-02T14:55:29: Maximum log likelihood ....: 47439.266 + 2019-04-02T14:55:29: Observed events (Nobs) ...: 7607.000 + 2019-04-02T14:55:29: Predicted events (Npred) ..: 7606.425 (Nobs - Npred = 0.575097306655152) + 2019-04-02T14:55:29: === GModels === + 2019-04-02T14:55:29: Number of models ..........: 2 + 2019-04-02T14:55:29: Number of parameters ......: 10 + 2019-04-02T14:55:29: === GModelSky === + 2019-04-02T14:55:29: Name ......................: Crab + 2019-04-02T14:55:29: Instruments ...............: all + 2019-04-02T14:55:29: Instrument scale factors ..: unity + 2019-04-02T14:55:29: Observation identifiers ...: all + 2019-04-02T14:55:29: Model type ................: PointSource + 2019-04-02T14:55:29: Model components ..........: "PointSource" * "PowerLaw" * "Constant" + 2019-04-02T14:55:29: Number of parameters ......: 6 + 2019-04-02T14:55:29: Number of spatial par's ...: 2 + 2019-04-02T14:55:29: RA .......................: 83.6331 [-360,360] deg (fixed,scale=1) + 2019-04-02T14:55:29: DEC ......................: 22.0145 [-90,90] deg (fixed,scale=1) + 2019-04-02T14:55:29: Number of spectral par's ..: 3 + 2019-04-02T14:55:29: Prefactor ................: 5.71422768206296e-16 +/- 7.28119011001326e-18 [1e-23,1e-13] ph/cm2/s/MeV (free,scale=1e-16,gradient) + 2019-04-02T14:55:29: Index ....................: -2.47772427704665 +/- 0.0108450088768338 [-0,-5] (free,scale=-1,gradient) + 2019-04-02T14:55:29: PivotEnergy ..............: 300000 [10000,1000000000] MeV (fixed,scale=1000000,gradient) + 2019-04-02T14:55:29: Number of temporal par's ..: 1 + 2019-04-02T14:55:29: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) + 2019-04-02T14:55:29: === GCTAModelIrfBackground === + 2019-04-02T14:55:29: Name ......................: CTABackgroundModel + 2019-04-02T14:55:29: Instruments ...............: CTAOnOff + 2019-04-02T14:55:29: Instrument scale factors ..: unity + 2019-04-02T14:55:29: Observation identifiers ...: all + 2019-04-02T14:55:29: Model type ................: "PowerLaw" * "Constant" + 2019-04-02T14:55:29: Number of parameters ......: 4 + 2019-04-02T14:55:29: Number of spectral par's ..: 3 + 2019-04-02T14:55:29: Prefactor ................: 0.925471278485926 +/- 0.0482291417226665 [0.001,1000] ph/cm2/s/MeV (free,scale=1,gradient) + 2019-04-02T14:55:29: Index ....................: -0.0649030558071282 +/- 0.0301870339200633 [-5,5] (free,scale=1,gradient) + 2019-04-02T14:55:29: PivotEnergy ..............: 1000000 [10000,1000000000] MeV (fixed,scale=1000000,gradient) + 2019-04-02T14:55:29: Number of temporal par's ..: 1 + 2019-04-02T14:55:29: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) :ref:`ctlike` has a hidden parameter called ``statistic`` that sets the statistic used for the fit. By default, :ref:`ctlike` will use ``CSTAT`` @@ -257,7 +254,70 @@ Alternatively, you can use ``WSTAT`` for an On/Off analysis, which treats the number of background counts in each energy bin as a nuisance parameter that is derived from the On and Off counts by profiling the likelihood function. In this case, the only assumption is that the background rate spectrum is the same -in the On and Off regions. +in the On and Off regions. Below the results for a :ref:`ctlike` run with +the ``statistic=wstat`` option. + +.. code-block:: none + + 2019-04-02T15:56:29: +=================================+ + 2019-04-02T15:56:29: | Maximum likelihood optimisation | + 2019-04-02T15:56:29: +=================================+ + 2019-04-02T15:56:29: Parameter "Prefactor" has zero curvature. Fix parameter. + 2019-04-02T15:56:29: Parameter "Index" has zero curvature. Fix parameter. + 2019-04-02T15:56:29: >Iteration 0: -logL=13.699, Lambda=1.0e-03 + 2019-04-02T15:56:29: >Iteration 1: -logL=13.645, Lambda=1.0e-03, delta=0.054, step=1.0e+00, max(|grad|)=0.226348 [Index:3] + 2019-04-02T15:56:29: >Iteration 2: -logL=13.645, Lambda=1.0e-04, delta=0.000, step=1.0e+00, max(|grad|)=0.001120 [Index:3] + 2019-04-02T15:56:29: Free parameter "Prefactor" after convergence was reached with frozen parameter. + 2019-04-02T15:56:29: Free parameter "Index" after convergence was reached with frozen parameter. + 2019-04-02T15:56:29: + 2019-04-02T15:56:29: +=========================================+ + 2019-04-02T15:56:29: | Maximum likelihood optimisation results | + 2019-04-02T15:56:29: +=========================================+ + 2019-04-02T15:56:29: === GOptimizerLM === + 2019-04-02T15:56:29: Optimized function value ..: 13.645 + 2019-04-02T15:56:29: Absolute precision ........: 0.005 + 2019-04-02T15:56:29: Acceptable value decrease .: 2 + 2019-04-02T15:56:29: Optimization status .......: converged + 2019-04-02T15:56:29: Number of parameters ......: 10 + 2019-04-02T15:56:29: Number of free parameters .: 4 + 2019-04-02T15:56:29: Number of iterations ......: 2 + 2019-04-02T15:56:29: Lambda ....................: 1e-05 + 2019-04-02T15:56:29: Maximum log likelihood ....: -13.645 + 2019-04-02T15:56:29: Observed events (Nobs) ...: 7607.000 + 2019-04-02T15:56:29: Predicted events (Npred) ..: 7606.133 (Nobs - Npred = 0.866926153597888) + 2019-04-02T15:56:29: === GModels === + 2019-04-02T15:56:29: Number of models ..........: 2 + 2019-04-02T15:56:29: Number of parameters ......: 10 + 2019-04-02T15:56:29: === GModelSky === + 2019-04-02T15:56:29: Name ......................: Crab + 2019-04-02T15:56:29: Instruments ...............: all + 2019-04-02T15:56:29: Instrument scale factors ..: unity + 2019-04-02T15:56:29: Observation identifiers ...: all + 2019-04-02T15:56:29: Model type ................: PointSource + 2019-04-02T15:56:29: Model components ..........: "PointSource" * "PowerLaw" * "Constant" + 2019-04-02T15:56:29: Number of parameters ......: 6 + 2019-04-02T15:56:29: Number of spatial par's ...: 2 + 2019-04-02T15:56:29: RA .......................: 83.6331 [-360,360] deg (fixed,scale=1) + 2019-04-02T15:56:29: DEC ......................: 22.0145 [-90,90] deg (fixed,scale=1) + 2019-04-02T15:56:29: Number of spectral par's ..: 3 + 2019-04-02T15:56:29: Prefactor ................: 5.71398803734648e-16 +/- 7.28140878478654e-18 [1e-23,1e-13] ph/cm2/s/MeV (free,scale=1e-16,gradient) + 2019-04-02T15:56:29: Index ....................: -2.47775827196727 +/- 0.0108569325078945 [-0,-5] (free,scale=-1,gradient) + 2019-04-02T15:56:29: PivotEnergy ..............: 300000 [10000,1000000000] MeV (fixed,scale=1000000,gradient) + 2019-04-02T15:56:29: Number of temporal par's ..: 1 + 2019-04-02T15:56:29: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) + 2019-04-02T15:56:29: === GCTAModelIrfBackground === + 2019-04-02T15:56:29: Name ......................: CTABackgroundModel + 2019-04-02T15:56:29: Instruments ...............: CTAOnOff + 2019-04-02T15:56:29: Instrument scale factors ..: unity + 2019-04-02T15:56:29: Observation identifiers ...: all + 2019-04-02T15:56:29: Model type ................: "PowerLaw" * "Constant" + 2019-04-02T15:56:29: Number of parameters ......: 4 + 2019-04-02T15:56:29: Number of spectral par's ..: 3 + 2019-04-02T15:56:29: Prefactor ................: 1 +/- 0 [0.001,1000] ph/cm2/s/MeV (free,scale=1,gradient) + 2019-04-02T15:56:29: Index ....................: 0 +/- 0 [-5,5] (free,scale=1,gradient) + 2019-04-02T15:56:29: PivotEnergy ..............: 1000000 [10000,1000000000] MeV (fixed,scale=1000000,gradient) + 2019-04-02T15:56:29: Number of temporal par's ..: 1 + 2019-04-02T15:56:29: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) .. warning:: Beware that the profiling may yield unphysical results (negative background @@ -270,11 +330,6 @@ in the On and Off regions. `XSPEC manual Appendix B `_ for more information. -Finally, you can also use ``CHI2`` as fit statistic which is a classical chi -square, i.e., a Gaussian signal and Gaussian background. As for ``CSTAT``, a -spectral model for the signal and a spectral model for the background are -jointly fit to the On and Off spectra. - .. note:: Many scripts can also be used in On/Off mode, including :ref:`ctbutterfly` and :ref:`csspec` that were used earlier. It is diff --git a/doc/source/users/tutorials/quickstart/residual_map.jpg b/doc/source/users/tutorials/quickstart/residual_map.jpg deleted file mode 100644 index 8d32ee4..0000000 --- a/doc/source/users/tutorials/quickstart/residual_map.jpg +++ /dev/null diff --git a/doc/source/users/tutorials/quickstart/residual_map.png b/doc/source/users/tutorials/quickstart/residual_map.png new file mode 100644 index 0000000..a752c41 --- /dev/null +++ b/doc/source/users/tutorials/quickstart/residual_map.png diff --git a/doc/source/users/tutorials/quickstart/residual_map.rst b/doc/source/users/tutorials/quickstart/residual_map.rst index d28e5f0..60ba068 100755 --- a/doc/source/users/tutorials/quickstart/residual_map.rst +++ b/doc/source/users/tutorials/quickstart/residual_map.rst @@ -34,11 +34,11 @@ subtracted. You run the :ref:`csresmap` script as follows: The :ref:`csresmap` script produces the FITS file ``resmap.fits`` that contains the residual counts map. The image is displayed below using -`ds9 `_ with a linear color scale and with a small Gaussian +`ds9 `_ with a linear color scale and with a 3 pixel Gaussian kernel smoothing applied. Obviously, there are no significant residuals, which indicates that the model fit was satisfactory. -.. figure:: residual_map.jpg +.. figure:: residual_map.png :width: 400px :align: center diff --git a/doc/source/users/tutorials/quickstart/residual_spectrum.png b/doc/source/users/tutorials/quickstart/residual_spectrum.png index 2fba09a..6c89ed1 100644 --- a/doc/source/users/tutorials/quickstart/residual_spectrum.png +++ b/doc/source/users/tutorials/quickstart/residual_spectrum.png diff --git a/doc/source/users/tutorials/quickstart/response.rst b/doc/source/users/tutorials/quickstart/response.rst index 57d26ce..82ae480 100755 --- a/doc/source/users/tutorials/quickstart/response.rst +++ b/doc/source/users/tutorials/quickstart/response.rst @@ -52,18 +52,18 @@ You run :ref:`ctpsfcube` as follows: Calibration database [prod2] Instrument response function [South_0.5h] Input counts cube file to extract PSF cube definition [NONE] + 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] Second coordinate of image center in degrees (DEC or galactic b) (-90-90) [22.51] - 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) [1.0] Size of the X axis in pixels [10] Size of the Y axis in pixels [10] Algorithm for defining energy bins (FILE|LIN|LOG) [LOG] Lower energy limit (TeV) [0.1] Upper energy limit (TeV) [100.0] - Number of energy bins (1-200) [20] - Output PSF cube file [psfcube.fits] + Number of energy bins (1-1000) [20] + Output PSF cube file [psfcube.fits] This produces the FITS file ``psfcube.fits`` that contains the point spread function as function of sky position and energy. You may have noted in the diff --git a/doc/source/users/tutorials/quickstart/selecting.rst b/doc/source/users/tutorials/quickstart/selecting.rst index ba8a0d1..4cf347d 100755 --- a/doc/source/users/tutorials/quickstart/selecting.rst +++ b/doc/source/users/tutorials/quickstart/selecting.rst @@ -33,43 +33,43 @@ event selection: .. code-block:: none - 2018-01-24T14:07:07: +===================+ - 2018-01-24T14:07:07: | Input observation | - 2018-01-24T14:07:07: +===================+ - 2018-01-24T14:07:07: === GObservations === - 2018-01-24T14:07:07: Number of observations ....: 1 - 2018-01-24T14:07:07: Number of models ..........: 0 - 2018-01-24T14:07:07: Number of observed events .: 202210 - 2018-01-24T14:07:07: Number of predicted events : 0 - 2018-01-24T14:07:07: - 2018-01-24T14:07:07: +=================+ - 2018-01-24T14:07:07: | Event selection | - 2018-01-24T14:07:07: +=================+ - 2018-01-24T14:07:07: === CTA observation (id=000001) === - 2018-01-24T14:07:07: Input filename ............: events.fits - 2018-01-24T14:07:07: Event extension name ......: EVENTS - 2018-01-24T14:07:07: GTI extension name ........: GTI - 2018-01-24T14:07:07: === Event selection === - 2018-01-24T14:07:07: Time range (MJD) ..........: 58849.0077451852 - 58849.0285785185 days - 2018-01-24T14:07:07: Time range (UTC) ..........: 2020-01-01T00:10:00 - 2020-01-01T00:40:00 - 2018-01-24T14:07:07: Time range (MET) ..........: 631109469.184 - 631111269.184 seconds - 2018-01-24T14:07:07: Selected energy range .....: 0.1 - 100 TeV - 2018-01-24T14:07:07: Requested RoI .............: Centre(RA,DEC)=(83.63, 22.51) deg, Radius=3 deg - 2018-01-24T14:07:07: RoI of data ...............: Centre(RA,DEC)=(83.63, 22.51) deg, Radius=5 deg - 2018-01-24T14:07:07: Selected RoI ..............: Centre(RA,DEC)=(83.63, 22.51) deg, Radius=3 deg - 2018-01-24T14:07:07: cfitsio selection .........: TIME >= 631109469.18399978 && TIME <= 631111269.18400002 && ENERGY >= 0.10000000 && ENERGY <= 100.00000000 && ANGSEP(83.630000,22.510000,RA,DEC) <= 3.000000 - 2018-01-24T14:07:07: FITS filename .............: /var/tmp/tmp.0.aFKvXb[EVENTS][TIME >= 631109469.18399978 && TIME <= 631111269.18400002 && ENERGY >= 0.10000000 && ENERGY <= 100.00000000 && ANGSEP(83.630000,22.510000,RA,DEC) <= 3.000000] - 2018-01-24T14:07:07: - 2018-01-24T14:07:07: +====================+ - 2018-01-24T14:07:07: | Output observation | - 2018-01-24T14:07:07: +====================+ - 2018-01-24T14:07:07: === GObservations === - 2018-01-24T14:07:07: Number of observations ....: 1 - 2018-01-24T14:07:07: Number of models ..........: 0 - 2018-01-24T14:07:07: Number of observed events .: 21974 - 2018-01-24T14:07:07: Number of predicted events : 0 + 2019-04-02T12:52:50: +===================+ + 2019-04-02T12:52:50: | Input observation | + 2019-04-02T12:52:50: +===================+ + 2019-04-02T12:52:50: === GObservations === + 2019-04-02T12:52:50: Number of observations ....: 1 + 2019-04-02T12:52:50: Number of models ..........: 0 + 2019-04-02T12:52:50: Number of observed events .: 206048 + 2019-04-02T12:52:50: Number of predicted events : 0 + 2019-04-02T12:52:50: + 2019-04-02T12:52:50: +=================+ + 2019-04-02T12:52:50: | Event selection | + 2019-04-02T12:52:50: +=================+ + 2019-04-02T12:52:50: === CTA observation (id=000001) === + 2019-04-02T12:52:50: Input filename ............: events.fits + 2019-04-02T12:52:50: Event extension name ......: EVENTS + 2019-04-02T12:52:50: GTI extension name ........: GTI + 2019-04-02T12:52:50: === Event selection === + 2019-04-02T12:52:50: Time range (MJD) ..........: 58849.0077451852 - 58849.0285785185 days + 2019-04-02T12:52:50: Time range (UTC) ..........: 2020-01-01T00:10:00 - 2020-01-01T00:40:00 + 2019-04-02T12:52:50: Time range (MET) ..........: 631109469.184 - 631111269.184 seconds + 2019-04-02T12:52:50: Selected energy range .....: 0.1 - 100 TeV + 2019-04-02T12:52:50: Requested RoI .............: Centre(RA,DEC)=(83.63, 22.51) deg, Radius=3 deg + 2019-04-02T12:52:50: RoI of data ...............: Centre(RA,DEC)=(83.63, 22.51) deg, Radius=5 deg + 2019-04-02T12:52:50: Selected RoI ..............: Centre(RA,DEC)=(83.63, 22.51) deg, Radius=3 deg + 2019-04-02T12:52:50: cfitsio selection .........: TIME >= 631109469.18399978 && TIME <= 631111269.18400002 && ENERGY >= 0.10000000 && ENERGY <= 100.00000000 && ANGSEP(83.630000,22.510000,RA,DEC) <= 3.000000 + 2019-04-02T12:52:50: FITS filename .............: /var/tmp/tmp.0.vWKyrl[EVENTS][TIME >= 631109469.18399978 && TIME <= 631111269.18400002 && ENERGY >= 0.10000000 && ENERGY <= 100.00000000 && ANGSEP(83.630000,22.510000,RA,DEC) <= 3.000000] + 2019-04-02T12:52:51: + 2019-04-02T12:52:51: +====================+ + 2019-04-02T12:52:51: | Output observation | + 2019-04-02T12:52:51: +====================+ + 2019-04-02T12:52:51: === GObservations === + 2019-04-02T12:52:51: Number of observations ....: 1 + 2019-04-02T12:52:51: Number of models ..........: 0 + 2019-04-02T12:52:51: Number of observed events .: 22708 + 2019-04-02T12:52:51: Number of predicted events : 0 -From the 202210 events in the input event list 21974 events were selected. +From the 206048 events in the input event list 22708 events were selected. .. note:: Event selection by :ref:`ctselect` in a given dimension can be avoided @@ -82,8 +82,8 @@ From the 202210 events in the input event list 21974 events were selected. $ ctselect Input event list or observation definition XML file [events.fits] - Radius of ROI around pointing or specified RA/DEC (degrees) (0-180.) [3.] NONE - Start time (UTC string, JD, MJD or time in seconds) [2020-01-01T00:10:00] - Stop time (UTC string, JD, MJD or time in seconds) [2020-01-01T00:40:00] + Radius of ROI around pointing or specified RA/DEC (degrees) (0-180) [3.0] NONE + Start time (UTC string, JD, MJD or MET in seconds) [2020-01-01T00:10:00] + Stop time (UTC string, JD, MJD or MET in seconds) [2020-01-01T00:40:00] Lower energy limit (TeV) [0.1] NONE - Output event list or observation definition XML file [selected_events.fits] + Output event list or observation definition XML file [selected_events.fits] diff --git a/doc/source/users/tutorials/quickstart/simulating.rst b/doc/source/users/tutorials/quickstart/simulating.rst index 9675157..468ae03 100755 --- a/doc/source/users/tutorials/quickstart/simulating.rst +++ b/doc/source/users/tutorials/quickstart/simulating.rst @@ -207,33 +207,33 @@ from this file are shown here: .. code-block:: none - 2018-01-24T14:01:30: === CTA observation === - 2018-01-24T14:01:30: Simulation cone ...........: RA=83.63 deg, Dec=22.51 deg, radius=5.5 deg - 2018-01-24T14:01:30: Time interval .............: 6.31109e+08 - 6.31112e+08 s - 2018-01-24T14:01:30: Photon energy range .......: 30 GeV - 70.310187347763 GeV - 2018-01-24T14:01:30: Event energy range ........: 30 GeV - 70.310187347763 GeV - 2018-01-24T14:01:30: Simulation area ..........: 1.90195e+09 cm2 - 2018-01-24T14:01:30: Use model ................: Crab - 2018-01-24T14:01:30: Normalization ............: 1 [Crab] - 2018-01-24T14:01:30: Flux .....................: 2.50006e-09 [Crab] photons/cm2/s - 2018-01-24T14:01:30: Normalized flux ..........: 2.50006e-09 [Crab] photons/cm2/s - 2018-01-24T14:01:30: Photon rate ..............: 4.75499 photons/s [Crab] - 2018-01-24T14:01:30: MC source photons ........: 17169 [Crab] - 2018-01-24T14:01:30: MC source events .........: 3870 [Crab] - 2018-01-24T14:01:30: MC source events .........: 3870 (all source models) - 2018-01-24T14:01:30: Photon energy range .......: 70.310187347763 GeV - 164.784081495918 GeV + 2019-04-02T12:47:42: === CTA observation === + 2019-04-02T12:47:42: Simulation cone ...........: RA=83.63 deg, Dec=22.51 deg, radius=5.5 deg + 2019-04-02T12:47:42: Time interval .............: 6.31109e+08 - 6.31112e+08 s + 2019-04-02T12:47:42: Photon energy range .......: 30 GeV - 70.310187347763 GeV + 2019-04-02T12:47:42: Event energy range ........: 30 GeV - 70.310187347763 GeV + 2019-04-02T12:47:42: Simulation area ..........: 1.90195e+09 cm2 + 2019-04-02T12:47:42: Use model ................: Crab + 2019-04-02T12:47:42: Normalization ............: 1 [Crab] + 2019-04-02T12:47:42: Flux .....................: 2.50006e-09 [Crab] photons/cm2/s + 2019-04-02T12:47:42: Normalized flux ..........: 2.50006e-09 [Crab] photons/cm2/s + 2019-04-02T12:47:42: Photon rate ..............: 4.75499 photons/s [Crab] + 2019-04-02T12:47:42: MC source photons ........: 17169 [Crab] + 2019-04-02T12:47:42: MC source events .........: 3870 [Crab] + 2019-04-02T12:47:42: MC source events .........: 3870 (all source models) + 2019-04-02T12:47:42: Photon energy range .......: 70.310187347763 GeV - 164.784081495918 GeV ... - 2018-01-24T14:01:30: MC source photons .........: 46368 [Crab] - 2018-01-24T14:01:30: MC source events ..........: 12749 [Crab] - 2018-01-24T14:01:38: MC events outside ROI .....: 0 - 2018-01-24T14:01:38: MC background events ......: 189461 - 2018-01-24T14:01:38: MC identifier 1 ...........: Crab - 2018-01-24T14:01:38: MC identifier 2 ...........: CTABackgroundModel - 2018-01-24T14:01:38: MC events .................: 202210 (all models) + 2019-04-02T12:47:42: MC source photons .........: 46368 [Crab] + 2019-04-02T12:47:42: MC source events ..........: 12749 [Crab] + 2019-04-02T12:47:50: MC events outside ROI .....: 0 + 2019-04-02T12:47:50: MC background events ......: 193299 + 2019-04-02T12:47:50: MC identifier 1 ...........: Crab + 2019-04-02T12:47:50: MC identifier 2 ...........: CTABackgroundModel + 2019-04-02T12:47:50: MC events .................: 206048 (all models) Each line starts with the UTC time at which the line has been written. In this run, 46368 Crab photons have been thrown. 12749 of these photons have been -registered by CTA as events. In the same time interval, 189461 background +registered by CTA as events. In the same time interval, 193299 background events have been registred by CTA. .. note:: diff --git a/doc/source/users/tutorials/quickstart/skymap.jpg b/doc/source/users/tutorials/quickstart/skymap.jpg deleted file mode 100755 index 013cb4b..0000000 --- a/doc/source/users/tutorials/quickstart/skymap.jpg +++ /dev/null diff --git a/doc/source/users/tutorials/quickstart/skymap.png b/doc/source/users/tutorials/quickstart/skymap.png new file mode 100644 index 0000000..915c5b7 --- /dev/null +++ b/doc/source/users/tutorials/quickstart/skymap.png diff --git a/doc/source/users/tutorials/quickstart/skymap.rst b/doc/source/users/tutorials/quickstart/skymap.rst index fc8284b..6efaad1 100755 --- a/doc/source/users/tutorials/quickstart/skymap.rst +++ b/doc/source/users/tutorials/quickstart/skymap.rst @@ -15,10 +15,10 @@ this with the :ref:`ctskymap` tool as follows: $ ctskymap Input event list or observation definition XML file [events.fits] selected_events.fits + 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] Second coordinate of image center in degrees (DEC or galactic b) (-90-90) [22.01] - 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] Size of the X axis in pixels [200] Size of the Y axis in pixels [200] @@ -34,10 +34,10 @@ nebula (Right Ascension 83.63 deg, Declination 22.01 deg) and consists of an area of 4 deg x 4 deg. Below an image of the sky map, displayed using -`ds9 `_ in logarithmic color scale with a small Gaussian +`ds9 `_ in logarithmic color scale with a 3 pixel Gaussian kernel smoothing applied: -.. figure:: skymap.jpg +.. figure:: skymap.png :width: 400px :align: center @@ -57,26 +57,26 @@ example: .. code-block:: bash $ ctskymap - Input event list or observation definition XML file [selected_events.fits] - First coordinate of image center in degrees (RA or galactic l) (0-360) [83.63] - Second coordinate of image center in degrees (DEC or galactic b) (-90-90) [22.01] - 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] - Size of the X axis in pixels [200] - Size of the Y axis in pixels [200] - Lower energy limit (TeV) [0.1] - Upper energy limit (TeV) [100.0] - Background subtraction method (NONE|IRF|RING) [NONE] IRF - Calibration database [prod2] - Instrument response function [South_0.5h] - Output skymap file [skymap.fits] skymap_subtracted.fits + Input event list or observation definition XML file [selected_events.fits] + 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] + Second coordinate of image center in degrees (DEC or galactic b) (-90-90) [22.01] + Image scale (in degrees/pixel) [0.02] + Size of the X axis in pixels [200] + Size of the Y axis in pixels [200] + Lower energy limit (TeV) [0.1] + Upper energy limit (TeV) [100.0] + Background subtraction method (NONE|IRF|RING) [NONE] IRF + Calibration database [prod2] + Instrument response function [South_0.5h] + Output skymap file [skymap.fits] skymap_subtracted.fits The resulting sky map, displayed using -`ds9 `_ in logarithmic color scale with a small Gaussian +`ds9 `_ in logarithmic color scale with a 3 pixel Gaussian kernel smoothing applied is shown below: -.. figure:: skymap_subtracted.jpg +.. figure:: skymap_subtracted.png :width: 400px :align: center @@ -89,50 +89,41 @@ the background subtraction method set to ``RING``. For each position in the map the number of background counts is estimated from a ring, scaled according to the background model in the :ref:`instrument response functions `, and subtracted from a circular -region centred on the trial source region. You need to avoid regions -with significant gamma-ray emission when deriving the background -rates, by providing an exclusion region (`ds9 `_ -region file, or FITS map) through the hidden ``inexclusion`` parameter. +region centred on the trial source region. Regions with with significant +gamma-ray emission need to be excluded from the ring background estimate. +You can do this by iteratively determining an exclusion map that contains +all pixels above a given significance threshold. We chose here 3 iterations +and a threshold of 5 sigma. .. code-block:: bash - $ ctskymap inexclusion=exclusion.reg + $ ctskymap Input event list or observation definition XML file [selected_events.fits] + 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] Second coordinate of image center in degrees (DEC or galactic b) (-90-90) [22.01] - 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] Size of the X axis in pixels [200] Size of the Y axis in pixels [200] Lower energy limit (TeV) [0.1] Upper energy limit (TeV) [100.0] Background subtraction method (NONE|IRF|RING) [IRF] RING - Source region radius for estimating on-counts (degrees) [0.2] - Inner ring radius (degrees) [0.6] - Outer ring radius (degrees) [0.8] + Source region radius for estimating on-counts (degrees) [0.1] 0.05 + Inner background ring radius (degrees) [0.6] + Outer background ring radius (degrees) [0.8] + Number of iterations for exclusion regions computation (0-100) [0] 3 + Significance threshold for exclusion regions computation [5.0] Calibration database [prod2] Instrument response function [South_0.5h] Output skymap file [skymap_subtracted.fits] skymap_ring.fits -The exclusion file used contained just a circle centered on the Crab -nebula. - -.. code-block:: bash - - $ more exclusion.reg - fk5 - circle(83.63,22.01,0.2) - The resulting sky map, displayed using `ds9 `_ in logarithmic color scale is shown below. -.. figure:: skymap_ring.jpg +.. figure:: skymap_ring.png :width: 400px :align: center *Sky map of the selected events with background subtraction from the ring method* - -Note that the map thus obtained is smoothed to the scale of the source -region size. diff --git a/doc/source/users/tutorials/quickstart/skymap_ring.jpg b/doc/source/users/tutorials/quickstart/skymap_ring.jpg deleted file mode 100644 index 5a61f77..0000000 --- a/doc/source/users/tutorials/quickstart/skymap_ring.jpg +++ /dev/null diff --git a/doc/source/users/tutorials/quickstart/skymap_ring.png b/doc/source/users/tutorials/quickstart/skymap_ring.png new file mode 100644 index 0000000..fc5cda1 --- /dev/null +++ b/doc/source/users/tutorials/quickstart/skymap_ring.png diff --git a/doc/source/users/tutorials/quickstart/skymap_subtracted.jpg b/doc/source/users/tutorials/quickstart/skymap_subtracted.jpg deleted file mode 100755 index 50a395c..0000000 --- a/doc/source/users/tutorials/quickstart/skymap_subtracted.jpg +++ /dev/null diff --git a/doc/source/users/tutorials/quickstart/skymap_subtracted.png b/doc/source/users/tutorials/quickstart/skymap_subtracted.png new file mode 100644 index 0000000..734f276 --- /dev/null +++ b/doc/source/users/tutorials/quickstart/skymap_subtracted.png diff --git a/doc/source/users/tutorials/quickstart/spectrum.jpg b/doc/source/users/tutorials/quickstart/spectrum.jpg deleted file mode 100755 index 8b8b9c5..0000000 --- a/doc/source/users/tutorials/quickstart/spectrum.jpg +++ /dev/null diff --git a/doc/source/users/tutorials/quickstart/spectrum.png b/doc/source/users/tutorials/quickstart/spectrum.png new file mode 100644 index 0000000..ce0cc4d --- /dev/null +++ b/doc/source/users/tutorials/quickstart/spectrum.png diff --git a/doc/source/users/tutorials/quickstart/spectrum.rst b/doc/source/users/tutorials/quickstart/spectrum.rst index ddbfb89..4334aa7 100755 --- a/doc/source/users/tutorials/quickstart/spectrum.rst +++ b/doc/source/users/tutorials/quickstart/spectrum.rst @@ -53,7 +53,7 @@ The :ref:`csspec` script creates the FITS file ``spectrum.fits`` that contains the results of the SED fitting. A graphical display of the fit results is shown below. -.. figure:: spectrum.jpg +.. figure:: spectrum.png :width: 600px :align: center diff --git a/doc/source/users/tutorials/quickstart/unbinned.rst b/doc/source/users/tutorials/quickstart/unbinned.rst index 7128516..43af7a8 100755 --- a/doc/source/users/tutorials/quickstart/unbinned.rst +++ b/doc/source/users/tutorials/quickstart/unbinned.rst @@ -56,54 +56,61 @@ the same events for both analyses. .. code-block:: none - 2018-01-24T14:42:49: +=========================================+ - 2018-01-24T14:42:49: | Maximum likelihood optimisation results | - 2018-01-24T14:42:49: +=========================================+ - 2018-01-24T14:42:49: === GOptimizerLM === - 2018-01-24T14:42:49: Optimized function value ..: 139993.152 - 2018-01-24T14:42:49: Absolute precision ........: 0.005 - 2018-01-24T14:42:49: Acceptable value decrease .: 2 - 2018-01-24T14:42:49: Optimization status .......: converged - 2018-01-24T14:42:49: Number of parameters ......: 10 - 2018-01-24T14:42:49: Number of free parameters .: 4 - 2018-01-24T14:42:49: Number of iterations ......: 2 - 2018-01-24T14:42:49: Lambda ....................: 1e-05 - 2018-01-24T14:42:49: Maximum log likelihood ....: -139993.152 - 2018-01-24T14:42:49: Observed events (Nobs) ...: 21974.000 - 2018-01-24T14:42:49: Predicted events (Npred) ..: 21973.998 (Nobs - Npred = 0.00192004605924012) - 2018-01-24T14:42:49: === GModels === - 2018-01-24T14:42:49: Number of models ..........: 2 - 2018-01-24T14:42:49: Number of parameters ......: 10 - 2018-01-24T14:42:49: === GModelSky === - 2018-01-24T14:42:49: Name ......................: Crab - 2018-01-24T14:42:49: Instruments ...............: all - 2018-01-24T14:42:49: Instrument scale factors ..: unity - 2018-01-24T14:42:49: Observation identifiers ...: all - 2018-01-24T14:42:49: Model type ................: PointSource - 2018-01-24T14:42:49: Model components ..........: "PointSource" * "PowerLaw" * "Constant" - 2018-01-24T14:42:49: Number of parameters ......: 6 - 2018-01-24T14:42:49: Number of spatial par's ...: 2 - 2018-01-24T14:42:49: RA .......................: 83.6331 [-360,360] deg (fixed,scale=1) - 2018-01-24T14:42:49: DEC ......................: 22.0145 [-90,90] deg (fixed,scale=1) - 2018-01-24T14:42:49: Number of spectral par's ..: 3 - 2018-01-24T14:42:49: Prefactor ................: 5.84749807744844e-16 +/- 1.02060247557839e-17 [1e-23,1e-13] ph/cm2/s/MeV (free,scale=1e-16,gradient) - 2018-01-24T14:42:49: Index ....................: -2.48831941472059 +/- 0.0150002799717933 [-0,-5] (free,scale=-1,gradient) - 2018-01-24T14:42:49: PivotEnergy ..............: 300000 [10000,1000000000] MeV (fixed,scale=1000000,gradient) - 2018-01-24T14:42:49: Number of temporal par's ..: 1 - 2018-01-24T14:42:49: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) - 2018-01-24T14:42:49: === GCTAModelIrfBackground === - 2018-01-24T14:42:49: Name ......................: CTABackgroundModel - 2018-01-24T14:42:49: Instruments ...............: CTA - 2018-01-24T14:42:49: Instrument scale factors ..: unity - 2018-01-24T14:42:49: Observation identifiers ...: all - 2018-01-24T14:42:49: Model type ................: "PowerLaw" * "Constant" - 2018-01-24T14:42:49: Number of parameters ......: 4 - 2018-01-24T14:42:49: Number of spectral par's ..: 3 - 2018-01-24T14:42:49: Prefactor ................: 1.00311308345777 +/- 0.0134769768119516 [0.001,1000] ph/cm2/s/MeV (free,scale=1,gradient) - 2018-01-24T14:42:49: Index ....................: 0.00710999567642649 +/- 0.00818686036300859 [-5,5] (free,scale=1,gradient) - 2018-01-24T14:42:49: PivotEnergy ..............: 1000000 [10000,1000000000] MeV (fixed,scale=1000000,gradient) - 2018-01-24T14:42:49: Number of temporal par's ..: 1 - 2018-01-24T14:42:49: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) + 2019-04-02T14:15:10: +=================================+ + 2019-04-02T14:15:10: | Maximum likelihood optimisation | + 2019-04-02T14:15:10: +=================================+ + 2019-04-02T14:15:10: >Iteration 0: -logL=143782.837, Lambda=1.0e-03 + 2019-04-02T14:15:10: >Iteration 1: -logL=143779.346, Lambda=1.0e-03, delta=3.491, step=1.0e+00, max(|grad|)=5.346881 [Index:7] + 2019-04-02T14:15:10: >Iteration 2: -logL=143779.343, Lambda=1.0e-04, delta=0.003, step=1.0e+00, max(|grad|)=-0.055672 [Index:3] + 2019-04-02T14:15:10: + 2019-04-02T14:15:10: +=========================================+ + 2019-04-02T14:15:10: | Maximum likelihood optimisation results | + 2019-04-02T14:15:10: +=========================================+ + 2019-04-02T14:15:10: === GOptimizerLM === + 2019-04-02T14:15:10: Optimized function value ..: 143779.343 + 2019-04-02T14:15:10: Absolute precision ........: 0.005 + 2019-04-02T14:15:10: Acceptable value decrease .: 2 + 2019-04-02T14:15:10: Optimization status .......: converged + 2019-04-02T14:15:10: Number of parameters ......: 10 + 2019-04-02T14:15:10: Number of free parameters .: 4 + 2019-04-02T14:15:10: Number of iterations ......: 2 + 2019-04-02T14:15:10: Lambda ....................: 1e-05 + 2019-04-02T14:15:10: Maximum log likelihood ....: -143779.343 + 2019-04-02T14:15:10: Observed events (Nobs) ...: 22708.000 + 2019-04-02T14:15:10: Predicted events (Npred) ..: 22707.995 (Nobs - Npred = 0.00519125738355797) + 2019-04-02T14:15:10: === GModels === + 2019-04-02T14:15:10: Number of models ..........: 2 + 2019-04-02T14:15:10: Number of parameters ......: 10 + 2019-04-02T14:15:10: === GModelSky === + 2019-04-02T14:15:10: Name ......................: Crab + 2019-04-02T14:15:10: Instruments ...............: all + 2019-04-02T14:15:10: Instrument scale factors ..: unity + 2019-04-02T14:15:10: Observation identifiers ...: all + 2019-04-02T14:15:10: Model type ................: PointSource + 2019-04-02T14:15:10: Model components ..........: "PointSource" * "PowerLaw" * "Constant" + 2019-04-02T14:15:10: Number of parameters ......: 6 + 2019-04-02T14:15:10: Number of spatial par's ...: 2 + 2019-04-02T14:15:10: RA .......................: 83.6331 [-360,360] deg (fixed,scale=1) + 2019-04-02T14:15:10: DEC ......................: 22.0145 [-90,90] deg (fixed,scale=1) + 2019-04-02T14:15:10: Number of spectral par's ..: 3 + 2019-04-02T14:15:10: Prefactor ................: 5.88338676901458e-16 +/- 1.02452856089807e-17 [1e-23,1e-13] ph/cm2/s/MeV (free,scale=1e-16,gradient) + 2019-04-02T14:15:10: Index ....................: -2.49375950219757 +/- 0.0149889370322137 [-0,-5] (free,scale=-1,gradient) + 2019-04-02T14:15:10: PivotEnergy ..............: 300000 [10000,1000000000] MeV (fixed,scale=1000000,gradient) + 2019-04-02T14:15:10: Number of temporal par's ..: 1 + 2019-04-02T14:15:10: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) + 2019-04-02T14:15:10: === GCTAModelIrfBackground === + 2019-04-02T14:15:10: Name ......................: CTABackgroundModel + 2019-04-02T14:15:10: Instruments ...............: CTA + 2019-04-02T14:15:10: Instrument scale factors ..: unity + 2019-04-02T14:15:10: Observation identifiers ...: all + 2019-04-02T14:15:10: Model type ................: "PowerLaw" * "Constant" + 2019-04-02T14:15:10: Number of parameters ......: 4 + 2019-04-02T14:15:10: Number of spectral par's ..: 3 + 2019-04-02T14:15:10: Prefactor ................: 1.0018169793538 +/- 0.0133053833141539 [0.001,1000] ph/cm2/s/MeV (free,scale=1,gradient) + 2019-04-02T14:15:10: Index ....................: -0.00708154249642314 +/- 0.00805278228449961 [-5,5] (free,scale=1,gradient) + 2019-04-02T14:15:10: PivotEnergy ..............: 1000000 [10000,1000000000] MeV (fixed,scale=1000000,gradient) + 2019-04-02T14:15:10: Number of temporal par's ..: 1 + 2019-04-02T14:15:10: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) .. note:: Many tools or scripts can also be used in unbinned mode, including