{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "\n",
    "<div class=\"alert alert-info\">\n",
    "\n",
    "**This is a fixed-text formatted version of a Jupyter notebook**\n",
    "\n",
    "- Try online[![Binder](https://static.mybinder.org/badge.svg)](https://mybinder.org/v2/gh/gammapy/gammapy-webpage/v0.19?urlpath=lab/tree/tutorials/api/catalog.ipynb)\n",
    "- You may download all the notebooks in the documentation as a\n",
    "[tar file](../../_downloads/notebooks-0.19.tar).\n",
    "- **Source files:**\n",
    "[catalog.ipynb](../../_static/notebooks/catalog.ipynb) |\n",
    "[catalog.py](../../_static/notebooks/catalog.py)\n",
    "</div>\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Source catalogs\n",
    "\n",
    "`~gammapy.catalog` provides convenient access to common gamma-ray source catalogs.\n",
    "This module is mostly independent from the rest of Gammapy.\n",
    "Typically you use it to compare new analyses against catalog results, e.g. overplot the spectral model, or compare the source position.\n",
    "\n",
    "Moreover as creating a source model and flux points for a given catalog from the FITS table is tedious, `~gammapy.catalog` has this already implemented. So you can create initial source models for your analyses.\n",
    "This is very common for Fermi-LAT, to start with a catalog model.\n",
    "For TeV analysis, especially in crowded Galactic regions, using the HGPS, gamma-cat or 2HWC catalog in this way can also be useful.\n",
    "\n",
    "In this tutorial you will learn how to:\n",
    "\n",
    "- List available catalogs\n",
    "- Load a catalog\n",
    "- Access the source catalog table data\n",
    "- Select a catalog subset or a single source\n",
    "- Get source spectral and spatial models\n",
    "- Get flux points (if available)\n",
    "- Get lightcurves (if available)\n",
    "- Access the source catalog table data\n",
    "- Pretty-print the source information\n",
    "\n",
    "In this tutorial we will show examples using the following catalogs:\n",
    "\n",
    "- `~gammapy.catalog.SourceCatalogHGPS`\n",
    "- `~gammapy.catalog.SourceCatalogGammaCat`\n",
    "- `~gammapy.catalog.SourceCatalog3FHL`\n",
    "- `~gammapy.catalog.SourceCatalog4FGL`\n",
    "\n",
    "All catalog and source classes work the same, as long as some information is available. E.g. trying to access a lightcurve from a catalog and source that doesn't have that information will return ``None``.\n",
    "\n",
    "Further information is available at `~gammapy.catalog`."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:34.423421Z",
     "iopub.status.busy": "2021-11-22T21:09:34.422535Z",
     "iopub.status.idle": "2021-11-22T21:09:35.105228Z",
     "shell.execute_reply": "2021-11-22T21:09:35.105399Z"
    }
   },
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "import astropy.units as u\n",
    "from gammapy.catalog import CATALOG_REGISTRY"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## List available catalogs\n",
    "\n",
    "`~gammapy.catalog` contains a catalog registry ``CATALOG_REGISTRY``, which maps catalog names (e.g. \"3fhl\") to catalog classes (e.g. ``SourceCatalog3FHL``). "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:35.108216Z",
     "iopub.status.busy": "2021-11-22T21:09:35.107906Z",
     "iopub.status.idle": "2021-11-22T21:09:35.109869Z",
     "shell.execute_reply": "2021-11-22T21:09:35.110050Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[gammapy.catalog.gammacat.SourceCatalogGammaCat,\n",
       " gammapy.catalog.hess.SourceCatalogHGPS,\n",
       " gammapy.catalog.hawc.SourceCatalog2HWC,\n",
       " gammapy.catalog.fermi.SourceCatalog3FGL,\n",
       " gammapy.catalog.fermi.SourceCatalog4FGL,\n",
       " gammapy.catalog.fermi.SourceCatalog2FHL,\n",
       " gammapy.catalog.fermi.SourceCatalog3FHL,\n",
       " gammapy.catalog.hawc.SourceCatalog3HWC]"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "CATALOG_REGISTRY"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Load catalogs\n",
    "\n",
    "If you have run `gammapy download datasets` or `gammapy download tutorials`,\n",
    "you have a copy of the catalogs as FITS files in `$GAMMAPY_DATA/catalogs`,\n",
    "and that is the default location where `~gammapy.catalog` loads from.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:35.112051Z",
     "iopub.status.busy": "2021-11-22T21:09:35.111745Z",
     "iopub.status.idle": "2021-11-22T21:09:35.231747Z",
     "shell.execute_reply": "2021-11-22T21:09:35.232078Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "2HWC.ecsv\r\n",
      "2HWC.yaml\r\n",
      "3HWC.ecsv\r\n",
      "3HWC.yaml\r\n",
      "README.rst\r\n",
      "\u001b[1m\u001b[34mfermi\u001b[m\u001b[m\r\n",
      "\u001b[1m\u001b[34mgammacat\u001b[m\u001b[m\r\n",
      "hgps_catalog_v1.fits.gz\r\n",
      "make_2hwc.py\r\n",
      "make_3hwc.py\r\n"
     ]
    }
   ],
   "source": [
    "!ls -1 $GAMMAPY_DATA/catalogs"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:35.236270Z",
     "iopub.status.busy": "2021-11-22T21:09:35.235741Z",
     "iopub.status.idle": "2021-11-22T21:09:35.361370Z",
     "shell.execute_reply": "2021-11-22T21:09:35.362197Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\u001b[1m\u001b[34mExtended_archive_v15\u001b[m\u001b[m\r\n",
      "\u001b[1m\u001b[34mExtended_archive_v18\u001b[m\u001b[m\r\n",
      "\u001b[1m\u001b[34mLAT_extended_sources_8years\u001b[m\u001b[m\r\n",
      "README.rst\r\n",
      "gll_psc_v16.fit.gz\r\n",
      "gll_psc_v20.fit.gz\r\n",
      "gll_psc_v27.fit.gz\r\n",
      "gll_psch_v08.fit.gz\r\n",
      "gll_psch_v09.fit.gz\r\n",
      "gll_psch_v13.fit.gz\r\n"
     ]
    }
   ],
   "source": [
    "!ls -1 $GAMMAPY_DATA/catalogs/fermi"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "So a catalog can be loaded directly from its corresponding class"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:35.372496Z",
     "iopub.status.busy": "2021-11-22T21:09:35.371714Z",
     "iopub.status.idle": "2021-11-22T21:09:35.634077Z",
     "shell.execute_reply": "2021-11-22T21:09:35.634316Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Number of sources : 5788\n"
     ]
    }
   ],
   "source": [
    "from gammapy.catalog import SourceCatalog4FGL\n",
    "\n",
    "catalog = SourceCatalog4FGL()\n",
    "print(\"Number of sources :\", len(catalog.table))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Note that it loads the default catalog from `$GAMMAPY_DATA/catalogs`, you could pass a different `filename` when creating the catalog.\n",
    "For example here we load an older version of 4FGL catalog:  "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:35.636285Z",
     "iopub.status.busy": "2021-11-22T21:09:35.636000Z",
     "iopub.status.idle": "2021-11-22T21:09:36.026181Z",
     "shell.execute_reply": "2021-11-22T21:09:36.026375Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Number of sources : 5066\n"
     ]
    }
   ],
   "source": [
    "catalog = SourceCatalog4FGL(\"$GAMMAPY_DATA/catalogs/fermi/gll_psc_v20.fit.gz\")\n",
    "print(\"Number of sources :\", len(catalog.table))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "\n",
    "Alternatively you can load a catalog by name via `CATALOG_REGISTRY.get_cls(name)()` (note the `()` to instantiate a catalog object from the catalog class - only this will load the catalog and be useful), or by importing the catalog class (e.g. `SourceCatalog3FGL`) directly. The two ways are equivalent, the result will be the same.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:36.028254Z",
     "iopub.status.busy": "2021-11-22T21:09:36.027950Z",
     "iopub.status.idle": "2021-11-22T21:09:36.185495Z",
     "shell.execute_reply": "2021-11-22T21:09:36.185684Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<gammapy.catalog.fermi.SourceCatalog3FGL at 0x1283dc070>"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# FITS file is loaded\n",
    "catalog = CATALOG_REGISTRY.get_cls(\"3fgl\")()\n",
    "catalog"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:36.187617Z",
     "iopub.status.busy": "2021-11-22T21:09:36.187336Z",
     "iopub.status.idle": "2021-11-22T21:09:36.583191Z",
     "shell.execute_reply": "2021-11-22T21:09:36.583361Z"
    }
   },
   "outputs": [],
   "source": [
    "# Let's load the source catalogs we will use throughout this tutorial\n",
    "catalog_gammacat = CATALOG_REGISTRY.get_cls(\"gamma-cat\")()\n",
    "catalog_3fhl = CATALOG_REGISTRY.get_cls(\"3fhl\")()\n",
    "catalog_4fgl = CATALOG_REGISTRY.get_cls(\"4fgl\")()\n",
    "catalog_hgps = CATALOG_REGISTRY.get_cls(\"hgps\")()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Catalog table\n",
    "\n",
    "Source catalogs are given as `FITS` files that contain one or multiple tables.\n",
    "\n",
    "However, you can also access the underlying `astropy.table.Table` for a catalog,\n",
    "and the row data as a Python `dict`. This can be useful if you want to do something\n",
    "that is not pre-scripted by the `~gammapy.catalog` classes, such as e.g. selecting\n",
    "sources by sky position or association class, or accessing special source information.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:36.585491Z",
     "iopub.status.busy": "2021-11-22T21:09:36.585167Z",
     "iopub.status.idle": "2021-11-22T21:09:36.586605Z",
     "shell.execute_reply": "2021-11-22T21:09:36.586759Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "astropy.table.table.Table"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "type(catalog_3fhl.table)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:36.588613Z",
     "iopub.status.busy": "2021-11-22T21:09:36.588288Z",
     "iopub.status.idle": "2021-11-22T21:09:36.589774Z",
     "shell.execute_reply": "2021-11-22T21:09:36.589932Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "1556"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "len(catalog_3fhl.table)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:36.592260Z",
     "iopub.status.busy": "2021-11-22T21:09:36.591960Z",
     "iopub.status.idle": "2021-11-22T21:09:36.594067Z",
     "shell.execute_reply": "2021-11-22T21:09:36.594298Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div><i>Table length=3</i>\n",
       "<table id=\"table4970136528\" class=\"table-striped table-bordered table-condensed\">\n",
       "<thead><tr><th>Source_Name</th><th>RAJ2000</th><th>DEJ2000</th></tr></thead>\n",
       "<thead><tr><th></th><th>deg</th><th>deg</th></tr></thead>\n",
       "<thead><tr><th>bytes18</th><th>float32</th><th>float32</th></tr></thead>\n",
       "<tr><td>3FHL J0001.2-0748</td><td>0.3107</td><td>-7.8075</td></tr>\n",
       "<tr><td>3FHL J0001.9-4155</td><td>0.4849</td><td>-41.9303</td></tr>\n",
       "<tr><td>3FHL J0002.1-6728</td><td>0.5283</td><td>-67.4825</td></tr>\n",
       "</table></div>"
      ],
      "text/plain": [
       "<Table length=3>\n",
       "   Source_Name     RAJ2000  DEJ2000 \n",
       "                     deg      deg   \n",
       "     bytes18       float32  float32 \n",
       "------------------ -------- --------\n",
       "3FHL J0001.2-0748    0.3107  -7.8075\n",
       "3FHL J0001.9-4155    0.4849 -41.9303\n",
       "3FHL J0002.1-6728    0.5283 -67.4825"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "catalog_3fhl.table[:3][[\"Source_Name\", \"RAJ2000\", \"DEJ2000\"]]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Note that the catalogs object include a helper property that gives directly the sources positions as a `SkyCoord` object (we will show an usage example in the following)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:36.596479Z",
     "iopub.status.busy": "2021-11-22T21:09:36.596123Z",
     "iopub.status.idle": "2021-11-22T21:09:36.598706Z",
     "shell.execute_reply": "2021-11-22T21:09:36.598925Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<SkyCoord (ICRS): (ra, dec) in deg\n",
       "    [(0.31067517,  -7.8075185), (0.4848653 , -41.93026  ),\n",
       "     (0.52826166, -67.48248  )]>"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "catalog_3fhl.positions[:3]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Source object\n",
    "\n",
    "### Select a source\n",
    "\n",
    "The catalog entries for a single source are represented by a `SourceCatalogObject`.\n",
    "In order to select a source object index into the catalog using `[]`, with a catalog table row index (zero-based, first row is `[0]`), or a source name. If a name is given, catalog table columns with source names and association names (\"ASSOC1\" in the example below) are searched top to bottom. There is no name resolution web query.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:36.601726Z",
     "iopub.status.busy": "2021-11-22T21:09:36.601411Z",
     "iopub.status.idle": "2021-11-22T21:09:36.602859Z",
     "shell.execute_reply": "2021-11-22T21:09:36.603090Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<gammapy.catalog.fermi.SourceCatalogObject4FGL at 0x1283e46a0>"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "source = catalog_4fgl[49]\n",
    "source"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:36.604871Z",
     "iopub.status.busy": "2021-11-22T21:09:36.604582Z",
     "iopub.status.idle": "2021-11-22T21:09:36.605931Z",
     "shell.execute_reply": "2021-11-22T21:09:36.606221Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(49, '4FGL J0010.8-2154')"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "source.row_index, source.name"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:36.688774Z",
     "iopub.status.busy": "2021-11-22T21:09:36.663590Z",
     "iopub.status.idle": "2021-11-22T21:09:36.690172Z",
     "shell.execute_reply": "2021-11-22T21:09:36.690414Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(49, '4FGL J0010.8-2154')"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "source = catalog_4fgl[\"4FGL J0010.8-2154\"]\n",
    "source.row_index, source.name"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:36.692225Z",
     "iopub.status.busy": "2021-11-22T21:09:36.691902Z",
     "iopub.status.idle": "2021-11-22T21:09:36.693321Z",
     "shell.execute_reply": "2021-11-22T21:09:36.693503Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'PKS 0008-222                '"
      ]
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "source.data[\"ASSOC1\"]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:36.700914Z",
     "iopub.status.busy": "2021-11-22T21:09:36.700618Z",
     "iopub.status.idle": "2021-11-22T21:09:36.737010Z",
     "shell.execute_reply": "2021-11-22T21:09:36.737234Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(49, '4FGL J0010.8-2154')"
      ]
     },
     "execution_count": 17,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "source = catalog_4fgl[\"PKS 0008-222\"]\n",
    "source.row_index, source.name"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Note that you can also do a `for source in catalog` loop, to find or process\n",
    "sources of interest.\n",
    "\n",
    "###  Source information\n",
    "\n",
    "The source objects have a `data` property that contains the information of the catalog row corresponding to the source."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:36.739172Z",
     "iopub.status.busy": "2021-11-22T21:09:36.738844Z",
     "iopub.status.idle": "2021-11-22T21:09:36.740217Z",
     "shell.execute_reply": "2021-11-22T21:09:36.740462Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "276.04663"
      ]
     },
     "execution_count": 18,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "source.data[\"Npred\"]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:36.742516Z",
     "iopub.status.busy": "2021-11-22T21:09:36.742177Z",
     "iopub.status.idle": "2021-11-22T21:09:36.743704Z",
     "shell.execute_reply": "2021-11-22T21:09:36.743879Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(<Quantity 60.28118 deg>, <Quantity -79.400505 deg>)"
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "source.data[\"GLON\"], source.data[\"GLAT\"]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "As for the catalog object, the source object has a `position` property."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:36.752480Z",
     "iopub.status.busy": "2021-11-22T21:09:36.752174Z",
     "iopub.status.idle": "2021-11-22T21:09:36.753483Z",
     "shell.execute_reply": "2021-11-22T21:09:36.753705Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<SkyCoord (Galactic): (l, b) in deg\n",
       "    (60.28120079, -79.40051035)>"
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "source.position.galactic"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Select a catalog subset\n",
    "\n",
    "The catalog objects support selection using boolean arrays (of the same length), so one can create a new catalog as a subset of the main catalog that verify a set of conditions.\n",
    "\n",
    "In the next example we selection only few of the brightest sources brightest sources in the 100 to 200 GeV energy band."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:37.156805Z",
     "iopub.status.busy": "2021-11-22T21:09:37.055322Z",
     "iopub.status.idle": "2021-11-22T21:09:37.680753Z",
     "shell.execute_reply": "2021-11-22T21:09:37.680937Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "352     3FHL J0534.5+2201    2.99e-10 1 / (cm2 s)\n",
      "553     3FHL J0851.9-4620e   1.24e-10 1 / (cm2 s)\n",
      "654     3FHL J1036.3-5833e   1.57e-10 1 / (cm2 s)\n",
      "691     3FHL J1104.4+3812    3.34e-10 1 / (cm2 s)\n",
      "1111    3FHL J1653.8+3945    1.27e-10 1 / (cm2 s)\n",
      "1219    3FHL J1824.5-1351e   1.77e-10 1 / (cm2 s)\n",
      "1361    3FHL J2028.6+4110e   1.75e-10 1 / (cm2 s)\n"
     ]
    }
   ],
   "source": [
    "mask_bright = np.zeros(len(catalog_3fhl.table), dtype=bool)\n",
    "for k, source in enumerate(catalog_3fhl):\n",
    "    flux = (\n",
    "        source.spectral_model()\n",
    "        .integral(100 * u.GeV, 200 * u.GeV)\n",
    "        .to(\"cm-2 s-1\")\n",
    "    )\n",
    "    if flux > 1e-10 * u.Unit(\"cm-2 s-1\"):\n",
    "        mask_bright[k] = True\n",
    "        print(f\"{source.row_index:<7d} {source.name:20s} {flux:.3g}\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:37.683453Z",
     "iopub.status.busy": "2021-11-22T21:09:37.683143Z",
     "iopub.status.idle": "2021-11-22T21:09:37.687401Z",
     "shell.execute_reply": "2021-11-22T21:09:37.687560Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<gammapy.catalog.fermi.SourceCatalog3FHL at 0x1283ee460>"
      ]
     },
     "execution_count": 22,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "catalog_3fhl_bright = catalog_3fhl[mask_bright]\n",
    "catalog_3fhl_bright"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:37.689630Z",
     "iopub.status.busy": "2021-11-22T21:09:37.689308Z",
     "iopub.status.idle": "2021-11-22T21:09:37.690683Z",
     "shell.execute_reply": "2021-11-22T21:09:37.690906Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "&lt;Column name=&apos;Source_Name&apos; dtype=&apos;bytes18&apos; length=7&gt;\n",
       "<table>\n",
       "<tr><td>3FHL J0534.5+2201</td></tr>\n",
       "<tr><td>3FHL J0851.9-4620e</td></tr>\n",
       "<tr><td>3FHL J1036.3-5833e</td></tr>\n",
       "<tr><td>3FHL J1104.4+3812</td></tr>\n",
       "<tr><td>3FHL J1653.8+3945</td></tr>\n",
       "<tr><td>3FHL J1824.5-1351e</td></tr>\n",
       "<tr><td>3FHL J2028.6+4110e</td></tr>\n",
       "</table>"
      ],
      "text/plain": [
       "<Column name='Source_Name' dtype='bytes18' length=7>\n",
       "3FHL J0534.5+2201 \n",
       "3FHL J0851.9-4620e\n",
       "3FHL J1036.3-5833e\n",
       "3FHL J1104.4+3812 \n",
       "3FHL J1653.8+3945 \n",
       "3FHL J1824.5-1351e\n",
       "3FHL J2028.6+4110e"
      ]
     },
     "execution_count": 23,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "catalog_3fhl_bright.table[\"Source_Name\"]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Similarly we can select only sources within a region of interest. Here for example we use the `position` property of the catalog object to select sources whitin 5 degrees from \"PKS 0008-222\":\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:37.717173Z",
     "iopub.status.busy": "2021-11-22T21:09:37.716875Z",
     "iopub.status.idle": "2021-11-22T21:09:37.740577Z",
     "shell.execute_reply": "2021-11-22T21:09:37.740747Z"
    }
   },
   "outputs": [],
   "source": [
    "source = catalog_4fgl[\"PKS 0008-222\"]\n",
    "mask_roi = source.position.separation(catalog_4fgl.positions) < 5 * u.deg"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:37.743031Z",
     "iopub.status.busy": "2021-11-22T21:09:37.742736Z",
     "iopub.status.idle": "2021-11-22T21:09:37.758566Z",
     "shell.execute_reply": "2021-11-22T21:09:37.758785Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Number of sources : 11\n"
     ]
    }
   ],
   "source": [
    "catalog_4fgl_roi = catalog_4fgl[mask_roi]\n",
    "print(\"Number of sources :\", len(catalog_4fgl_roi.table))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Source models\n",
    "\n",
    "The `~gammapy.catalog.SourceCatalogObject` classes have a `sky_model()` model\n",
    "which creates a `gammapy.modeling.models.SkyModel` object, with model parameter\n",
    "values and parameter errors from the catalog filled in.\n",
    "\n",
    "In most cases, the `spectral_model()` method provides the `gammapy.modeling.models.SpectralModel`\n",
    "part of the sky model, and the `spatial_model()` method the `gammapy.modeling.models.SpatialModel`\n",
    "part individually.\n",
    "\n",
    "We use the `gammapy.catalog.SourceCatalog3FHL` for the examples in this section."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:37.801447Z",
     "iopub.status.busy": "2021-11-22T21:09:37.801160Z",
     "iopub.status.idle": "2021-11-22T21:09:37.802289Z",
     "shell.execute_reply": "2021-11-22T21:09:37.802460Z"
    }
   },
   "outputs": [],
   "source": [
    "source = catalog_4fgl[\"PKS 2155-304\"]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:37.806359Z",
     "iopub.status.busy": "2021-11-22T21:09:37.806067Z",
     "iopub.status.idle": "2021-11-22T21:09:37.815412Z",
     "shell.execute_reply": "2021-11-22T21:09:37.815590Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "SkyModel(spatial_model=<gammapy.modeling.models.spatial.PointSpatialModel object at 0x128724430>, spectral_model=<gammapy.modeling.models.spectral.LogParabolaSpectralModel object at 0x128724940>)temporal_model=None)"
      ]
     },
     "execution_count": 27,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "model = source.sky_model()\n",
    "model"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:37.817399Z",
     "iopub.status.busy": "2021-11-22T21:09:37.817113Z",
     "iopub.status.idle": "2021-11-22T21:09:37.818567Z",
     "shell.execute_reply": "2021-11-22T21:09:37.818339Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "SkyModel\n",
      "\n",
      "  Name                      : 4FGL J2158.8-3013\n",
      "  Datasets names            : None\n",
      "  Spectral model type       : LogParabolaSpectralModel\n",
      "  Spatial  model type       : PointSpatialModel\n",
      "  Temporal model type       : \n",
      "  Parameters:\n",
      "    amplitude               :   1.44e-11   +/- 1.7e-13 1 / (cm2 MeV s)\n",
      "    reference    (frozen)   :   1118.643       MeV         \n",
      "    alpha                   :      1.767   +/-    0.01             \n",
      "    beta                    :      0.041   +/-    0.00             \n",
      "    lon_0                   :    329.714   +/-    0.00 deg         \n",
      "    lat_0                   :    -30.225   +/-    0.00 deg         \n",
      "\n",
      "\n"
     ]
    }
   ],
   "source": [
    "print(model)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:37.820899Z",
     "iopub.status.busy": "2021-11-22T21:09:37.820591Z",
     "iopub.status.idle": "2021-11-22T21:09:37.821863Z",
     "shell.execute_reply": "2021-11-22T21:09:37.822011Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "PointSpatialModel\n",
      "\n",
      "  type   name    value    unit   error      min        max    frozen link\n",
      "------- ----- ----------- ---- --------- ---------- --------- ------ ----\n",
      "spatial lon_0  3.2971e+02  deg 3.735e-03        nan       nan  False     \n",
      "spatial lat_0 -3.0225e+01  deg 3.227e-03 -9.000e+01 9.000e+01  False     \n"
     ]
    }
   ],
   "source": [
    "print(model.spatial_model)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:37.824330Z",
     "iopub.status.busy": "2021-11-22T21:09:37.824030Z",
     "iopub.status.idle": "2021-11-22T21:09:37.825332Z",
     "shell.execute_reply": "2021-11-22T21:09:37.825506Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "LogParabolaSpectralModel\n",
      "\n",
      "  type      name     value         unit        error   min max frozen link\n",
      "-------- --------- ---------- -------------- --------- --- --- ------ ----\n",
      "spectral amplitude 1.4365e-11 cm-2 MeV-1 s-1 1.655e-13 nan nan  False     \n",
      "spectral reference 1.1186e+03            MeV 0.000e+00 nan nan   True     \n",
      "spectral     alpha 1.7668e+00                1.163e-02 nan nan  False     \n",
      "spectral      beta 4.0969e-02                4.120e-03 nan nan  False     \n"
     ]
    }
   ],
   "source": [
    "print(model.spectral_model)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:37.878288Z",
     "iopub.status.busy": "2021-11-22T21:09:37.840360Z",
     "iopub.status.idle": "2021-11-22T21:09:38.135879Z",
     "shell.execute_reply": "2021-11-22T21:09:38.136112Z"
    },
    "nbsphinx-thumbnail": {
     "tooltip": "Access and explore thew most common gamma-ray source catalogs."
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAaAAAAEKCAYAAABUsYHRAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAyf0lEQVR4nO3dd3xcV5n/8c8zvahXS1a15RY73ekhToEQsqEsvS6QUBZCSWCBQFjKQgiwu6ETCAuEpSSUDQvhF5KQkLLpsRMnseOmLlm9j0bTNHN+f8xIcRxbGtuauSPpeb9e85Lmzsydr3Mz8+ice+45YoxBKaWUyjab1QGUUkotT1qAlFJKWUILkFJKKUtoAVJKKWUJLUBKKaUsoQVIKaWUJRxWB1gsysrKTENDg9UxlFJqUdm2bduQMab8UI9pAUpTQ0MDW7dutTqGUkotKiLScbjHtAtOKaWUJbQAKaWUsoQWIKWUUpbQAqSUUsoSWoCUUkpZQguQUkopS2gBUkopZQm9DkgptewZY5ienmZ6epp4PD57SyQSszdjDIlEYvb5BxIRRASbzTb70263z/50OByzP202/bt/hhYgpdSSlkgkiEQiRKPRF91isdjsLR6PZy2PzWbD6XTicDhwuVw4nU5cLhculwu3243L5cLhWB5fzcvjX6mUWvJisRihUIhwOPyiWywWszrai8wUxEgkQjAYPORz7HY7Ho8Ht9uN1+vF4/Hg9Xpxu91ZTptZWoCUUovOzJf31NQUU1NThEIhpqenrY61YOLxOMFg8CUFymaz4fV68Xq9+Hw+/H4/Xq8XEbEo6bHRAqSUymkzX8aTk5OzX8rZ7DLLJYlE4iWFSURmi1FeXh55eXk4nU4LU6ZPC5BSKqfE43ECgQCBQIDJyUlCodBLTvqrFxhjZovSwMAAAG63m7y8PPLz88nPz8flclmc8tC0ACmlLDXzBToxMcHExARTU1NacI7RzDmm4eFhIFmQCgoKKCgoID8/H7vdbnHCJC1ASqmsm56eZnx8nPHxcSYmJpZtl1q2RCIRBgcHGRwcRETIy8ujsLCQoqIiSwc2aAFSSmVFJBJhbGyMsbExJicnF2SfxhgCUcNIKM5wKMFYOM54JMF4OMFkNEEwZpiMJghPm9lbLGGYjsN0whBPNbQMIIDdBnYRnDZw2QWXQ/DYBZ9T8Lts+J1CvttGQepW7LFT4rVR4rVT4JJFMRjAGDPbxdnd3Y3H46GoqIji4mJ8Pl9Ws2gBUkplTCQSYXR0lNHRUaampo5qH3FjGAjG6QnE6Q1M0zsZpz84TX8wzlAwTjTx0td47EKeW/A7beS5hBKvHY9D8DgEl11w2MBhE+xy8HtBPGGIJSAaN0RSRWsqlmAwGKctliAQMUTiL+0idNmhzGenwmdnRZ6DFXl2qvLsrCxwUOG3Y8/R4hQOh+nr66Ovrw+Xy0VxcTHFxcX4/f6Mv7cWIKXUgorFYoyOjjIyMnLY61wOJxBN0D4Wo210mvbxGJ3j0+yfmH5RkfE7hUq/nboCB5ur3JT67JR67ZR6bRR5ki0TjyOzsw1Epg3jkTij4QSjoQTDoThDU3EGp+IMBOPs7QwxFXuhSDltsDLfQV2hg8YiJ/VFDlYVOcl359asCNFolP7+fvr7+3G73ZSUlFBSUoLH48nI+2kBUkods0QiwdjYGMPDwwQCgbQGEcQShrbRGHuGY+wbidE8EqM/+MK5oBKPjbpCB5uafNQWOFiZ76A630F+DnR1uR1ChcNBxWEaCTNdgz2BafYHpumemKZrYpodA1Ee7AzPPq/Sb2d1sZO1pU7WlTppLHLiPLhZZpFIJEJvby+BQIB169Zl5D20ACmljtrU1BRDQ0OMjIzMO5AgEjfsGYqyczDK84NRmkdisy2bMq+NphInr1jlZVWxk4YiJ4U51jo4EiJCgVsocLtYX/biIdATkQRtYzFaR2O0jMbYNxLlke5kUXLaoKnEycZyF8eVu1hf6sLtyI2ClAlagJRSRyQejzMyMsLQ0NCc53USxtA+Ns32vgjb+6PsGY4ynUhOwd9Y7ODi1T7Wl7lYX+qk2Jsbw4KzocBt48RKNydWvjD6bDQUZ89wjN3DUXYNxrhtV5A/7ArisMH6UhcnVLo4sdLNqmIHthw9l3Q0tAAppdISCoUYHBxkeHh4dlbog03FEjzTH2Vbb4SneiOMR5LPayh0cGmTj00VLjaUufA5F2/rJhOKvXbOrLFzZk3yXMtULMHuoRjPDUR4tj/Kb3ZM8psdkxS6bZy8wsWpVR5OWrH4/ztqAVJKHZYxhvHxcQYGBggEAod8zlg4zhM9EZ7YH+a5/ijTBvKcwkkr3Jxc5ebEShfFnuXTwlkIPqeNU6rcnFKVbCWNh+Ns74/ydF+Erb0R7u8I47DBpnIXp6/0cMZKN0WL8L+xFiCl1EskEgmGhoYYGBggEom85PHxcJzH9kd4pCvM84NREiRPqF+6xsdp1R7WlTqx2469q8hut+N0OmeXLzhwbZ2ZtXYOXIvnYDNr+BhjZtf3mVnrZ2b9n+npaWKxWE5PZlrosbOl3suWei9xY9g7HOOJ/WGe7Ilw01MT/NdTsKHcxVk1bs6u8VC4SIqR6JQX6dm8ebPZunWr1TGUyqhYLMbAwACDg4MvGVQQnk7wxP4ID3aGeKY/SsLAynw7Z9d6OHOlh/pCxxGPTjtw2YGZ28zaOC6XK6uj3YwxxGKx2fWCZqazmVnWIRdnazDG0DUxzSPdYR7tjtA9MY1N4IQKFy+r83JGjRvvMQ5Jz8vLO6ZRcCKyzRiz+ZCPaQFKjxYgtZRFIhH6+voYHh5+0RDqhDHsGopxX3uIR7vChOOGcp+Nl9V5Oac2/aIjIng8Hnw+3+xyAh6PJ2cnyTyUWCxGOBwmFArNLgGRaxOldozHeKgzzEOdYQam4njswpk1bs5v8LKx3HVUAxgyWYC0C06pZSwcDtPb28vo6OiLvkhHw3HuawtxT1uI/mAcr0M4t87Dlnov68uc836RuVyu2eUBfD4fPp9v0S9FPdMVmJ+fP7stkUgQCoVmZ6OenJwkGo1alrG+0En98U7evimP3cMx7m8P8UhXmPs7wlT67VzU6OXCBm/OjDrUFlCatAWklpJQKDRbeGYYY3huIMqdLVNs7YkQN7Cx3MlFjT7OXOmZ83oUl8tFQUHB7BIAi6lls9BisdjsXGuBQOCQ59CyKRI3PN4d5p62EDsHo9gETqt2c8lqH8dXzN/NqS0gpdSCCIfD9PT0vKjwhGIJ7usIcWfzFPsDcfJdwj+s8fGKVT6q8w/9FWG328nPz5+d4n+pLRV9LJxO5+wUNpDs3pxZaiIQCGT9XJLbLpxX7+W8ei89gWnuaQvx97YpHt8foTrPzquafFzQ4MVrwZBubQGlSVtAajGbmVZlZn0YgIHgNHc0T3Fva4ipaUNTsZNXNfk4u9aD6xDTwbhcLoqKiigsLCQ/P9/y6XAWo5mZqMfHxxkbG7Osuy4aNzzSHebO5in2jcTwOYQLG71cusZHpf/Ff3RoC0gpdVRisRi9vb0MDQ3NnuNpGY3xpz1BHu0Kg8BZNR4uW+NjbelLu82snKp/KRKR2VZjbW0twWCQsbExRkdHs9pV57IL59d7Ob/ey97hKHc0T/HX5inu2DfFmTUeXrvOT1NJ5pf11gKk1BKUSCTo6+ujv79/9jqY5wai3LY7yHMDUXwO4dXr/Fza5KPM9+IT0jNT8peUlGjRyTC/34/f72flypUEg8HZWcRjsVjWMqwtdbG21MW7Tohzx74p7m6Z4pHuMJvKXbx+g5+zM7gsg3bBpUm74NRiYIxheHiYnp4eYrEYxhi29Ub4w64g+0ZiFHtsXLbGxytW+/Af0Odvs9koLi6mtLT0RaO8VPbNdNMNDw8zNjZ22GmPMmUqluBvrSFu3xtkNJxgbZmba19zIlvWlh/V/rQLTqllYGJigu7u7tlrU7b2Rvjtzknaxqap8Nn5wCkFXNDgfdH5Hb/fT1lZGSUlJYt+mPRScWA3XboTvy4kn9PGa9f5eVWTj/s7Qvx57xTNA5NHXYDmogVIqUUuEonQ1dXF+Pg4xhie6otw645JWsemWeG3c+VpBZxX58WRmhrHbrdTUlJCeXk5Xq/X4vRqLna7nfLycsrLywkGgwwODjI6OpqVVpHLLly8ysdrN5XTtKY+I++hBUipRSqRSNDb20t/fz/GGHYORvn1cwH2DMeoTBWeLXXe2TnZPB4PFRUVlJaWamtnEZo5X1RTU8PQ0BCDg4NZGUVntwmuDK0wqwVIqUVoZGSE7u5uYrEY7WMxfvlsgO39UUq8Nj54agEXNrzQ4ikoKKCyspKCggKLU6uF4HA4WLFiBZWVlYyNjdHf33/ES5/nCi1ASi0i4XCYzs5OAoEAg1Nxbt0R4IGOMH6n8E8n5HNJkw+3PTk7dElJCZWVldrNtkSJCMXFxRQXFzM5OUlfXx/j4+NWxzoiWoCUWgRmhlX39fUxFYvzx91Bbt8TxACvWefn9ev95Lls2Gw2ysrKqKysXNbT4Sw3eXl5NDU1EQqF6Ovre8ncfrlKC5BSOW5iYoLOzk5C4TD3tYf4zY5JxsIJXlbn4e2b8qnwJ9fGqaiooLKyEodDP9bLldfrpbGxkerqanp7exkZGcnpQqT/pyqVo6anp+nu7mZ4eJg9w1F++vQELaPTrCt18pmzi1hb6tLCow7J7XbT0NBAVVXV7Nx/uViI9P9YpXLQyMgIXV1dDE1G+OWzAe7vCFPisfHx0wt5WZ1ntvCsWLFCC486LLfbTWNjI1VVVezfv5+xsTGrI72I/p+rVA6JxWJ0dnYyPDrK31pC/HpHgOi04R/X+3nDBj8+p53S0lKqqqr0HI9Km8fjYfXq1QSDQfbv308gELA6EqAFSKmcMTw8TFdXF/uGwvxo2zgto9McX+HifScXUFPgoLCwkJqaGjwej9VR1SLl9/tZu3Yt4+PjdHd3Ew6HLc2jBUgpi8ViMTo6OugbHuXWHZPcsW+KAo+Nq88o5Jxaz+zFhzpHm1oohYWFFBQUMDQ0RE9PD9PT05bk0AKklIVmzvU80RXkx0+NMzSV4OJVXt55fD5FfjfV1dWUlZVZHVMtQSJCeXk5JSUl9Pb2MjAwkPWBClqAlLJAPB6no6ODzr5hfv7MBA90hKnJt3PdBSVsKHdTUVFBVVUVdrt9/p0pdQzsdjs1NTWUlZXR1dXFxMRE1t5bC5BSWTYxMUF7ezsPtgX4yVMTTEYTvHGDnzduyKOkqIC6ujo9z6OyzuPxsGbNGkZHR+nu7s7KPHNagJTKEmMM3d3dNHf18V9PT/BwV5hVRQ6+cF4xa8p91NTUUFJSYnVMtcwVFxdTWFhIT08PAwMDGX2vZVWARGQVcC1QaIx548H3rU2nlrJwOExraysPtoxy49YJgtEEb9uYx+vW+6mqrGDlypXa3aZyhs1mm/2DaGRkJHPvk6kdi0itiNwnIrtEZKeIfPwY9vUzERkQkR2HeOwSEdkjIs0ics1c+zHGtBpjrjjcfaUyYXBwkG3P7uQ/Huzj6w+PUeyx8Y2Xl/KuU8vZdNwG6urqtPionOTzJVvmmZLJFtA08EljzFMikg9sE5G/GWOen3mCiFQAIWNM4IBtTcaY5oP2dTPwfeC/D9woInbgB8ArgG7gSRH5M2AHrj9oH5cbYzLbnlTqADMDDR7d1893Hh9naCrO69f7ecvGfOpqqlmxYgUiMv+OlFqi5ixAIvKJNPYRNMb8+OCNxpheoDf1e0BEdgErgecPeNoW4EMicqkxJiwi7wf+Ebj0oH09KCINh3jv04FmY0xrKu+twGuNMdcDl6WRXamMCAaD7G1u4dfbR7htV5Byv52vXFDCKXVFNDQ06BIJSjF/F9yngDwgf47bJ+d7k1TxOBl4/MDtxpjfA3cCt4rIO4DLgTcfQf6VQNcB97tT2w6Xo1REfgScLCKfPfj+YV7zahG5abGts6Gs09/fzwPbdvLpu/r4w64g59V7uOHiMi46sZH169dr8VEqZb4uuF8aY/5trieIiH+ex/OA/wGuMsa8ZIC5MeabqZbLjcBqY8zkPJletPtDbDvslVTGmGHgnw/afPD9g19zO3D75s2b338EudQyFI/HaW9v5/Zne/nxtglE4BNnFvKKdaXa6lHqEOYsQMaYT8+3g7meIyJOksXn18aY2w7znJcBm4A/Al8EPjLfex6gG6g94H4N0HMEr1dqQUxNTfH83mZ++NgQf28Psa7UydVnFnH8qpVUV1fruR6lDiGtQQgishl4GVANhIAdwD3GmMOOz5PkJ+6nwC5jzA2Hec7JwE+AfwDagF+JyFeNMZ9PM/+TwBoRaQT2A28F3p7ma5VaEENDQzz0XAv/8cgo+wNx3rDBzz+dXMrqVY3k5eVZHU+pnDXnOSAReY+IPAV8FvACe4AB4FzgbyLyCxGpO8zLzwHeBVwoIttTt0sPeo4PeJMxpsUYkwDeDXQcIsctwKPAOhHpFpErAIwx0yRbTHcBu4DfGWN2pvUvV+oYJRIJ2tvbufmB3XzmnmECUcO/nlfMR86rZ9PG47T4KDWP+VpAfuAcY0zoUA+KyEnAGqDz4MeMMQ9x6HM0Bz7n4YPux0i2iA5+3tvm2McdwB1zvY9SCy0ajfL8nn1875F+7m0LsbHcxSfPLuGkdY06m4FSaZrvHNAP5nl8+4KmUWoRCAQCPPLsXr7+f0O0jU3zhg1+3ru5gjVNq3G73VbHU2rRSGsmBBH5pogUiIhTRO4VkSEReWemwymVa/r7+7nlgef4xJ39DATjfO7cIq66cDUbj9ugxUepI5TuVDwXp4ZQX0Zy5NlaktcIKbUsJBIJWlpbueHu3Vz/0CgVfjs3vLKCN5+7kdraWh3lptRRSHcqHmfq56XALcaYEf3AqeUiGo3y7K59fOOBXp7sibCl3sPV565gw9ombfUodQzSLUC3i8hukkOwPywi5YC1i4krlQXBYJAHn97NVx4YpH8yzhUn5fPOM2qpr6/HZsvYXL5KLQtpFSBjzDUi8g1gwhgTF5Ep4LWZjaaUtUZGRrjt0d3c8NgYDoEvbSnhklObKC8vtzqaUktC2rNhG2NGD/g9CAQzkkipHNDd3c1ND7byq2cD1BU6+Pz55Zx1wjr8/jlnnlJKHYFltSCdUvNJJBLsaW7hm3/v4r72EGfVeLjm/GqOW9eE0+mcfwdKqbRpAVIqJRaLsXXHHv7t3h52DcV483F+PnxeA/X19TrKTakMSHcuOGdqloIDt5UZY4YyE0up7AqFQty3bRdfum+AsVCcT5xZxFvPXkNFRYXV0ZRasuabC+4CEekGekTk7oMWhbs7o8mUypKJiQl+98Az/MtdfUSnDV+9qJz3XHi8Fh+lMmy+caTfBF5pjCkHbiI5AemZqce0T0ItekNDQ/z4b8/ybw8MU+q1c8OrqnjNOSdQUFBgdTSllrz5uuBcM7NLG2P+kFpW+zYRuYY5Fn5TajHYv38/37pnH3/YFeTEShdffPlKTtiwFodDT40qlQ3zfdJiIrLCGNMHYIzZKSIXAX8BVmc8nVIZYIyhubWN6+5u4/6OMBc2ePnsyxtoWt2ogw2UyqL5CtA1QCXQN7PBGNMtIucDV2YullKZkUgkeHbXXr5wVxfPDkR5y8Y8PnZhEytXrrQ6mlLLznzLMdxzmO1jwHWZCKRUpkxPT/P4s7v4/N09dIxPc+Vphbx3y3rKysqsjqbUspTuMOzLgK8A9anXCGCMMXqmVi0KkUiE+7c9z7/e2894OMHnXlbCm845jsLCQqujKbVspXu29dvA64HnjDE6+EAtKlNTU/z18ef58n0DAHz1onIuPeM4nVZHKYulW4C6gB1afNRiEwgEuO3hnVz/0Aj5LuGrL1/B+acep8soKJUD0i1AnwbuEJEHgMjMRmPMDRlJpdQCGBsb4zcPPs8Nj41SlefguourOOvEDTqnm1I5It0CdB0wCXgAV+biKLUwhoeH+dl9u7hx6zhrSp189ZU1nLJxHXa73epoSqmUdAtQiTHm4owmUWqBDAwMcOO9u/n5MwFOrHTx1Uvq2biuSReQUyrHpPuJvEdEtACpnNfT08N/3rmLnz8T4MyVbr5xWSOb1q/R4qNUDkq3BXQl8GkRiQAxdBi2ykHd3d188+59/HnvFBc2eLn2latY3dhgdSyl1GGkuyR3fqaDKHUs2js6uP7uFu5qCXFpk4/PXLya+ro6q2MppeaQVr+EiPyjiBQecL9IRF6XsVRKHYGW1ja+/Ndk8XndOj+fu2SNFh+lFoF0O8a/aIwZn7mTmornixlJpFSajDE0t7TypTtbua89xFuOy+PTr1xLTU2N1dGUUmlI9xzQoQqVzlmvLGOMYW9zC1++q4NHusO84/g8PnrROlasWGF1NKVUmtJtAW0VkRtEZLWIrBKRbwHbMhlMqcMxxrBnXwtfuLOdR7rDvPuEfD728vVafJRaZNItQB8FosBvgd8BIXQ5BmUBYwy79zbzxbvaeXx/hPeelM+HL1pHZWWl1dGUUkco3VFwQZJrAyllmUQiwe69zXz5b508vj/CFSfl84ELtPgotVjN2QISkZtE5PjDPOYXkctF5B2ZiabUCxKJBHv2afFRaimZrwX0Q+BfU0VoBzBIcj64NUAB8DPg1xlNqJa9RCLB3n3NfPnuF4rP+89fq8VHqUVuvhVRtwNvFpE8YDNQRfL8zy5jzJ7Mx1PLnTGGvfua+co9XTy2P8J7T8znfVvW6IADpZaAdM8BTQL3ZzaKUi82U3yu/3sXD3eFedcJ+bx/SxNVVVVWR1NKLQCdoVHlpORFpi38+31dPNAR5m2b8vjgeauorq62OppSaoFoAVI5xxhDa2sr37q/i3vaQrxxg58PndeoMxwotcTMNwrusyJycrbCKAXQ1tbG9x7s4o7mKV691seV59VTp3O7KbXkzHcOqA34uIicCDwD/BW42xgzmvFkallqb2/nJ4908b97gly8ysvHz6uloaHB6lhKqQyYbxTcrcCtAKmW0CXAbSJiB+4B7jTGPJHxlGpZ6Ozs5JePdXLLjkm21Hv45Pk1rF69GhGxOppSKgPSnlDUGPM08DRwvYgUAK8A3gdoAVLHbP/+/fz+yU5+uj3A6dVuPr2lmjVNTVp8lFrCjmpGa2PMBPA/qZtSx6S3t5c/bWvnB0+Oc0KFi2svqGL9urW6jLZSS5x+wpWlBgYGuGt7Gzc8NkZTiZPPn1/BcevXYrfbrY6mlMowLUDKMsPDw9z3TCtff3iMqjwHXzi/nBM3rsfpdFodTSmVBfMNw35eRK4VkdXZCqSWh7GxMR55roXr/m+EApeNL24p5ZRN63C73VZHU0plyXwtoLcBecDdIvK4iFwlInopujomgUCAJ3fs48sPDiMifGFLCacfvxafz2d1NKVUFs1ZgIwxzxhjPmuMWQ18HKgHHhORv4vI+7OSUC0pU1NTPPP8Xr7y4AiTUcPnX1bM2cc3UVBQYHU0pVSWpX0OyBjzmDHmauCfgGLg+xlLpZakcDjMjl17ue7/hukJTPPps4s47/hGSkpKrI6mlLJAWsOwReQ0kt1xbwDagZuA32cullpqYrEYu/fs5T8fGWLXUIyrzyjkFcfX6po+Si1jcxYgEfka8BZglOSMCOcYY7qzEUwtHfF4nL179/KjJ4Z5PLWmz2UnVFFbW2t1NKWUheZrAUWAVxlj9mYjjFp6EokEzc3N3LJ9iDtbpnjNWh9vObmCxsZGq6MppSw23yCELxtj9oqIT0T+VUR+AiAia0TksuxEVIuVMYa2tjb+385BfvXcJOfUenjfaWU0NTXpLAdKqbQHIfycZGvorNT9buCrGUmkloyuri7+b28/P3xynI3lLq4+q5R1a9ficBzVDFBKqSUm3QK02hjzTSAGYIwJATpLpDqs3t5etjX38M1HxliR7+Cac4pZv7ZJLzRVSs1KtwBFRcQLGIDUzAiRjKVSi9rw8DA7W7q47qFRXHbh2nOL2bh2FXl5eVZHU0rlkPmm4rk79euXgDuBWhH5NXAv8OnMRlOL0cTEBHta2rj+4VECEcPnzi3mpDW1eq2PUuol5uuMLwcwxtwtItuAM0l2vX3cGDOU6XBqcZmammLvvma+9dgYbaPTfOacIk5bXUlVVZXV0ZRSOWi+AlQoIq8/xPbzRARjzG2ZCKUWn2g0SnNzMzdvH+fJnghXnJTPljWlupy2Uuqw5i1AwGUcesCBAbQAKeLxOM3Nzdy+e5y/7Jvi0iYfr9tYTJOuaKqUmsN8BajDGHN5VpKoRckYQ0tLC4+2jfHTpyc4tcrN5ScX0tTUpMOtlVJzmm8UnP75qubU0dHBrv2j/OejY9QWOLj6jEKaVq/C6/VaHU0plePmK0DvykoKtSj19fXR1jPI11LDra85p5imhlqKioqsjqaUWgTmK0Bfn28HIvKXBcqiFpHR0VHau7r55iOjjITifOacItbXlrNixQqroymlFon5OunPFZE/z/G4AMctYB61CASDQdra2vjJUxPsGopx1RmFnFRbRH19vdXRlFKLyHwF6LVp7CO6EEHU4hCNRmlpaeH2vZPc2xbiDRv8XLi6gNWrV+sEo0qpIzJnATLGPJCtICr3zQy3frJrkl9sD3B6tZu3bcpn1apVuFwuq+MppRYZ/ZNVpWVmaYWWgQD/+dgYtYUOPnZGIQ319TrHm1LqqGgBUmnZv38/PUOjfP3hURwC15xTRO2KCsrKyqyOppRapNIuQCLiFZF1mQyjctPQ0BA9fX185/Fx+ibj/MvZxTRWFFJXV2d1NKXUIpZWARKRVwPbSc6IjYicNM/oOLVETE5O0tnZya07JtnWG+Hykwo4qdrP6tWrdZodpdQxSbcF9CXgdGAMwBizHWjIRCCVO2ZGvD3cGeK23UFescrLJU0+Vq1ahdPptDqeUmqRS7cATRtjxjOaROWURCJBc3MzLcMhvv/kOOtKnVxxUgG1tbU66EAptSDSnS1yh4i8HbCLyBrgY8AjmYulrNbW1sbAeJBvPDyGzyl86qwiKstLqaiosDqaUmqJSLcF9FFgI8lluG8BJoCrMpRJWay3t5fh0VG+/dgYw1NxPnV2EdUleTrTgVJqQaXVAjLGTAHXpm5qCRsbG6Onp4dbd0yyvT/KB08t4LgKr850oJRacHMWIBG5neTCc4dkjHnNgidSlgmFQrS1tfHE/jC37Q5yUaOXi1f5aGhowO12Wx1PKbXEzNcC+o/Uz9cDK4Bfpe6/DWjPUCZlgXg8TktLC93jUb73xDirix287+QCKisrdXkFpVRGpDUXnIh8xRhz3gEP3S4iD2Y0mcqq1tZWxoIhvvHIGA4bfOrsYkoK81m5cqXV0ZRSS1S6nfrlIrJq5o6INALlmYmksm3//v2Mj49z49YJeiamuerMIqoK3KxatUovNlVKZUy6w7CvBu4XkdbU/QbggxlJpLJqdHSUvr4+/to8xcNdYd6xKY8TK900NjbqxaZKqYxKdxTcnanrf9anNu02xkQyF0tlQzgcpr29nT3DUW5+JsDmKjevW++nqqqKgoICq+MppZa4dFtAAKeSbPk4gBNFBGPMf2cklcq4mUEHo1Mx/uPRMcp8dj52eiFFhYVUV1dbHU8ptQykVYBE5JfAapITksZTmw2gBWiRam9vJxgK8e3Hx5mMJPjaRaUU+ZNdb0oplQ3ptoA2A8cZYw57TZBaPHp7exkbG+MPz0/y7ECUD20uoLHISWNjIw7HkTSKlVLq6KU7Cm4HyeuA1CI3MTFBT08P2/si/P75IOfXe7iowUtVVRX5+flWx1NKLSPp/rlbBjwvIk+QnA8O0JkQFptoNEpbWxvDoTjfeXyM2gIHHzilkIKCAj3vo5TKunQL0JcyGUJlnjGG1tZWItEY33psjGgC/uXsIvwep573UUpZIt1h2A9kOojKrK6uLoLBILfunGTXUIyrzihkZb6DhoYGvd5HKWWJ+SYjDTD3ZKR6scgiMDIywuDgIE/3RWZXNn1ZnZeKigoKCwutjqeUWqbmmwsuH0BE/g3oA34JCPAOQM9YLwLhcJiOjg6Gp+J89/Ex6gsdvPekAnw+HzU1NVbHU0otY+mOgnulMeaHxpiAMWbCGHMj8IZMBlPHLpFI0NLSQmw6zrceHyMah0+eVYTXadd53pRSlku3AMVF5B0iYhcRm4i8gxcuSFU5qqOjg3A4zO+eT573+cCpBazMd1BXV6fr+yilLJduAXo78GagP3V7U2qbylGDg4OMjIzw3ECE/9kV5IIGL1vqvZSUlFBaWmp1PKWUSnsUXDvw2sxGUQtlamqKrq4uxsNxvv34OFX5dt53cj4ul4u6ujqr4ymlFJBGC0hEXikiV4hI/UHbL89cLHW04vE4ra2txBMJvvfkOMFogk+emTzv09jYiN1utzqiUkoB8xQgEfkacC1wPPB3EfnoAQ9/JJPB1NHp6OggEonw//ZN8XRflHefmE9DkZMVK1aQl5dndTyllJo1Xwvo1cCFxpirSC7H8CoR+VbqMR1ClWMGBwcZHR2ldTTGr54NcFq1m0tW+/D7k2v8KKVULpmvADmMMdMAxpgxkgWpQER+D7gynE0dgZnzPqFYghseG6PAY+PKzYXY7cmuNx1yrZTKNfMVoBYR2TJzxxgTN8ZcAewBNmQ0mUrbzHkfYww/3R6gbzLOVWcUke+2UVtbq0OulVI5ab4C9CbgiYM3GmM+D9RmJJE6Yp2dnUQiER7uCnFfe4g3bPCzsdxFUVERZWVlVsdTSqlDmrMAGWNCxpiQJL1TRL4AICJ1wMqsJFRzGh4eZmRkhIFgnB9tm2BNiZM3HZeH0+mkvr5+/h0opZRF0r0Q9YfAWcDbUvcDwA8ykkilLRwO09nZSdwYvvvEGMbAVWcU4rAJ9fX1urqpUiqnpVuAzjDGXAmEAYwxo+ggBEslEglaW1tJJBL8cXeQXUMx3ndKASvyHJSXl+ss10qpnJduAYqJiJ3U0gwiUg4kMpZqgYnIKhH5qYj8Ya5ti0l3dzehUIh9I1F+u3OSc2s9bKnz4Ha7dZZrpdSikG4B+i7wR6BCRK4DHgK+lrFUBxCRn4nIgIjsOGj7JSKyR0SaReSaufZhjGlNjd6bc9tiMTY2xuDgIKHpBN9+fJwSr40PnFKAzWajsbERmy3dw6qUUtZJdy64X4vINuAikhegvs4YsyujyV5wM/B94L9nNqRaYz8AXgF0A0+KyJ8BO3D9Qa+/3BgzkJ2omReNRuno6ADg5u0B+ifjfPn8EvwuG5WVlfj9fosTKqVUetI+S22M2Q3szmCWw73vgyLScNDm04FmY0wrgIjcCrzWGHM9cFmWI2ZVe3s709PTPLE/zD1tIV63Ljnk2uv1Ul1dbXU8pZRK22Ltq1kJdB1wv5s5hoWLSKmI/Ag4WUQ+e7hth3jdB0Rkq4hsHRwcXMD4R6e3t5dAIMBoOM6NW8dpLHLw1k15iIjOdqCUWnQW6zjdQ33TmsM92RgzDPzzfNsO8bqbgJsANm/efNj9Z0MwGKS3txdjDDdunSA8bfj4GUU4bUJ1dTVer9fKeEopdcQWawuomxfPxFAD9FiUJePi8ThtbW0YY/hbW4htvRHeeUI+tQUO/H4/lZWVVkdUSqkjtlgL0JPAGhFpFBEX8FbgzxZnypiuri4ikQi9k9PcvD3ACRUuXtXkw2az0dDQoF1vSqlFKecLkIjcAjwKrBORbhG5IjVD90eAu4BdwO+MMTutzJkpo6OjDA8PE08Yvvv4OA4bfOS0QmyS7HrzeDxWR1RKqaOS8+eAjDFvO8z2O4A7shwnqw4ccv2nPUH2jsS46oxCSn128vLytOtNKbWo5XwLaDlrb29Pnv8Zi/HbnZOcXePh3FoPNptNJxpVSi16WoByVF9fH4FAgFg82fWW507OdiDa9aaUWiK0AOWgqakpenqSg/pu3TlJ58Q0H95cQL7bpl1vSqklQwtQjkkkErNDrncPRfnTniAvb/RyapV2vSmllhYtQDlm//79hMNhItOG7z85TpnPxntOzAfQrjel1JKiBSiHTExMMDCQnDf1V88F6J2Mc+VphXidNvx+PxUVFRYnVEqphaMFKEfE43Ha29sB2DEQ4Y7mKV7V5OP4CjciyRVO9YJTpdRSogUoR3R2dhKLxQhNJ/jB1glW5Nl55/F5AFRVVelcb0qpJUcLUA4YHR1lZGQEgF89O8lgMM5HTivE47Dh9XpZsWKFxQmVUmrhaQGyWCwWo7OzE4DnBiLc2TLFP6zxsaHMhYjoXG9KqSVLC5DFZhaYC00n+GGq6+3tm5Kj3ioqKvD5fBYnVEqpzNACNA8RebWI3DQ+Pr7g+x4cHGRiYgJ4oevtys2FuB2C2+3WFU6VUkuaFqB5GGNuN8Z8oLCwcEH3G4lE6O7uBmDnYJQ7W6a4dI2P48pdANTX12Oz6eFRSi1d+g1nkfb2dhKJBJFpww+fHGeF3847Ul1vZWVl5OfnW5xQKaUySwuQBfr7+5mcnATgNzsC9AXjfGhzAW6H4HQ6qampsTihUkplnhagLAuHw7MTje4ZjvL/9k3xytVeNlW4AaitrcVut1sZUSmlskILUBYZY2a73qJxww+eHKfUZ+Ndxye72woLCykuLrY4pVJKZYcWoCzq7+8nGAwC8PvnJ9kfiPPPpybnerPb7dTV1VmcUCmlskcLUJaEQqHZrre2sRj/uyfIBQ1eTl6R7Hqrrq7G5XJZGVEppbJKC1AWzHS9GWOIJ5JdbwUuG+9OLbPg8/l0pmul1LKjBSgLBgYGmJqaAuDPe4O0jU3zvlMKyHfZZme6Vkqp5UYLUBaEw2EAegLT/HbnJGesdHNWTXJhOZ1uRym1XGkBypKEMfxo2zguu/D+kwsAcLlcOt2OUmrZ0gKUJfe2hdg5GOOfTsyn2Ju8zqeurk6n21FKLVv67ZcFw8EY//1sgE3lLi5qSC4sV1RUxELPL6eUUouJFqAs+M5DvcTihg+eWoCIYLPZqK2ttTqWUkpZSgtQht21s48HWid488Y8qvMdgF7zo5RSAA6rAyx1G1YU8LqNJbxmrRPQa36UUmqGtoAyrK7Ux9XnVeOwJZfVrqur0yW2lVIKLUBZVVZWht/vtzqGUkrlBC1A81ioJbkdDgcrV65coFRKKbX4aQGax0Ityb1y5UocDj3lppRSM/QbMQvy8vIoLS21OoZSSuUUbQFlgRYfpZR6KS1ASimlLKEFSCmllCW0ACmllLKEFiCllFKW0AKklFLKElqAlFJKWUILkFJKKUtoAVJKKWUJLUBKKaUsIcYYqzMsCiIyCHQcsKkQOHiG0kNtKwOGMhhtLofKk639pPOaY3nOsW7X43Jsz5vr8SM5NvqZObLXLMbPTL0xpvyQjxhj9HYUN+CmNLdtzaWM2dpPOq85lucc63Y9Lsf2vLkeP5Jjo5+ZI3vNUvvMaBfc0bs9zW1WWqg8R7OfdF5zLM9ZqO1WyPXjks7z5nr8SI5BLh0XyP1js6Q+M9oFl2EistUYs9nqHOrF9LjkLj02uSkTx0VbQJl3k9UB1CHpccldemxy04IfF20BKaWUsoS2gJRSSllCC5BSSilLaAFSSillCS1AFhKRVSLyUxH5g9VZljsR8YvIL0TkJyLyDqvzqBfo5yQ3icjrUp+XP4nIxUezDy1AR0lEfiYiAyKy46Dtl4jIHhFpFpFr5tqHMabVGHNFZpMuX0d4jF4P/MEY837gNVkPu8wcybHRz0n2HOFx+d/U5+U9wFuO5v20AB29m4FLDtwgInbgB8CrgOOAt4nIcSJyvIj85aBbRfYjLzs3k+YxAmqArtTT4lnMuFzdTPrHRmXPzRz5cfl86vEj5ji6jMoY86CINBy0+XSg2RjTCiAitwKvNcZcD1yW5YjL3pEcI6CbZBHajv5hlnFHeGyez3K8ZetIjouI7AK+DvzVGPPU0byfftAW1kpe+Csakl9qKw/3ZBEpFZEfASeLyGczHU4Bhz9GtwFvEJEbyb3pYZaLQx4b/ZxY7nCfmY8CLwfeKCL/fDQ71hbQwpJDbDvslb7GmGHgqA6cOmqHPEbGmCDw3myHUS9yuGOjnxNrHe64fBf47rHsWFtAC6sbqD3gfg3QY1EWdWh6jHKXHpvclLHjogVoYT0JrBGRRhFxAW8F/mxxJvVieoxylx6b3JSx46IF6CiJyC3Ao8A6EekWkSuMMdPAR4C7gF3A74wxO63MuZzpMcpdemxyU7aPi05GqpRSyhLaAlJKKWUJLUBKKaUsoQVIKaWUJbQAKaWUsoQWIKWUUpbQAqSUUsoSWoCUOkYiEheR7Qfc5lyGI1sOyFWdup8nIjeKSIuIPC0i20Tk/fPs434ReeVB264SkR+KyOrU/icz+e9QS5fOBafUsQsZY05ayB2KiCN1AeCxODjXfwGtwBpjTEJEyoHL59nHLSSvfL/rgG1vBT5ljGkBTtICpI6WXoiq1DESkUljTN4htrcDvwBeDTiBNxljdouIH/gecDzJPwK/ZIz5k4i8B/gHwAP4SS7hcTOwnuQV6A3AlcCJwCZjzNWp93k/sMEY84nD5RKR1cDfgCZjTOIw/45PAW8G3MAfjTFfFJFSYDdQY4yJpKbqfxCoN6kvj8P9+5Waj3bBKXXsvAd1wR24OuSQMeYU4EbgX1LbrgX+bow5DbgA+PdUUQI4C3i3MeZC4MPAqDHmBOArwKmp59wKvEZEnKn77wV+Pk/GjcAzcxSfi4E1JNd+OQk4VUTOS81E/QQvLFL2VuC3Rv9yVQtAu+CUOnZzdcHdlvq5jeSy3wAXkywgMwXJA9Slfv+bMWYk9fu5wHcAjDE7ROTZ1O9BEfk7cFlqUTCnMea5IwksItcCbwIqjDHVqUwXA0+nnpJHsiA9yAvdcH9K/Zyv206ptGgBUiqzIqmfcV74vAnwBmPMngOfKCJnAMEDN82x3/8CPkeye2y+1g8kVxU9UURsxpiEMeY64LoDzt8IcL0x5seHeO3/AjeIyCmA92hXv1TqYNoFp1T23QV8VEQEQEROPszzHiJ5TgYROY7kOSMAjDGPk1yj5e0kWyhzMsY0A1uBr4qIPbVPDy8UubuAy0Vk5pzRShGpSL12Ergf+Fk676VUurQFpNSx84rI9gPu32mMmWso9leAbwPPpopQO8kBBwf7IfCLVNfb08CzwPgBj/8OOMkYM5pmzvcB/w40i8gIEAI+A2CMuVtENgCPpuriJPBOYCD12ltIdie+Nc33UmpeOgpOqRyVaqk4jTHh1Ci2e4G1xpho6vG/AN8yxtx7mNdnZXSajoJTR0u74JTKXT7gIRF5Bvgj8CFjTFREikRkL8nBD4csPikTB16IutBmLkQF+jOxf7X0aQtIKaWUJbQFpJRSyhJagJRSSllCC5BSSilLaAFSSillCS1ASimlLKEFSCmllCX+P20P0mC8/tSAAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "energy_bounds = (100 * u.MeV, 100 * u.GeV)\n",
    "opts = dict(sed_type=\"e2dnde\", yunits=u.Unit(\"TeV cm-2 s-1\"))\n",
    "model.spectral_model.plot(energy_bounds, **opts)\n",
    "model.spectral_model.plot_error(energy_bounds, **opts);"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "You can create initial source models for your analyses using the `.to_models()` method of the catalog objects. Here for example we create a `Models` object from the 4FGL catalog subset we previously defined:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:38.140640Z",
     "iopub.status.busy": "2021-11-22T21:09:38.140274Z",
     "iopub.status.idle": "2021-11-22T21:09:38.209798Z",
     "shell.execute_reply": "2021-11-22T21:09:38.210044Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<gammapy.modeling.models.core.Models at 0x12872b940>"
      ]
     },
     "execution_count": 32,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "models_4fgl_roi = catalog_4fgl_roi.to_models()\n",
    "models_4fgl_roi"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Specificities of the HGPS catalog\n",
    "\n",
    "Using the `.to_models()` method for the `gammapy.catalog.SourceCatalogHGPS` will return only the models components of the sources retained in the main catalog, several candidate objects appears only in the Gaussian components table (see section 4.9 of the HGPS paper, https://arxiv.org/abs/1804.02432). To access these components you can do the following:\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:38.212232Z",
     "iopub.status.busy": "2021-11-22T21:09:38.211951Z",
     "iopub.status.idle": "2021-11-22T21:09:38.213157Z",
     "shell.execute_reply": "2021-11-22T21:09:38.213335Z"
    }
   },
   "outputs": [],
   "source": [
    "discarded_ind = np.where(\n",
    "    [\n",
    "        \"Discarded\" in _\n",
    "        for _ in catalog_hgps.table_components[\"Component_Class\"]\n",
    "    ]\n",
    ")[0]\n",
    "discarded_table = catalog_hgps.table_components[discarded_ind]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "There is no spectral model available for these components but you can access their spatial models:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:38.217929Z",
     "iopub.status.busy": "2021-11-22T21:09:38.217642Z",
     "iopub.status.idle": "2021-11-22T21:09:38.218796Z",
     "shell.execute_reply": "2021-11-22T21:09:38.219019Z"
    }
   },
   "outputs": [],
   "source": [
    "discarded_spatial = [\n",
    "    catalog_hgps.gaussian_component(idx).spatial_model()\n",
    "    for idx in discarded_ind\n",
    "]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "In addition to the source components the HGPS catalog include a large scale diffuse component built by fitting a gaussian model in a sliding window along the Galactic plane. Information on this model can be accessed via the propoerties `.table_large_scale_component` and `.large_scale_component` of `gammapy.catalog.SourceCatalogHGPS`."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:38.221462Z",
     "iopub.status.busy": "2021-11-22T21:09:38.221170Z",
     "iopub.status.idle": "2021-11-22T21:09:38.222577Z",
     "shell.execute_reply": "2021-11-22T21:09:38.222803Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div><i>Table length=5</i>\n",
       "<table id=\"table4978974384\" class=\"table-striped table-bordered table-condensed\">\n",
       "<thead><tr><th>GLON</th><th>GLAT</th><th>GLAT_Err</th><th>Surface_Brightness</th><th>Surface_Brightness_Err</th><th>Width</th><th>Width_Err</th></tr></thead>\n",
       "<thead><tr><th>deg</th><th>deg</th><th>deg</th><th>1 / (cm2 s sr)</th><th>1 / (cm2 s sr)</th><th>deg</th><th>deg</th></tr></thead>\n",
       "<thead><tr><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th></tr></thead>\n",
       "<tr><td>270.000000</td><td>0.205357</td><td>0.251932</td><td>6.149827e-10</td><td>4.064108e-10</td><td>0.269385</td><td>0.137990</td></tr>\n",
       "<tr><td>272.959198</td><td>-0.120154</td><td>0.058201</td><td>1.426735e-09</td><td>7.346488e-10</td><td>0.088742</td><td>0.041882</td></tr>\n",
       "<tr><td>275.918365</td><td>-0.095232</td><td>0.089881</td><td>1.193710e-09</td><td>6.117877e-10</td><td>0.167219</td><td>0.111797</td></tr>\n",
       "<tr><td>278.877563</td><td>-0.257642</td><td>0.065071</td><td>1.506986e-09</td><td>5.230542e-10</td><td>0.156525</td><td>0.056130</td></tr>\n",
       "<tr><td>281.836731</td><td>-0.283487</td><td>0.066442</td><td>1.636973e-09</td><td>4.336444e-10</td><td>0.205192</td><td>0.049676</td></tr>\n",
       "</table></div>"
      ],
      "text/plain": [
       "<Table length=5>\n",
       "   GLON       GLAT   GLAT_Err Surface_Brightness Surface_Brightness_Err  Width   Width_Err\n",
       "   deg        deg      deg      1 / (cm2 s sr)       1 / (cm2 s sr)       deg       deg   \n",
       " float32    float32  float32       float32              float32         float32   float32 \n",
       "---------- --------- -------- ------------------ ---------------------- -------- ---------\n",
       "270.000000  0.205357 0.251932       6.149827e-10           4.064108e-10 0.269385  0.137990\n",
       "272.959198 -0.120154 0.058201       1.426735e-09           7.346488e-10 0.088742  0.041882\n",
       "275.918365 -0.095232 0.089881       1.193710e-09           6.117877e-10 0.167219  0.111797\n",
       "278.877563 -0.257642 0.065071       1.506986e-09           5.230542e-10 0.156525  0.056130\n",
       "281.836731 -0.283487 0.066442       1.636973e-09           4.336444e-10 0.205192  0.049676"
      ]
     },
     "execution_count": 35,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# here we show the 5 first elements of the table\n",
    "catalog_hgps.table_large_scale_component[:5]\n",
    "# you can also try :\n",
    "# help(catalog_hgps.large_scale_component)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Flux points\n",
    "\n",
    "The flux points are available via the `flux_points` property as a `gammapy.spectrum.FluxPoints` object."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:38.274184Z",
     "iopub.status.busy": "2021-11-22T21:09:38.236173Z",
     "iopub.status.idle": "2021-11-22T21:09:38.276761Z",
     "shell.execute_reply": "2021-11-22T21:09:38.277006Z"
    }
   },
   "outputs": [],
   "source": [
    "source = catalog_4fgl[\"PKS 2155-304\"]\n",
    "flux_points = source.flux_points"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 37,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:38.278968Z",
     "iopub.status.busy": "2021-11-22T21:09:38.278643Z",
     "iopub.status.idle": "2021-11-22T21:09:38.280026Z",
     "shell.execute_reply": "2021-11-22T21:09:38.280253Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<gammapy.estimators.points.core.FluxPoints at 0x128af9d00>"
      ]
     },
     "execution_count": 37,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "flux_points"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:38.282577Z",
     "iopub.status.busy": "2021-11-22T21:09:38.282295Z",
     "iopub.status.idle": "2021-11-22T21:09:38.286942Z",
     "shell.execute_reply": "2021-11-22T21:09:38.287113Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div><i>Table length=7</i>\n",
       "<table id=\"table4974795840\" class=\"table-striped table-bordered table-condensed\">\n",
       "<thead><tr><th>e_ref</th><th>e_min</th><th>e_max</th><th>flux</th><th>flux_errp</th><th>flux_errn</th><th>flux_ul</th><th>sqrt_ts</th><th>is_ul</th></tr></thead>\n",
       "<thead><tr><th>MeV</th><th>MeV</th><th>MeV</th><th>1 / (cm2 s)</th><th>1 / (cm2 s)</th><th>1 / (cm2 s)</th><th>1 / (cm2 s)</th><th></th><th></th></tr></thead>\n",
       "<thead><tr><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float32</th><th>bool</th></tr></thead>\n",
       "<tr><td>70.71067811865478</td><td>49.99999999999999</td><td>100.00000000000004</td><td>8.811208118686409e-08</td><td>3.097255785178277e-08</td><td>3.039940565940924e-08</td><td>nan</td><td>2.5337915</td><td>False</td></tr>\n",
       "<tr><td>173.20508075688775</td><td>100.00000000000004</td><td>299.99999999999994</td><td>6.859372092549165e-08</td><td>3.0207172319052233e-09</td><td>3.0207172319052233e-09</td><td>nan</td><td>25.057798</td><td>False</td></tr>\n",
       "<tr><td>547.722557505166</td><td>299.99999999999994</td><td>999.9999999999998</td><td>3.348723609519766e-08</td><td>6.225928661507396e-10</td><td>6.225928661507396e-10</td><td>nan</td><td>88.745026</td><td>False</td></tr>\n",
       "<tr><td>1732.0508075688763</td><td>999.9999999999998</td><td>2999.9999999999977</td><td>1.274911376469845e-08</td><td>2.198041054723987e-10</td><td>2.198041054723987e-10</td><td>nan</td><td>128.2931</td><td>False</td></tr>\n",
       "<tr><td>5477.225575051666</td><td>2999.9999999999977</td><td>10000.00000000001</td><td>5.410327741373067e-09</td><td>1.2490801448716837e-10</td><td>1.2490801448716837e-10</td><td>nan</td><td>120.772316</td><td>False</td></tr>\n",
       "<tr><td>17320.50807568877</td><td>10000.00000000001</td><td>30000.000000000007</td><td>1.7687512565700558e-09</td><td>6.802798602212334e-11</td><td>6.802798602212334e-11</td><td>nan</td><td>82.056656</td><td>False</td></tr>\n",
       "<tr><td>94868.32980505146</td><td>30000.000000000007</td><td>299999.9999999999</td><td>6.965684695714458e-10</td><td>4.1411239021238444e-11</td><td>4.1411239021238444e-11</td><td>nan</td><td>54.02705</td><td>False</td></tr>\n",
       "</table></div>"
      ],
      "text/plain": [
       "<Table length=7>\n",
       "      e_ref              e_min              e_max                 flux          ...       flux_errn          flux_ul    sqrt_ts   is_ul\n",
       "       MeV                MeV                MeV              1 / (cm2 s)       ...      1 / (cm2 s)       1 / (cm2 s)                 \n",
       "     float64            float64            float64              float64         ...        float64           float64    float32    bool\n",
       "------------------ ------------------ ------------------ ---------------------- ... ---------------------- ----------- ---------- -----\n",
       " 70.71067811865478  49.99999999999999 100.00000000000004  8.811208118686409e-08 ...  3.039940565940924e-08         nan  2.5337915 False\n",
       "173.20508075688775 100.00000000000004 299.99999999999994  6.859372092549165e-08 ... 3.0207172319052233e-09         nan  25.057798 False\n",
       "  547.722557505166 299.99999999999994  999.9999999999998  3.348723609519766e-08 ...  6.225928661507396e-10         nan  88.745026 False\n",
       "1732.0508075688763  999.9999999999998 2999.9999999999977  1.274911376469845e-08 ...  2.198041054723987e-10         nan   128.2931 False\n",
       " 5477.225575051666 2999.9999999999977  10000.00000000001  5.410327741373067e-09 ... 1.2490801448716837e-10         nan 120.772316 False\n",
       " 17320.50807568877  10000.00000000001 30000.000000000007 1.7687512565700558e-09 ...  6.802798602212334e-11         nan  82.056656 False\n",
       " 94868.32980505146 30000.000000000007  299999.9999999999  6.965684695714458e-10 ... 4.1411239021238444e-11         nan   54.02705 False"
      ]
     },
     "execution_count": 38,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "flux_points.to_table(sed_type=\"flux\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 39,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:38.289055Z",
     "iopub.status.busy": "2021-11-22T21:09:38.288758Z",
     "iopub.status.idle": "2021-11-22T21:09:38.487702Z",
     "shell.execute_reply": "2021-11-22T21:09:38.487947Z"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAaAAAAEKCAYAAABUsYHRAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAa5UlEQVR4nO3de7RdZXnv8e8vITSGYDiEzTiUkOYKgiik7OChSEFrECkBBeTajiAZ4WLBWmy5iK0XbLH2jBQZIBcJJ0etCQHlQCiCKLIRG2MSJCEQkAAh7IN2S0ICMVwS8vSPNTesLtZl7r3XXHNdfp8x9iDvO29PXjc8PnO+c76KCMzMzBptWN4BmJlZZ3ICMjOzXDgBmZlZLpyAzMwsF05AZmaWCycgMzPLxU55B9Aq9thjj5gwYULeYZiZtZQVK1a8GBFd5bY5AaU0YcIEli9fnncYZmYtRdJzlbb5FpyZmeXCCcjMzHLhBGRmZrlwAjIzs1w4AZmZWS6cgMzMLBdOQGZmlgsnIDNreqfesIRTb1iSdxhWZ05AZmaWCycgMzPLhROQmTW19Ru2srJ3E0uf3ciMuT2s37A175CsTvwtODOru3o+r1nZu4nXtu0A4Km+LRx9VQ8HjdutLue+5dzD6nIeGxxXQGbW1PqTT6W2tS5XQGZWd/WsLGbM7eGpvi0ADBNM7hrtyqVNuAIys6Y2b9Z0Ro4o/Kdqctdo5s2annNEVi+ugMysqY0fO+qtZz6ufNqLKyAzM8tFR1VAkiYBlwNjIuLk0na+0ZlZJa582lPmFZCk4ZJ+JemuIZzjZkl9klaX2XaMpCclrZV0abXzRMQzETG7UtvMzBqnEbfg/hpYU26DpD0l7VrSN6XMrvOBY8ocPxy4FvgYcABwuqQDJL1P0l0lP3sO9S9iZmb1k2kCkjQO+HPgpgq7HAncIWlksv8c4OrSnSLiQWBjmeMPBdYmlcwbwELghIh4NCKOK/npq8ffyczM6iPrCugq4GKg7JtjEXErcA+wUNKZwNnAKQM4/97A80Xt3qSvLEljJV0PTJN0WWm7wjEzJd24efPmAYRlZma1ZDYJQdJxQF9ErJB0VKX9IuLrkhYC1wGTI2LLQC5T7pRVrrUBOK+ku7RdesxiYHF3d/ecAcRlZmY1ZFkBHQ4cL2kdhVtjH5b03dKdJB0BHAjcDnxxgNfoBfYpao8DXhhUtGZm1lCZJaCIuCwixkXEBOA04P6I+IvifSRNA74FnAB8Cthd0lcHcJllwFRJEyXtnFznzrr8BczMLFN5v4g6CvhkRDwdETuAWcBzpTtJWgAsAfaT1CtpNkBEbAcuAO6lMNNuUUQ81rDozVqUVxi1ZqCIio9MrEh3d3csX7487zDMhmz9hq0cfVUPr23bwdQ9C99WGz92VN5hWZuStCIiustt66gvIZi1miyqlCzX1+nnLxdko//3oV3GN+9bcGbWYF5fx5qFKyCzJpbF/9P1+jrWLFwBmXUYr69jzcIVkFmHGT92FE9c8bG8wzBzBWRm1grWb9jKyt5NLH12IzPm9rB+w9a8QxoyV0Bm1rFa6V2oRsxebDRXQGZmLaAdZy/6RdSU/CKqmeWp3OzF+y46Mueoaqv2IqorIDOzFtCOsxf9DMjMrAWMHzvqrWc+7fLelisgMzPLhROQmZnlwgnIzMxy4QRkZma58CQEM7MW0S6TD/q5AjIzs1w4AZmZWS6cgMzMLBdOQGZmlgsnIDMzy4UTkJmZ5cIJyMzMcuH3gMyGaP2GrRx9VQ+vbdvB1D0LXykeP3ZU3mGZNT0nIOs49V4FM8uVKtvtxUOzYjUTkKRu4AjgD4FXgdXAjyNiY8axmbWEdlyp0qwRKiYgSWcBnwGeBVYATwIjgQ8Cl0haDfx9RKxvQJxmdVPvqqLcSpWuXMxqq1YB7QIcHhGvltso6WBgKtAyCUjSJOByYExEnFzazjc6a1XzZk1n9v9dxjO/+z2TunZpi5UqzRpBEZHNiaWRwIPAH1BIdLdFxBcHea6bgeOAvog4sGTbMcA3gOHATRHxtRTnu6044ZS2y+nu7o7ly5cPJnwzs44laUVEdJfbVvUZUJJEjuOdz4D+PSIeq3Hd14EPR8QWSSOAhyT9MCJ+UXT+PYFXI+KVor4pEbG25FzzgWuAb5fENxy4FpgB9ALLJN1JIRldWXKOsyOir0bMZmbWINWeAX0JmAk8ACwF+ig8A9oX+FqSnD4XEavKHR+F0mpL0hyR/JSWW0cC50s6NiJekzQH+ARwbMm5HpQ0ocxlDgXWRsQzScwLgRMi4koKiXPIJM0EZk6ZMqUepzMzs0S1F1GXRcQhEfG5iPheRPw4Iu6KiLkRMRM4E9i52sklDZf0CIXkdV9ELC3eHhG3AvcACyWdCZwNnDKA+PcGni9q9yZ9leIZK+l6YJqky0rb5Y6JiMURcc6YMWMGEJaZmdVSsQKKiH8v7ZM0DBgdES8nt7Oq3tKKiDeBgyXtBtwu6cCIWF2yz9eTyuU6YHJEbClzqkpU7rJV4tkAnFfSXdo2M7MGqPkpHknfk/RuSbsAjwNPSvq7gVwkIjZRuJV3TJnzHwEcCNwODHSSQi+wT1F7HPDCAM9hZmY5SPMtuAMi4mXg48DdwHjgL2sdJKkrqXyQ9C7gI8ATJftMA74FnAB8Cthd0lcHEP8yYKqkiZJ2Bk4D7hzA8WZmlpM0CWhEMovt48AdEbGNKre5iuwF/FTSKgqJ4r6IuKtkn1HAJyPi6YjYAcwCnis9kaQFwBJgP0m9kmYDRMR24ALgXmANsCjF7DwzM2sCab4FdwOwDlgJPCjpj4CXax2UzI6bVmOfn5e0t1GoiEr3O73KOe6mUJmZmVkLqVkBRcTVEbF3RBybTK1eD3wo+9DMzKydDfhr2EkS2p5BLGZm1kG8IJ2ZmeXCCcjMzHJRNQEl7/9MLtP//uxCMjOzTlAxAUk6hcJ7O9+X9Jik4m/Mz886MDMza2/VKqDPA4dExMEUXhL9jqQTk23lPoFjZmaWWrVZcMMj4jcAEfFLSR8C7pI0jnQvopqZmVVUrQJ6pfj5T5KMjqLw2Zz3ZhyXmZm1uWoV0PmU3GqLiFeSFUgHsmSCmZnZO1RbjmFlcVvSu4v2/2GWQZmZWfur+SUESecCX6GwHHf/s58AJmUYl5mZtbk0n+L5W+C9EfFi1sGYmVnnSPMlhKeBrVkHYmZmnSVNBXQZ8B+SlgKv93dGxGcyi8rMzNpe2vWA7gceBXZkG46ZmXWKNAloe0RclHkkZmbWUdI8A/qppHMk7SVp9/6fzCMzM7O2lqYCOiP552VFfZ6GbWZmQ1IzAUXExEYEYmZmnaXmLThJfyVpt6L2/5D06Uyjso526g1LOPWGJXmHYWYZS/MMaE5EbOpvRMRLwJzMIrKOtn7DVlb2bmLpsxuZMbeH9Rv8CppZu0rzDGiYJEVEAEgaDuycbVjWCrKoUlb2buK1bYXZ/k/1beHoq3o4aNxudb3GLeceVtfzmdngpElA9wKLJF1PYfLBecA9mUZlHas/+VRqm1n7SJOALgHO4e3lGX4E3JRlUNYasqgkZszt4am+LQAME0zuGu2KxaxN1XwGFBE7IuL6iDg5Ik6KiBsi4s1GBGedZ96s6UzdczTDJSZ3jWberOl5h2RmGalYAUlaDNwI3BMR20q2TQLOAtZFxM2ZRmgdZfzYUdx30ZF5h2FmDVDtFtwc4CLgKkkbgd8BI4GJwFrgmoi4I/sQzcysHVVbEfW3wMXAxZImAHtRWJTu1xHhubFmZjYkaSYhEBHrgHWZRmJmZh0lzYuoZmZmdecEZGZmuaiYgCTdKOkTknZtZEBmZtYZqlVANwMHAXdL+omkSyQd1KC4zMyszVWbBfcL4BfAlySNBY4GPifpfcCvKLwftKgxYZqZWbtJOwtuA7Ag+UHSIcAxGcZlZmZtLlUCKhURK4AVdY7FzMw6iGfBmZlZLpyAzMwsFzVvwUk6sUz3ZuDRiOirf0hmZtYJ0jwDmg0cBvw0aR9FYXbcvpK+EhHfySg2MzNrY2luwe0A9k/WAjoJOAB4HfgAhcXqWoakSZLmSbqtXNvMzBonTQKaEBH/WdTuA/aNiI3AtgrHIGkfST+VtEbSY5L+erBBSrpZUp+k1WW2HSPpSUlrJV1a7TwR8UxEzK7UNjOzxklzC+5nku4Cbk3aJwEPStoF2FTluO3A5yLi4eRzPisk3RcRj/fvIGlP4NWIeKWob0pErC0513zgGuDbxZ2ShgPXAjOAXmCZpDuB4cCVJec428+szMyaR80EFBGflnQS8EFAFJLA9yMigA9VOe43wG+SP78iaQ2wN/B40W5HAudLOjYiXpM0B/gEcGzJuR5M1iQqdSiwNiKeAZC0EDghIq4Ejqv1dzMzs/xUTUCShgGrIuJA4PuDvUiSPKYBS4v7I+JWSROBhZJuBc6mUM2ktTfwfFG7l8KzqUpxjAX+EZgm6TIKS46/1U4SV+kxM4GZU6ZMGUBYZmZWS9UEFBE7JK2UND4i1g/mApJGU0hen42Il8tc4+tJ5XIdMDkitgzk9OXCrrRz8kmh80q6S9ulxywGFnd3d88ZQFxmZlZDmmdAewGPSfol8Pv+zog4vtaBkkZQSD7/FhE/qLDPEcCBwO3AF4ELUsTUrxfYp6g9DnhhAMebmVlO0iSgLw/mxJIEzAPWRMTcCvtMA74F/DnwLPBdSV+NiC+kvMwyYGpyG+//A6cBZwwmXjMza6ya07AjogdYB4xI/rwMeDjFuQ8H/hL4sKRHkp9jS/YZBXwyIp6OiB3ALOC50hNJWgAsAfaT1CtpdhLbdgoV073AGmBRRDyWIraGOPWGJZx6w5K8wzAza0ppPsUzBzgH2B2YTOHB//XAn1U7LiIeovwzmuJ9fl7S3kahIird7/Qq57gbuLvadczMrPmkeRH1ryhUMy8DRMRTwJ5ZBmVmZu0vTQJ6PSLe6G9I2okqM83MzMzSSJOAeiR9HniXpBkUvoiwONuwzMys3aVJQJcCvwMeBc6l8Lwl7Sw1MzOzstJ8imcHhYkB75gcYGZmNlgVKyBJiyXNTF4mLd02SdJXJJ2dbXhmZtauqlVAc4CLgKskbaRwG24kMBFYC1wTEXdkH6KZmbWjigkoIn4LXAxcnHxMdC/gVeDXEbG1MeGZmVm7SvMpHiJiHYWvIZiZmdVFmllw1uL8SSAza0ZOQGZmlotUCUjSuyTtl3UwZmbWOWomoGRF0EeAe5L2wZLuzDguq5P1G7aysncTS5/dyIy5Pazf4PkjZtYc0kxC+BJwKPAAQEQ8ksyKs4zU83nNyt5NvLZtBwBP9W3h6Kt6OGjcbnU59y3nHlaX85hZZ0pzC257RGzOPBLLRH/yqdQ2M8tLmgpotaQzgOGSpgKfAf4j27A6Wz0rixlze3iqbwsAwwSTu0a7cjGzppCmAroQeC/wOrCAwrpAn80wJqujebOmM3JE4X/myV2jmTdres4RmZkVpPkY6Vbg8uTHWsz4saPeeubjysfMmknFBCRpMVUWnouI4zOJyMzMOkK1Cuh/J/88EfifwHeT9un4szxmZjZE1T5G2gMg6YqI+NOiTYslPZh5ZFY3vvVmZs0ozSSELkmT+huSJgJd2YVkZmadIM007L8BHpD0TNKeQGFpbjMzs0FLMwvunuT9n/ckXU9ExOvZhmVmZu0u1XpAwCEUKp+dgIMkERHfziwqMzNrezUTkKTvAJMpfJD0zaQ7ACcgMzMbtDQVUDdwQERUfCfIzMxsoNLMgltN4T0gMzOzuklTAe0BPC7plxS+Bwf4SwhmZjY0adcDMjMzq6s007B7GhGImZl1lmofI32F6h8jfXcmEZmZWUeo9i24XQEkfQX4LfAdQMCZwK4Nic7MzNpWmllwH42Ib0bEKxHxckRcB5yUdWBmZtbe0iSgNyWdKWm4pGGSzuTtF1KtgvUbtrKydxNLn93IjLk9rN+wNe+QzMyaimq9XyppAvAN4HAKz4R+Dnw2ItZlHVwz2f2P9o8Zn7859f4rezfx2rYdb7VHjhj21sqktXj5BDNrF5JWRER3uW1pZsGtA06od1Dtrjj5lGubmXW6qglI0keBccCPI+K5ov6zIyJ9OdAGJnXtMqDKZMbcHp7q2wLAMMHkrtGubMzMilR8BiTpn4DLgfcB90u6sGjzBVkH1urmzZrOyBGF4Z3cNZp5s6bnHJGZWXOpVgHNBKZFxHZJXwK+J2lSRPwNhenYVsX4saPeeubjysfM7J2qzYLbKSK2A0TEJgoJ6d2SbgV2bkBsZmbWxqoloKclHdnfiIg3I2I28CSwf+aRmZlZW6uWgD4J/LK0MyK+AOyTWURmZtYRKiagiHg1Il5VwV9I+gcASeOBvRsWoZmZtaU0X0L4JnAYcHrSfgW4NrOIzMysI6RZD+gDEfHHkn4FEBEvSfIkBDMzG5I0FdA2ScNJlmaQ1AX4tX4zMxuSNAnoauB2YE9J/wg8BPxTplGZmVnbS/MtuH+TtAL4MwovoH48ItZkHpmZmbW1NM+AiIgngCcyjsXMzDpImltwZmZmdecEZGZmuXACMjOzXDgBmZlZLpyAzMwsF05AZmaWCycgMzPLhROQmZnlwgnIzMxy4QRkZma5cAIyM7NcOAGZmVkunIDMzCwXTkBmZpYLJyAzM8uFE5CZmeXCCcjMzHLhBGRmZrlwAjIzs1w4AZmZWS6cgMzMLBdOQGZmlgsnIDMzy4UTkJmZ5cIJyMzMcuEEZGZmuXACMjOzXDgBmZlZLpyAzMwsFx2RgCRNkjRP0m3V+szMrHGaPgFJullSn6TVJf3HSHpS0lpJl1Y7R0Q8ExGza/WZmVnj7JR3ACnMB64Bvt3fIWk4cC0wA+gFlkm6ExgOXFly/NkR0deYUM3MLK2mT0AR8aCkCSXdhwJrI+IZAEkLgRMi4krguHpdW9I5wDkA48ePr9dpzcyMFrgFV8HewPNF7d6kryxJYyVdD0yTdFmlvlIRcWNEdEdEd1dXVx3DNzOzpq+AKlCZvqi0c0RsAM6r1WdmZo3TqhVQL7BPUXsc8EJOsZiZ2SC0agJaBkyVNFHSzsBpwJ05x2RmZgPQ9AlI0gJgCbCfpF5JsyNiO3ABcC+wBlgUEY/lGaeZmQ1M0z8DiojTK/TfDdzd4HDMzKxOmr4CMjOz9tT0FVAru+Xcw/IOwcysabkCMjOzXDgBmZlZLpyAzMwsF05AZmaWCycgMzPLhRNQDZJmSrpx8+bNeYdiZtZWnIBqiIjFEXHOmDFj8g7FzKytOAGZmVkunIDMzCwXiqi4jI4VkbQZeKrC5jFAuYdE5frT9O0BvDiIMAejUuz1PjbNvtX2GcoYt+r4DvT4WvvWa3zL9bfqGLfK73C5vkaN8VB/h6dGRPlnGBHhnxQ/wI0D3VauP00fsLwZ/l71PDbNvlmNcauOb73HuF7j205j3Cq/w3mOcZa/w74Fl97iQWwr15+2r1GGcu2BHJtm36zGuFXHd6DH19q3XuNbrr9Vx7hVfofTXj8Lmf0O+xZcE5K0PCK6846jXXl8s+cxzl47jLEroOZ0Y94BtDmPb/Y8xtlr+TF2BWRmZrlwBWRmZrlwAjIzs1w4AZmZWS6cgFqApI9L+pakOyQdnXc87UbS/pKul3SbpPPzjqddSdpF0gpJx+UdS7uRdJSknyW/x0flHU9aTkA5kXSzpD5Jq0v6j5H0pKS1ki4FiIj/FxFzgLOAU3MIt+UMcHzXRMR5wClAS09rbaSBjHHiEmBRY6NsXQMc3wC2ACOB3kbHOlhOQPmZDxxT3CFpOHAt8DHgAOB0SQcU7fKFZLvVNp8BjK+k44GHgJ80NsyWNp+UYyzpI8DjwH82OsgWNp/0v8M/i4iPUUjyX25wnIPmBJSTiHgQ2FjSfSiwNiKeiYg3gIXACSr4Z+CHEfFwo2NtRQMZ32T/OyPiT4AzGxtp6xrgGH8I+F/AGcAcSf5vTw0DGd+I2JFsfwn4gwaGOSQ75R2A/Td7A88XtXuBDwAXAh8BxkiaEhHX5xFcGyg7vsk98xMp/It7d+PDaitlxzgiLgCQdBbwYtF/MG1gKv0Onwh8FNgNuCaHuAbFCai5qExfRMTVwNWNDqYNVRrfB4AHGhtK2yo7xm/9IWJ+40JpS5V+h38A/KDRwQyVy+Dm0gvsU9QeB7yQUyztyOObPY9xttpqfJ2AmssyYKqkiZJ2Bk4D7sw5pnbi8c2exzhbbTW+TkA5kbQAWALsJ6lX0uyI2A5cANwLrAEWRcRjecbZqjy+2fMYZ6sTxtcfIzUzs1y4AjIzs1w4AZmZWS6cgMzMLBdOQGZmlgsnIDMzy4UTkJmZ5cIJyGyIJL0p6ZGin0trH5W9orj+MGmvk/Szkn0eKf3cf5nzPCtpv5K+qyRdLOkISY/XOodZOf4WnNnQvRoRB9fzhJJ2Sl46HIpyce0qaZ+IeF7S/inPs5DCG/dfTmIbBpwMHB4Rz0k6FrhriLFaB3IFZJaRpOL4sqSHJT0q6T1J/y7JYmPLJP1K0glJ/1mSbpW0GPiRpFGSFklaJekWSUsldUuaLelfi64zR9LclGEt4u1FDU8HFhSdZ7ikf0niWiXp3GTTAgoJqN+fAusi4rlBDYxZwgnIbOjeVXILrnjV2hcj4o+B64C/TfouB+6PiOkU1sn5F0m7JNsOA2ZFxIeBTwMvRcT7gSuAQ5J9FgLHSxqRtD8F/J+Usd5GYekJgJnA4qJts4HNSVzTKazbMzEiVgE7JB2U7HcaRYnLbLB8C85s6Krdguv/RP4K3v4P/9EUEkh/QhoJjE/+fF9E9C9C9kHgGwARsVrSquTPv5d0P3CcpDXAiIh4NGWsG4GXJJ1G4VtiW4u2HQ28X9LJSXsMMBV4lqQKkvQYhQXm/iHl9cwqcgIyy9bryT/f5O1/3wScFBFPFu8o6QPA74u7qpz3JuDzwBOkr3763UJhWeezSvoFXBgR95Y5ZgHwI6AHWBURfQO8ptk7+BacWePdC1woSQCSplXY7yHglGSfA4D39W+IiKUU1oU5g4HfDrsd+HoSR2lc5/ff2pO0b/+twYh4GtgAfG0Q1zMrywnIbOhKnwF9rcb+VwAjgFXJ9OUrKuz3TaArufV2CbAK2Fy0fRHw84h4aSDBRsQrEfHPEfFGyaabgMeBh5O4buC/3yVZALyHQgIzGzIvx2DWpCQNp/B85zVJk4GfAPv2Jw5JdwH/GhE/qXD8logY3YA4JwB3RcSBWV/L2osrILPmNQp4SNJKClXH+RHxhqTdJP2awuSHsskn8XLxi6hZkHQEhZl0L2Z1DWtfroDMzCwXroDMzCwXTkBmZpYLJyAzM8uFE5CZmeXCCcjMzHLhBGRmZrn4L+IYUO1dGapPAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "flux_points.plot(sed_type=\"e2dnde\");"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Lightcurves\n",
    "\n",
    "The Fermi catalogs contain lightcurves for each source. It is available via the `source.lightcurve()` method as a `~gammapy.estimators.LightCurve` object."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 40,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:38.550988Z",
     "iopub.status.busy": "2021-11-22T21:09:38.538721Z",
     "iopub.status.idle": "2021-11-22T21:09:38.555595Z",
     "shell.execute_reply": "2021-11-22T21:09:38.555792Z"
    }
   },
   "outputs": [],
   "source": [
    "lightcurve = catalog_4fgl[\"4FGL J0349.8-2103\"].lightcurve()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 41,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:38.557809Z",
     "iopub.status.busy": "2021-11-22T21:09:38.557497Z",
     "iopub.status.idle": "2021-11-22T21:09:38.558935Z",
     "shell.execute_reply": "2021-11-22T21:09:38.559099Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<gammapy.estimators.points.core.FluxPoints at 0x128ec3910>"
      ]
     },
     "execution_count": 41,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "lightcurve"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 42,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:38.560767Z",
     "iopub.status.busy": "2021-11-22T21:09:38.560481Z",
     "iopub.status.idle": "2021-11-22T21:09:38.626852Z",
     "shell.execute_reply": "2021-11-22T21:09:38.627089Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div><i>Table length=10</i>\n",
       "<table id=\"table4981051552\" class=\"table-striped table-bordered table-condensed\">\n",
       "<thead><tr><th>time_min</th><th>time_max</th><th>e_ref [1]</th><th>e_min [1]</th><th>e_max [1]</th><th>flux [1]</th><th>flux_errp [1]</th><th>flux_errn [1]</th><th>flux_ul [1]</th><th>ts [1]</th><th>sqrt_ts [1]</th><th>is_ul [1]</th></tr></thead>\n",
       "<thead><tr><th></th><th></th><th>MeV</th><th>MeV</th><th>MeV</th><th>1 / (cm2 s)</th><th>1 / (cm2 s)</th><th>1 / (cm2 s)</th><th>1 / (cm2 s)</th><th></th><th></th><th></th></tr></thead>\n",
       "<thead><tr><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float32</th><th>float32</th><th>bool</th></tr></thead>\n",
       "<tr><td>54682.65603794185</td><td>55045.301668796295</td><td>3872.9833462074166</td><td>49.99999999999999</td><td>299999.9999999999</td><td>8.182966837466665e-08</td><td>4.0442960091979785e-09</td><td>4.0442960091979785e-09</td><td>8.991825950488419e-08</td><td>866.5454</td><td>29.437143</td><td>False</td></tr>\n",
       "<tr><td>55045.301668796295</td><td>55410.57944657407</td><td>3872.9833462074166</td><td>49.99999999999999</td><td>299999.9999999999</td><td>3.4780882174345606e-08</td><td>3.377941926174799e-09</td><td>3.377941926174799e-09</td><td>4.1536765138516785e-08</td><td>177.31786</td><td>13.316075</td><td>False</td></tr>\n",
       "<tr><td>55410.57944657407</td><td>55775.85722435185</td><td>3872.9833462074166</td><td>49.99999999999999</td><td>299999.9999999999</td><td>1.4417943283717705e-08</td><td>2.6934492414198985e-09</td><td>2.530577081216734e-09</td><td>1.980484221064671e-08</td><td>45.33303</td><td>6.7329807</td><td>False</td></tr>\n",
       "<tr><td>55775.85722435185</td><td>56141.13500212963</td><td>3872.9833462074166</td><td>49.99999999999999</td><td>299999.9999999999</td><td>1.8312720229118895e-08</td><td>2.9073525809053535e-09</td><td>2.697997159017973e-09</td><td>2.4127425390929602e-08</td><td>69.82312</td><td>8.356023</td><td>False</td></tr>\n",
       "<tr><td>56141.13500212963</td><td>56506.412779907405</td><td>3872.9833462074166</td><td>49.99999999999999</td><td>299999.9999999999</td><td>1.929730331085011e-08</td><td>3.0470421741313203e-09</td><td>2.900109263848094e-09</td><td>2.5391386770934332e-08</td><td>63.285576</td><td>7.955223</td><td>False</td></tr>\n",
       "<tr><td>56506.412779907405</td><td>56871.690557685186</td><td>3872.9833462074166</td><td>49.99999999999999</td><td>299999.9999999999</td><td>2.279863409171412e-08</td><td>3.0816984519788093e-09</td><td>2.9492861486346555e-09</td><td>2.896203099567174e-08</td><td>98.43889</td><td>9.921638</td><td>False</td></tr>\n",
       "<tr><td>56871.690557685186</td><td>57236.96833546296</td><td>3872.9833462074166</td><td>49.99999999999999</td><td>299999.9999999999</td><td>3.4101141466180707e-08</td><td>3.5666940512157908e-09</td><td>3.438366702468443e-09</td><td>4.123452868043387e-08</td><td>164.55066</td><td>12.82773</td><td>False</td></tr>\n",
       "<tr><td>57236.96833546296</td><td>57602.24611324074</td><td>3872.9833462074166</td><td>49.99999999999999</td><td>299999.9999999999</td><td>1.2460412257553344e-08</td><td>2.889906092207184e-09</td><td>2.7565447702215806e-09</td><td>1.824022533014613e-08</td><td>26.10482</td><td>5.1092877</td><td>False</td></tr>\n",
       "<tr><td>57602.24611324074</td><td>57967.523891018514</td><td>3872.9833462074166</td><td>49.99999999999999</td><td>299999.9999999999</td><td>1.5963367516746985e-08</td><td>2.667339016326764e-09</td><td>2.5512609802547104e-09</td><td>2.1298045993489723e-08</td><td>57.259422</td><td>7.5669956</td><td>False</td></tr>\n",
       "<tr><td>57967.523891018514</td><td>58332.801668796295</td><td>3872.9833462074166</td><td>49.99999999999999</td><td>299999.9999999999</td><td>2.1374123804207557e-08</td><td>2.9909645871128987e-09</td><td>2.8739282065259886e-09</td><td>2.7356053422522567e-08</td><td>84.4538</td><td>9.189875</td><td>False</td></tr>\n",
       "</table></div>"
      ],
      "text/plain": [
       "<Table length=10>\n",
       "     time_min           time_max          e_ref [1]          e_min [1]     ...      flux_ul [1]         ts [1]  sqrt_ts [1] is_ul [1]\n",
       "                                             MeV                MeV        ...      1 / (cm2 s)                                      \n",
       "     float64            float64            float64            float64      ...        float64          float32    float32      bool  \n",
       "------------------ ------------------ ------------------ ----------------- ... ---------------------- --------- ----------- ---------\n",
       " 54682.65603794185 55045.301668796295 3872.9833462074166 49.99999999999999 ...  8.991825950488419e-08  866.5454   29.437143     False\n",
       "55045.301668796295  55410.57944657407 3872.9833462074166 49.99999999999999 ... 4.1536765138516785e-08 177.31786   13.316075     False\n",
       " 55410.57944657407  55775.85722435185 3872.9833462074166 49.99999999999999 ...  1.980484221064671e-08  45.33303   6.7329807     False\n",
       " 55775.85722435185  56141.13500212963 3872.9833462074166 49.99999999999999 ... 2.4127425390929602e-08  69.82312    8.356023     False\n",
       " 56141.13500212963 56506.412779907405 3872.9833462074166 49.99999999999999 ... 2.5391386770934332e-08 63.285576    7.955223     False\n",
       "56506.412779907405 56871.690557685186 3872.9833462074166 49.99999999999999 ...  2.896203099567174e-08  98.43889    9.921638     False\n",
       "56871.690557685186  57236.96833546296 3872.9833462074166 49.99999999999999 ...  4.123452868043387e-08 164.55066    12.82773     False\n",
       " 57236.96833546296  57602.24611324074 3872.9833462074166 49.99999999999999 ...  1.824022533014613e-08  26.10482   5.1092877     False\n",
       " 57602.24611324074 57967.523891018514 3872.9833462074166 49.99999999999999 ... 2.1298045993489723e-08 57.259422   7.5669956     False\n",
       "57967.523891018514 58332.801668796295 3872.9833462074166 49.99999999999999 ... 2.7356053422522567e-08   84.4538    9.189875     False"
      ]
     },
     "execution_count": 42,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "lightcurve.to_table(format=\"lightcurve\", sed_type=\"flux\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 43,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:38.650763Z",
     "iopub.status.busy": "2021-11-22T21:09:38.650462Z",
     "iopub.status.idle": "2021-11-22T21:09:38.759996Z",
     "shell.execute_reply": "2021-11-22T21:09:38.760235Z"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZ8AAAE6CAYAAADJM1/JAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAstUlEQVR4nO3deZwU5bX/8c9hE3EBHTEqI7K6oA4oEMRERSOIu0YUyI0ieoO5GrlqNJrfTaImmpioBOKuQfAmueCKCwhqwhIVEgG3qJgoogiaAENAEYdl5vz+qJqhGWaGmab7qV6+79erXtNV1d3nOd09ffqpeqrK3B0REZGQmiXdABERKT4qPiIiEpyKj4iIBKfiIyIiwan4iIhIcCo+IiISXIukG5DL9tprL+/UqVPSzRARySsLFy5c5e7tG7qPik8DOnXqxIIFC5JuhohIXjGzj7Z3H212ExGR4FR8REQkOBUfEREJTvt86mBmpwOnd+vWLemmSA7btGkTy5Yto6KiIummiCSidevWlJaW0rJlyyY/1nRi0fr16dPHNeBA6rNkyRJ22203SkpKMLOkmyMSlLtTXl7O559/TufOnbdaZ2YL3b1PQ4/XZjeRNFVUVKjwSNEyM0pKStLu+av4iOyAphaeoffNY+h987LUGpGwduSHl4qPSB7r1KkThx9+OL169aJPn7q3crg7o0ePplu3bpSVlfHqq6/WrJsxYwYHHXQQ3bp145ZbbslImx566CG6d+9O9+7deeihh2qW33nnnXTr1g0zY9WqVXU+dvbs2ZgZ48ePr1n22muvYWbcdttt9cacPXs2/fv332rZ5s2b+cpXvsKnn37a6LZffPHF9OzZk7KyMoYMGcK6deualOOSJUvo168f3bt3Z+jQoWzcuLHRseuzcOFCDj/8cLp168bo0aOp3lUyceJE2rdvT69evejVqxe//e1v63y8mXH++efXzG/evJn27dtz2mmn1Rvziy++oKSkhLVr1261/KyzzuKRRx7Z4ZyA6IOpqe6pd+/enoTz7p3r5907N5HY0njvvPNOkx+T6ff2gAMO8JUrVzZ4n2nTpvngwYO9qqrK582b51/96lfd3X3z5s3epUsXX7x4sW/YsMHLysr87bffbnTs4447zpcsWbLVsvLycu/cubOXl5f76tWrvXPnzr569Wp3d3/11Vd9yZIlDbZ51qxZfvjhh/vAgQNrlv3gBz/wnj17+q233lpvWyorK720tHSr9kyfPt1POOGERufj7r527dqa21deeaX/4he/2OY+DeV47rnn+qRJk9zd/ZJLLvG777670bGvv/56nzBhwjbL+/bt63PnzvWqqiofPHiwP/vss+7uPmHCBL/sssu2+7y77LKL9+rVy9evX+/u7s8++6z37NnTTz311AYfN2zYMJ84cWLN/Jo1a7ykpMS/+OKLre5X1/8BsMC38/2qno9IIEvL1/PGsjX8dclqBo6Zw9Ly9UHiPvXUU1xwwQWYGUcddRRr1qzh008/5ZVXXqFbt2506dKFVq1aMWzYMJ566ikAFi9ezODBg+nduzfHHHMM7777bqNiPffccwwcOJA999yTPfbYg4EDBzJjxgwAjjjiCBpzuqqOHTtSUVHBv/71L9ydGTNmcPLJJ9esr6ttzZo149xzz+Xhhx+uud/kyZMZPnx4E14p2H333YHoR/mXX35Z52al+nJ0d2bOnMmQIUMAGDFiBE8++SQQ9SQuuugi+vbtyxFHHFHzOm/Pp59+ymeffUb//v0xMy644IKa52yKk08+mWnTpgEwadKkrV6X+to2fPhwJk+eXHO/KVOmMHjwYNq0adPk+HVR8ckxSX1BSXZU7+MZet88Bo2dQ8WmKgDeW7GOQWPn7PA+IDNj0KBB9O7dm/vvv7/O+yxfvpz999+/Zr60tJTly5fXuxxg1KhR3HHHHSxcuJDbbruNSy+9tFHtaeg5m2LIkCE8+uijzJ07lyOPPJKddtqpZl19bUv9stywYQPPPvss55xzTpNjjxw5kn322Yd3332Xyy+/fJv19eVYXl5Ou3btaNGixVbLAW6++WZOOOEE5s+fz6xZs7jmmmv44osvttuW5cuXU1pauk2sao8//njNJsKPP/643ucZNmwYkydPpqKigjfffJN+/frVrKuvbYMHD2bhwoWUl5cD6RXzhug4nyxK50vljWVrtvmC6lnarknP8fAl/bd/Jwmu+n2tbz4dL7/8Mvvttx8rVqxg4MCBHHzwwRx77LFb3cfrOJzCzOpdvm7dOubOncu5555bs3zDhg0ATJgwgXHjxgHw/vvvc8opp9CqVSs6d+7MlClT6n3OpjrvvPMYOnQo7777LsOHD2fu3LkADbatb9++rFu3jr///e8sWrSIo446ij322KPJsSdMmEBlZSWXX345Dz/8MCNHjtxqfVNfT4Dnn3+ep59+uma/VUVFBUuXLmXz5s01+2P++c9/0qpVK8aOHQvAn/70pwaf8/TTT2f48OHstNNO3HvvvYwYMYKZM2fWmVNZWRkffvghkyZN4pRTTtlqXX1tO+SQQzjjjDN47LHHOOecc3j99dcZNGhQva9bU6n45JhsfEFJclJ/CAwcM4f3VkQ7sJsZdG2/6w7/UNhvv/0A2HvvvTn77LN55ZVXtik+paWlW/0qXrZsGfvttx8bN26sc3lVVRXt2rXj9ddf3ybeyJEja76MBwwYwMSJE7falFZaWsrs2bO3es4BAwY0Oa999tmHli1b8sILLzBu3Lia4tNQ22DLL/xFixbV+Su9srKS3r17A3DGGWfw05/+tM7nad68OUOHDuXWW2/dpvjUl+Nee+3FmjVr2Lx5My1atKh5PSEqWI8//jgHHXTQNrGqc7nhhhvo1KkTF154Yc26jRs3smzZsq1iVT9nSUlJzfLvfOc7XHvttXXmUu2MM87g6quvZvbs2TW9me21bfjw4dx00024O2eeeWZaB5PWa3s7hYp5SmLAwYm3z/YDrp3qB1w71TtfN9VPvH128DZI4zR1wMFHq77wg370rB9wbfS+frTqi+0/qAHr1q3zzz77rOZ2//79ffr06dvcb+rUqVsNOOjbt6+7u2/atMk7d+7sH3zwQc2Ag7feesvd3fv37++PPPKIu7tXVVX566+/vs3z1jfgoFOnTr569WpfvXq1d+rUycvLy7e6z/YGHFTvCH/55Zd9ypQp7h7tjK8ecNBQ295++23v1q2bt2/f3tetW1f/i1eHqqoqf++992puf//73/fvf//729yvoRyHDBmy1YCDu+66y93df/jDH/pll13mVVVV7h4NvqitvgEHffr08Xnz5tUMOJg2bZq7u3/yySc193niiSe8X79+dea1yy67uLv7xx9/7GPHjnX3rV/nhtpWWVnpHTp08MMOO8xnzZpV5/OnO+Ag8S/4XJ6SKD6Z/oKS7El6tNvixYu9rKzMy8rKvEePHn7TTTfVrLvnnnv8nnvucffoi/TSSy/1Ll26+GGHHebz58+vud+0adO8e/fu3qVLl60e/8EHH/hJJ53kZWVlfsghh/iNN964Tfy6io+7+/jx471r167etWtXf/DBB2uWjxs3zjt06ODNmzf3fffd1y+++OJtHpv6pZgqtfhsr21lZWU+dOjQ+l62elVWVvrRRx/thx12mB966KH+rW99q2b02/z587dqb305Ll682Pv27etdu3b1IUOGeEVFhbu7r1+/3keNGlXz3PXlWFfxmT9/vh966KHepUuXrYrEdddd5z169PCysjIfMGCAL1q0qM68qotPqtTXeXttGz16tO+7775eWVlZ5/OnW3x0ep0GJHV6nep9Rdp3k9sWLVrEIYcc0qTH6L2VQlPX/0FjTq+jfT45SF9MhUvvrUhEQ63rYGanm9n9tY/uFRGRzFDxqYO7P+Puo9q2bZt0U0RECpKKj8gO0D5TKWY78vlX8RFJU+vWrSkvL1cBkqLkHl3Pp3Xr1mk9XgMORNJUWlrKsmXLWLlyZdJNEUlE9ZVM06HiI5Kmli1bbnMFRxFpHG12ExGR4FR8REQkOBUfEREJTsVHRESCU/EREZHgVHxERCQ4FR8REQlOxUdERIJT8RERkeBUfEREJDgVHxERCU7Fpw66mJyISHap+NRBF5MTEckuFR8REQlOxUdERIJT8RERkeBUfEREJDgVHxERCU7FR0REglPxERGR4FR8REQkOBUfEREJTsVHRESCU/EREZHgVHxERCQ4FR8REQlOxUdERIJT8RERkeBUfEREJDgVHxERCa5FyGBm1gzoCewHfAm87e7/CtkGERFJXpDiY2ZdgWuBE4H3gJVAa+BAM1sP3Ac85O5VIdqzPWZ2OnB6t27dkm6KiEhBMnfPfhCzScA9wIteK6CZ7Q18C/i3uz+U9cY0QZ8+fXzBggVJN0NEJK+Y2UJ379PQfYL0fNx9eAPrVgBjQ7RDRERyQ+h9PnsDX2PLPp+3gAW5srlNRETCCLXP53jgOmBP4DVgBdE+n7OArmb2GHC7u38Woj0iIpKsUD2fU4DvuPvS2ivMrAVwGjAQeDxQe0REJEGhis9t9Q2pdvfNwJOB2iEiIjkg1EGmb5jZC2Z2kZm1DRRTRERyVKji0wG4DTgG+IeZPWlmQ81s50DxRUQkhwQpPu5e6e7PuftIYH9gAtFggyVm9ocQbRARkdwR/Nxu7r4ReAdYBHwG9AjdBhERSVaw4mNmHc3sGjN7FZgKNAfOdPcjQrVBRERyQ6jjfOYS7fd5FBjl7jpnjYhIEQs11PqHwJ9rn9dNRESKU6gBB3Pc3c3sQDP7k5m9BWBmZWb2oxBtEBGR3BF6wMEDRL2gTQDu/iYwLHAbREQkYaGLTxt3f6XWss2B2yAiIgkLUnzMrGN8c1V8YTmPlw8BPg3RBhERyR2hBhw8CRwJfI/oqqUHm9lyYAnw7UBtEBGRHBGq+BiAuy8GTjSzXYBm7v55oPgiIpJDQhWfDmb2m9oLzQwAdx8dqB0iIpIDQhWfL4GFgWKJiEiOC1V8yt39oUCxREQkx4Uaar0xUBwREckDoc5wcFSIOCIikh+CX1JBRERExacOZna6md2/du3apJsiIlKQghYfM+tqZjvFtweY2WgzaxeyDY3h7s+4+6i2bdsm3ZSght43j6H3zUu6GSJSBEL3fB4HKs2sGzAe6Az8X+A2iEgC9ONGUoUuPlXuvhk4Gxjr7lcC+wZug4iIJCx08dlkZsOBEUSX0gZoGbgNIiKSsNDFZyTQH7jZ3ZeYWWfg94HbICIiCQt1hgMA3P0dYHTK/BLglpBtEBGR5GmotQCwtHw9byxbw1+XrGbgmDksLV+fdJNEpIAF7flIGOmMKHpj2RoqNlUB8N6KdQwaO4eepe2a9BwPX9K/yXFFpDip5yMANYWnvnmRHaGetdRm7p79IGbNgf8ESoEZ7v5yyrofuftNWW9EGvr06eMLFixIuhlBDBwzh/dWrAOgmUHX9rvywlXHJdwqyUU72rMGaN2yWZN71qDedb4ws4Xu3qeh+4Tq+dwHHAeUA78xszEp674ZqA3SgPEj+tK6ZfRx6Np+V8aP6Jtwi6SQqGcttYXa5/NVdy8DMLM7gbvN7AlgOPEltiVZHUva1PwS1a9LaUg6n4+6etb6nBW3UD2fVtU33H2zu48CXgdmArsGaoOIJEQ9a6ktVM9ngZkNdvcZ1Qvc/adm9glwT6A2iEhC1LOW2kJdTO7bqYUnZflv3V2n1xERKTJBj/OJR72dCnRKje3uY+p7jIiIFJ7QB5k+A1QAfwM03EVEpEiFLj6l1aPeRESkeIUuPtPNbJC7Px84rjSCdgSLSCihi89fgClm1gzYRHSMj7v77oHbISIiCQpdfG4nup7P3zzEeX1EJGeoZy2pQp9Y9D3gLRUeEZHiFrrn8ykw28ymAxuqF2qotYhIcQldfJbEUytSTrkjIiLFJfRltG8MGU9ERHJT0H0+ZvaCmbVLmd/DzJ4L2QYREUle6AEH7d19TfWMu/8b2DtwG0REJGGhi0+lmXWsnjGzAwCNfBMRKTKhBxz8D/CSmc2J548FRgVug4iIJCz0gIMZZnYkcBTR2Q2udPdVIdsgIiLJC1J8zKyTu38IEBebqbXWG9DB3ZeFaI+IiCQrVM/n1vh8bk8BC4GVQGugG3A88A3gekDFR0SkCAQpPu5+rpn1AP4DuAjYF1gPLAKeBW5294oQbRERkeQF2+fj7u8QDTgQEZEiF3qotYiIiIqPSLEZet88ht43L+lmSJELUnzMLPTxRDvEzE43s/vXrl2bdFNERApSqJ7PX8zsSTP7rpl1ChQzbe7+jLuPatu2bdJNEREpSKFGu/WJT6VzMjDWzDoALwHTgTnuvqHBJxARkYISbJ+Pu3/k7ve6+1nA0cAzwInAi2Y2LVQ7REQkeYnsi3H3TcDMeCLuCYmISJHIidFu7r486TaIiBSSXB/VmBPFR0REiouKj4iIBJd48TGz+5Nug4iIhBXqINM965lKgFNCtEFEYGn5et5Ytoa/LlnNwDFzWFq+PukmSZEKNdptJfAR0QXkqnk8v3egNogUlHR2Jr+xbA0Vm6oAeG/FOgaNnUPP0nZNeo6HL+nf5LgitYUqPh8A33D3pbVXmNnHgdogUvSqC0998yKhhCo+Y4E9gG2KD/CrQG0QKSjp9EAGjpnDeyvWAdDMoGv7XdWTkUQE2efj7ne5+xv1rLsjRBtEcklSx2CMH9GX1i2jf/uu7Xdl/Ii+wdsgAoF6Pmb2dXd/qYH1uwMd3f2tEO0RKVYdS9rU7ONRj0eSFGqo9TlmNtfMfmJmp5rZV83sWDO7yMx+B0wFdg7UFhGRgpYPoxpDndX6SjPbAxgCnAvsC3wJLALua6hXJCJSzAp1VGOwE4u6+7+BB+JJJCdU/2NrE5QUknwY1ZhXVxgVESk2hTqqMfHT64iISGblw6hG9XxERApMPoxqDNrzMbM2ZvZjM3sgnu9uZqeFbIOIiCQv9Ga3CcAGoLoULwNuCtwGERFJWOjNbl3dfaiZDQdw9y/NzLb3IBHJnFzdDCPFJXTPZ6OZ7Ux0RmvMrCtRT0ikaOTDAYAi2Ra653M9MAPY38z+AHwNuDBwG0QyJqkDAEE9GMlvQYuPu79gZq8CRxFdy+e/3X1VyDaIJC0fDgAUybZQJxY9staiT+O/Hc2so7u/GqIdIplWqAcAimRbqJ7P7fHf1kAf4A2ink8Z8Ffg64HaIZK48SP6MmjsHCo2VeXsAYAi2Rbqej7Hu/vxRJfSPtLd+7h7b+AI4P0QbRCpLakd/9UHAPbrvCcvXHUcHUvaBIkr4SV13aZ8EHrAwcHu/rfqGXd/y8x6BW6DFKBCPfOvSLpy/bMZuvgsMrPfAr8nGm79baLLKogEpx3/IskJXXxGAv8F/Hc8/2fgnsBtkAKkHf8i+SXoQabuXuHuv3b3s+Pp1+5eEbINItXy4cy/IoUqaM/HzL4G3AAckBrb3buEbIcI5MeZf0UKVejNbuOBK4GFQGXg2CIikiNCF5+17j49cEwREckxoYvPLDO7FXiClBOK6gwHUj1UWpu/RIpD6OLTL/7bJ2WZAycEboeIiCQo9IlFjw8ZT0REclOoE4te1dB6dx8Toh0iIpIbQvV8dov/HgT0BZ6O508nOtBUpKho31bhqz53YMWmKgaOmcP4EX11Hr8UQYqPu98IYGbPE51Y9PN4/gbg0RBtEBFJly4amHmhL6PdEdiYMr8R6BS4DSIiWadzBzYs9Gi33wGvmNkUolFuZwMPBW6DiEiT6NyBmRd6tNvNZjYdOCZeNNLdXwvZBpFU+jKQbNFFAxsWuudTfUCpDioVkYKmcwc2LPQ+HxGR4HRF0dyj4iMiIsGp+IiISHAqPpK46oPx/rpkNQPHzGFp+fqkmyQiWRZ8wIEUrnS3qWfiYDzt0BXJL+r5SOJ0MJ5I8VHPRzIm3d6HDsYTKT7q+Ujixo/oS+uW0UdRB+OJFAf1fCRxOhhPpPio5yMiIsGp+IiISHAqPiIiElxR7fMxs47AncAq4B/ufkvCTRKRAqZ9mPXLm56PmT1oZivM7K1ayweb2d/N7H0zu247T3MgMM3dLwJ6ZK2xIiLSoLwpPsBEYHDqAjNrDtwFnExUTIabWQ8zO9zMptaa9gZeA4aZ2UxgVuD2i4hILG82u7n7n82sU63FXwXed/cPAMxsMnCmu/8COK32c5jZ1cD18XM9BkzIcrNFRKQO+dTzqUsH4OOU+WXxsvrMAEab2b3Ah3XdwcxGmdkCM1uwcuXKjDVURES2yPfiY3Us8/ru7O5vufsQd/+uu19dz33ud/c+7t6nffv2GWuoiCRDZ03PTXmz2a0ey4D9U+ZLgU8SaouIZFk6Z07XWdNzU773fOYD3c2ss5m1AoYBTyfcJhHJITprem7Km56PmU0CBgB7mdkyooED483se8BzQHPgQXd/O8FmikgWpdMD0VnTc1PeFB93H17P8meBZwM3R0TyxPgRfRk0dg4Vm6p01vQckjfFR0QkHTprem5S8ZGcoC8FkeKS7wMOREQkD6n41MHMTjez+9euXZt0U0RECpKKTx3c/Rl3H9W2bdukmyIiUpBUfEREJDgVHxERCU7FR0REglPxERGR4FR8REQkOBUfEREJTsVHRESCU/EREZHgVHxERCQ4FZ866PQ6IiLZpeJTB51eR0Qku1R8REQkOBUfEREJTsVHRESC05VMRaTg6Uq5uUc9HxERCU7FR0REglPxERGR4FR8REQkOBUfEREJTsVHRESCU/Gpg87tJiKSXSo+ddC53UREskvFR0REgjN3T7oNOcvMVgIfBQy5F7AqYLxciF1scZOMrZyLI3Yu5HyAu7dv6I4qPjnEzBa4e59iil1scZOMrZyLI3a+5KzNbiIiEpyKj4iIBKfik1vuL8LYxRY3ydjKuThi50XO2ucjIiLBqecjIiLBqfiIiEhwKj45xMys2GIr58KPm2Rs5Zy7sVV8EmRm/c1sjJkNAfCAO+CSiq2cCz9ukrGVc/7krOKTEDP7BvAAsBT4rpndamZ7FXJs5Vz4cZOMrZzzK2cVn+T0BP7o7mOBC4EOwKlmtksBx1bOhR83ydjKOY9yVvEJxMxONrNzzawkXrQYqDKzEndfBjwN9AO6F0ps5Vz4cZOMrZzDxc1GbBWfLLLITmY2EfgpMBi4w8z6AB8DrYGD4rs/CuwC9Igfu0PvTVKxlXPhx00ytnIunJxVfLIo3vlmQAtgsLtfDMwCHnD3V4H1wNFm1tHdK4GXgfPjx1blY2zlXPhxk4ytnAsnZxWf7OsMtAM2mVkzd38gvn0JcCfQEbgqvu9ewMwCiK2cCz9ukrGVcyHk7O6aMjwRFfVmKfMvAqNT5o8E3gNaAfsCv4vfsFeArvkYWzkXflzlrJwzGTvtxmna6o3qCowFhgGtU5Z/Jf7bH1gBtExZ9wfguPj2TkCnfIqtnAs/rnJWztmMrc1uO8jMyoh2tK0BTgEmmtl+ZrYTcJOZ9XD3ecCfiHbUtTWzNkBb4F0Ad9/g7h/mS2zlXPhxlbNyznrsdKqkpq1+MXwTmBLfbgFMAK4k6o5ayv32INo+OhlYBNxH9GvB8i22ci78uMpZOWc953QbXKwT0Tj2kUBJPH8kMJ64ywkcH88fW8djjejArCPyKbZyLvy4ylk5h4idOmmzWyPF491vBn4LDADGmtnZwD+BSuKx7u4+C/iM6A3GzNqY2RVmdqBH3nD31/IhtnIu/LjKWTmHiF2nHa1exTQB9xBXfOAkonMatQJ+QHQA1iHxuuOAv6U87ox8ja2cCz+uclbOIWLXntTzaSQz25vo3EXrANz9OWAu8CvgbmBX4L/NrCVQCvwxvo27P52PsZVz4cdNMrZyLo6c65XpalZoE1uPdf8dcGfK/B5Evxz2B3YHxgHPA+8AR+drbOVc+HGVs3IOEbvBdmXzyfNtIjqKdwzR6SFKU5bvTjQKZG+i4YhdU9bdCQxLme+cT7GVc+HHVc7KOUTspk7a7BYzs/2BJ4GWwKHAJDPrGq++Efiau68geqPGmdl+8boOxGPdAdx9Sb7EVs6FHzfJ2Mq5OHJOW4gKlw8T0Ad4KWX+jvgNawc0r3XfMcBE4E1gClDCjo2zTyS2ci78uMpZOYeInVZ7QwbLpQk4ELgM+ArR2PV9iLaH9ozXHwo8CJxOyjbTeJ0BXYBj8im2ci78uMpZOYeInYmpKDe7mdmlwFSi8xX9Avge8G/gC+BAM2vh7m8D/wAGuHuVmbUws1Fm1hfA3T9w9xfzJbZyLvy4ylk5h4idMUlVvSQn4Cbgm/HtnsByotEe3wZ+DfSO13UAPgJ2jedPA3bOx9jKufDjKmflHCJ2pqai6/lYdLK8zkRH8OLubxCdSuI+d/89sAG4xMw6AJ2APwMe33equ3+Zb7GVc+HHTTK2ci6OnDMu6eoXciLeoQZcD0yvte5d4GtEwxH/B3iO6HoV5+RzbOVc+HGVs3IOETvTU+INyEpS0c60kUCHeL5Z6hsX314OnJQy/yPgqpT5A/IptnIu/LjKWTmHiB1qSrwBGU8IegGLgQ+AX9exvkX89wJgPtA9nn8QOC0fYyvnwo+rnJVziNghp0Lc57MC+A1wFtDRzI4FMLNmAO6+Of77v8B04AdmtoBo2OGiPI2tnAs/bpKxlXNx5BxW0tVvRyei00X8iGjER9t4WWugDXAF8Ls6HlPdhW0O7AUMzKfYyrnw4ypn5RwidpJTXvd8zOx4YBbQHfgvotOF4+4V7r6e6FdBczO7MPVx7l4V36xy91Xu/kK+xFbOhR83ydjKuThyzglJV78dmYBvAbfHt5sRjewYnrJ+5/g+k+P5LsAu+RxbORd+XOWsnEPETnrK654PUXdzpZm19uiXwNXAVbblGhRfApOAVmb2OTCW6KR7+RxbORd+3CRjK+fiyDl5SVe/NH8tVI91PxmYCeyWsm4mcG3KL4k7gGXApfkcWzkXflzlrJxDxM6VKfEGNOJNOpI6upls2eH2LNElYKvfzJOAR9kyHPG8uh6fy7GVc+HHVc7KOUTsXJ4Sb0ADb9hhwBKibaDHpyy36jctnj8EeBEYGs9fAfwoH2Mr58KPq5yVc4jY+TDl5D4fM2sB7AfcCjwCHBWfpwiPVJlZFzObSnRVvl8CR5vZy8BFRG9kXsVWzoUfVzkr5xCx80bS1a96ItqJdgqwZzy/W/y3DPgDcAZbuqlHAk8BN6Q83oBj8ym2ci78uMpZOYeInY9T4g2IX/STiHaoTQUeIz5VeMr6q4HbgcPi+Z1J2QZKrav05UNs5Vz4cZWzcg4RO1+nXNnsdgRwtbufBjwBnGlm30xZ/3ugPdG2UYDW7v6FmbUyM3P3yjyMrZwLP26SsZVzceSct3Kl+Hyd6I0BmAa8AJxv0XUrcPd/El0e9hwz+yvRG4m7b/T4Z0MexlbOhR83ydjKuThyzluJFh8zax7fvI/oTWnu7muJLn70CTAk5e6DiLaZTnf3U/M1tnIu/LhJxg4dNyVe0eScK7HzXbDiY2bHmNmhqcvcvdLMDJhHdPrw78Wr1gCriU6ah5mVAP8EDnL3G9KIfZyZlYWObWZ7x39bhIwbP/5rZnZw6rJAOX/DzHpbfAbegHHPMrMh8e3moeLGj7/AzEalLguU80gzG2NmnULGjR9/nplNBwaEjG1mZ8RTdS8jZM6JfYcVpIZ2CGViIrqG+IvAH4mGHI4E9onX7ZRyv+OA14AD4/mx7OARvUTd4Llx/D8Rnbgv67GBHsAbwCcpy6oPIGuR5ZyPIDpCehYwB/gOW86U2zKLOfcmOgniHKJt3lew5brx2c65lOg69e+yZYRR8wBxewLPxK/111Pf62zFjp97J+DO+LN9AtAq23FTnrNf/P6+SHTNmX4hXu/4PZ4dTxOILjmwV6DPV2LfYYU8hej5fAP4s7ufSDSW/TDgUgB332Bme5nZFUQf5EeBMWY2BTgRWLiDsbsCr7v7McBPiD5E12QzdvyLfwjwAPC+mf00XlVzLQ4zK8lGzma2C/Bj4El3Px4YRzSks20ce1M2YptZK+CHwFPufhzwv0QXxNoYx81azrF1RNvQ5wO3xMuqsvxatyR6fVe4+/Hu/lJ1D9cjWYnt8bcasCtwqrvPTF0Xx83K/5SZtQNuA2bE/1NPABfGsSuz/D5/A3jJ3QcQnW6mJVHPIsTnK8nvsMKVjYoG7E80mgPgp8DDKevGAy8R7aA7GFgJ/DhlfTdgBPGvuTRip/66vxCYWb2cqFfwJNGvxQMzGbtW3APivz2AtUBJPN8c6E/U/c5Wzj2AnVPmFwE949tHZzJ2rbitU27fSnQMQ0+ia5KUZeO1ZssxE4cQ9UA6Aq+z5dLDRwCrsvhaX0B09civAFfFn/VvACXA4Vn8fB0BTIlvXxu/1t8D+hL94MpmzqmXcT4HuJstPdxsfr6uA96Kb99M9EPjFKKzPB+c5Zx/QsDvsGKZMvtkUbfzI2AK0YiP5sChRDvfRgLHEx1sdXP8j2rA/imPtx2IPYRoc9N3Upa1AV4mPrUFsBvRprcxcezSHY1dV9xa6x8EHkuZbwZ0zHDOo+pY1zLOdxrQJZOxG8qZ6HiGGcDlwG+JekRGXBAyHbf6fQRuTom/IP5StCy81t+ptfxJ4B9EI5muIPrl++Ns5hwvf4to89OtRD+mfhz/nxmwbzY/X8THpADDgScy/dluIOcHiL70X4xjjwMeinPeL1vvM9F32EtEP2az9h1WbFPGNrvFQwpHAle4+9nA+vgfYhXwM6LNPz8k+od5C2o2E3xsZs3MzDx+95oQ0+K/hxCdmG8J0MvMDozvspHodOT/Gcf7nGgESguiA7yWpRN7e3Hj56x+bS8DTkgZ8LCbuy81s5YZzLlnauw4101E+7xau/sH8cNbxbFbZCHn6kEVd7r7YHe/g+iLuAfRF9LybLzW8XPtC+xpZl2IDvbrDqyOP19Ls/j5gugKlNe4+/nuPpao4HYl+mGTlZxjNxKdcHKyR5vebgY+A05y90/NrHmWPl8GVMV3fxroHd8XgHRf70bmfBXwD3c/xt0nEe33agP0cfdPsvA+HxTn9DbR5tyMf4cVsx0qPvG2fiDa9gm0Y8u1Jq4BDgAGEm36utzdB7n7H4mKwi4pj61K4wPTqvox7r4IGAb8nGiUybnx8s1EmyR2NbNr4oeuIPpluC6d2I2MW+XRuZtaenQ9ju8Dk8zsV8DlFg3H3JSlnKtSHnIU0U5pzOyXwHAza+bum7OQ8+b47htSHrqWaGfwR9Vty3Tc2Io415eIenoXAzUH+GXr8xWve8vdn0p56OdEO6E/Tid2E+I+CvwdOCte1J7oB9+r8frKLH2+3N09/sLeSNTzO7w61yzn/Dmwi5mNjBe1IPq+WZRO3EbErhkm7e5TgSsz+R1W7NIuPmZ2OfCKmf3MzM6IF78O7GZmbdz9Q6Ivg6OItoljZrvGO+Z+RdS1zVhsd//A3RcDfwEOMLNj4+UfE/1K/LaZ3QVMBl6yWBbiHhPft5lHvQ+AVkT7JUqAX3saRzM3MXZ1L6QzcK6Z/QXYnehXclUdT5+puNW/IHczs+8B9wN/yeJrfWx89wqi97i7R72Px4FxFvcumxI3jZybxX93jnO+l2jYLVnK+ZiUhwwFysxsHPA80S/28hA5x1+0lUSFp1V9z5vhuC2JetM3m9nPgIeJRpWuD5RzZaa+w4T09vkQ7dv5K1E3dBjRiI4yYDBwF9A3vt8uwCtE3WKIfkn8gXgHeIZiLwC+lrJ+H6Kdk2NqPa4jcHa6sZsQ9/ZajzuGaAdliJxrx55A9E/UPVRcol+i1ft9Dgv5HsfrWqQTMwOv9Wiio9qznXPtuLsT7fPpETJntlxr5kyga8j3mWgU2VUJ/U+dDfzfjsTWFE3Vx55sV7ypqDK+PQQoc/efxPOXARe4ez8zGwssBx5x94/M7B5gtrs/3KhATY/9XeC/3L1nyv2/Gn9IVhH1Nm539/IE4t7i7p8llPMvgUqPNy8Girsn0S/CDe5eETjumHTe4wzELgF+AVT4lp5uqLhpfa4zEHtPoh78qsBx9wJ+6e5r0kg5E7F/7tGmP8mA7W52s2jn9M+Bn5vZSfHizUD1Jg/c/S6gjZmdRzT6pj0w3sxuI/qVMj+dxjUy9r3xfVOPMH+P6JfN/8T3adI/aAbjNrnwZCi2ufvaphSeDMdtdOHJYNx0flxk6n3+vCmFJ6nPdQZjW1MLT4biejqFJ4OxVXgyqaFuEdEL/wZwD9FO3PnER3IT7ew8P+W+pwHTUuZHEI0E+ko6XbImxj6V6HxJ1fMTiXaEluRL3GLMWa+1ci7UnDU14r3Zzht3TK03Zxxwa3z7dKJjeqoPJj2S6PiZjBxYlWbs6vmd8y1uMeas11o5F2rOmrY/NbjPx8zaEI1o2ezRSI/hREdXX+fRUOIJREMO/0g0JPIzd//Pep+wCZKKrZwLP26SsZVzceQs29fgPh93X+/uG3zL0OCTgI99y3DdK4gONBsK/D2Tb1pSsZVz4cdNMrZyDhc36diyHY3pHhGdJqcZ0VmLu8bLDmPLOZ2ydg6jpGIr58KPq5yVc4jYmuqeGnuQaRXR8RuriE618QzRUfvVZ/Hd2MjnSUdSsZVz4cdNMrZyLo6cpT5N+OVwFNEb+BJwccgKmVRs5Vz4cZWzctaUzNSUg0xLgfOJDubbsL37Z1JSsZVz4cdNMrZyLo6cpW6NLj4iIiKZEuJKpiIiIltR8RERkeBUfEREJDgVHxERCU7FR0REglPxEckCMysxs9fj6Z9mtjy+vc7M7s5CvIlmtiS+Lg1m9l0zuyCN5xlqZu+b2dRMt1EklYZai2SZmd0ArHP327IYYyIw1d0fy8BzDQCudvfTdvS5ROqjno9IQGY2oLpXYWY3mNlDZva8mX1oZt80s1+Z2d/MbIaZtYzv19vM5pjZQjN7zsz2bUScG8zs6vj2aDN7x8zeNLPJ8bI9zezJeNlfzKwsm3mL1KbiI5KsrkQXMTsT+D0wy90PB74ETo0L0B3AEHfvDTwI3NzEGNcBR7h7GfDdeNmNwGvxsv8H/O8OZyLSBC2SboBIkZvu7pvM7G9EZ16eES//G9AJOIjo7MsvmBnxfT5tYow3gT+Y2ZNEV+YE+DpwDoC7z4z3UbV197XppyLSeCo+IsnaAODRhc02+ZadsFVE/58GvO3u/XcgxqnAscAZwI/N7ND4eWvTDmAJRpvdRHLb34H2ZtYfwMxaxsWjUcysGbC/u88CfgC0A3YF/gz8R3yfAcAqd/8soy0XaYB6PiI5zN03mtkQ4Ddm1pbof3Ys8HYjn6I58Pv4sQb82t3XxCPwJpjZm8B6YETGGy/SAA21FikAGmot+Uab3UQKw1rgZ9UHmabLzIYCdwP/zkirROqhno+IiASnno+IiASn4iMiIsGp+IiISHAqPiIiEpyKj4iIBPf/AQjajijY39dQAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "lightcurve.plot();"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Pretty-print source information\n",
    "\n",
    "A source object has a nice string representation that you can print.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 44,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:38.764730Z",
     "iopub.status.busy": "2021-11-22T21:09:38.763129Z",
     "iopub.status.idle": "2021-11-22T21:09:38.766167Z",
     "shell.execute_reply": "2021-11-22T21:09:38.766357Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "*** Basic info ***\n",
      "\n",
      "Catalog row index (zero-based) : 18\n",
      "Source name          : HESS J1514-591\n",
      "Analysis reference   : HGPS\n",
      "Source class         : PWN\n",
      "Identified object    : MSH 15-52\n",
      "Gamma-Cat id         : 79\n",
      "\n",
      "\n",
      "*** Info from map analysis ***\n",
      "\n",
      "RA                   :  228.499 deg = 15h14m00s\n",
      "DEC                  :  -59.161 deg = -59d09m41s\n",
      "GLON                 :  320.315 +/- 0.008 deg\n",
      "GLAT                 :   -1.188 +/- 0.007 deg\n",
      "Position Error (68%) : 0.020 deg\n",
      "Position Error (95%) : 0.033 deg\n",
      "ROI number           : 13\n",
      "Spatial model        : 3-Gaussian\n",
      "Spatial components   : HGPSC 023, HGPSC 024, HGPSC 025\n",
      "TS                   : 1763.4\n",
      "sqrt(TS)             : 42.0\n",
      "Size                 : 0.145 +/- 0.026 (UL: 0.000) deg\n",
      "R70                  : 0.215 deg\n",
      "RSpec                : 0.215 deg\n",
      "Total model excess   : 3502.8\n",
      "Excess in RSpec      : 2440.5\n",
      "Model Excess in RSpec : 2414.3\n",
      "Background in RSpec  : 1052.5\n",
      "Livetime             : 41.4 hours\n",
      "Energy threshold     : 0.61 TeV\n",
      "Source flux (>1 TeV) : (6.434 +/- 0.211) x 10^-12 cm^-2 s^-1 = (28.47 +/- 0.94) % Crab\n",
      "\n",
      "Fluxes in RSpec (> 1 TeV):\n",
      "Map measurement                : 4.552 x 10^-12 cm^-2 s^-1 = 20.14 % Crab\n",
      "Source model                   : 4.505 x 10^-12 cm^-2 s^-1 = 19.94 % Crab\n",
      "Other component model          : 0.000 x 10^-12 cm^-2 s^-1 =  0.00 % Crab\n",
      "Large scale component model    : 0.000 x 10^-12 cm^-2 s^-1 =  0.00 % Crab\n",
      "Total model                    : 4.505 x 10^-12 cm^-2 s^-1 = 19.94 % Crab\n",
      "Containment in RSpec                :  70.0 %\n",
      "Contamination in RSpec              :   0.0 %\n",
      "Flux correction (RSpec -> Total)    : 142.8 %\n",
      "Flux correction (Total -> RSpec)    :  70.0 %\n",
      "\n",
      "*** Info from spectral analysis ***\n",
      "\n",
      "Livetime             : 13.7 hours\n",
      "Energy range:        : 0.38 to 61.90 TeV\n",
      "Background           : 1825.9\n",
      "Excess               : 2061.1\n",
      "Spectral model       : ECPL\n",
      "TS ECPL over PL      : 10.2\n",
      "Best-fit model flux(> 1 TeV) : (5.720 +/- 0.417) x 10^-12 cm^-2 s^-1  = (25.31 +/- 1.84) % Crab\n",
      "Best-fit model energy flux(1 to 10 TeV) : (20.779 +/- 1.878) x 10^-12 erg cm^-2 s^-1\n",
      "Pivot energy         : 1.54 TeV\n",
      "Flux at pivot energy : (2.579 +/- 0.083) x 10^-12 cm^-2 s^-1 TeV^-1  = (11.41 +/- 0.24) % Crab\n",
      "PL   Flux(> 1 TeV)   : (5.437 +/- 0.186) x 10^-12 cm^-2 s^-1  = (24.06 +/- 0.82) % Crab\n",
      "PL   Flux(@ 1 TeV)   : (6.868 +/- 0.241) x 10^-12 cm^-2 s^-1 TeV^-1  = (30.39 +/- 0.69) % Crab\n",
      "PL   Index           : 2.26 +/- 0.03\n",
      "ECPL   Flux(@ 1 TeV) : (6.860 +/- 0.252) x 10^-12 cm^-2 s^-1 TeV^-1  = (30.35 +/- 0.72) % Crab\n",
      "ECPL   Flux(> 1 TeV) : (5.720 +/- 0.417) x 10^-12 cm^-2 s^-1  = (25.31 +/- 1.84) % Crab\n",
      "ECPL Index           : 2.05 +/- 0.06\n",
      "ECPL Lambda          : 0.052 +/- 0.014 TeV^-1\n",
      "ECPL E_cut           : 19.20 +/- 5.01 TeV\n",
      "\n",
      "*** Flux points info ***\n",
      "\n",
      "Number of flux points: 6\n",
      "Flux points table: \n",
      "\n",
      "e_ref  e_min  e_max        dnde         dnde_errn       dnde_errp        dnde_ul     is_ul\n",
      " TeV    TeV    TeV   1 / (cm2 s TeV) 1 / (cm2 s TeV) 1 / (cm2 s TeV) 1 / (cm2 s TeV)      \n",
      "------ ------ ------ --------------- --------------- --------------- --------------- -----\n",
      " 0.562  0.383  0.825       2.439e-11       1.419e-12       1.509e-12       2.732e-11 False\n",
      " 1.212  0.825  1.778       4.439e-12       2.489e-13       2.654e-13       4.970e-12 False\n",
      " 2.738  1.778  4.217       7.295e-13       4.788e-14       4.898e-14       8.302e-13 False\n",
      " 6.190  4.217  9.085       1.305e-13       1.220e-14       1.282e-14       1.571e-13 False\n",
      "13.991  9.085 21.544       1.994e-14       2.723e-15       2.858e-15       2.588e-14 False\n",
      "31.623 21.544 46.416       9.474e-16       3.480e-16       4.329e-16       1.919e-15 False\n",
      "\n",
      "*** Gaussian component info ***\n",
      "\n",
      "Number of components: 3\n",
      "Spatial components   : HGPSC 023, HGPSC 024, HGPSC 025\n",
      "\n",
      "Component HGPSC 023:\n",
      "GLON                 :  320.303 +/- 0.005 deg\n",
      "GLAT                 :   -1.124 +/- 0.007 deg\n",
      "Size                 : 0.057 +/- 0.005 deg\n",
      "Flux (>1 TeV)        : (2.01 +/- 0.23) x 10^-12 cm^-2 s^-1 = (8.9 +/- 1.0) % Crab\n",
      "\n",
      "Component HGPSC 024:\n",
      "GLON                 :  320.298 +/- 0.020 deg\n",
      "GLAT                 :   -1.168 +/- 0.021 deg\n",
      "Size                 : 0.206 +/- 0.020 deg\n",
      "Flux (>1 TeV)        : (2.54 +/- 0.29) x 10^-12 cm^-2 s^-1 = (11.2 +/- 1.3) % Crab\n",
      "\n",
      "Component HGPSC 025:\n",
      "GLON                 :  320.351 +/- 0.005 deg\n",
      "GLAT                 :   -1.284 +/- 0.007 deg\n",
      "Size                 : 0.055 +/- 0.005 deg\n",
      "Flux (>1 TeV)        : (1.88 +/- 0.22) x 10^-12 cm^-2 s^-1 = (8.3 +/- 1.0) % Crab\n",
      "\n",
      "\n",
      "*** Source associations info ***\n",
      "\n",
      "  Source_Name    Association_Catalog    Association_Name   Separation\n",
      "                                                              deg    \n",
      "---------------- ------------------- --------------------- ----------\n",
      "  HESS J1514-591                2FHL    2FHL J1514.0-5915e   0.098903\n",
      "  HESS J1514-591                3FGL     3FGL J1513.9-5908   0.026914\n",
      "  HESS J1514-591                3FGL    3FGL J1514.0-5915e   0.094834\n",
      "  HESS J1514-591                COMP            G320.4-1.2   0.070483\n",
      "  HESS J1514-591                 PSR              B1509-58   0.026891\n",
      "\n",
      "*** Source identification info ***\n",
      "\n",
      "Source_Name: HESS J1514-591\n",
      "Identified_Object: MSH 15-52\n",
      "Class: PWN\n",
      "Evidence: Morphology\n",
      "Reference: 2005A%26A...435L..17A\n",
      "Distance_Reference: SNRCat\n",
      "Distance: 5.199999809265137 kpc\n",
      "Distance_Min: 3.799999952316284 kpc\n",
      "Distance_Max: 6.599999904632568 kpc\n",
      "\n"
     ]
    }
   ],
   "source": [
    "source = catalog_hgps[\"MSH 15-52\"]\n",
    "print(source)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "You can also call `source.info()` instead and pass as an option what information to print. The options available depend on the catalog, you can learn about them using `help()`"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 45,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:38.768331Z",
     "iopub.status.busy": "2021-11-22T21:09:38.768010Z",
     "iopub.status.idle": "2021-11-22T21:09:38.769270Z",
     "shell.execute_reply": "2021-11-22T21:09:38.769466Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Help on method info in module gammapy.catalog.hess:\n",
      "\n",
      "info(info='all') method of gammapy.catalog.hess.SourceCatalogObjectHGPS instance\n",
      "    Info string.\n",
      "    \n",
      "    Parameters\n",
      "    ----------\n",
      "    info : {'all', 'basic', 'map', 'spec', 'flux_points', 'components', 'associations', 'id'}\n",
      "        Comma separated list of options\n",
      "\n"
     ]
    }
   ],
   "source": [
    "help(source.info)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 46,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-11-22T21:09:38.771362Z",
     "iopub.status.busy": "2021-11-22T21:09:38.771029Z",
     "iopub.status.idle": "2021-11-22T21:09:38.772335Z",
     "shell.execute_reply": "2021-11-22T21:09:38.772511Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "*** Source associations info ***\n",
      "\n",
      "  Source_Name    Association_Catalog    Association_Name   Separation\n",
      "                                                              deg    \n",
      "---------------- ------------------- --------------------- ----------\n",
      "  HESS J1514-591                2FHL    2FHL J1514.0-5915e   0.098903\n",
      "  HESS J1514-591                3FGL     3FGL J1513.9-5908   0.026914\n",
      "  HESS J1514-591                3FGL    3FGL J1514.0-5915e   0.094834\n",
      "  HESS J1514-591                COMP            G320.4-1.2   0.070483\n",
      "  HESS J1514-591                 PSR              B1509-58   0.026891\n",
      "\n"
     ]
    }
   ],
   "source": [
    "print(source.info(\"associations\"))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.9.0"
  },
  "latex_envs": {
   "LaTeX_envs_menu_present": true,
   "autoclose": false,
   "autocomplete": true,
   "bibliofile": "biblio.bib",
   "cite_by": "apalike",
   "current_citInitial": 1.0,
   "eqLabelWithNumbers": true,
   "eqNumInitial": 1.0,
   "hotkeys": {
    "equation": "Ctrl-E",
    "itemize": "Ctrl-I"
   },
   "labels_anchors": false,
   "latex_user_defs": false,
   "report_style_numbering": false,
   "user_envs_cfg": false
  },
  "nbsphinx": {
   "orphan": true
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}