Fermi_LAT.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "%matplotlib inline\n",
    "from astropy.io import fits"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Filename: photons_100GeV.fits\n",
      "No.    Name      Ver    Type      Cards   Dimensions   Format\n",
      "  0  PRIMARY       1 PrimaryHDU      31   ()      \n",
      "  1  EVENTS        1 BinTableHDU    284   21917R x 50C   [E, E, E, E, E, E, E, E, E, D, J, J, I, 3I, 32X, 32X, I, D, E, E, E, E, E, E, E, E, E, E, E, E, E, E, E, E, E, E, E, E, E, E, E, E, E, E, E, E, E, E, E, E]   \n",
      "  2  GTI           1 BinTableHDU     46   53027R x 2C   [D, D]   \n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "XTENSION= 'BINTABLE'           / binary table extension                         \n",
       "BITPIX  =                    8 / 8-bit bytes                                    \n",
       "NAXIS   =                    2 / 2-dimensional binary table                     \n",
       "NAXIS1  =                  206 / width of table in bytes                        \n",
       "NAXIS2  =                21917 / number of rows in table                        \n",
       "PCOUNT  =                    0 / size of special data area                      \n",
       "GCOUNT  =                    1 / one data group (required keyword)              \n",
       "TFIELDS =                   50 / number of fields in each row                   \n",
       "TTYPE1  = 'ENERGY  '           / energy of event                                \n",
       "TFORM1  = 'E       '           / data format of field: 4-byte REAL              \n",
       "TTYPE2  = 'RA      '           / right ascension (J2000) of event               \n",
       "TFORM2  = 'E       '           / data format of field: 4-byte REAL              \n",
       "TTYPE3  = 'DEC     '           / declination (J2000) of event                   \n",
       "TFORM3  = 'E       '           / data format of field: 4-byte REAL              \n",
       "TTYPE4  = 'L       '           / Galactic longitude of event                    \n",
       "TFORM4  = 'E       '           / data format of field: 4-byte REAL              \n",
       "TTYPE5  = 'B       '           / Galactic latitude of event                     \n",
       "TFORM5  = 'E       '           / data format of field: 4-byte REAL              \n",
       "TTYPE6  = 'THETA   '           / inclination angle of event in instrument coordi\n",
       "TFORM6  = 'E       '           / data format of field: 4-byte REAL              \n",
       "TTYPE7  = 'PHI     '           / azimuthal angle of event in instrument coordina\n",
       "TFORM7  = 'E       '           / data format of field: 4-byte REAL              \n",
       "TTYPE8  = 'ZENITH_ANGLE'       / zenith angle of event                          \n",
       "TFORM8  = 'E       '           / data format of field: 4-byte REAL              \n",
       "TTYPE9  = 'EARTH_AZIMUTH_ANGLE' / Earth azimuth (from north to east) of event   \n",
       "TFORM9  = 'E       '           / data format of field: 4-byte REAL              \n",
       "TTYPE10 = 'TIME    '           / Mission Elapsed Time                           \n",
       "TFORM10 = 'D       '           / data format of field: 8-byte DOUBLE            \n",
       "TTYPE11 = 'EVENT_ID'           / ID number of original event                    \n",
       "TFORM11 = 'J       '           / data format of field: 4-byte signed INTEGER    \n",
       "TTYPE12 = 'RUN_ID  '           / Run number of original event                   \n",
       "TFORM12 = 'J       '           / data format of field: 4-byte signed INTEGER    \n",
       "TTYPE13 = 'RECON_VERSION'      / version of event reconstruction software       \n",
       "TFORM13 = 'I       '           / data format of field: 2-byte signed INTEGER    \n",
       "TTYPE14 = 'CALIB_VERSION'      / versions of calibration tables for the ACD, CAL\n",
       "TFORM14 = '3I      '           / data format of field: 2-byte signed INTEGER    \n",
       "TTYPE15 = 'EVENT_CLASS'        / event class: 0, 1, 2,...                       \n",
       "TFORM15 = '32X     '           / data format of field: 32 bit array             \n",
       "TTYPE16 = 'EVENT_TYPE'         / event types encoded as bit array               \n",
       "TFORM16 = '32X     '           / data format of field: 32 bit array             \n",
       "TTYPE17 = 'CONVERSION_TYPE'    / type of conversion: 0=Front converting, 1=Back \n",
       "TFORM17 = 'I       '           / data format of field: 2-byte signed INTEGER    \n",
       "TTYPE18 = 'LIVETIME'           / Accumulated livetime since mission start       \n",
       "TFORM18 = 'D       '           / data format of field: 8-byte DOUBLE            \n",
       "TTYPE19 = 'DIFRSP0 '           / Diffuse response component                     \n",
       "TFORM19 = 'E       '           / data format                                    \n",
       "TTYPE20 = 'DIFRSP1 '           / Diffuse response component                     \n",
       "TFORM20 = 'E       '           / data format                                    \n",
       "TTYPE21 = 'DIFRSP2 '           / Diffuse response component                     \n",
       "TFORM21 = 'E       '           / data format                                    \n",
       "TTYPE22 = 'DIFRSP3 '           / Diffuse response component                     \n",
       "TFORM22 = 'E       '           / data format                                    \n",
       "TTYPE23 = 'DIFRSP4 '           / Diffuse response component                     \n",
       "TFORM23 = 'E       '           / data format                                    \n",
       "CHECKSUM= '7VRA8UR97URA7UR9'   / HDU checksum updated 2020-05-06T21:50:48       \n",
       "DATASUM = '4283111551'         / data unit checksum updated 2020-05-06T21:50:48 \n",
       "TELESCOP= 'GLAST   '           / name of telescope generating data              \n",
       "INSTRUME= 'LAT     '           / name of instrument generating data             \n",
       "EQUINOX =                2000. / equinox for ra and dec                         \n",
       "RADECSYS= 'FK5     '           / world coord. system for this file (FK5 or FK4) \n",
       "DATE    = '2020-05-06T21:50:48.0000' / file creation date (YYYY-MM-DDThh:mm:ss U\n",
       "DATE-OBS= '2008-10-27T00:00:00.0000' / start date and time of the observation (U\n",
       "DATE-END= '2018-04-30T05:19:56.0000' / end date and time of the observation (UTC\n",
       "OBSERVER= 'Peter Michelson'    / GLAST/LAT PI                                   \n",
       "ORIGIN  = 'FSSC    '           / name of organization making file               \n",
       "EXTNAME = 'EVENTS  '           / name of this binary table extension            \n",
       "HDUCLASS= 'OGIP    '           / format conforms to OGIP standard               \n",
       "HDUCLAS1= 'EVENTS  '           / extension contains events                      \n",
       "HDUCLAS2= 'ALL     '           / extension contains all events detected         \n",
       "TSTART  =           246758401. / mission time of the start of the observation   \n",
       "TSTOP   =           546758401. / mission time of the end of the observation     \n",
       "MJDREFI =               51910. / Integer part of MJD corresponding to SC clock s\n",
       "MJDREFF =  0.00074287037037037 / Fractional part of MJD corresponding to SC cloc\n",
       "TIMEUNIT= 's       '           / units for the time related keywords            \n",
       "TIMEZERO=                   0. / clock correction                               \n",
       "TIMESYS = 'TT      '           / type of time system that is used               \n",
       "TIMEREF = 'LOCAL   '           / reference frame used for times                 \n",
       "CLOCKAPP=                    F / whether a clock drift correction has been appli\n",
       "GPS_OUT =                    F / whether GPS time was unavailable at any time du\n",
       "PASS_VER= 'P8R3    '           / IRF pass version corresponding to a specific se\n",
       "NDIFRSP =                    5 / Number of diffuse response labels              \n",
       "DIFRSP0 = 'NONE    '           / Diffuse response label for component 0         \n",
       "DIFRSP1 = 'NONE    '           / Diffuse response label for component 1         \n",
       "DIFRSP2 = 'NONE    '           / Diffuse response label for component 2         \n",
       "DIFRSP3 = 'NONE    '           / Diffuse response label for component 3         \n",
       "DIFRSP4 = 'NONE    '           / Diffuse response label for component 4         \n",
       "TUNIT1  = 'MeV     '           / physical unit of field                         \n",
       "TLMIN1  =                   0. / minimum value                                  \n",
       "TLMAX1  =            10000000. / maximum value                                  \n",
       "TUNIT2  = 'deg     '           / physical unit of field                         \n",
       "TLMIN2  =                   0. / minimum value                                  \n",
       "TLMAX2  =                 360. / maximum value                                  \n",
       "TUNIT3  = 'deg     '           / physical unit of field                         \n",
       "TLMIN3  =                 -90. / minimum value                                  \n",
       "TLMAX3  =                  90. / maximum value                                  \n",
       "TUNIT4  = 'deg     '           / physical unit of field                         \n",
       "TLMIN4  =                   0. / minimum value                                  \n",
       "TLMAX4  =                 360. / maximum value                                  \n",
       "TUNIT5  = 'deg     '           / physical unit of field                         \n",
       "TLMIN5  =                 -90. / minimum value                                  \n",
       "TLMAX5  =                  90. / maximum value                                  \n",
       "TUNIT6  = 'deg     '           / physical unit of field                         \n",
       "TLMIN6  =                   0. / minimum value                                  \n",
       "TLMAX6  =                 180. / maximum value                                  \n",
       "TUNIT7  = 'deg     '           / physical unit of field                         \n",
       "TLMIN7  =                   0. / minimum value                                  \n",
       "TLMAX7  =                 360. / maximum value                                  \n",
       "TUNIT8  = 'deg     '           / physical unit of field                         \n",
       "TLMIN8  =                   0. / minimum value                                  \n",
       "TLMAX8  =                 180. / maximum value                                  \n",
       "TUNIT9  = 'deg     '           / physical unit of field                         \n",
       "TLMIN9  =                   0. / minimum value                                  \n",
       "TLMAX9  =                 360. / maximum value                                  \n",
       "TUNIT10 = 's       '           / physical unit of field                         \n",
       "TLMIN10 =                   0. / minimum value                                  \n",
       "TLMAX10 =         10000000000. / maximum value                                  \n",
       "TLMIN11 =                    0 / minimum value                                  \n",
       "TLMAX11 =           2147483647 / maximum value                                  \n",
       "TLMIN12 =                    0 / minimum value                                  \n",
       "TLMAX12 =           2147483647 / maximum value                                  \n",
       "TLMIN13 =                    0 / minimum value                                  \n",
       "TLMAX13 =                32767 / maximum value                                  \n",
       "TLMIN15 =                    0 / minimum value                                  \n",
       "TLMAX15 =                32767 / maximum value                                  \n",
       "TLMIN16 =                    0 / minimum value                                  \n",
       "TLMAX16 =                32767 / maximum value                                  \n",
       "TLMIN17 =                    0 / minimum value                                  \n",
       "TLMAX17 =                32767 / maximum value                                  \n",
       "TUNIT18 = 's       '           / physical unit of field                         \n",
       "TLMIN18 =                   0. / minimum value                                  \n",
       "TLMAX18 =         10000000000. / maximum value                                  \n",
       "TLMIN19 =                   0. / minimum value                                  \n",
       "TLMAX19 =              1.0E+38 / maximum value                                  \n",
       "TLMIN20 =                   0. / minimum value                                  \n",
       "TLMAX20 =              1.0E+38 / maximum value                                  \n",
       "TLMIN21 =                   0. / minimum value                                  \n",
       "TLMAX21 =              1.0E+38 / maximum value                                  \n",
       "TLMIN22 =                   0. / minimum value                                  \n",
       "TLMAX22 =              1.0E+38 / maximum value                                  \n",
       "TLMIN23 =                   0. / minimum value                                  \n",
       "TLMAX23 =              1.0E+38 / maximum value                                  \n",
       "EXTVER  =                    1 / auto assigned by template parser               \n",
       "TTYPE24 = 'CalEnergyRaw'       / label for field                                \n",
       "TFORM24 = 'E       '           / format of field                                \n",
       "TTYPE25 = 'PtAlt   '           / label for field                                \n",
       "TFORM25 = 'E       '           / format of field                                \n",
       "TTYPE26 = 'PtDecx  '           / label for field                                \n",
       "TFORM26 = 'E       '           / format of field                                \n",
       "TTYPE27 = 'PtDecz  '           / label for field                                \n",
       "TFORM27 = 'E       '           / format of field                                \n",
       "TTYPE28 = 'PtLat   '           / label for field                                \n",
       "TFORM28 = 'E       '           / format of field                                \n",
       "TTYPE29 = 'PtLon   '           / label for field                                \n",
       "TFORM29 = 'E       '           / format of field                                \n",
       "TTYPE30 = 'PtMagLat'           / label for field                                \n",
       "TFORM30 = 'E       '           / format of field                                \n",
       "TTYPE31 = 'PtPosx  '           / label for field                                \n",
       "TFORM31 = 'E       '           / format of field                                \n",
       "TTYPE32 = 'PtPosy  '           / label for field                                \n",
       "TFORM32 = 'E       '           / format of field                                \n",
       "TTYPE33 = 'PtPosz  '           / label for field                                \n",
       "TFORM33 = 'E       '           / format of field                                \n",
       "TTYPE34 = 'PtRax   '           / label for field                                \n",
       "TFORM34 = 'E       '           / format of field                                \n",
       "TTYPE35 = 'PtRaz   '           / label for field                                \n",
       "TFORM35 = 'E       '           / format of field                                \n",
       "TTYPE36 = 'PtSCzenith'         / label for field                                \n",
       "TFORM36 = 'E       '           / format of field                                \n",
       "TTYPE37 = 'Tkr1FirstLayer'     / label for field                                \n",
       "TFORM37 = 'E       '           / format of field                                \n",
       "TTYPE38 = 'WP8BestSXX'         / label for field                                \n",
       "TFORM38 = 'E       '           / format of field                                \n",
       "TTYPE39 = 'WP8BestSXY'         / label for field                                \n",
       "TFORM39 = 'E       '           / format of field                                \n",
       "TTYPE40 = 'WP8BestSYY'         / label for field                                \n",
       "TFORM40 = 'E       '           / format of field                                \n",
       "TTYPE41 = 'WP8BestXDir'        / label for field                                \n",
       "TFORM41 = 'E       '           / format of field                                \n",
       "TTYPE42 = 'WP8BestYDir'        / label for field                                \n",
       "TFORM42 = 'E       '           / format of field                                \n",
       "TTYPE43 = 'WP8BestZDir'        / label for field                                \n",
       "TFORM43 = 'E       '           / format of field                                \n",
       "TTYPE44 = 'WP8CTAllBkProb'     / label for field                                \n",
       "TFORM44 = 'E       '           / format of field                                \n",
       "TTYPE45 = 'WP8CTAllProb'       / label for field                                \n",
       "TFORM45 = 'E       '           / format of field                                \n",
       "TTYPE46 = 'WP8CTBestEnergyProb' / label for field                               \n",
       "TFORM46 = 'E       '           / format of field                                \n",
       "TTYPE47 = 'WP8CTCalTkrBkProb'  / label for field                                \n",
       "TFORM47 = 'E       '           / format of field                                \n",
       "TTYPE48 = 'WP8CTCalTkrProb'    / label for field                                \n",
       "TFORM48 = 'E       '           / format of field                                \n",
       "TTYPE49 = 'WP8CTPSFCore'       / label for field                                \n",
       "TFORM49 = 'E       '           / format of field                                \n",
       "TTYPE50 = 'WP8CTPSFTail'       / label for field                                \n",
       "TFORM50 = 'E       '           / format of field                                \n",
       "HISTORY Input merit file: /lustre/ki/pfs/fermi_scratch/glastmp/P305-FITS.mergeRu\n",
       "HISTORY n/0.0/IN/r0239557414_v301_merit.root                                    \n",
       "HISTORY Filter string: (((FswGamState==0)&&((GltGemSummary&0x60)==0)&&(TkrNumTra\n",
       "HISTORY cks > 0))&&((-log10(1.0-WP8CTAllProb)) >= (0.010000)*(EvtJointLogEnergy \n",
       "HISTORY < 1.250000) + ((EvtJointLogEnergy <= 1.750000)*((0.010000)*(1.0)+(0.0000\n",
       "HISTORY 00)*(pow((EvtJointLogEnergy-1.250000)/0.500000,1))+(0.018669)*(pow((EvtJ\n",
       "HISTORY ointLogEnergy-1.250000)/0.500000,2)))+((EvtJointLogEnergy > 1.750000)*(E\n",
       "HISTORY vtJointLogEnergy <= 2.250000))*((0.028669)*(1.0)+(0.037338)*(pow((EvtJoi\n",
       "HISTORY ntLogEnergy-1.750000)/0.500000,1))+(-0.017111)*(pow((EvtJointLogEnergy-1\n",
       "HISTORY .750000)/0.500000,2)))+((EvtJointLogEnergy > 2.250000)*(EvtJointLogEnerg\n",
       "HISTORY y <= 2.750000))*((0.048897)*(1.0)+(0.003117)*(pow((EvtJointLogEnergy-2.2\n",
       "HISTORY 50000)/0.500000,1))+(0.001967)*(pow((EvtJointLogEnergy-2.250000)/0.50000\n",
       "HISTORY 0,2)))+((EvtJointLogEnergy > 2.750000)*(EvtJointLogEnergy <= 3.250000))*\n",
       "HISTORY ((0.053980)*(1.0)+(0.007050)*(pow((EvtJointLogEnergy-2.750000)/0.500000,\n",
       "HISTORY 1))+(-0.003525)*(pow((EvtJointLogEnergy-2.750000)/0.500000,2)))+((EvtJoi\n",
       "HISTORY ntLogEnergy > 3.250000)*(EvtJointLogEnergy <= 3.750000))*((0.057505)*(1.\n",
       "HISTORY 0)+(0.000000)*(pow((EvtJointLogEnergy-3.250000)/0.500000,1))+(0.121963)*\n",
       "HISTORY (pow((EvtJointLogEnergy-3.250000)/0.500000,2)))+((EvtJointLogEnergy > 3.\n",
       "HISTORY 750000)*(EvtJointLogEnergy <= 4.250000))*((0.179468)*(1.0)+(0.243925)*(p\n",
       "HISTORY ow((EvtJointLogEnergy-3.750000)/0.500000,1))+(0.493075)*(pow((EvtJointLo\n",
       "HISTORY gEnergy-3.750000)/0.500000,2)))+((EvtJointLogEnergy > 4.250000)*(EvtJoin\n",
       "HISTORY tLogEnergy <= 4.750000))*((0.916468)*(1.0)+(1.230076)*(pow((EvtJointLogE\n",
       "HISTORY nergy-4.250000)/0.500000,1))+(-0.501532)*(pow((EvtJointLogEnergy-4.25000\n",
       "HISTORY 0)/0.500000,2)))+(EvtJointLogEnergy > 4.750000)*((1.645012)*(1.0)+(0.227\n",
       "HISTORY 011)*(pow((EvtJointLogEnergy-4.750000)/0.500000,1))+(0.029483)*(pow((Evt\n",
       "HISTORY JointLogEnergy-4.750000)/0.500000,2))))*(EvtJointLogEnergy >= 1.250000 &\n",
       "HISTORY & EvtJointLogEnergy <= 5.750000) + (2.216967)*(EvtJointLogEnergy > 5.750\n",
       "HISTORY 000))) || (((-log10(1.0-WP8CTCalTkrProb)) >=  -0.035931+ (   0.27029)*po\n",
       "HISTORY w(EvtJointLogEnergy,1)+(  -0.23349)*pow(EvtJointLogEnergy,2)+ (   0.0786\n",
       "HISTORY 7)*pow(EvtJointLogEnergy,3)+ (-0.005295)*pow(EvtJointLogEnergy,4))&&((Fs\n",
       "HISTORY wGamState==0)&&((GltGemSummary&0x60)==0)&&(TkrNumTracks > 0)&&EvtCalCsIR\n",
       "HISTORY Ln>4 && ((log10(max(CalTrackAngle,1E-4))) <= (0.529795)*(EvtJointLogEner\n",
       "HISTORY gy < 3.000000)+ ((1.0)*((0.529795)*(1.0)+(-1.379791)*(pow((EvtJointLogEn\n",
       "HISTORY ergy-3.000000)/0.916667,1))+(0.583401)*(pow((EvtJointLogEnergy-3.000000)\n",
       "HISTORY /0.916667,2))+(-0.075555)*(pow((EvtJointLogEnergy-3.000000)/0.916667,3))\n",
       "HISTORY ))*(EvtJointLogEnergy>= 3.000000 && EvtJointLogEnergy <= 5.750000) + (-0\n",
       "HISTORY .398962)*(EvtJointLogEnergy > 5.750000)))&&(WP8CTPSFTail > 0.05 && WP8CT\n",
       "HISTORY BestEnergyProb > 0.1)) && (EvtElapsedTime >= 239557417)  && (EvtElapsedT\n",
       "HISTORY ime <= 239558070)                                                       \n",
       "HISTORY CFITSIO used the following filtering expression to create this table:   \n",
       "HISTORY /lustre/ki/pfs/fermi_scratch/glastmp/P305-FITS.mergeRun/0.0/foo.fit[EVEN\n",
       "HISTORY TS][gtifilter(\"/lustre/ki/pfs/fermi_scratch/glastmp/P305-FITS.mergeRun/0\n",
       "HISTORY .0/OUT/gll_xe_p305_r0239557414_v305.fit_tempgti\")]                      \n",
       "HISTORY CFITSIO used the following filtering expression to create this table:   \n",
       "HISTORY /lustre/ki/pfs/fermi_scratch/glastmp/P305-FITS.mergeRun/0.0/OUT/gll_xe_p\n",
       "HISTORY 305_r0239557414_v305.fit[EVENTS][((EVENT_CLASS&o200010) != o0) && gtifil\n",
       "HISTORY ter()]                                                                  \n",
       "HISTORY Filter string: ((EVENT_CLASS&o200010) != o0) && gtifilter()             \n",
       "HISTORY CFITSIO used the following filtering expression to create this table:   \n",
       "HISTORY /lustre/ki/pfs/fermi_scratch/glastmp/P305-FITS.mergeRun/0.0/OUT/foo.fits\n",
       "HISTORY [EVENTS][gtifilter(\"/lustre/ki/pfs/fermi_scratch/glastmp/P305-FITS.merge\n",
       "HISTORY Run/0.0/OUT/gll_ev_p305_r0239557414_v305.fit_tempgti\")]                 \n",
       "HISTORY CFITSIO used the following filtering expression to create this table:   \n",
       "HISTORY /Volumes/Seagate Backup Plus Drive/weekly/extended/lat_extend_weekly_w00\n",
       "HISTORY 9_p305_v001.fits[EVENTS][((EVENT_CLASS&o4000) != o0) && ((EVENT_TYPE&o3)\n",
       "HISTORY  != o0) && 1000 < ENERGY && ENERGY <= 10000000 && 246758401 < TIME && TI\n",
       "HISTORY ME <= 1000000000000000 && 0 < ZENITH_ANGLE && ZENITH_ANGLE <= 100 && gti\n",
       "HISTORY filter()]                                                               \n",
       "HISTORY Filter string: ((EVENT_CLASS&o4000) != o0) && ((EVENT_TYPE&o3) != o0) &&\n",
       "HISTORY  1000 < ENERGY && ENERGY <= 10000000 && 246758401 < TIME && TIME <= 1000\n",
       "HISTORY 000000000000 && 0 < ZENITH_ANGLE && ZENITH_ANGLE <= 100 && gtifilter()  \n",
       "HISTORY CFITSIO used the following filtering expression to create this table:   \n",
       "HISTORY f.fits[EVENTS][gtifilter(\"ph_2048_Jan2020.fits_tempgti\")]               \n",
       "HISTORY CFITSIO used the following filtering expression to create this table:   \n",
       "HISTORY ph_2048_Jan2020.fits[EVENTS][((EVENT_CLASS&o4000) != o0) && ((EVENT_TYPE\n",
       "HISTORY &o3) != o0) && 100000 < ENERGY && ENERGY <= 10000000 && 246758401 < TIME\n",
       "HISTORY  && TIME <= 546758401 && 0 < ZENITH_ANGLE && ZENITH_ANGLE <= 100 && gtif\n",
       "HISTORY ilter()]                                                                \n",
       "DSTYP1  = 'BIT_MASK(EVENT_CLASS,2048,P8R3)'                                     \n",
       "DSUNI1  = 'DIMENSIONLESS'                                                       \n",
       "DSVAL1  = '1:1     '                                                            \n",
       "DSTYP2  = 'BIT_MASK(EVENT_TYPE,3,P8R3)'                                         \n",
       "DSUNI2  = 'DIMENSIONLESS'                                                       \n",
       "DSVAL2  = '1:1     '                                                            \n",
       "DSTYP3  = 'TIME    '                                                            \n",
       "DSUNI3  = 's       '                                                            \n",
       "DSVAL3  = 'TABLE   '                                                            \n",
       "DSREF3  = ':GTI    '                                                            \n",
       "DSTYP4  = 'ENERGY  '                                                            \n",
       "DSUNI4  = 'MeV     '                                                            \n",
       "DSVAL4  = '100000:10000000'                                                     \n",
       "DSTYP5  = 'ZENITH_ANGLE'                                                        \n",
       "DSUNI5  = 'deg     '                                                            \n",
       "DSVAL5  = '0:100   '                                                            \n",
       "NDSKEYS =                    5                                                  \n",
       "HISTORY Filter string: ((EVENT_CLASS&o4000) != o0) && ((EVENT_TYPE&o3) != o0) &&\n",
       "HISTORY  100000 < ENERGY && ENERGY <= 10000000 && 246758401 < TIME && TIME <= 54\n",
       "HISTORY 6758401 && 0 < ZENITH_ANGLE && ZENITH_ANGLE <= 100 && gtifilter()       "
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "photon_list=fits.open('photons_100GeV.fits')\n",
    "photon_list.info()\n",
    "photon_list[1].header"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "There are overall 21917 photons\n"
     ]
    }
   ],
   "source": [
    "photons=photon_list[1].data\n",
    "print('There are overall',len(photons),'photons')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnMAAAG/CAYAAAAzclFcAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3de5QlZX3u8e+jSBA07agIRhiGcUQDJvEsR4MaA5qjicGBqECOx+NdBpN4N5oRNeIlCV6ix6AYx/stR6Mx4mQMXkCRgxIDJwpiUFAGxaigYHNTFPidP2q3NpvePbVn9qWr+/tZa6/eu97aVb+NFjzrrXrfN1WFJEmSuulW0y5AkiRJO84wJ0mS1GGGOUmSpA4zzEmSJHWYYU6SJKnDDHOSJEkdZpiTJEnqsGUT5pI8OUkl+aNp1yJJkjQpyyLMJdkPOAY4a9q1SJIkTdLEw1ySfZKcmOSLSa7r9aatGbDvvkk+kmQ2yVVJPppkdd8+twLeATwTuH7sP0CSJGkJmUbP3DrgaOBK4IxBOyXZHTgNuBfwRODxwD2AzybZY96uzwPOrKpzxlaxJEnSErXLFM75+araCyDJ04CHD9jvGGAtcM+quqi3/7nAhcCxwOuTHAQcCTx47FVLkiQtQRPvmauqm1ruejhw1lyQ6333YuBM4Ijept8F9gMuTLINOBjYnORPRlexJEnS0rWUB0AcBHx1ge3nAwcCVNVbququVbWmqtbQDIDYWFVvmVyZkiRJ0zON26xt3ZHmubp+VwCrduSASTYCGwH22GOP+97rXvfa8eokSZIm5JxzzvlhVe25UNtSDnMAtcC2DNy56tBFD1a1GdgMsH79+jr77LN3qjhJkqRJSHLJoLalfJv1SpreuX6rWLjHTpIkacVZymHufJrn5vodCHxtwrVIkiQtSUs5zH0cODjJ2rkNvcmFH9Rr2yFJNiTZPDs7u9MFSpIkTdtUwlySI5McCdy3t+kRvW2HzNvtbcA24OQkRyQ5HDgZ+A7w1h09d1VtqaqNMzMzO3oISZKkJWNaAyA+3Pf5pN7f04FDAarq2iQPBd4AvI9m4MOpwHOq6poJ1SlJkrSkTSXMVdXAEal9+30beMyYy5EkSeqspfzM3Fj4zJwkSVpOVlyY85k5SZK0nKy4MCdJkrScGOYkSZI6zDAnSZLUYUt9bdaRS7IB2LBu3bqJnG/Npq0TOc+2Ew6byHkkSdLSsuJ65hwAIUmSlpMVF+YkSZKWE8OcJElShxnmJEmSOswwJ0mS1GErLsy5nJckSVpOVlyYczSrJElaTlZcmJMkSVpODHOSJEkdZpiTJEnqMMOcJElShxnmJEmSOmyXaRcwaUk2ABvWrVs37VJGas2mrRM717YTDpvYuSRJ0uJWXM+cU5NIkqTlZMWFOUmSpOXEMCdJktRhhjlJkqQOM8xJkiR1mGFOkiSpwwxzkiRJHWaYkyRJ6rAVF+aSbEiyeXZ2dtqlSJIk7bQVF+acNFiSJC0nKy7MSZIkLSeGOUmSpA4zzEmSJHWYYU6SJKnDDHOSJEkdZpiTJEnqsF2mXYC6Z82mrRM717YTDpvYuSRJ6iJ75iRJkjrMMCdJktRhKy7MuZyXJElaTlZcmHM5L0mStJysuDAnSZK0nBjmJEmSOswwJ0mS1GGGOUmSpA4zzEmSJHWYYU6SJKnDDHOSJEkdZpiTJEnqMMOcJElShxnmJEmSOmyXaRcgLWbNpq0TOc+2Ew6byHkkSRo1e+YkSZI6zDAnSZLUYSsuzCXZkGTz7OzstEuRJEnaaSsuzFXVlqraODMzM+1SJEmSdtqKC3OSJEnLiWFOkiSpwwxzkiRJHWaYkyRJ6jDDnCRJUocZ5iRJkjrMMCdJktRhhjlJkqQOM8xJkiR1mGFOkiSpwwxzkiRJHbbLtAuQloI1m7ZO7FzbTjhsYueSJC1/9sxJkiR1mGFOkiSpwwxzkiRJHWaYkyRJ6jDDnCRJUocZ5iRJkjpsWUxNkuRU4M5AAVcDz6yqL0+3KkmSpPFbFmEOeHRVzQIkeRTwbuA+U61IkiRpAqZymzXJPklOTPLFJNclqSRrBuy7b5KPJJlNclWSjyZZPX+fuSDX86tjLF2SJGlJaRXmkuya5GVJLuiFrxv7XjcMed51wNHAlcAZi5x3d+A04F7AE4HHA/cAPptkj759P5DkUuCVwP8ash5JkqROanub9bXAnwH/CnwUuH4nz/v5qtoLIMnTgIcP2O8YYC1wz6q6qLf/ucCFwLHA6+d2rKrHzTveqwHXTJIkScte2zB3JPCyqvqrUZy0qm5quevhwFlzQa733YuTnAkcwbwwN887gJOS3KmqfrTz1UqSJC1dbcPc7YAvjrOQAQ4CTl5g+/nAUQBJVgG7VdX3em2PAS4DrphIhdKQ1mzaOrFzbTvBDmpJWu7ahrktwO/SPL82SXekea6u3xXAqt77VcCHkuwG3EQT5B5ZVdX/pSQbgY0Aq1ev7m+WJEnqnLZh7kTgvUluAj7BAr1eVfWtURY2/9ALbEvfee/X6kBVm4HNAOvXr1/ouJIkSZ3SNszN3WI9HnjZgH1uvdPV3NKVNL1z/VaxcI+dJEnSitI2zD2FhXvIxu18mufm+h0IfG3CtUiSJC05rcJcVb17zHUM8nHgdUnWzt3G7U0u/CBg044cMMkGYMO6detGVaMkSdLUTGUFCIAkRyY5Erhvb9MjetsOmbfb24BtwMlJjkhyOM3o1u8Ab92R81bVlqraODMzsxPVS5IkLQ0De+aSvBN4ZW9et3du5zhVVU8d8twf7vt8Uu/v6cChvYNem+ShwBuA99EMfDgVeE5VXTPk+SRJkpadxW6zPgR4Y+/9Q1n8mbmhn6erqmx/L6iqb9PMHSdJkqQ+A8NcVe0/7/2aiVQjSZKkoUztmblpSbIhyebZ2dlplyJJkrTTWoe5JLdO8uQkm5Ns7f19UpJxzC83Ng6AkCRJy0mrMJdkP5o5394B/AFwl97fdwJf7bVLkiRpwtr2zL0J+FXgd6pqdVXdr6pWAw8GZmiW+5IkSdKEtQ1zDwVeVFVfmL+xqs4Ejuu1S5IkacLahrlrgMsGtF0GXDeacsbPARCSJGk5abs26/uBpwP/ukDbscB7R1bRmFXVFmDL+vXrj5l2LdK4rdm0dSLn2XbCYRM5jyTplhZbAeIp8z5eCByV5Dzgn4AfAHsBRwK3Z+GQJ0mSpDFbrGfu7Qts2wc4aIHtbwb+fiQVSZIkqbXFwtz+i7RJkiRpCVhsOa9LJlmIJEmShudyXpIkSR224sKcy3lJkqTlZMWFOUmSpOXEMCdJktRh2w1zSW6d5LeS7DmJgiRJktRem565As4G/tuYa5EkSdKQthvmquom4DvAHuMvR5IkScNo+8zcW4HnJNl1nMVMglOTSJKk5WSxFSDmuz1wd+BbSU4Bvkdz+3VOVdXLRl3cOFTVFmDL+vXrj5l2LZIkSTurbZg7bt77pyzQXkAnwpyk0VuzaevEzrXthMMmdi5J6oJWYa6qnMJEkiRpCTKkSZIkdVjrMJfG4Ulel+RdSfbrbT8kya+Nr0RJkiQN0uo2a5JVwCeA3wauohkQcSJwCXAMcAXwrDHVKEmSpAHa9sy9FtgXeBBwZyDz2j4D/N6I65IkSVILbcPcEcCLq+qL3HxKEoBv0wS9TnCeOUmStJy0DXO3A747oG03bt5Tt6RV1Zaq2jgzMzPtUiRJknZa2zD3deDhA9oOAc4bTTmSJEkaRttJg98MvDnJLPAPvW13SPJk4BnAxnEUJ0mSpMW1nTT4bUnuDrwceEVv86eBm4DXVNUHxlSfJEmSFtG2Z46q2pTkLTS3W/cEfgR8uqq+Na7iJEmStLjWYQ6gqi4B3jamWiRJkjSk1mEuya2BJwAPAO5GM7r1C8D7qurG8ZQnSZKkxbRdAWI/4JPAAcClwA+A3wCeBvxFkj/o9dpJ0lit2bR1YufadsJhEzuXJO2otlOTvAn4VeB3qmp1Vd2vqlYDDwZmaJb2kiRJ0oS1DXMPBV5UVV+Yv7GqzgSO67VLkiRpwtqGuWuAywa0XQZcN5pyxs/lvCRJ0nLSNsy9H3j6gLZjgfeOppzxczkvSZK0nLQdzXoRcFSS84B/ohkAsRdwJHB74F+TPGVu56p656gLlSRJ0i0Ns5wXwD7AQQu0nzTvfQGGOUmSpAloG+b2H2sVkiRJ2iFt12Z1DjlJkqQlqO0ACEmSJC1BhjlJkqQOM8xJkiR1mGFOkiSpwwxzkiRJHdZqNGuSWwG3qqob5m37feDewGlV9R9jqk+SpmbNpq0TO9e2Ew6b2LkkLS9t55n7P8D1wBMAkjydX04U/PMkh1XVZ8ZQnyRJkhbR9jbrwcAn5n1+AfB2YAb4KPDiEdclSZKkFtqGubsA3wVIso5mRYg3VdXVwLuA3xhPeZIkSVpM2zB3FXCn3vtDgR9W1bm9zzcCu424LkmSJLXQ9pm5LwCbktwAPIeb33JdB1w66sLGJckGYMO6deumXYokSdJOa9sz90LgjsDHaXrhjp/X9sfAF0db1vhU1Zaq2jgzMzPtUiRJknZaq565qroQOCDJnarqR33Nzwa+P/LKJEmStF1tb7MCsECQo6rOG105kiRJGkbrMJdkLXA0sJpbDnioqnrqKAuTJEnS9rVdAeII4MM0z9hdRjOB8Hw14rokSZLUQtueuVcBnwMeV1WXj68cSZIkDaNtmFsLPN8gJ0mStLS0nZrkAn45abAkSZKWiGHmmTuuNwhCkiRJS0Tb26zH0/TM/WeSC4Er+tqrqg4ZZWGSJEnavrZh7kbg6+MsRJIkScNruwLEoWOuQ5IkSTug7TNzkiRJWoJah7kkd03yuiT/nuSbSb6U5DVJ9h5ngZIkSRqs7QoQBwBnAKuAM4GLgL2BZwNPSPLgqrpwbFVK0jK3ZtPWiZxn2wmHTeQ8kian7QCIVwNXAb9dVdvmNibZD/hUr/3RI69OkiRJi2p7m/UhwEvnBzmAqrqEZtqSh4y2LEmSJLXRNsztClw9oO3qXrskSZImrG2Y+zLwzCQ32z9JgD/ttUuSJGnC2j4z9wrgX2hWgPgQ8D2aARBHAfcApvZEbZLdgA8C9wSuB34A/ElVfWtaNUmSJE1Kq565qjoFeCTNLdUXA28GXgJcAzyyqj41tgrbeUtV/XpV3QfYArx9yvVIkiRNROt55qrqlKpaD9we2Be4fVXdv6o+OexJk+yT5MQkX0xyXZJKsmbAvvsm+UiS2SRXJfloktXz6vppXw1nAWuHrUmSJKmLWoW5JO9Msj9AVV1XVd+tqut6bfsleeeQ510HHA1cSTN/3aDz7g6cBtwLeCLweJrbup9NsseArz0TOHnIeiRJkjqpbc/ck4A9B7TdmSZoDePzVbVXVf0h8OFF9juGppftj6rqY1V1MnA4sB9wbP/OSV4EHAC8aMh6JEmSOmmYtVlrwPa9gZ8Mc9KquqnlrocDZ1XVRfO+ezHNKhRHzN8xyZ8DjwEeMddrKEmStNwNHM2a5FHAo+ZtenmSH/btdlvgwcA5Y6gN4CAWvmV6Ps1IWgCSPA94LPDfq+rHY6pFkiRpyVlsapLVNEENml65+9BM/THf9cAXGN9tzTvSPFfX7wqadWJJsg/wt8C3aJ6lA7ihN1jjZpJsBDYCrF69ur9ZkiSpcwaGuap6I/BGgCQX0zy39pVJFTa/lAW25ReNVZfO/7zogao2A5sB1q9fP+i2sSRJUme0mjS4qvYfdyEDXEnTO9dvFQv32EmSJK0obVeAIMldgecDh9AErB8BnwNeX1XfH0t1zbNxBy2w/UDga2M6pyQtW2s2bZ3YubadMLXFgaQVpe08cwcAXwGeRbPqw5eAa4FnA19Oco8x1fdx4OAkv5gEuDe58IN6bUNLsiHJ5tnZ2ZEUKEmSNE1tpyZ5NTALHFBVD6mqx1bVQ2jmdJvttQ8lyZFJjgTu29v0iN62Q+bt9jZgG3BykiOSHE4zuvU7wFuHPSdAVW2pqo0zMzM78nVJkqQlpe1t1ocAT6+qbfM3VtUlSY4HTtqBc/dPFjx3jNOBQ3vHvzbJQ4E3AO+jGehwKvCcqrpmB84pSZK0rLQNc7sCVw9ou7rXPpSqajsC9ds0kwFLkiSpT9vbrF8GnpnkZvunmdTtT3vtneAzc5IkaTlp2zP3CuBfgP9M8iHgezTLeB1Fs/B9Z4YsVdUWYMv69euPmXYtkiRJO6vtPHOnJHkk8CrgxTTPrhXNMl6PrKpPja9ESZIkDdJ6nrmqOgU4Jcnu9CbtdUF7SZKk6Wod5ub0ApwhTpIkaQkYZgWItcDRwGpgt77mqqqnjrKwcUmyAdiwbt26aZciSZK001qFuSRH0MwLdyvgMuD6vl06s2i9AyAkSdJy0rZn7lU067A+rqouH185kiRJGkbbMLcWeL5BTpIkaWlpO2nwBcCdxlmIJEmShtc2zL0QOK43CEKSJElLxMDbrEk+37fpTjQrQFwIXNHXVlV1yKiLGwdHs0qSpOVksZ65m4Ab572+DnwBuLxv+429fTuhqrZU1caZmZlplyJJkrTTBvbMVdWhE6xDkiRJO2DoFSAkSWpjzaatEzvXthMOm9i5pKVmmBUg7gA8F3gAcDfguzS3Xf93Vf14POVJkiRpMa1Gsyb5LeBC4EU0S3l9rff3OOAbSX5jbBVKkiRpoLY9c38H/AhYX1WXzG1MsgY4BTgROHTEtUmSJGk72s4zdz/gpfODHEBVbQNeBtx/xHWNTZINSTbPzs5OuxRJkqSd1jbM/Qi4fkDbT3vtneDUJJIkaTlpG+beArwgyW7zNya5LfDnwJtHXZgkSZK2r+0zc7sD+wHfTvIJ4AfAXsAfAj8B9kjyit6+VVUvG3mlkiRJuoW2Ye64ee+fsED7i+e9L5rn6CRJkjRmrcJcVbW9HStJkqQJMqRJkiR1mGFOkiSpwwxzkiRJHbbiwpyTBkuSpOVkxYU5Jw2WJEnLycAwl+SjSdb13j8hyZ0mV5YkSZLaWKxn7gjgjr337wLuPv5yJEmSNIzFwtwPgAf03odmMmBJkiQtIYuFuX8E3pDkRpogd1aSGwe8bphMuZIkSZpvsRUgngucCRxIszzXu4HvTqAmSZIktTQwzFVVAR8GSPIk4I1V9ZUJ1SVJkqQW2q7Nuv+4C5EkSdLwWs8zl+SuSV6X5N+TfDPJl5K8Jsne4yxQkiRJg7UKc0kOAL4CPAu4BvgScC3wbODLSe4xtgolSZI0UKvbrMCrgVng/lW1bW5jkv2AT/XaHz3y6sYgyQZgw7p166ZdiiRJ0k5LM85hOzslPwaeXlUfXKDtscBJVbVqDPWNzfr16+vss88e+3nWbNo69nNIkiZj2wmHTbsErVBJzqmq9Qu1tX1mblfg6gFtV/faJUmSNGFtw9yXgWcmudn+SQL8aa9dkiRJE9b2mblXAP8C/GeSDwHfA/YGjgLuAdjvLEmSNAVt55k7JckjgVcBL+aXa7WeAzyyqj41vhIlSZI0SNueOarqFOCUJLsDq4Arq+q6sVUmSZKk7Wod5ub0ApwhTpIkaQlovQKEJEmSlh7DnCRJUocZ5iRJkjrMMCdJktRh2w1zSXZN8v+SPHwSBUmSJKm97Y5mraqfJdkfuGEC9UiStGRNcr1t14FVW21vs34asGdOkiRpiWk7z9yJwPuT7AJ8jGY5r5q/Q1V9a8S1jUWSDcCGdevWTbsUSZKknda2Z+50YF/geb333wAu7Ht1QlVtqaqNMzMz0y5FkiRpp7XtmXvyWKuQJEnSDmkV5qrqPeMuRJIkScMbap65JLdKcu8khyTZY1xFSZIkqZ3WYS7JnwHfB84FTgPu2dv+sSTPGk95kiRJWkyrMJfkGOCNNCNZjwYyr/kM4DGjL02SJEnb07Zn7nnA31bVRuCf+9ouoNdLJ0mSpMlqG+b2Bz45oO1a4A6jKUeSJEnDaBvmfgisGdB2T+C7I6lGkiRJQ2kb5rYAf5lk7bxtleTOwHNpnqWTJEnShLUNcy8Brge+CnyGZimvvwP+E7gReMVYqpMkSdKiWoW5qvoRsB74G+A2wDdpJhx+E/CAqpodW4WSJEkaqO1yXlTV1cArey9JkiQtAa3DHECSXwXuDdwNuBQ4v6quGkdhkiRJ2r7WYS7JXwLPB27HLycNvjrJa6vqVeMoTpIkSYtrFeaSvBx4KfB24IPAD4C9gMcCL0+yS1UdP64iJUmStLC2PXPH0KwA8YJ5284HTksyC2wEjh9xbZIkSdqOtlOTzDB4BYhTeu2SJEmasLZh7t+A+w1ou1+vXZIkSRM2MMwludXcC3gW8JQkL0iyJslte39fCDwFeMakCh5Q63FJvp7kpiR/NM1aJEmSJmmxZ+ZuoFnpYU6AE3ov+rafu51jjdupwIeAd0yxBkmSpIlbLIC9gpuHuZFJsg/wFzSrSvwWcFtg/6ratsC++wJvAB5GExw/Azynqr49t09V/Vtv33GUK0mStGQNDHNjnmpkHXA0cA5wBvDwhXZKsjtwGs26sE+kCZevAj6b5Der6tox1ihJkrTkTevW6Oerai+AJE9jQJijmRJlLXDPqrqot/+5wIXAscDrJ1CrJEnSkjXMChC/DhwJ7Avs1tdcVfXEtseqqpta7no4cNZckOt99+IkZwJHYJiTJEkrXNsVIJ4AvJPmNudlwM/6dhnLs3XAQcDJC2w/HzhqTOeUJEnqjLbzzL2UJlTtWVV3q6r9+15rx1TfHYErF9h+BbBq7kOSlyS5FHgA8PYklybZu/9LSTYmOTvJ2ZdffvmYSpYkSZqctmFub+CkqvrxOIsZYKFev5sNW62qV1XVPlX1K1V1597779/iQFWbq2p9Va3fc889x1awJEnSpLQNc2cCvz7OQga4kqZ3rt8qFu6xkyRJWlHaDoB4BvDRJD8CPsUCQWqIQQ3DOJ/mubl+BwJfG8P5JEmSOqVtz9ylwH8A76cZAPHzvlf/gIhR+ThwcJJfPJOXZA3woF7b0JJsSLJ5dnZ2JAVKkiRNU9ueubcBfwx8DLiAEYS3JEf23t639/cRSS4HLq+q0+ed9xnAyUleQvP83CuB7wBv3ZHzVtUWYMv69euP2eHiJUmSloi2Ye4I4AVV9cYRnvvDfZ9P6v09HTgUoKquTfJQmuW83kcz8OFUmuW8rhlhLZIkSZ3UNsxdy4ifUauqVgup9tZgfcwozy1JkrRctH1m7l3A/xxnIZIkSRpe2565S4DHJvk0cAoLj2Z95ygLG5ckG4AN69atm3YpkiQNtGbT1omda9sJh03sXBq9tmHuLb2/+wG/t0B70Sz3teQ5AEKSJC0nbcPc/mOtQpIkSTukVZirqkvGXYgkSZKG13YAhCRJkpagVj1zSS5m4QXvf6Gq1i7WvlQ4AEKSJC0nbZ+ZO51bhrk7AQ8ErgFOG2VR4+QACEmStJy0fWbuSQttT3IHmqlKPjPCmiRJktTSTj0zV1U/Bl4L/OVoypEkSdIwRjEA4qfAPiM4jiRJkobU9pm5W0iyC3Bv4Hjg/FEVJEmSpPbajma9icGjWa8COrMOiKNZJUm6uUktHeayYePRtmfuFdwyzP2UZs3Wf62q2ZFWNUaOZpUkSctJ29Gsx4+5DkmSJO0AV4CQJEnqsIE9c0mGmm6kql6x8+VIkiRpGIvdZj2+xffnP0dnmJMkSZqwxW6z3mY7r/sBnwICXDTeMiVJkrSQgWGuqm5c6AWsBd4P/BtwILCx97cTkmxIsnl2tjMDcCVJkgZqPQAiyb5J3k4zQfBDgT8H7lFVb++FvE6oqi1VtXFmZmbapUiSJO207U5NkuQuwItpeuB+SvNs3Buq6tox1yZJkqTtWGw06wzwF8AzaZ6LeyPw6qq6ckK1SZIkaTsW65m7GJihGeTwKuB7wKokqxbauaq+NfryJEmStJjFwtwden9/H3h4i2PdeufLkSRJ0jAWC3NPnlgVkiRJ2iEDw1xVvWeShUiSJGl4K25tVueZkyRJy8mKC3POMydJkpaTFRfmJEmSlhPDnCRJUocZ5iRJkjrMMCdJktRhhjlJkqQOM8xJkiR1mGFOkiSpwwxzkiRJHWaYkyRJ6rAVF+ZczkuSJC0nKy7MuZyXJElaTlZcmJMkSVpODHOSJEkdZpiTJEnqMMOcJElShxnmJEmSOswwJ0mS1GGGOUmSpA4zzEmSJHWYYU6SJKnDDHOSJEkdZpiTJEnqMMOcJElShxnmJEmSOswwJ0mS1GErLswl2ZBk8+zs7LRLkSRJ2mkrLsxV1Zaq2jgzMzPtUiRJknbaigtzkiRJy4lhTpIkqcMMc5IkSR1mmJMkSeoww5wkSVKHGeYkSZI6zDAnSZLUYYY5SZKkDjPMSZIkdZhhTpIkqcMMc5IkSR1mmJMkSeoww5wkSVKHGeYkSZI6zDAnSZLUYYY5SZKkDjPMSZIkdZhhTpIkqcM6H+aS3D3J/03yjST/kWT9tGuSJEmalM6HOeDvgXdX1QHAC4EPJMmUa5IkSZqIiYe5JPskOTHJF5Ncl6SSrBmw775JPpJkNslVST6aZPW89j2Bg4H3AFTVp3tN9x3zz5AkSVoSptEztw44GrgSOGPQTkl2B04D7gU8EXg8cA/gs0n26O22Gvivqvr5vK9e0tsuSZK07O0yhXN+vqr2AkjyNODhA/Y7BlgL3LOqLurtfy5wIXAs8PoB3/MWqyRJWjEm3jNXVTe13PVw4Ky5INf77sXAmcARvU3fBn4tyW3mfW+/3nZJkqRlbykPgDgI+OoC288HDgSoqsuBLwFPAkjyMJqeuXMWOmCSjUnOTnL25ZdfPo6aJUmSJmoph7k70jxX1+8KYNW8z08HnpzkG8BrgcdVVS10wKraXFXrq2r9nnvuOfKCJUmSJm0az8wNY6FQdrNn4qrqQuCBkylHkiRpaVnKPXNX0vTO9VvFwj12kiRJK85SDnPn0zw31+9A4GsTrkWSJGlJWsph7uPAwUnWzm3oTS78oF7bDkmyIcnm2dnZnS5QkiRp2qYS5pIcmeRIfrlSwyN62w6Zt9vbgG3AyUmOSHI4cDLwHeCtO3ruqtpSVRtnZmZ29BCSJElLxrQGQHy47/NJvb+nA4cCVNW1SR4KvAF4H83Ah1OB54RTEQUAAApRSURBVFTVNROqU5IkaUmbSpirqlarNFTVt4HHjLkcSZKkzlrKz8yNhc/MSZKk5WTFhTmfmZMkScvJigtzkiRJy4lhTpIkqcMMc5IkSR224sKcAyAkSdJysuLCnAMgJEnScrLiwpwkSdJyYpiTJEnqMMOcJElShxnmJEmSOmzFhTlHs0qSpOVkxYU5R7NKkqTlZMWFOUmSpOXEMCdJktRhhjlJkqQOM8xJkiR1mGFOkiSpw1JV065hopJsADYAfwxcBgw7R8nMkN+5M/DDIc+xkg37z3cpmGbNkzj3qM8xiuPtzDF25LvDfMdrfnhdu+6nXe+4zz+O40/zuh/3NQ+Tue73q6o9F2ypqhX7AjaP+zvA2dP+nV167cj/JtN+TbPmSZx71OcYxfF25hjjvu695qfz/4mVVO+4zz+O40/zul8J/61f6bdZt0zoO2qvi/98p1nzJM496nOM4ng7cwyv+6Wna/98p13vuM8/juNP87pf9tf8irvNOmlJzq6q9dOuQ9JkeM1LK8+0r/uV3jM3CZunXYCkifKal1aeqV739sxJkiR1mD1zkiRJHWaYkyRJ6rBdpl3ASpZkG3A98JPepjdV1dunV5GkcUuyK/Aa4DDgZ8AlVfWH061K0rgk+TXgE/M27QHsD9ylqq4YxTkMc9P3x1X15WkXIWli/hrYFbhnVd2U5K7TLkjS+FTVfwH3mfucZBPwwFEFOfA261CS7JPkxCRfTHJdkkqyZsC++yb5SJLZJFcl+WiS1ZOtWNLOGuV1n2R3YCOwqapuAqiq703id0hqb8z/vX8K8I5R1muYG8464GjgSuCMQTv1/oV9GnAv4InA44F7AJ9Nskff7u9Ncl6S9ya523jKlrQTRnndr+sdZ1OSf0/yhSRHjLN4STtkHP+9J8nvArcHto6yWG+zDufzVbUXQJKnAQ8fsN8xwFqa2ygX9fY/F7gQOBZ4fW+/Q6rqkiS7AMcBHwEeMMb6JQ1vlNf9bYDVwDer6rgk9wLOSPLVqvrmmH+HpPZG/d/7OU8F3lNVN4yyWHvmhjB3W6SFw4Gz5v6H7X33YuBM4Ih52y7p/b0BeAPw20luM7qKJe2sEV/3lwAFvK/XfgHwZeC/jaxgSTtt1P+9B0jyq8BjGPEtVjDMjctBwFcX2H4+cCBAkj2S3GFe2+OAr1bVzydQn6TR2+51X1U/BD4J/AFAb/DDvYHzJlSjpNHa7nU/z2OBc6rqwlEX4W3W8bgjzX32flcAq3rv9wL+KcmtgQDfAY6aTHmSxqDNdQ/wJ8A7kvwVcBPw51X19QnUJ2n02l730NxiPXEcRRjmxmehddLyi8aqb+GtFWm5WfS6B6iqbcDvTaQaSZOw3eseoKruP64CvM06HlfSpPV+q1g4wUvqPq97aeVZEte9YW48zqe5j97vQOBrE65F0mR43Usrz5K47g1z4/Fx4OAka+c29CYbfFCvTdLy43UvrTxL4rpP1UK3ejVIkiN7b38PeDrwp8DlwOVVdXpvnz2Ar9CsufoSmvvpr6SZKPA3q+qaSdctacd53UsrT5eue8PckJIM+gd2elUdOm+/1TRzxz2M5kHIU4Hn9B5+ltQhXvfSytOl694wJ0mS1GE+MydJktRhhjlJkqQOM8xJkiR1mGFOkiSpwwxzkiRJHWaYkyRJ6jDDnCRJUocZ5iQtGUmelKQGvH487fqmIckDknwwyaVJfpbkqiT/nuSVSe465LGe2/tnef9F9jkjycVJsvPVS5qEXaZdgCQt4Cjg0r5tN0yjkGlK8nzgtcBnaZYK+hZwO+CBwEZgPfCIIQ75AeA1wOOBLy1wvv1p1pR8ZTmjvNQZhjlJS9GXq+qiaRaQ5Feq6vopnv8hNEHujVX13L7mTyT5G5rQ21pVXZbkFOB/JHleVf28b5cn0CxH9N4drVvS5HmbVVLnzLsde3CSD/RuPf5Xkr9LslvfvrsneXXv1uHPen9fnORW8/Y5tHe8Ryd5W5LLgR/Ma39skguS/DTJeUkOT/K5JJ/rte/dO/azF6j1+CTXJVk15M/8C+CHvb+3UFXXVtW7h/2twHuAO7Nwj97/As6sqm8OWaukKTLMSVqKbp1kl77XQv++eh/wTeDRwFuAPwNeNNeYZBfgk8DTgDfSBJi3Ay+l6fXqdyJNz9TjgSf1jvEwmtuTFwCPAV4H/G/ggLkvVdX3gY8Bx84/WJJbA08F/rGqrmz743t1HwJ8uqp+NsR32vzWLcCVvd84//sPBNbRhD1JHeJtVklL0QULbNsKPLJv2z9U1ct67z+T5LeBxwJz2x4L/A5wSFV9vrft1N6z/S9L8uqqumze8b5UVU/rO8fLga8Bj5p7jizJecA5wDfm7XcS8NkkD66qM3rbDgP2Af5+u7/45u4E7AZ8u7+hF9p+oarmniVs9Vur6vokHwSekuQOVTU3sOQJwE+BfxyyVklTZs+cpKXoUcD9+l7PWWC/rX2fzwNWz/v8B8AlwBfm9/IBnwJuAxzc9/1/nv+h17O2Hvin+QMCqur/ARfP37eqPkcT+ub3zh0LnFtVZw38pQtbcCRpkr2Bn89/zQt3w/zW9wC/AhzdO+7c+5OranbIWiVNmT1zkpair7YcAHFF3+fraULKnLsA+9EEn4Xcqe/z9/o+35kmCF3GLf1ggW1vAV7Xe3budjQB6xkDzr2YH9L0kq1eYPv9eu83AsfMa2v9W6vq35JcQNMbtxnYAKzCW6xSJxnmJC1nP6LpQTt6QPu2vs/903H8kCYc3WWB7+7FLW+Dvhf4G5rn7VYBP6F53m4oVXVDks8DD0uy69xzc71bqmcDJOm/5Tzsb30v8NdJ1tI8P/d9ml48SR1jmJO0nJ1CM2jhmqpa6Dm8RVXVjUnOBh6T5Ph5z8zdF9ifvjBXVVcl+QDN7dXb0TzTd9UO1v4a4NPAq4H+qUkWMuxvfR/wKprb148A/q6qbtzBWiVNkWFO0lJ0nyR3XmD72fMe+G/jA8CTaQYC/C3wFWBX4O7A4cAfVdV12znGy2h6rP45yWaaW6/H0/Rk3bTA/ifxy+fmbjHwIcmTgHcBD+k9Z7egqjo1ySbghCS/SdOTdjHNwIgDgP8BXMsvexOH+q1VdWmS02huAwdvsUqdZZiTtBR9eMD2PWlufbZSVT9P8vvAJppnzPanCUDfpBk8sd1pP6rq00keRxPq/hm4CHg+8JfALQYLVNW5Sb4BXNUbKNFvj97fhZ656z/Wa5KcCTwb+Gua3/9T4OvAh4C/n+tN28Hf+h7gv9NM0nze9uqRtDTFFVskaThJ9qEJdX9VVa/sazuAZmqVY6rqHQt89x+AO1TVH06kWEnLnmFOkhaR5LbA64HP0PQKrgVeSDMA4qCq+l5vv31oJt19ee/vuqr6yQLHuxQ4uqq+MJlfIGm58zarJC3uRmBv4E0003tcC5wBHDUX5HqeRnPr9RvA/1woyAFU1T7jLVfSSmPPnCRJUoe5AoQkSVKHGeYkSZI6zDAnSZLUYYY5SZKkDjPMSZIkdZhhTpIkqcP+PzJ2mZcuZZChAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 720x504 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "energy=photons['ENERGY']\n",
    "\n",
    "plt.figure(figsize=(10,7))\n",
    "plt.hist(energy,bins=np.logspace(5,7,20))\n",
    "plt.xscale('log')\n",
    "plt.yscale('log')\n",
    "plt.ylim(0.5,1e4)\n",
    "plt.tick_params(axis='both', which='major', labelsize=16)\n",
    "plt.xlabel('Energy, GeV',fontsize=16)\n",
    "plt.ylabel('Number of photons per bin',fontsize=16)\n",
    "\n",
    "plt.savefig('Spectrum.png',format='png',dpi=100)\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [],
   "source": [
    "ra=photons['RA']\n",
    "dec=photons['DEC']\n",
    "l=photons['L']\n",
    "b=photons['B']\n",
    "time=photons['TIME']\n",
    "\n",
    "l=(l>180)*(l-360)+(l<180)*l"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA1sAAAG3CAYAAABCJU/VAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOy9edg1RXnn/ylAeV0GFFwABRGNZtxAIW5x14gTN1zym2AmODqKOqOJyxh1ohNM1OjgGJxkzGiMicaJOlHjMhlFUWNc4wouSTSAaBTEIEpcEaR+f5xTz9NPPVXVVdXVffqc5/u5rvd637eXqrurq6tP39+67zLWWoQQQgghhBBCtGW/VRsghBBCCCGEEJuIPraEEEIIIYQQYgT0sSWEEEIIIYQQI6CPLSGEEEIIIYQYAX1sCSGEEEIIIcQIHLBqA+bM1c2Bdh/XWrUZQgghhBBCiJnyPb5zibX2+qF9+thKsI9rcSdz31WbIYQQQgghhJgpZ9k3fzW2T9MIhRBCCCGEEGIE9LElhBBi4zngZjflgJvddNVmCCGE2GPoY0sIIYQQQgghRkAxW0IIITaeK8/7yqpNEGKjccqxnjUhdiJlSwghhBBCCCFGQMqWEEKISZDnO80c22eONu1l5nw/5miTEC0Y+txJ2RJCCCGEEEKIEZCyJYQQYhLm5Pn2MxPOwbY52OAzR5v2Mrofm8OcVcpNpqbdh94jKVtCCCGEEEIIMQJStsRoyGsjhKhl7PFD45JYFXo3CtD9XxWraHcpW0IIIYQQQggxAvrYEkIIIYQQQogR0DRCMRqSyIUQMF1A8qZNz9q062lNTvvMsQ3nZIuYDtcXLz/6EAD2f9+nV2nOSkg9j3N6VlvbImVLCCGEEEIIIUbAWGtXbUMVxpifB34LOA7YB5wL/IG19jWdY/YBvwP8O+A6wNnAs6y1f5NTx0HmEHsnc9/WposAc/JoCCHEXkdjspg73eUb1E/bsolt667pa488HIBrXbT4/jn4dR8D4MwLz9k69rYvfxIAR735IiBPjTzLvvnT1toTQvvWUtkyxtwOOAu4GvB44BHAJ4E/NsY8qXPoHy/3/1fgQcBFwJnGmOOmtVgIIYQQQgix11hLZcsY8yLgPwOHWGu/39n+ccBaa+9ijDmWhZL1WGvtnyz3HwB8EfiStfYhffX4ypY8fUJsM/bzMMfnbY42CSHEXsJfkLyLxmZRSqv3+sYpW8DVgSuAH3nbv8v2NT1kecyb3E5r7ZXAG4ETjTEHTmCnEEIIIYQQYo+yrtkI/xR4EvA/jDEvBH4I/BJwX+BXl8fcGviKtfaH3rlfZPGxdvPlv7ORx2R6fnrf47f+feAFl+7Yp/uxWsZu/yHl13iqcs5Rn2uDFMLVoHZPo/ZZD+Z8fzYx1mnTmeI+reXHlrX2C8aYewF/CfzH5eYrgCdaa9+4/P8hwHcCp1/a2b8LY8ypwKkA+7hmK5OFEEIIIYQQe4y1/NgyxvwM8BYWytQTWUwnfCjwv4wxP7bW/m/AAKGANJMq21r7KuBVsIjZamm3KKeb+eXKFdpRSomHdFO8qXO6jhob5mD3XqF1W8+p782FkIdd7ZNmzu2jPj5P/Pui+1NOrG/Psc/X2rSWH1vAi1goWQ+y1l6x3PY+Y8yhwMuNMW9goWAdFTj3usu/Lw3sE0IIIYQQQogmrOvH1m2BczofWo5PAI8CbsBC9XqYMeaaXtzWrYCfsFiXS8yEOXowhlByHXvxmsdi0/qRyEP3ezdqk81i7PupsbMOtddwYm04x5lBtfWsazbCbwLHGWOu7m2/E/BjFqrVO1isw/VLbucy9fu/Bd5jrb18IluFEEIIIYQQe5B1Vbb+APgL4J3GmFewiNl6CHAy8HvW2p8AZxtj3gScYYy5GvAVFhkMbwr8ymrMzmPTPEwlWd72YqyTaMde7gt6HtKofcQqWJfsdHO2Lca6tG0LcsavvTjGrcu1ruXHlrX2zcaYXwSeBbwa2AecB/wn4JWdQx8DvBB4AXAd4BzgAdbaz0xrsRBCCCGEEGKvYaxVwr0YB5lD7J3MfVdthhAiQMqLN7WHby96FIUQ/ewl9aWEFmPm0Lbts0H3bsHY77c5vj9rbDrLvvnT1toTQvvWNWZLCCGEEEIIIWaNPraEEEIIIYQQYgQ0jTDBpk0jnItUK2l+s5hLvxJCCCGEWAWaRiiEEEIIIYQQE7OW2QhFHXNRHuZih2jDGIHJNWW1Kk+IFkjxbcemtOVY455j3dsnl03pD2In/n2d030eaouULSGEEEIIIYQYAcVsJZhLzNYmfd3Phamv46f3PX7r3/u/79OT1DkGe9mbWss6PTPrZKvYe8T6517ut1L2xRD8d3rrpVTW9dlU6nchhBBCCCGEWAOkbCVooWyt61f9OjFWG+d4fIQQm0vu2LKX1IU5vNPmYMOq2UttsJeudQg57bRX2nIVY7KULSGEEEIIIYSYGClbCeYSsxViDO/EXvLOCrGp7BXPpRiXMTLnjTX7QH19vugeib2ClC0hhBBCCCGEmBgpWwnmrGyVIM/S+KxDfJf6wfrT+h76/fbyow8B6jJmlijjUjqE2Dz0/Im9jJQtIYQQQgghhJgYKVsJNkXZEnubvrVpQvtEmFV7blddf5c52dKSOTwX69a2U9m7afVMhX89c+jjYzMkM59bF3MOa2LOsc+PNcNinfpiyGYpW0IIIYQQQggxMVK2EqxK2VI8w95mlfdpU/tI6Lo29VrH4vw3HAfAP97ztQCceMSxqzRHzJC9+EzVXPOQdmrRxpuobG1q3xtbRdr0dgOtsyWEEEIIIYQQG4s+toQQQgghhBBiBDSNMMGY0whbSZxjLDzZZdOk5b3AXKYFbOJUlRhzmvo5tS1z6W9d5mjTEFo/S3NOmKDgezGnBBWlTPksXXbKXQA4+HUfG7WevcLQe6dphEIIIYQQQggxMVK2EkyVIOOOZ/8UgE8ct//odYkF8naKqVi18rTuqL3aMae2nMtC8HNS9UKMMXtljtc5lE2+tlLUFtPRHT/efe7pUraEEEIIIYQQYkqkbCU44dh99hNnHsX9fvWxABx4waXAdDFWLdLD1p6/6cjzM08ufNZdt/59xEs+ukJLNgfFwGwmug9CiBQaI9rjK/J/9aG3bf17/8P/UcqWEEIIIYQQQkyJlK0EBx10Y3vCnZ7M9290dUAZX+bCGN4al/0IdmdAWtdFJNfRq5W6DzWsYxvMnZbPwxzvyybPCsiJHxxzfF3H7HLryJyfr01lThlpxXR0Z+P83YufLmVLCCGEEEIIIaZEylaCltkIczwPc8nOJITYxq1lApuvbstD2g615d5A97kNmv2xk5LfjGPYP/bafpuI1tkSQgghhBBCiImRspVgqnW2WrDuXoN1t3/dUfsLIdaNTRu3prqeTY5LbEnrGGJR18fX5TnfWGXLGPOLxpi/McZ83xjzL8aYTxlj7tPZf11jzKuNMZcYY35gjDnLGHPbVdoshBBCCCGE2BscsGoDajHGPAH4g+Wf32Hx4XgccM3lfgO8A7gp8BTgO8BzgA8YY46z1n59FXa3Zl2++PtYd/v7mLsncY42rSOb8jyOjdqpjNT4EYv1HTrmTK2y1NQz5+uqYSqb5njtc2SompWT/VPEWYcxKLuO0UoeEWPM0cAZwDOttWd0dp3Z+fdDgLsB97HWfmB53seArwC/AfzaJMYKIYQQQggh9iRrGbNljPlt4BnAodbaH0eO+WPgAdbaG3nbXwvcy1p7k7561ilmS4zPWN6Pscp1WfQO/cg3e8vvsyG1f1O8dZtyHSLMXNZ6mrvipOdgmzm2xRxtSrFu9u51tF5YmBzbNjFm627APwC/bIw5zxhzpTHmXGPMf+occ2vgC4FzvwgcZYy59hSGCiGEEEIIIfYm6/qxdQTwM8DpwIuB+wPvBf7AGPPry2MOYRGn5XPp8u/rhgo2xpy6TLTxqSu4vK3VQgghhBBCiD3Duk4j/DKLj61HWGvf2tn+LuD2wOHAl4FPWWtP9s59PPAq4Chr7T+l6jl432H2rkeesvX/OUqbOeRKs3NP4tCSvTAtTgxjFf2gZmHzFnZuUiDyOqP2yaPmORF5zKlt9TyMz6raeBN/b27iNMJvL/9+r7f9PcANWXxsXcpC3fJxilZI9RJCCCGEEEKIJqxlNkIWcVd3Dmw3y7+vWh5z/8AxtwK+Zq39fl8l9vKfDEoqMAU5NuTa11rliZ07J49GyMYxbKpJ3dyqzqmD8OfwXLRgFfZPlf7aMbUqNlV/Wtc+uG72rpoW7TX2+2hO77sc5pQAquX9nUPbz0k1nMqGWPu3/p3Tutyq+s+NH7euytZfLv8+0dt+IvB1a+03WayxdSNjzD3dTmPMQcCDl/uEEEIIIYQQYjTWNWbLAO8DjgV+EzgfeCTweOAx1to/NcbsB3wYOBJ4JtuLGt8OOLYvXguU+n2u+CmcYwsHdrftJea8kGKNbesUk1SC75GDedknpqNvFkArVa+kvLmkyq8ldzmLFHN8HqdO/d+qvjmOwVOxl699L5GK2VrLaYTWWmuMOQn4XeD5LOKw/gH4FWvtny+PucoY8yDgpcArgH3Ax4B753xoCSGEEEIIIcQQ1lLZmgqnbI3hlVgX9UUeGdGH+sj4tFT3UuVowdwwLe1uNfav2yLr61L/FKw6++fcF9bOraf176g59L29FKPqM8frKIkJ28RshEIIIYQQQggxa6RsJYgpW1PFmeQQm1cvj09d2aXlz6Fd+hiSrWfdrm/Odoo85qzYzGnsH1L/XlxbZ1XZcTfxXdyCOffTdWGd2mfOtrayTcqWEEIIIYQQQkxMkbK1zPC3n7X2ys62E4HbAO+31n62vYmrI5aNcA5ewS+/4GAAjjn57OLza5QNxxy9Ei0o8T6ONd997KyKmxaPM2cFZBNsmJq5xIGMxVjZ/XIUG58596uxs9+tapyask4h1olVv++mULZKsxG+AbgcOAXAGPNEFpn+AK4wxjzQWntWtaVCCCGEEEIIsSGUKltfBZ5lrX3j8v/nsVjv6hnAq4DDrLX3HsPQVTBmNsIUl51yFwAOft3HqsuQN609JesitcgopHtYzzqtzbUu93nV3sccfBudmgT5itK63A+Hs/fyow8B2ilnm640zplQH9zUdlrFOOtYh/dDa+Zs9xxtK7GpZczWDYBvABhjbg7cFPgDa+33gD8BbltYnhBCCCGEEEJsJKXK1sXAE6y1bzPGPA54kbX2Bst99wXebq299jimTk8sZivFkNXX3bF/9aG3AXDiEcdm19MCp6jBMFVN7C3m6I2aik2LG9sLGcKmXjcxJ+615j3R946Zozoylmq46uuakk0bc8am9bpVm9pOc2aObR6KwW0Zs/VR4NnGmCuBpwL/r7Pv5sDXC8sTQgghhBBCiI2kdBrhbwCHAO8A9gGndfb9W0ByiBBCCCGEEEJQuaixMeZQa+23vW23Bb5prf3nVsatmppphJvCHGXbubCKAPo534+Wtq1bcoLW7JWkBKuuP5ch08Jb1CfKGDp1co7tH5uGqqluO5nj2DlHmzah3la0tr/lNEIA/A+t5bbP15QlhBBCCCGEEJtIsbJljLk98DzgHsB1gDtaaz9jjHkR8DfW2ne3N3M1rIOy5ZJaHPqRbwJ1ge1z8k6sylu3SYrKqhZNFuvDnJ75lqzaw5tTd81zl5OAo0Ua66HlpMr1meo9tKl9vYZNa4vadO5TK9ZTsanXNUdCbdss9bsx5m4s4rJ+Fvhz7/yrgCeWmSuEEEIIIYQQm0lp6vcPA98GTgL2B34CnLBUth4OnGGtPWoUS1fAwfsOs3c98pQiD9ymehI29bqmosbzndPWofSjNeXo/k7Purd5TZr4qcfOmkWNa2gVK9R6YeIhxFSpVim0S45pcU4LNBtgHoTu/1z60VRjjlgw5B7OOWbrDsDDrbXWGON/pV0CXL/GQCGEEEIIIYTYNEqVrUuBx1lr32qM2R+4gm1l698CL7fWHjaSrZPTImZrag+c86oceMGlW9tida+bZz2l4tSybou4zsmWqVi1F7vGezp2jM26MPW1ThWXWBKLVFNuzgLFQ0jdlzHG2aE2rSNjtWPrdhpi55zv2aoyh9bGkeWWvylsokrcLGYL+DDw1OWHlsN9rf0H4P0V9gkhhBBCCCHExlGqbB0LfAS4AHgzi6yEvw8cCxwP/Jy19kvtzVwN65CNcGzm5F2Zky0tmeN1deedw/CYsJbnjsUcbYqR4xWMrc+TOnaV93IuMRcp/JkDQ9qre/4Ya9VNtbZPlzk+O3Po22OVt07UZMprmWlzanW9y16833uRZsqWtfYcFinfLwZ+EzDAk5e777lJH1pCCCGEEEIIMYTidba2TjRmH3AI8F1r7Q+bWjUT+pQtP4MU5M97nuN81a6aMcRzW8Kc4nHmXP7U9q76vrSqeyqv9lSqTKssXLn1dBli56pjt0psaL0OVqmduftr+qdPKrYtFsszZwVyju/VoYwdG5Rb317KsjfVb7u9oHS2fBfMvQ+2zEa4hbX2x8CF1VYJIYQQQgghxAZTrGwZY27PIlbrHsB1gDsusxG+CPgba+2725u5GsaM2dK83nrOf8NxANzk1Ys8LavKmgXr5blqnYGsz+M6h2sewljXofXPwrRUnsaOI6spv1V2ur3UJ0oZOx5uFeSqn6FjVq0wrwvrNtNFzJNmMVvGmLsBHwN+Fvhz7/yrgCfWGimEEEIIIYQQm0RpNsIPA98GTgL2B37C9jpbDwfOsNYeNYqlK2AdsxFOpZhtiqdmU64jxaoyOm1i/MQcGaIITR1L81cfehsAD7z7SZPW27qudYvJnAutVck5ZBosUY+mft7mkIExt56aGSM1MVXrhmZCrI/y2DJm6w7Aw6211hjjf6VdAly/xkAhhBBCCCGE2DRKFzX+MXDNyL7DgcuGmSOEEEIIIYQQm0HpNMJ3sEiKce/lpiuA4621nzXGvAe4xFr7qPZmroa5TSNsNUUwJd+XlnvZKXfZ+vfBr/tYsS0lKNh39Yw9NXATp31tSiKRmhTUfc/sutzvWLkl115j5xzTrG9Kv52qnpL3a+upjX31DWXqtPR99Y5Zzhz7fcvpomOPza3br1XCoaF02+nd557ebBrh84CPAOcAbwYs8GhjzMuA44Gfq7JWCCGEEEIIITaMmtTvdwBOZ5H6fX8WWQg/BDzdWvvZ5hbm2fRu4ETghdba53a2X3dp60nANVhkUnyatfbzOeXOTdnqMlXSg01dTHcuXpG9Rszjum5ewzl5RKfyYvv1zele1TB1uzlyFuZskTygZkHn3PJWzaZdzxxZ1fMxFpuYsGnVyvcqFMYapvq913RRY2vtZ4D7GmP2AYcA37XW/nCgjdUYY04Gjg1sN8A7gJsCTwG+AzwH+IAx5jhr7dcnNVQIIYQQQgixp8hWtowxVwfeBPyetfZvRrUqE2PMdYB/AJ7GYt2vLWXLGPNQ4G3Afay1H1huOxj4CvB6a+2v9ZU/lbI1Zw+S7xFYlRd9irpidc8hXmKMcucQL9Pi2BylYNMoWcg0p5w5jTktGBpL1aLuWFxXzUyCmrT9c4q5mGM/a60ATsW6P99T2zSHGTw+7p114AWXTlovjPO7Zg6q4Rz6epNFja21PwHuV3LOBPw34IvW2jcE9j0EuNB9aAFYay8D3gk8dCL7hBBCCCGEEHuU0myE7wb+2lr74vFMyrblbsBZwLHW2i8t1/3qKlsfBy6z1p7onfcbwEuAf2Wt/X6qjpiyNcSzNLZaMlacSQuvpqNkIcKxMznNwRvSghrPd0m5694+U1NyP8a6d6Hy/XrWgZwxzdHaez16/3//jRd/32f4rPapMj3W2OJvH8umOal5JUwVUzKVyjP0PuyV986QMan0vDFsGavukqyvNeWX7s+lZczWM4C3GWO+z2KK3kUsMhJuYa29qsrKAowxVwNeCbzUWvulyGGHABcEtl+6/Pu6wK6PLWPMqcCpAPuiS4oJIYQQQgghRJpSZct9SMVOstba4qQbpRhjngs8Fri1tfZHy22+svWPwKestSd75z4eeBVwlLX2n1L1+MrW1Bns9op3J0XXo+EUsTHmOa+rR7QVuX17Ez2Wdzz7pwD8zg2+AMAD734SkFYIVr12z5CyWpc3lZd8lR5XR03MQslsAz+Wo6/e3HpCduaWW1NPi3NaM5Xatm5jwzpnHJ5DrFANc3geYkzdHzaJlsrWbxP/0JoEY8xRwG8CjwMONMYc2Nl94DJpxvdYKFiHBIq47vLv74xqqBBCCCGEEGJPU7zO1qoxxtwL+EDPYbcHfg24v7X2xt75fwrc21p7k766fGVrVR6Zdc8+5DMHb9QqvTc13vK+Y3PrXiemukdOMb34hIXf5oiXfHSU+qZiqvGiJEaopQq2ivEj9Cx26w8p8ENmQJRkK5s6s1mruOA5sg5Z11ow9+tYVT8aq945PRc1tlx2yl0AOPh1H5uk/rFjcWPnDn0umq6zNQPOBu4d2P4B4PXAHwPnslhj6zHGmHtaaz8IYIw5CHgwizTxQgghhBBCCDEapTFb70/svgq4DPg08MfW2osH2lZEIGZrP+DDwJHAM9le1Ph2LDIYJuO1YLp1tnzm5AWZippYghoPfkn5OWWl5v8PrafUlpZ1zjEuY2jmptJzW9vSsvw5rS02h/EqJy4k1hdaK0OrzkA71Zgwd3VkLzKGl791v62pZ+oxZg5jWg4t7/cqGWPW2NjKeyjOvck6W0sMcEvgXsBNgH3Lv+8F/GvgpsDzgC8YY25VWHZTllkRHwS8F3gF8JfAT1lMIez90BJCCCGEEEKIIZQqWw8GzgAeaa39bGf78cD/AZ7OQtl6D/Ala+3D2po7LUOUrVUpA2PFa8wpbmxq9agFJfEsq1ZNpi5nKPKwlzPnsaYvM2ZOdshYjFXKhpQa5radd8adAbjZUz/ecxXx8nPWGfTbINXHWyhzY2XcbB17MUfFXQxjLCW2NXutH7WaPeGPey3KgrpZRX2/tcaM2SpVtl4AnNb90AKw1n4aeD7wAmvt14HTgXsUWyqEEEIIIYQQG0Lpx9YtgEsi+/4ZuPny3+cB16o1SgghhBBCCCHWndJphH8P/J219hGBfX8J/Ky19l8bYx4O/L619kbtTJ2eVSXIqKG1xF2zgPOcpw/2LWhZk2Sj1IYhrNtCmWNQMt0y1lc05bA9rQOSY9NOaqYGllCSVCN2biubppoqral7q2WObdv3rkztW5f7nTudbE73JafOsW2o+V1Ywjqk/E8d23pR49cbYz4HvAX4FnAD4BHAbYBHLY+7H/C3hWULIYQQQgghxMZQvKixMeYXWMRnHQ9cDbgC+BTwW9bas5bH7AN+aq29oq2503LwvsPsXY88pddrkOPVnFr1SdkypPyxAurn5GGtUcH6PNFDPd/rkPq25P6uSxKQlsS8gqtU26ZSAMe6H365LRffhPxkHaHritmSk4xibA+yb0urpBpT2Z3DnGxxTP0+mrNCEDu39vwh7MX30RyZ032o+c3SdFFja+17gfcu17G6HnDJMs1695gfl5YrhBBCCCGEEJtEsbK1l6iJ2RriuXrme98OwO8ec7uiOlvRTfXp0gj7lMzjLmHVnp/a62jhSV/1tZcwdtzac87/3Na/3XOw6lT/U/bxFrGFsTK7DCmnpdKVYtU2lpSf00dKZkCkjhujj0zFWCnmS+pusfzAUBscOX1lVffzwmfdFYAjXvLRpjYNeYamYpXPUov3XasYpBrWYRwai5ap3zHG3N4Y81ZjzCXGmCuNMXdYbn+RMeYBQ40VQgghhBBCiE2gNBvh3YCzgPOXfz8ZOMFa+xljzAuA21hrTxrF0hXglK2WnvVUxquWma5qMgflxH2tQ5aeKbMGrsrLLG/UNpuirg5lTvbH5rs7hsRmrlKpi11HyQKgKfv9faGFQGPH+Isa16pIue0/NHZ1TkrWkHOGMHYs6xxjMjcl899YZbRQpcb6LTeEOb2fxqKlsvVi4Ezg1sDTvX2fAe5Qbp4QQgghhBBCbB6lytYPgYdba99tjNmfRSZCp2zdAzjTWnuNkWydnLmssxXyCLSYU57yPo/qKX7/jbf+eeGZRwHb88JrmCrTWc4xq5qP7mLHYHf82Byzcg1hDh6yVfa5McusUd7nSEo9arGGYCulrmXbjr02YmsFYmo1pIZ1m0GwTs9qizjUdbjOOVDy+2OqmUEpm/yxbF2UzJbK1o+Ba0b2HQ5cVlieEEIIIYQQQmwkpcrWO4DrAPdebroCON5a+1ljzHtYpIF/VLSANaNlNsKar+SxMyJN7SVsfT193ubuvtg5/vZQeTle61h9OdRcR+zcmvrGYsrYuXVkXb3zYzFkXCqJTfHLzRmXasb1VNyET03Mhc/YytacmZN6vyltOoQ5ZRacEy3j/6ek5e/aOWR6jNWdmiGUU0bLdbaeB3wEOAd4M2CBRxtjXsZikeOfKyxPCCGEEEIIITaS4nW2lqneTwfuAewPXAV8CHi6tfazzS1cIb6yFfOejeXFWaUHYIw5sy3WpBrKOngdp4oTyMmKVhOv5pfVirHbZVPjKYbWm6O8+gy5Rjf25K7117XJt9Fl6oP+saxGGe9mIfTraUHOdYQyFnb3d+2dM2N7yaeKs0z1kdg9WvcxZx3eq1NREx8193YreVfmXtMqM2KORUtlC2vtZ4D7GmP2AYcA37XW/nCgjUIIIYQQQgixURQrW3uJg/cdZu965CmzXIdgVXFXsfoBzrzwHAAeePeTkraM7XEd27OUU/4qPYot51cPtaXFsX3nrJsHfwgtMnjVnu/wVe9UvEwre3PL8M/pKkKxdahSsVt9ipmvMuXaWUpJ/ONY96OGkliqOXimY9TEn7a6niFK7KrYRHWnhjlfW4t38JxjtlJ9cCwbBilbxpj3F9RlrbWrz5UuhBBCCCGEECsmZxrhfiwSYThuCRwGXABcDNwQOBq4CPhSW/OEEEIIIYQQYj0pTf1+EvBy4JestZ/obL8T8Cbg1621b29u5YoYsqjxqlK9zyGphsOftpNDKhh/1VL82NNCxjp3SGr5IaQCw0uY81SMGCVTGNalj49xH2pSnKfKaT3e5qSU99mfSa8AACAASURBVOv1E3zkPHd903NSKYn9Y887484A3PL3Lw7aFyLHthRTTQdfx7GgS5/9U0317tLSlk2ZNjon5tBOY4/9qw7HaEXLRY1/B3he90MLwFr7t8BpwAuqLBRCCCGEEEKIDaNU2fox8DBr7bsC+34ReIu19hoN7VspTtnK9W52960TU6fUHurJz/UOpmyoqackGUUOfV7AsRIalFCTAn5Iea28XevmEesjJ1lEi2suUV9qyEne4NebYkhK4pgC1VVk+5JqxGxvZX+rAPSaZ9Qx52co5347Wl9HrP+kKBn/HHNu/3VdTLpvaZt1e3+sk71+kiFo03/m0AYtla2vAE+I7HsCizguIYQQQgghhNjzlCpbvwz8b+DvgTeznSDjkcDPAr9irX3TCHauhCExW1MR+5q/8Fl33fr3ES/56Oj1h2zIOcd5N3I9yLX2lcQstPCMjK0WDjl2rLadYyxSi7iDHKV0bC9661g9v5ySdOEt+/aUcTljqEh+2aHyUox5zd36v/3zhwGrXVB+CFN7rafqg33HtbZhqpiwsWM+h6hYq2QO6stUxNp9DjNTZpn6vYu19o3GmEuA5wPPAa4GXAF8EjjRWvu+ocYKIYQQQgghxCZQvaixMWY/4HrAJdbaq5paNROmXtS4Bakv9r55ymPb5ChRZ1rP6/XLb6XuOMbO1lPiFSqxadVZsvx6xqprlRmQWqqTrc99zvmfA+B3j7kdMDyTZItYxpxYmL5jc54hv4wufcqlyxLYVYxi46y/EHIId+yXX3AwAMecfPauY9Y53qRVTNKQWL2UTatSnHKIXUcoU2WrWD8RJ6Zcr+K90RdLqntaT8lY00zZ6rL8wPpW7flCCCGEEEIIsclUK1t7gTFitlaZPWkMW8bKXjb1nN2x58h36fNwryIj1arm08fObWXLkPJbe4dbPg9TrmGW6zVNxSi1iDUradOSc0rG3SHrCKXU+pr4R7fNj8cK2TiG2j02XcXm0I98s9gW/5pzlMVYGUPqHToWxI5tHYu0rnEyq2Kq2KOp4uA2jVWo0C2zEQohhBBCCCGEyGAtlS1jzCOBk4ETgBsAXwPeCrzIWvu9znHXBU4HTgKuAXwMeJq19vM59bRQtqbyFjnm4IWvoWU8Syuvf41tLWgVl9UiVmFIZr6a9irJohmqZ8jaZVMpsauKCSuJV+uSa2drT2KfRx92x1ml+mCuEhCKgfFt8unWUzIm56qDofL7nr+x1rEJ2RSzLVVG33jRWr31KXlWQ+cMiZNx41xOpuAW7/iSc+cYI9RSIRpazpC6a2Kva8aPVV3f2HXnjAmt3sU1bbmJytZ/Bn4K/BfgAcAfAk8C3rtM3IExxgDvWO5/CvAIFtkTP2CMufEqjBZCCCGEEELsHdZV2bq+tfafvW2nAK8F7mutfb8x5qHA24D7WGs/sDzmYBYLM7/eWvtrffWsg7I1FX2erNC+lvV2GWv+c+yYFt6P7nXE1hZbV4Z41nM8+1PHq5XEy4zdF1vQqp6+dVO6jNEuJWpAyLa+tabG7nuh6/KfFde2Tl37weEGgKPefNEum1rE7LTy7OY+Q7XeZtceQ2K3WqsKpfXVljek7lWV0YohazX5lNyHFupeq3pi55SeN6TOknNzVc9VKo0t6V7Hu889fXxlyxjzX40xjzfG7GtVZgz/Q2vJJ5d/32j590OAC92H1vK8y4B3Ag8d10IhhBBCCCHEXqeZsmWMcWtt/TPwMmvtS5oUnF//E1lMJ/w5a+2njDEfBy6z1p7oHfcbwEuAf2Wt/X6qzDGyEc6BIXPAW5ffIs6kxrap5ji38A7DbjWsZM0en1XEFbXoc2PZ1Nf3cuI0Yue2OjZHgcgtI0SLZ3TofRnDKx+K74o9+zlKUd8aWgDfv9HVgWFrHrl6/LJy7RyDmnqnVq5zniVfIU31kZy6WjwXqWd4DAVzbPWlG3vrVNkhauRY54htWrxjxmKsd0yLMa3LKOtsBbgpcC3gbsA9GpbbizHmRsBvA2dZaz+13HwIcEHgcBdRfV1g18eWMeZU4FSAfVyzua1CCCGEEEKIvUGzjy1r7VeX//w74FWtyu3DGHNt4O3AlcBjuruAkGxnUuVZa1/F0v6DzCHrF9AmhBBCCCGEmAXFH1vGmGsAPw9cH/g68DFr7ZWtDcu0ZR+LjIPHAPe01n69s/tSFuqWz3WXf3+nt4J/dU1+eqfjowtP5jB2MGLJOX0pdmuk1BA1i0fGyLmOHPtj155T/pAg8popJCHctbkH7UBXfqK8WCKOnH4VC24dmk46ds9qbBpKX3mh/bmpjkPntpyG1S1/yDSgnL4XmwJWk2QmVWZf2+Ys1uuzo6177N+ashcoZ6ueZXlbU8+W/9+/U5Y73z+mpK3dM7v/+3a38dbz7NUTe85D22qmUZc8q885/3MAPPu0ZWILzw7fvlaEbDr/DccBcMzJZwO7E6z4ywcMqauWVFl940XOlMAhttac223jy91zu+LpfWMlOqqxYZVTHV2yGT9RUMm7N2dJmBpK3qc1+OWkkj01D68pOdgYc3sWHzdHsK0QfdcY83vA70750WWMuRrwFuCOwP0Ca2d9Ebh/4NRbAV/ri9cSQgghhBBCiCEUJcgwxnwUuCHwH4HPA4cDvww8Afgs8ABr7Y9GsNO3Yz/gjSwyDj7QWvu+wDEnAX8J3Mta+8HltoNYpH7/c2vtU/rqOXjfYfauR55SnZazy5yDNlPelhbe+FQZQ1MCp47LLbe03pwg71Sgvt8OLpg4FUhcokq6bWf92WsAeODdT4qWGyt/jPsy9JwW5abuXUwhGGJHqJyUVzBmm19uTbKFkAd/rPsd857GbAvh7HUp23PSfreYSRBSb3PGAId/jLP/2t/4CbDzPsTub8k1++fGys45J/e8WlvG8hyP9QwNmdUw5NzS8lqUP0VZq+S8M+689e+bPfXjWedMpYbVJlTyyXnvjaE85Rybk6SlZnmPFrR6Vlsuanw74DestWdaay+01n7aWvtMFmrRkcDzC8ur5X8CvwS8FPiBMebOnT9uweJ3AB8DXm+M+WVjzInLbQb4bxPZKYQQQgghhNijlCpbXwSeYa19d2Dfo4EXWWtvtPvMthhjLgBuEtn9fGvtacvjDmHxQXYSsI/Fx9fTrbXn5NTjUr/HvshDsUmxr+B1UwqmUpFaXGtJWtuaeKxYWbXlxs7JocRjfOaFi26eo2zVUNPXW6aQrVEYa8rN8XaVKB0t+kTOc5dS6vqusUTxaJ2WNxZrWKPQpRTlGKlzYjaUtG2XmAJYMhvAL6PGS9u1eUjMbcmY0LIPhoi1S0lZJTM7WsQl1pSfOjY2PoXuw6rimFrNfOljlTFbNf1oLJs2RamcC91+lVrUuDdma6kIfdJaeylwBvAkYNfHFvBPhGOLm2OtPTrzuEuBxy7/CCGEEEIIIcRk9Cpby8WKLYs1qz7JIhPh54DftNaevTxmP+C1wDWstY8c0+ApGbKocZ+3qMR7OlRRaeG5yonDcvNtP//rfwjAiUccu+NcR8q2sbLc+IylPLW2oaUtQzyjJeWPFRs2xLapKLFtSBbNoV7aGi//GPPpQ4rKxSccGKynRNnyt4fqrFFJXAxVX71d/Jit0EyImJrnH5eqs+YZKpmdkbsf8rKNjqGopNQwdw9D98G3uyZuc+qxbFUxXK3Kn+od3KKNQ+evarZSzm+isd/1UxOK7Y1lEmx9HSVjwtBFjW8J3AE4fvnnWsC/AR5gjDkPuBC4GXA58KBM+4UQQgghhBBioymK2do6yZhjWHx43aHz9yGAtdbu39TCFTJE2RqbmCI0RA3r0hcHkuON98sqiSVI1dNnQ9cLkuMpjtnk8Of+l8SBpNYG8uusaVvfI95l7BibPttSXuy+c7vHtojzGvuaazyHqYx9U3l0a9o2x+Pqk3Pf+1SiVF+P1ZfzrMbUpZBNNbFIfmbBbrZDRywGrGTMbOU9H6L8xWLPUu+lmtkZsTK6jBGXk4pta1lPyTtmKvUlR02fk9qWqqfkeR5STw01GaDHoGRWUasZSC1+G+XESvaVUXJOl6HK1i6stecD5wN/4bYZY27C4sNLCCGEEEIIIfY8VcrWXqFlzNYQajzGKXWnJK4hdsxY86Brjh2iLuSoeiXlOy92LANWqp4cYh6enHYakhFpaH8eMz5jbLWqxnOcc9/9PgN5GeUg7PnOjf/JYai3vE/dzokVKlHgh3g+/Wc2dR0+ITUstq5WTnxXSzU3hK80pTIPDlEp/LYNrRc2RF0oUcNiY3KNB7zm2KGqYQ1jqyMt++fUCle3rrnENg1VVKaOu5oq3q7EhrHGzBJarrMlhBBCCCGEECIDKVsJDt53mL3rkac0VXX8MkLltJjDnqqzJh6ghJr4MZ8hKlU3JiK1Hk7I5i59MR6p8vo84imbQtc3xFM5RA1L2VRzbM39jpVbohqmMunVxOO0jpMpPScnLtE/t8sYsS8l8VElY04qQ2lMmRkytoXa1q/PkVKEfNtS8ZuOWN/u1lsTI5RSFH361hIbErdbwljvyCHjU45N/v6UCl0Tm9KCVjE2fnklauHYSsTYmY1bze7xqZmBNHZMVd84NZUqWfIbuOS91wopW0IIIYQQQggxMVK2EvgxWznxGX3xBmPHl7Qoq1teTBEo8Q7mtFeJTWOoCylvtjsmlT0uVm6LuIfU+VOpniUqUgk519xyXn2OKpIqy9+W0yf6bGvlgSu5hzUe9Rb3IaQgO/rUqVA7xTIU1qg9fhxTSbld21IxYJAXb1ejbNQ8o2MpKa1nUbQcV3Pqzel7fbaMdV0197nVzIRcct7bU9mQKuP8NxwHwDEnn71j+5D45hQl74shTDXjYo60noVT0h5StoQQQgghhBBiYoo/towxtzfGvNUYc4kx5kpjzB2W219kjHlAexOFEEIIIYQQYv0omkZojLkbcBaLNbbOAp4MnGCt/Ywx5gXAbay1J41i6QpwCTJ8SmTFFskJhtI3JSJnyoKfwjd1bqy+kgV4SxYDriFnKoaz4eITDgS2p1KWLJqcWijVPzZ1XX1JR3L6k5vC4N/DLmO3bWwh0KlsqZn2WrLAaMqGloHJQ+vtSwRQszB4N4V9bJzISZ7j99OcPhObCliTCv5LT7nh1jG3/P2Li+3vqyfUn/xrrpmm6ogluCi1N3bO0Klufbakpu/2nRsi1ratKEl2MWTKYW69rY4ZMsUq532UmqrX97tp6D0cOzlHqR0lv6NaTXmrCR2Y47TEknrOvPAcAE484tjic3OObTmN8MXAmcCtgad7+z4D3KGwPCGEEEIIIYTYSA4oPP4OwMOttdYY40tilwDXb2PWPLCX/2THV+yWt8XzUKa8p1se0OX/Ux6/Em9Li8DmLe9voh5nw8HuGM/GVGIJ3xOwf6CMqLc847r6UhSHyPHeuG0/Xf7/Whft7Oohj7F/rsPdn8u77RQ9O17WAZGFUd3/cxJk+GpDKBFArE1bKXR4qbNjisFQWnrVdlzz4NK2ybGxxnNccu0xb3zoPsf6nDunq774492u5zrx3LnFgLds8/pMiP29cp39316qGt1y3bH+tbs+fvQ7r9hVfkwVzkmi4tezo48vj42pLlvqXkAN85ONuO1b19fZt3Vfl//PSdISI8ezntMnh/TTEsVmq1/6fTFBSYp8/13ltl/w4KsBcMsL6lW+ULlDFit35ChPNeNSSImN/a5JtW1s35DZAaFzx07A0UcrlarGBv+cGgUtdGzsuWhBqN+WjGUPvPvyt919l++UghlnQ39LlCpbPwauGdl3OHDZIGuEEEIIIYQQYkMojdl6B3Ad4N7LTVcAx1trP2uMeQ9wibX2Ue3NXA1+6vep6JtD3bcvtD+nnpqFUrtxGs6bVZMKNRbjURLTk7KpJpVrXwrn1LkO5yH9/o2uXmxTSVxfTur6kvilVNySX0/JfYi1T05sjX8/ajyAJX08p5ySGLrYIps1CnZJ30gt6NxXj19XiJJYw9T99mOmbvbUj+/YHurHrk841cqpCU6d6tYTU09j/apbp9+GqXsXu+ZQX+mLby2JbfTvYe0i0znvn9zrSJGj2saOLfGw59yHGDH1MOfY0Dk5z3Ouba3iTWrio0pitoYsZxBjrHignH5cspD6WDG3pZQ8j1PHWoXq9rf3nd93zpCYtthYHCorFbNVOo3wecBHgHOANwMWeLQx5mXA8cDPFZYnhBBCCCGEEBtJ8aLGy1TvpwP3APYHrgI+BDzdWvvZ5hauEJeNsEWWnhyvoL9vrIw5OXO/+zzSIW9YX7a7UD0liybnxhPVeI5TlChCvk2tlYIYOZ6rlAc/tz1S5/g2hJQtv5xUhsoWHtESZXRI+SVKbIvnOrUwbo2nu0Tl9p9DX12CbSU3dt8dXVv71KMcFSa1aHLfdYT6Yt+iyc7mG37q8l31+MeWeMkdqYx2qWfHJ9aGOapnTh+PKeH+dYXuT07/baGc5ZTh9+U52Nb3jIZmTcTKDVGiPvuUxM7FrqNGXS25HkeJjalyW2RIzLnffeeW1plrSytaKa+QHjdqbKlR+mL7u7RUtrDWfga4rzFmH3AI8F1r7Q9LyxFCCCGEEEKITaZY2dpL+MrWEA/A2B52v54aL1HJekI1XgtHrSIY856mFLqYF7jES+uXkYq/6vOed8vpyyqWsiEnhifmsUwRsylkW64i263bt8WtYXbUmy+K2jQktsM/Nse7mYrVq/G41sQLOvoU8hQlnuOp4nFK1oBySvkPDjfAtnqUWrfovDPuDGyvjxXKuPm1Rx4O7F43z19PL1S+fy9Dar6LOXM2+ITWI4v1uZI4kJi6HsJXHksU2RJ1Z8i5IRticY855KhvOWsRxsqLlTs0NqUkC2Fff805p4aa3xI1at5UrCJuyTH1tffNLsqxKec6/O0hYscOeYeGbBmrTQcpW8aY9xfUZa2102eUEEIIIYQQQoiZkTONcD8WiTActwQOAy4ALgZuCBwNXAR8qa15q8VfZytnPulYK553yyotzx3rr4Hhf+1feML2ClA/ftxxANzk1fsDZZ5EX4U5wIvtOLhju/NEu8xjfhnOCw1wxEvCc9UPvIAd1xfySu6Ky+i9inh53f9/bekV8pWZrXZatlt3vZ+t85fHuPW8/DXZurh7tX/C+xuzOxbbEfL6X+n93z82Fefg7quvMsB2driuVx92t1vXplifc+WXZDJM4T+jh7p6es/M8+DXxJzFYgz99dZSNjlC6+n552epkv7/E7FCMe/+Aak+6JXh+u8PlmPAVn2BcdD1+aPfuVyj64+W2Qjvs9327ln90eFX7bDpvK0MhpfvKAt2K9dbipBn685+bXfY5thW3XaPT77C5MoLjWW7FHbv2BxFy8Ukuevqjv2+h9t/b4RU760xLFrzghKvc/fYmNIee04gPosh5ZWPZcXNvRboPKOBff7YvHWfI9cDbL1DcsYj/3nw3zX+cdE6O9tr1vTMOdZv2x1jv/cuyaFk1kHfNbWavdQXb5dai2qIMlQSdxeLkW0Ra9Wl5Bj/ms/6s9cAcOIRx/ba0Eq1ahmjBxkfW9bae7l/G2NOAl4O3Nla+4nO9jsBb1ruE0IIIYQQQog9T+k6W58HTrfWvi6w798Dz7DW3radeavFX2erJK5o6jm/NfNgfdzcXYDff8L/AuB3L/jFRRmPv9qOc2uyfoXic/z28pWPnLnzjpBH14/TcHEfX33cwhd76LuuAYRXvE/FXfXhzg1lK4uRiq2KZcUqiVtyOJuuddH2s++3c0y9CMWbxGIhLk54y1P9x+Hq7lt3KUTI0w1l2SFz4g9qYg1z4rscfZk3Q9tK5tPn9O1YFstY5rbuMVHb3n/jxd/3+frWObFrz7mHPr4a093miI0tqXX6HG78cMp/N47zkjssnit/ra9U/Gbf2JJas64kK2HsnoWeR/f8unEiNRbHMpP6fSeU5dK/v36Wy669rp1jcXzdvu7GfF89Tz0fqYyt/rkl8WL++TVrLMaygHbvS01W2b4Y5ZJ4r5o2CRGbjROrv3tOn425+2KUxNCNzZDfmWPE6qX2lWTBrpkNkvM7t+aaS85JxWztl13jgp8B/jmy71vAzQvLE0IIIYQQQoiNpDT1+1eAJwDvCux7Aos4ro2lJk7K0SJbT+4+n12e4Yji9OPjtjP4v+DxjwE6Htyj999xbk6GO99LcUOO32VzzBPmPFmXZ6gKfgxDd70ft8/PPOY80vCTXXbE2jTl6YtlKeuqRz4x1ahrv1++m+N9/hsWMXXHnHz2zv3Ac9/7dmD7Hvr3KpT5z2/DXTEFy/93PYxfWsagXOOihQfZXWvKfr/tXP+yF+8DvNi95T1xCoGvTByc81x4/TP13MXue/c+ubiJvuevu93ZvdUuCa+drybs/76dHuPQc7d1Hctj+mIzQ/jZIUOq57cjMYxbnvZAXKLrT9+/0VIxWMZSXfanRy6Ou9nVtgvynp3FyiLb17elErNbhfGfpS3VtbPN945/zcvC5frGt//Nj7aOufY3do6Rrl2OOfljO23qKNc/OPzAHTZd7qlG3Xu3rYjtVKd8Jb5ru7Pzso6aBnDBg3dmQQzFLzlF6ODXLdrUxZG5Zwy22/uoC/KyHnaJxT9eHnr+exTlQzvb/HFoq9xELKkb51Lxb90yusdc7rVxSJk7wHsWc2PEusf6fcO1Tyh+ye2LPYfdcon83nBlhOKXfPv93wklakNofOrLJBlqpz5FqzZWKHZ+KoOvs3urDb2ycuvus2PsWVElmUJ9+pTGbnk+fruljhkLNxb8wMtI6xj6uzxF6cfW84H/bYz5AvBmthNkPBL4WeBXCssTQgghhBBCiI2k6GPLWvtGY8wlLD66ngNcDbgC+CRworX2fe1NFEIIIYQQQoj1o3pRY2PMfsD1gEustVc1tWomuEWNHX0BmLB7wcy+4PvUMf6xOdOAQhJ97jSH0FQ6N13tFs+9bJfdsXpSAeF95FyrH8TsT6kL2eeX6y88GgrG9qf0uGkpoek5/hTA2PYuoUDwLm7KD+wOBPenpoXauC+xQRdXzl3++yLJ6Meecccd9Z154TnAztSrsdT4X37BIjG26zPQ6Ru/uWjv5xz9/wB437/cGoDfucEXFuX/wwO3zrnwzKOA7VTdboqhPx2lm9jF4e6Vs8X1je6x/hQCRyrxgD+9qyRIOpYAoia9dGiaYmzqp7uOCx68PWXvGhftt+M6/Ge2m9zEtWV3GYbQud1pnn6/9Nsy9Fx835sW5xJNOFvdFL5QMpvYNbtpqhBOZgHbz5nrD/50VdieWugS6nz8xYsEQjf/P0/YYWPXTh93Pd3Fjvv6kz9edc+PjUvuGXPPD2w/Q9f7zOJaXbu7lMr3+9XH7rLXtd0/3vO1O47xE9fAzmmIXVJjf2yh89C46CfEcHT7qY8/jTDWJ7u4vuEnOUkteO3jj6/dPu7X7U/bDk1Vv9BbXsTdd3dPu/0pllgilXgn9k73n6mchAapaXF+fb5NOdPzhyTvyJkiFisjZP9YCRlKbMpNvlSTSKtV4recPuLXWbLsQEnSMFeeG/ccbmxLJdDKaYOWCTK2sNZeZa391tw/tIwxRxpj3myMucwY8y/GmLcaY47qP1MIIYQQQggh6ilN/f7+nkOstfa+PcdMhjHmmsA5wOXAc1msNvkC4JrA7ay1P0id71K/53iXHX0LoTlC6kiMkGfGT6OeSv+b6/0IpeVN2QDhdM+xlLshb5Sv0FywtcDoTo8ibHs1/UQYqfTFPr7neFfygk75fhBoTlpbX01yXkmXSh/gKa984o7rcTaFFCG/XEdMfYPdqo7z6Lu2dB7RLr56EEutHFIvfJxn3SVBAHjxaa8EtpcSOPNn/yp47vO+dZutf7/hr38e2O2Nd33EbQ8pgDGPccijuyvxxvLeue2h8mOknr9Y8HXonNg+X6nr2vkLT/8wAO992d12nOOrMrDdLk6dCD1vDtdf/D7STSQBO/utK/fke30EgM88YpGo1vVjZ6tTNGFb1Tz/4usB8Mv/eucY+sa/95PbxMdB1wd/+sJt9eW5f/QnwHbiGD+NeEi59scjX+1xNp/3uRvj8+qHvArY7vNOzX38356y69hjbnjJjmvfd/Y1ge3+5hRngPse9MXkdfjPMGwnP3KJaNz1hJQtX1k6+N//E7D9PPuKY9fOUIIVSCv7fh90pMZ+X23ZXjB6W+Vxz8pW8p3bLZYZ6PaJbpldXNv57ebGHNitVMYScaTS9rtr95Xm0LW7vueeC/c8hJ5rH1ePueGPgcisg55EFqGkGjFCyxuE9oX2d4+JKSgpxaPk91qfUpOTOCtVVslCy46+5T9aLPPTiiHKXUhprFlY27chtdyO/9vRH99zFLSYbQAfOOs5UWWrNEHGfiw+WLocCtySRUr4LxeWNzaPB44BbmmtPRfAGPM54B9ZZE982QptE0IIIYQQQmww1TFbOwox5mbA24CnWWvPGlxgI4wx7wP2WWt/3tv+QQBr7T1T5/uLGpcQS6daMqc5NRc1FaPlU7LoZV8ZjpCa5KsfMQ9fd7u/7ZnLtOX32rfol90YHuc9dl6IkkVP/VgUPz6j63n1vbO+VzZHifDn1TuvKmx7VmPxK93YGofzSPoeVt+b3SWm6qQ8ibG58KH+FYt78z3h3W2OmLLVxcXDnPv/LVSxn/ngo4Ftr/9WDE+nbNe2rg1dnJe7z101xnmEY3FvvkrSvSZfLfJjI7re8q7CANsq5+m/8FBgpyrpvHHOft8T7drYedW7uHZy6uDb3nB3AD7/638IbLcfbCspTm1x+DF0AH/944X9j3vHqcC2YuM/ozn31NngK0MAdzzkq8C2cuNscQqnuz5nT9cGx21f/iQg3Df8epwC6OKvHnj3k4BtRQS2VQRnU6zeodz52U/cYUsN7r77tkKZvf49KsE9s66PuD7jlJvu8iIOX9Vz/NGdXrf1bzcbwJ3vnoe/yndRIgAAIABJREFU+tDbora463Dluxi2I078GrD9PnH9qktXYe/SVY3dM+iPxb4K1+2Drk5/zIyVBdvvENemThl174KQwphajB62lWWATxy3/459vqIfitn0YwAdoZTgMUU/FocFu2eXlKSfL1G2fBtT773cuKWcVPApBcr//ZQTf+XIeX/XxHn1/SatmeU1ZCHmEH553RkrsVT1sd+qOeWHbHr3uae3j9nqYq09D3gxcHqL8hpya+ALge1fBG41sS1CCCGEEEKIPUQTZQvAGHN/4C+ttddqUmADjDE/AV5mrX22t/0FwLOttclplL6yleMx8b/maxYE9OvL8QDkHBvyIHVJxbPkfPnHFKZUdjc/K5OfBSqUQc15rWMZnFJZ0fx6Xfnduf7+Aq+pzHMhZax7bCgewS/XzcV3ik3IaxjLKuUIxR3EPGOhLJH+vdtazPrVO72e3Xp8b6yLy3r2aU/YZYvz2DtVwve0O0+yi0OA3d5r5/331byu18qPv0r1W/+a/diXrXi1jrLlqyIOPxYj5HH3FaCQguCUmZNO/hAAn7j0JjvOKVGRhBiTlMLo9213rFMrYftZz1HQunV1yVHsYkqdK9PF1MHu58pXid1zCbvtd8c6JdZ/X8HuuEo3Zjq1ylftoZMtc6nIpZ59Z4PLJuvGNKfEh8r33y05M0ZimW3dO8HFasL2GOaUxdC43S0Tdisnsdk5JfFRXXJ+A4XKCpXn29ptk9x3cOo6UnHyqXJyKfnt6MhRwUrq6csoWGJbKlNln8oaUsP6YrW6fSalbJXGbAUxxhwCPB04r0V5jQmNxuGRGzDGnAqcCrCPa45lkxBCCCGEEGLDKc1G+BV2f7xcHXDpfR5hrX1HI9sGY4y5GHibtfYJ3vZXAL9krb1+6vy+mK2S7DNDPBA5cV6OlNcgtt5WKM7E4RQbP+tXyOsSuzb/2rsxEb6nbZda1VEVDnj81XaU72edCcVS3eEt5wLbnr6Y166b/SnmeQuplLHMQb46FvJUxrK8+eslQX/cWGqdsJJMSDFcW3fXmYmt5eIrdV38tW+cJ9SPR4DtmAQXt+Rwild3nR8fp7b5ylZIlfTvhx+HEFonx/VhP2bLrfnWvXY/o5m/TlI3O6SL7XPxK35slYuB2aF6LmNCfI+3vwaU84yHjnUeceeF7mbX82Na/GN91a27zR3rYsD8/7ciVm439qbPhpRtfjxUSlGpUR9jiu9UdOP5+pSmkArt4qtchkQ3Y8CP78vBj9UM1e3Gcz87ZBc/xtM9Z365rr7QPj9DpnuGYfsZ9BV+9zx21znzcWOCb7dfX+zaIB1b58px70yfUJyo/67Kmc3ix1ynYuicau/X58epwvYYGcuU6Ai9i2NZFEvW7szJPOzwZ+d08ZUtX/0KxRKVqDouC6eLTY6VlaMMpWKeSuLHYse0yJidomZNrtS5fUppSEFLrbNVqmx9kN0fWz8Gvgr8xTJ2a058kUXcls+tgL+b2BYhhBBCCCHEHqJZzNYcMcY8FXgpcAtr7fnLbUezSP3+bGvtf0+df63rHWlv/aCnJVdQh7JMLyXzVkPl+/StQwFxhSO2Mn13m19P6qs+ts5WiWfJb+vUqu9+eaFsh7F2T8Wi+Vn2Uvh1+pkFXdt3PX4nHnEsEFfmQm3rHxtbR6p7Xq5npkvfWitd/HsUW2MH+jN3+cop7G5Lhx+H0F1nxl97LRZb0D3Wb0vXf52Kl1pTJhaX2PVy9sUWdvtgTH1299/FjDjPfhd/jSznPfezLnbL99vFXx8Odnu6Xfv42Tq71+y8+/66Uf5Y073fTi3si63prpPkX3MsgyjE18zq9h8ffy0lZ7evVnZxSqy7Zj/7XVeV9J8DZ5NPtx6nfrr4QafEuWyKOeqC60cuFimkvPuKg4vDcX0vNBvAxSQ51bObkRTC6xk6/DGimzHPz0jprqOrtMPOccqN4+4ct5aYH6frZj/A7vXgHP66dLB7bSwff02f0D7/mQqt9eZwdrt6Q1lG/dkE/kyCUAxzLFttyv6+LMihPujH3Mbee91tvt197zaIzzbp/r6JjRP+2JxafzNHffFnfzhScVkt4v1T5Ko7OepeTkxbTMXLyWRdosiFYqhSdvbRl1Hc0Z2pdf4vP3fcbIQz5o+AC4C3G2Meaox5CPB24J+A3XMUhBBCCCGEEKIRpTFb5wMPs9aeE9h3G+Ad1tpjGto3GGPMUcDvAb/AIjHG+4CnWmsv6Dv3oINubE+405N7V7fuUpK5xGfI17fv/UjNt62JH/PjulLzufvUtm4b+J6rEsXGry/HExNb46Nr85kXLrq3m8vv5kOnYrZitqUy4zh8L1pOFp2SudMxD1BqNfacdorFbIXw76vvOU55XB3OFqe6OO9kd00aX5H1CXlPfc+ur1am1iVzpObg56zhEqPrjYVtT35ITe+LYUyt5eII9Ve/f/oe71Aspq8W9Smb3XKdYukrXM6D72L2uvb6ZYTKd8TiP3zPNey+R1uq51KtcipbN7bUbXMqlYtncvGILg4Pdis//pps3Xg+R3ddpZAtIXXMKR6+Mu7HK3bL91WRWNa67jG+UuCva5jKDOeUG9dHQoqKu7ZjTj57R3khpdThZzP11cod75jlfXTxVlsZSb2xp7uv753Yfe5y4ox9YiqbU+Sc6geduLEe9SUn1sY/NrQOVt9vi+557hyXmdbdD1dGd/zw41trshGmfnP1ZUEeoiaFsiHnZGZ25GQh9Km5zyVrwfbVk0Pqt1Lfb7ia+nLImanliI1xY8VsHQ0cGNm3D7hJYXmjY639GvCIVdshhBBCCCGE2FuUKltXAXey1n4ysO+JwAuttYc2tG+lxLIRprLc5K4yniovh5pVzGvKz1mR3NE3NzelpIRiamCnp9ePBdo1r/f9i5iIbgammJpTMr+3ry26x/hluGxBqfivWLxaSaafFCXxdiXeTYffHiGFK6bI+etipebIl9jU53kN2RRT6lLertj96PbnmP2pa3LluuyGbr2znL6XQ6yPhWLO+jJ0+eob7FYEYipSVw3wYyVjcSahzJux5zp1rB+zGmqTPgU25aF2+876s9cA23FGXfysdC7bnluHyXncu2pWLA7HqUgulscpB9BZqymhMPWRMy7F3o2peFqHf6xTQLrX0qcUpGIlS8aEkL0hm0PXnOPB73sfpeKnY3Gd3bpiswJysq7FFI/UuJ5SFv0YT3df3cwRpyal4shi8UWha+xbpwp2xxf3rccJ8Vk3qd9GY/wezMkRENsfOjZ2HaGYqpLy+yj5fRP6bVmSD8HRN8MsFQOY89tx0DpbxpinAU9b/tcC71wuFtzlGsAhwBv7yhNCCCGEEEKIvUDONMLzWcQ5ATwa+BTwz94xl7NIpf7qdqYJIYQQQgghxPpSOo3wT4DfcWnUN53cRY27DFmcrUYCzpkiVmPLEPqmqIQC9R2p4MPcRQpzEgHkSM0hG2KU3A//WH86Uwh/WkjN4tg+oekzjpx+HJP4c6Z4pKZf+TYMWayyZkpByTROR83U4tBC0bFyS6YO5Uz5KEmQ4dfjkzOFKyetdN+0ztR99ssNLQsQS86S6rf+VNbYVNPQtftLCHSTKzj8lOlu+qBfRjfRR+7i8V38clIJa2LXlpo+lRtsH0rc5IhN98uhZHmRnP6Us1xG37XmTJ/PSebg7PangqaSUQyZ9lXyO8G3ITQW+Il1YskourbGxrucaamx+xoKVXBt2rfUTWhb6yl0jr7piTlJQXLCPmL2h9LVt0hcUZPKPidZhz/O+vtD5+c8F64dXnzaInH5s09bTHtNJZ9pliDDWvuYkuOFEEIIIYQQYq/Sq2wZY34K3MVa+4llgozUCdZaW5rhcLY4ZatF6uzYcV36PCe1Hvxcr0TKi5OTwryvnpAXOBZwHvJuxjxXqVTdQxSgmPcjJ+A5R+nwiS3GGDonx7vWFygcUtv6FKdUfbEFHEtt6SNVT6y+IeTcuxz8/h8KynZpqt0CtX59ITWsRvGNpenPSctbot76xBIP9NkL6UWgY/3HJaiB3QsT+wpBKL23nzzDkROE7wdUh9S8vkD9VHv1JTnpElMIQs97LJ2+b1vKw54TvN73nqsZx3PUkZzzY2071ayQnEQc/vbUvlSCiZyFgkvtzhk//Gemu0h2bMHj1Hsi9g7OmUURS+udSqCV0zdiyRVq2jznN2PMlpJ32NBkHjXPUK7qFXq+HX2LD/ft821zx7h38olHHLtje+g6hipbvw18vfPv/HmHQgghhBBCCLFHKYrZ2mvkpn5PURNfUlJ+jpc/17bQPt+22P5YOaFjQx7L2Hz6mjZOeZl9SuYRD5mDXxI7VOPVHOpxjfWjVPkpJS52Ts11xPpPTnv5tow1B781sTi+Ic9DarHsnPTF/r4cb3PsHEco7sePk0qNRX0eytD+2Fjjp57v4ux23vdY+vXuNbVYbiI1pvWNnaHyh4x3qXdXrL6S911qe9+itjXKWc34FHp3xa4rJw4r9t4OKU85964vbXuJMuRTcv+79KUWL4kPLZnJkaonR3nto2Z2yTpQ83zUKKU1MVUlY1urGSklv+1SytZ+oY1CCCGEEEIIIYZRmo3w/YndVwGXAZ8G/thaG1/FdU3wY7YcOV47/1hHSUaZGs99iNw5tCVzv4d4bFILy/r1pezsy1gE/ZnyUnEafr0l1Mz1r4nzCtXXN0c+J+6uL/Yix6YusaxVJYpKjO69i2Ukypkn3sILmXq+HS08oSVe/5x6Ssay2PPX9fr3KSchr3/J+JRrYwklGTJzFBuHOzZnYV+fUD3+s+P6v7+QdE55qWvuUwhKVOKYspmyN1R+64xmsWNrVNscW/ueu5LY6Jo4n5zfEv6xOUqH37dzFMYcu2uej7Fmf/T17Vj9JWXnnJ869vw3HAfALZ57WbSslvYO/V0wxKaad1hf+SGFzt+X+g3TLBshYIBbAIcDXwEuBm4I3BS4aPn/XwSeZoy5p7X27wrLF0IIIYQQQoiNoPRj62XAGcDx1trPuo3GmOOB/wM8n4Wy9R7ghcDDGtm5UnJVnpxjc4h5Vrtf3Qf4X9eJ8v26+zzI3brcV3vLOcepLHsue5iLm0h5yLauORET5itarm6XkeyoCxbndhWRKyN2l3hVazwzNR6aVP/aVV4Du0vmcYfY5QH19u+IiYjU6Xufu0pk7N6F2FJEI/trYhpT8SY5a5/k0rXpwMS+bj3dtv3p8u/9I3bn3Oetc90BAdViqz7v2XF/X95VuSP2++TEbKWOjSk2JYpWSr2IjRM5Xv/YmBCcDbD8f46aG1Nsdt1DOvGbXjv1PS+lfPVxi154iwt22huKV/THsByVwW9D/5jUmLm/dz9K1jtzhN4Tfe/iUBvXxK/EtoeuvS+mJvU8un2HRo+Iz9hJUfM7Kud3jc+QuKJY/X3lhI4d+l51HHPy2YtzEscMud+xMobaH1PYc57VlA1+OQ4/s607d8dvU+/cM//sNcB2VsLUb+4QpTFbLwBO635oAVhrP83iQ+sF1tqvA6cD9ygsWwghhBBCCCE2htKYrR8BD7fWviuw798Ab7XWXsMYcw/gPdbafe1MnZ4h62z1zXdPkZOVK+Tx7JKaJz6EIfOvHaFza7Iq+uX2ZXjqnjMkBiDHJrc2w8988NHAtscpZUts/n5L+/pwyuLNnvrx4nNbxjyFyumLLai1qSaDV1+fTmVnGmJ3zvXkxFf22ZLj0S2xJSdjnl9fjnrhKzX+OJLy/g+5H2deeA6wvR5a6jocQ8a41KyDGsVg7ExwJbFCNapCyzigknpqnovuWOC/0/tshXg8VM5z7dNyhkotObE2OUq7f+7UM5Bqzm1xbOo6+tb1Sh0bKt8nZ/ZN33s79FyM/fum5ndfTRxZy2yEFwCPi+w7dbkf4HrAtwvLFkIIIYQQQoiNoVTZOhl4PfBF4C3At4AbAI8AbgM8ylr7RmPMK4DDrLUPb2/ydMTW2XLUeClC8RstlSe/zC5+3WPUm2NTScxC6HxHrvcLtr3KDt+zmLNORI032Lc1hxKvkSNnXZYhbdtXf235NZS0T8l6ZzXe+Fgb56iSNSpViiHPUioLWol32THEc1yjusUUtBCxjKSh/Tn9xqevbWs87DnntF5/ru+cHHI8yC6DmlP/W6gMkG73WPk1z1ALda3knNR15M6sST1DfWMbtMnWGNsfOiaHPqWm1Tg75NyS/hU7t0uL96h/L2tmbuWUH5s5EqqrZsxJ2RaL0ap5lmL1Arz73NPbZCO01r7BGHMJi/is/wJcDbgC+BRwf2vtWctDn852DLYQQgghhBBC7DmKlK0dJxqzH4vpgpdYa69qatVMqFG2+jyHJcrWkPiMEC28IDUemRb15RwbItcrMeVc9rHbckwvXc4861QZLWLkxurjNTa0eGZr6suJHcmpz2V+cxntWqtUJbb4x8TqLzmnq2h317fqHuMrXam1v2Ie/aHK8pDnokbBHqoa+dQo4UNsKmmvEmVoiG2x9sn5feCTM9MiVGaL913J+6lk5kDI3tJzUscOUbb6bAydP4eZQTW25MaOhpRMR8lMkSFjQ8qmlr9nU+WXzOxopmx1WX5gfav2fCGEEEIIIYTYZIo/towxjwZOBo4C/GyD1lp7sxaGCSGEEEIIIcQ6U5og43ks4rW+sPxzuX+MtfYxzaxbMbHU70Pk3CHBlLXkBsN3j3Opy11q4yFBrUOCKocGadac23J6QE3CBEdJe7nEAJCfWrpm+knKvhaBtal6hty71OKzsSkRrQOFc6YhtOjTLaYb5ZTvKJkanTP9y5/aGLMnVH5f/SlbSqbJ5bRtLCV7aCyomZaamv4YO8en1RS+kj7eZ1/rd0ssIUrOWJNDi2lSJQtr+6Seh5owA/+YUNKZlmNaqj/FbJsq0ViXnGmi/nFThV3EGPJ8h6ZT10x3znlfDJmO2mIKc4iad1cq9XupsvUfgJdba59WeJ4QQgghhBBC7ClKla3vAQ+11r5/PJPmw0EH3diecKcn4762+1LJdqnx4PvnDgkCHnpsn3c5hyFB3zXqYYvEA91tOanyWyqZNQxRWadMvVpyTN85Od7alu2d8lTGaKVGtwrM7yMVNF2i7PvnxBYdvvY3fgKEPflfesoNgbyFtf0F4F0SDD+gu1tXX5ILV1bXzlhqf7/MVLk5KbT7FtgO7athiIe69fOd49VukeQg533RYrwekoQiZmtu+UOu2TH2EjFDlPecdsoZr4Yo7zUqzBBKfmeOPXOn5F3Td2zN7+ecPuIv8wO7U7775YVsqmmflosafxA4ttgCIYQQQgghhNhjlCpbNwfeCrwU+H/ALrfkJqWB70v97sjxOk49Z7dkPnqOFyTn2D6v0FQesiG25ZQ31PsRo2YB3qEq3hjntop1ip1T4t30jxmqgrVU/FJz40vIVXpbxTrlnBPzlrf2Bpd4WmNxMSUe9lj7nHfGnbf+ffQ7rwDKVCNfWazxnucsXBs7J9UHa5TNsRli0xBldmi5YxPrRzlxamPT+h51jys9dqzxKFbnWM/OkFkHLdTbEGP89m313h57BlLLmK0vL//+k8h+W1GmEEIIIYQQQmwcpcrWaSw+qKJYa58/0KbZ4GK2+hbJa+0tyvHKxzwZNZmihlDjBS49Zgg1c9jnQq2CNkbcwVjxDWMxVezFVLSIzcuJQ/Bjh2C39zenXXIXzqwh9Fz0LWjapW/s7Ko8oUxssWP7yg+RGzs19B1TsvCxX+dUY+aQOMWa+KWctq+5h67vh+JEXHmxftVytkaKmndLSR8cWv7UY26NkjmknpLfjq3jilpQMn44Ws+omeq3Y9943r3OZosaW2tPKzleCCGEEEIIIfYqRcrWXiMWsxWLR+hSE7vVYv7tnLzxDmeTW7vrfr/62K19QzxIU13rkLnmNRl3hhKL95jatrGyDY3lAR/yrLYit9ycTJJDFOUcFT3Hpr5zupmjcjNGhfa5MfjiEw4EyrKoxhSu7nX0KWdzGG9TtgxRQUrqaVFeq/LHULdLskLmzDJpHcPqlxNTmEsU39A5MWWuxP6YzTnnjq381VBiW4tsz7W0fH5zZrzkxMqmfit0j615L+VQ8jz4/+++wz752mc0y0aIMeZwY8xLjTGfNMacZ4z5hDHmvxljDus/WwghhBBCCCH2BqUxW7cAPgRcF/gI8E3gMOCuwHeAu1tr/3EEO30b/hNwb+AY4HvAJ4HnWWvPCRz/eOAZwE2BC4Dfs9b+r5y6nLKVO/cb+r1drVSSlllVStbYKbE3h6nOKWFI+Tlewb5yh16fO79vvaKuTX1rDw2NC3DUxE2U1DMHhaGPqZS5nP0l8aGxc/zMXt1zclXW7rGt+kLXti45GUkduc9qap2t1D3M9fKH1LY+1WIVcQ4l5HrAu8f0lZE6pq+sHFrHIrWYvdLqd0FNO9W871rej1pK2tTRMvasJnNoTcxcyXheoxamxvHc38JdWii8NbFtLWaDQNt1tl4C/AtwC2vtva21J1tr7w3cArhsuX9s7s/iQ+u1wIOB/whcH/hbY8zx3QOXH1qvBN4CPAD4C+AVxpgnTWCnEEIIIYQQYg9Tqmx9F3iitfaNgX0nA6+w1l63oX0hG64HfNt2DDfGHMxCtXqntfaU5bYDgAuBd1lrH9059jXAQ4DDrbVXpOrqW2erZp2cOSo5c7JhDl6vVVPjwco5toVNqbnZY927FrGMNfWk6PPetfIk5tByfbua/tQqq1vMhhJvfE42sVh8QKgfx8orUY9WFd9XG98whrJfc06N/aF7F1NV/bKhLjZvyPqILfpGKn6zz44uLe/ZkL7TutxU+SX3I9eGVVxH37FDxx4/w2Zu/a0JXfuYmW9LaalsXZ3FtL0Q31vuHxVr7SXW+0K01l7GYg2wG3U234WF4vV6r4g/Aw4F7jamnUIIIYQQQoi9TekCxGcDTzHGvMtae5XbaIwxLKbznd3SuFyMMYcAt2HnYsu3Xv79Be/wLy7/vhXwgWS5B16dA4686S7v5tbX9fL/VybK8D0Ltd7l3GNSnreSWIiY3WPHA/TN704dm7rOdVIUa7zOKe/sVJ72Ft66nDJyrqPGlpL22T/j2c+tJ7SvhC1lfUAcQklsmH9MTdxB6jpT9saOdezvxf90xwLnAb32N34StDtk24GuHq+MEi/qVua2xDXXqCN91KpXQ2JTtt6VBWNzrL8Ofa772jK0PZZ98kqvX3X3xcaCkr6eQ9/zl6qz5neIo+T3Q6rcvnJy2qnFb5kU/jHnnXHnrX/HYp9LKFGhW9BKoXNjZsnv1yHvAJ+t5y7Q151tsXpK6xqL0o+t3wb+L/D3xpg3ARexSJDxS8DPAA9sa142vw8Y4IzONhex/B3v2Eu9/TswxpwKnAqw74CDGpoohBBCCCGE2EsUTSO01r4beBCLKYO/CfxP4LnA94EHWWvfU2qAMeZ+xhib8eevI+c/B3gU8GRr7bndXc7sEnusta+y1p5grT3h6vtfo/RyhBBCCCGEEAIYsKixMeaaLFLAf8da+8NqAxblHJVx6A+ttV/zzn0i8IfAc621L/T2PQl4BXCEtfaizvYbABez+Dj7n6kKD953mL3rkac0CS4eklo0Z+peXxk5lKSwXzeGBDHn0Hff3VQiiC/eWlKPIzS9pebeDelHLYN+xwrcL6l7SArZsWzyt3cpmU47JjXTIkMJDb5/o0Xob2yqXmoaXSpBRiyZhqMkwUer/uofk0oJ7agJ7o/V65dZyxySOPnktk+r6a9DzmmVWMLR10dCSyKM9c5v0XZT37saG2sSprW2oe+4oTat8j3dLTtEqq/ExrsW19Mt+93nnh5NkFE6jXCL5QdW9UeWV84/lJ5njPlVFh9T/93/0FriYrNuzWK6o+NWy7//rrROIYQQQgghhMilWtlaJcaYh7FYM+s11tpTI8dcjUXq9/9rrX1MZ/urgYexSP0ejqxb4qd+nzolbqvyc72YQ70fY3tZ+tTBVm3re/+mVggcJSl9x1Jfatq4dV8fkmZ41Z64UuaYXjhWT875LcaRkqQ8qf7qpy+uGQdjtriyc8oPlRdaGDpWRp+6nZMgyAX+u6D/mvThofJL1MKYbf727r6p3qN+valjp1ZUVqVct7ahJCnPqsbmVfbBOf7OHMKc262mvND+VOr3XmXLGHMV+XFP1lpbrZblYIy5B/AG4HPAnxpj7tzZfbm19rNLQ64wxjyPxSLG3wDOAu4DPBZ4St+HlhBCCCGEEEIMoVfZMsacRkGSCWvt8wfalGRpz29Fdn/VWnu0d/wTgGcANwG+BvyetfYVOXX5MVs+NWpSq6/vMeJMhnqz+2ysPcbn/DccB8AxJ/evNNBSbRt73n6NTUMYWx1pbVPNOWN5T8fw7Jba4JczpN+29pb3qcMlY6cjNbb1KVyhcmJp3Ev6U+y4FDkqUsn9yLm/Yy/8OUR9HlJ+7LgWdfbV3ze2tL52v55WMWFjxCOO/e4am6lVmDm21xx/H+RQ815t9Y4cpGxZa0/rO2ZKlvacVnD8K4FXjmWPEEIIIYQQQoRYy5itqfBjtnxqPEuhufm53viQpyEny96qvRBjxVaNxRCv3RhewiG21trUqu650Fq19cudaj5/DWP1wRJveUu693JIXGVJrNBU5MSLrlrlSZU/9VgzVXxLzbE5+O3U4r7n2NYtY4wYsFZ9b4zYz6E2zSFmbgh+u8QyoIZU3BbPyhxjt4aoYV1S2QiL1tkSQgghhBBCCJGHlK0ELZWtknPm6B2fik1p05qYC0fOdeQcW0NNvExNuVPP9V+3Z6rFc1ATU1pDjlcwx7YcD2suqeevGzMF6VkBsXJKvOZDPOFzVobGKKe0rKme69Yq7pzHo1bKgKOvnKmUopJ3y1hxPznkKO65NoXOnTr+aoyxocsYvxNS5aauJxWzJWVLCCGEEEIIIUZAylaCmnW2po7ZmRNziSWosSl1zlT3cJVzwceKBxjMDI1BAAAgAElEQVSzvjnMnW8R9xPy1n3tkYcDcMRLPjrUxObjS4kqNvY9SilMvm19sbAp1aIvNmwV8VH+viFxLV3Gtn9IPX1jcyv1YlPeyT5zVOFKYoSmij8dcs0XPuuuW/9uMX6PRc3YvE7PRc2YPDQeWMqWEEIIIYQQQkxMkbJljHkWcGNr7VMC+/4H8E/W2tMb2rdS+mK2Smg9D9pRoyLUrLky9bzqOXtOWnuxp8qotQ6sMhNcC4bejzEzg9WOPatWO3PaxJ37pafcEICbPfXjveWPFUtQUl5fNtnWcSCOuSsbY1ISf1Vj2xwU93VgjrN+5vA+7bMhZ52+3LJyzi09v8Vv3dbj4FQZElsqW48BPhfZd/ZyvxBCCCGEEELseUqVrR8Cv2it/evAvnsBf2WtvVYz61aMU7bGiGcJeQ1aZK0KsWr1JYch1x7LYgZ1noxcL0iXqdp4rGyENZ6fqTKQtbjmse0vqbt1bFtuxrxapo6fGJvcOKwUviIVouT+1pzTwks7pH/lUPJcO1pl/BuD0Lul5PkYO655SDv1jSNdxla7p44h7itjaDlzZqzfFLF6auIrS8qbWoXr0nKdrR8CN4rsuzFweWF5QgghhBBCCLGRlH5sfQh4pjHmwO7G5f+fsdwvhBBCCCGEEHue0mmExwIfBS4BXg98g4XS9e+AQ4Gft9aeM4KdK2HVCTJalbeqBAxTnTtkSmBrm/ZiQG8Lxprm13rqwhgMTTfbV+5Y0whbMNb9GDImhALQ1+E560u20Yqpgu5bTzMqOaem3iHvrJJ6WyfgWNV06rHLn/q3S03isRpW+Q6b6ndnTlmrStyTqjeVIOOAksqttecYY+4NvBR4Fgtl7Crgw8AjNulDSwghhBBCCCGGUL2osTHmGsB1ge9Ya3/U1KqZcPC+w+xd///27jxekqq++/jn64CMJE8QBmVRFAU1ghERVJaogCBRE4WIUaMZNAkIURIVfQiIZnhQUYkLJmJANMgSiRuKgiD7jrJJAgjILpvsmDgwDPB7/qhqpqanu29Vd63d3/fr1a++t+pU1ak6VXVvnd85p9Zb2OrazKp0oSZ3mDYNzT5ImQOhDMpT3UMPVzWAxThpu3zeZg0b9KWszt792hgFmWud2fVWNaz7XANhVD18+KgoVVX7XKauD4Ne5r1t1PKTROCb+HvXpRYcdQ2cUPZ66xq4ogxFhqNvQpWDSGXLadQAGYUiW1npA9ZUPmSZmZmZmZlNas7IlqTHgS0j4ueSngBGLRARMfYDXNsM67NVdc36ODWWZUQGigxrm0cba1rLME5NYhvbWY/KU55hgJselncalXntj6oZzTt8+FzbzpvHMkw6BHWRocb7l5mr3X6efneDomR5y6zItTpO35FJ7/3jnIPj6K2vd4z7I39N3mfrumcWWVfV0ZGxrv0znp18b3dbqXkZpqyh6+seanycv9tt1saIadkm7bP1/4DbMj+P1+7QzMzMzMxshozdZ2sW9Ee26hr9ZFitXVXa0D65qjx1vaakX57atWnf5ya3M9fxb/K8LWM0tDacM0Vq48cZWbBN/Yh6UaiLPvNvAOy47iYTr3PUPaH/b0tPG/qzFNGmfjhdcModydhlb3r1TivMG/b/xrB+o9k0eeQ97m1s/VH2MlWsoypdi6D17qULzr8LyPc/URnRtrx9tgq9Z0vSjenw74PmvUTSjUXWZ2ZmZmZmNq2KvmfrCWCLiPj5gHmbAz+LiHkl5q9RdUe2qpI33028X6bMviNzLT/uOibRRI3ZJPs8Tg1Q1X0VhuWh6ghFE9f7JOdLT5l9bsbN0yTbq1qRqF4VEcBB10X/uTxJVKGqPiplqCo6UnV/kLKOaRV5GrXssNE079t67Sd/HhQJyC476f21ivt0G/4Xq6tFUJsi8WXuc52Rszv22QqAdT97QaXbGdVnq1BkKzXs6Wxz4MEx1mdmZmZmZjZ18oxG+CHgQ+mvzwLuAR7tS/Y0YA3guIh4V9mZbEovspV3pKVBacpUVpRn0ppuyDeylnVX0xHNIsuUPYrmJHmcdP11R5GqqkUdp59GTxl9O9pQ893Thr6r/UZFzoaNpjgsb3ny10ubjagUGS2xTFVHMus+9waVw6j+K3nX14ZrZxJ1R+vrUldEqKp3Z7Xxf4uy1jvpaIQ3AqenP+8KXELywJW1BLgaOCLH+szMzMzMzKZe0T5b/w4cGBEzMRDGsPdsjVJlLeaoGo26nubH2c4kteWjohb900ct33T7+jZuZ5QyIhxdGc2orn6Ddfd9acN51FNmVL3q86pIjW6bjnG/qvtMzrXNQar6+1RkW0XKrM39QkctO0lka9h72uq8n7e5ZUKZ97KyZO9ZUH3/aRt8nUwa2XpSRLx3ksyZmZmZmZnNiqJDv+8j6V+GzPuypI+Wky0zMzMzM7NuK9qM8Brg8xHxtQHz/hr4SERsVGL+GjVJM8J+RZrh1R16LtJEIk8TvkleIDdqe9Ma/q6rWULV2+lCk6o2DRXcpuHP61q2Tea6V49apupBWvoHIBo0iEqe/PYv04Zh3Ietp8h+lLHdslR9LVV932j6em56+02qa5CkNiqrKXxV18c46ylz6PfnAL8aMu9G4LkF12dmZmZmZjaVika27gX+ISKOHTDv3cCXI2KNEvPXqLleajxJlGqcYXOrfmFjkx0+J4mk3LrLOkD5L6xrU+1/07VdbYhKVvUKhqaPbdc1GQUtI9IxbECAcbdXpMXAOEPk5x0oaNxo2zhlU+QYlrHdMiOMRYa9r6ulRRtadHRhiO6qoiN1a0N5d9U4/5eXMWjboPWXGdk6F/iopFWyE9Pf907n10rSOyWFpNuGzN9N0jWSlki6VtIedefRzMzMzMxmT9HI1ibABcC9wDHA7SQvOn43sADYOiKuqCCfw/LzdOAaIIDHI+LZffN3Aw4DDgJOA14H7Ae8PyK+Otf6V5u/dmy13sLaX05YRg1xkWHiJ40M5M3TpMbpC9Z0jdW0aFNEswvKeiFkG183UEa0viqjaiyrPJZF7rdzLT9X2rx5KdIPZJyWFll1l3eZZTnueVtG/+wipvVeOY42nIPjaNNrB+pS5N5WRn//uV4/kXc94xz/Mod+v0LStsA/A/uQRMaeAM4D3lrng1bqc8AVwJ3A9tkZklYCPgUcHREfSyefKWld4EBJR0TE0lpza2ZmZmZmM6NQZGu5BaWnAasDD0TEw6XmKt/2twZOBV4K7A9sn41sSXo1cA7w+og4NTN9W+AMYLuIOHPUNiYZjfC+rdcGymnD3jNuzU2VI1E1MeJLFdpQK1+3skdQK7LNaTnGdddeT8vxa8N+VBVxKmO0wDzR+knWX1bauZZpQ3/gprZjo41qYVP36JbTci7MUguUItH6njLSjlr25OsPLq3P1pMi4uGIuKOhB62VgcOBgyPi+iHJNk6/r+ybflX6PTVD1JuZmZmZWfuMFdlK+269CJjfPy8ijiohX3Ntf3/gPcBLIuIRSUeyYmRrP5JmhE+LiEcy01cClgKfiIgDB6x7d2B3gPmsutkf64211fbPNRJS2durK5pRJE8949Sq1nX8imhTBHBUW+amR1xsUxv8smtGq4wsT7reJkcSzK5r0Pp6I9wtOP+u0raXx6j+dm2sNS+jBcQ451PVrQHKGmlu1HqHGacGvA1/R8vo59OG/RhHkUjHief+AIAd192k+oxR3z2zDX0Z8663zaNDZhXJS2l9ttIBKU4EtuhNSr+zT2yFHrYkbU/SHHAuZ0fENpI2BD4G7Jx9iBq06gF5m1NEHE4SNeMPtMZ4bSzNzMzMzGzmFR2N8FBgO+BvSIZ53xl4CPhrYEvgHRFRqEpF0qokL0uey+KIuFXSSSQPUO/KzDsUeC1J08ElEfGwpD3T6etGxJ2Z7T0T+A3wgYj4yqgN5u2zlWfksUme1LvWDrfuWomyj23VowC1qdZmElXUopUdAZzFPmhN5anJfo9lXrOT7kee0Rrzrn+SkVbrLI+5omttigi2YZS3qo5HF/u9tfEeOkrT99e6/6+qapvT2E++tMgWsCNwAHBR+vtt6cPVWZK+CvwDsLDICiNiMcnw7XltBDwXeGDAvAeAQ4APsqxv1sYkoxVmlwe4ukg+zczMzMzMiiga2VpMMrrfeenPb4yIs9J5OwDHRcSCSnK6LA9bsGJfsX8ENgPeRvIAeH06iMYdwI8j4r2Z5Y8gicitExGPjtrWJKMRTsuT+iTa9J6RNkXxulQz5vN5POP2QauyxngaaxKbily2oS9mkb6qbbmO6xwBtYx+wE0fr1Gqup9P675Xvb4yruVJ89BmXf//Y1heen2LAS7+5t6lRbbuAp6e/nwLSdPBs9LfNyy4rrFExEX90yS9h6T54FmZdEslfRw4VNLtJC813o6kyeNecz1omZmZmZmZTaJoZOso4NcR8bF0tL9PAN8EHgN2BU6IiL+sJKej83UkfaMRZua9D9ibpOnhrcAXI+LQPOvtRbba9HQ9TBfyOMos9q2py7DRt6qq3WyyT0SZ/WUm2d4sacMx6NUuTvJew0FGRQez0wcpch204Rj2tCkvk2hz1HOcqNu0lEu/NvdLzDNa6jj3nGkty1lXZp+tA4B1058PBhYAbwdWBU4A9ho3k5OIiPeMmHcYcFh9uTEzMzMzMyv4sBURNwA3pD8vJYkY7V1BvszMzMzMzDptrJcaz4pxBsgYR1s6YY/TebnsMH7Xwutdy2+V2nws2py3qnXhJc2T5qVnWHO/udZTpmH5HtWssEhzxTI78w9a/1zLtLHJ27QMTlD1edum5vptuo/YbKmqq8NEzQglPUH+FwNHRBRtmmhmZmZmZjZ15oxsSVpE/octIuKACfPUGnVFturWxtqnubbfZB7G0fRxG6WNeWtDRKWMPDR5bOvadpteIDzXepu8x7VpSPM2XDttXv9c2+0pa9jqIi1Gxtl23u1Muuywe8E41/WoZdpwPU/CUbwVDSrvNrxwvF+RPE0U2YqIRQXzZmZmZmZmNvPcZ2uEYZGtstvM11Vj1XtCX+Xm+8febhnaUJttk6u670iby7ONNZXjDPE/iarug2W+/LLIvab/JcGw7F6Zdx3j5DE7r984UZFJ8pZdb5nbafO1PK68Ua9xzsG6+mln0/a06Z5W5vrb1EqmzddDm/M2St3n8SBlDv0OgKRNgBcB8/vnRcRR46zTzMzMzMxsmhR9qfHTgROBLXqT0u8nVxIR80rLXcMm6bPVpVrBqtvJdrWmZJg21ZB1zbSdC+Ooq49QG/rB1b2dItH7In1S+tdXdQSiqr5CZSrSr6gr131X8gmD81pmtLnuljt1/l3tUjlPYpJWAlm9aH+b+lK10ajI1lMKruvTJC8yfg3Jg9bOwHbAscCNwCsnyKeZmZmZmdnUKBrZugE4gOThainwioi4NJ33VeD3ImJhFRltwjiRrXFql4ct269XywrFahhmpRanbG08brOUpzbuaxVmIVJaVT+KqlsQzBX1v2OfrQBY65IlT04bNjJbnv5XdfVn6fE1W7029FEuu1xczvl06Tg19W68Qdsss99uVQZtr8zI1jrAjRHxOPAI8H8y874PvKng+szMzMzMzKZS0cjWjcDfR8SPJf0SOCoiDkrn7Ql8MiIWVJPV+lXxnq021HRU/e4Ns3G1KdIxjjbWIDcxwlkZ6yo76tO//DjvCBqWh0GtDh5auCUAC86/a2ie2tifpX+b03qvn/b9K6pN95ph6+3pWpn5XFvRNLboKHM0wvOALYEfA0cD/yRpfeAxYFfghPGzaWZmZmZmNj2KRrY2ANaNiHMlrQx8Bng7sCpwMrBXRNxXSU4bsNr8tWOr9RZOzVP3tHKtkc2l6+dIk/0z2mac9xblSZtnmSIRuTzbLaOVQZ7tVKGN0bY21Ja3aRTQqrShf0zTmr7+mlTknYFNRUqLjJZaltIiWxFxA3BD+vNSYO/0Y2ZmZmZmZhmFIlvLLSj1D64RMe7KWqqKPltVGdYfAfK9c8ZW1HRtWhtqaSdRdv6bLo+spkdg7Pq5UUSRY132OwPz9hsrMlLsqHtznlrhvNGdrp4XXc9/m40TIe+peoTMrvUpLXP9Xe+r3ISm97n00QglrS3pREkLM9PmkQz9nv08KGmtiXJvZmZmZmY2JfI0I/w74OXALn3TBXwNuCP9+e3AHiTv4TIzMzMzM5tpczYjlPRz4OSI+ERm2jzgUZKXGl+WTns/sDAiXlVhfmvV5QEy8oSlmw7Ddk2bj1ee8u4NRb3aURfmXibvNsddxzjbK6PD7aTrGZanIi+wrUvZTet6yh60o5fmxHN/AMCG334fABt88KKJ8jmJ/mM3STPCQeVQxTnRRMfwYZpublunNg+Y0Ibj1YY8NK2Me3HXm5C36Two+4Xzk77U+EXABQOmq+/369K0ZmZmZmZmMy9PZGsJ8LqIOK9v+lrAvRHxePr7q4FTI2J+VZmtW/8AGeM8kU+yTE8XIgZNqKJDbFnrMxumq9db08oeUrtN5dDLy5L11wDGq/luajjuPNtsw7FuQx7qMmxf+1s3WLWqPueqarXQtDZHiUeZNLJ1N/D8/okR8Zveg1bqecA942XRzMzMzMxsuuSJbH0LWDMidpgj3Wkkka53lJi/RjU99HtdNXFtrhEta1uTbGcWao/GeXmrraiJyIrLqpguHK829HWa1vteHmVEGifZ7iz0ccujq/kuU5uPQZvzVpYi+zhpZOvLwHaS/lnSCqMXSlpJ0heAbYBDcqzPzMzMzMxs6uV6qbGkvYHPkTQTPBW4NZ31HGAHYE1g34g4uKJ8NqLpyJYtb1gNwyzUrnRRV/rBVT0iXJdGAy07T01Fzdt0TPu14WXf09pqomxtHF3UbFa1/X+KUZGtPO/ZIiI+L+kyYB/grUBvEIxHgHOAz0XEGWVk1szMzMzMbBrkimwtt0Dyjq0FJEO/39s3SMZUqSKy1ebasLbXGtSlzWU0rWbhmOcZYWlaIw5tir5My7nWhZFu22Rayr0M/ltvs6TK1ivZdU4c2cpKH67uHjN/ZmZmZmZmM6FwZGuWzEqfra7V+E1r7f8kXFM5vnHeFD9s2TbrjS4Hy0ZX87XUPbN4LKt4p+I46xp0ny0jb2XnaZbN4vXRlKqPdZE+k02Ve/b6O/n6gycajbCVJD1L0jck3SVpiaSbJB00IN1ukq5J01wraY8m8mtmZmZmZrOlk5EtSesD5wM3kQxN/xtgfWDDiPh4Jt1uwGHAQcBpwOuA/YD3R8RX59pO05GtLtTQVFXTN2r9c6Vt8/Gy2dDfNwbaeV62+Zppc97qNuvRi7r6XORdZtByedY3V5qy89RGvq7r42Ndr1L7bLXEvwG3A9tGxNJ02tnZBOk7wT4FHB0RH0snnylpXeBASUdkljUzMzMzMytV5yJbkjYArgcWRsTRI9K9mmRY+tdHxKmZ6dsCZwDbRcSZo7aVN7JVV81SV2qwXJtSvl5/m15fm7pNc5lO877VoYwafWufpt4/10ZtGk2zzWZxn6dVV6/VJo2KbHWxz9bW6ffDkk5N+2I9IOkoSQsy6TZOv6/sW/6q9HujSnNpZmZmZmYzrYsPW+um398ArgPeQPKy5TcBp0jq7dMa6fcDfcvf3zd/OZJ2l3SJpEuWsqS8XJuZmZmZ2UxpvM+WpO2BU+dMCGdHxDYse0A8KyLen/58hqSHgOOAHYGfkLx0GaBQO8mIOBw4HJJmhHmWqSvE2uRLUItoenjersmzr0WaD1Zx7KalHAY1jWjzi33bZJKhdru6z21RZ5OeKs/Trp4H4+S7jKHgx13PsDJ8aOGWACw4/6451z/OedCm8u3S/XbcgcaGpSlj31/+veuf/Pmyt2643LymXnNQtjrz0vjDFnAB8OIc6Ran3/el3/0PaD9NvzcledjKRrDuzKTrRbTux8zMzMzMrCKNP2xFxGLgmgKL9PpcDYs6PdGXbmOWf9jq9dW6usA2W2taahj6lZ2nvEPuVpmHufTnYdztt7E826LJYzNJ5/o2XKs+r6o3rJzrPPbTWs5NDzJURLYMxrn279t6bQBW61tmtaMuBGBJ9gXnFURFBr1AvQzTEk3v349BeZ1kH8vY95+/bF7mt+lsJVPn/89d7LN1EXAX8Cd903u/X5x+XwjcC7yrL927SaJa51eVQTMzMzMzs84N/Q4gaVfgSJIXFn8f2JDknVq/IBnSPdJ0ewCHAp8meanxdsD+wF4R8ZW5tjPX0O9l1TZ3qcbNytc7j56sjUxrH61aHs65mFne9yKqPk51DcleZX8QszzKaO3h87V803ZMJ9mf7Dl68vUHT9dLjSPim5KeIBmF8L0kkapjgH0j8/QYEf8mKYC9gY8CtwIfiIhDG8i2mZmZmZnNkE5GtuqS96XGeYx6ci4zsjVoO1XWQvjFd9ZVs3ZdjMpTVceizPVOW23qKF3d167m26wJvl4S03Icpu2lxmZmZmZmZq3nyNYIq81fO7Zab2Hnn7aH6WLfgkk5T2ZmVhXfz4vp2vHqUn6byGubj0/VeXNky8zMzMzMrGaObI0wSZ+tNj/d2zJ1lVMb+/CMY1r2wwbzfWu2znGXt42jjNHbqhoFtq4WOz1NXju+ftvFkS0zMzMzM7OaObI1Ql3v2eqySWuBu3gMq6757uIxKWoW9rEKPm7Va0Nka9rKedr2pwllHMM2lUPZeZm2kU9v+NIWAGzwwYsay4MV48iWmZmZmZlZzRzZGqGM92y1oYZkEnPlvw21wJZP18/FadFUOVR1rfq8mtxDC7cEYLWjLmw4J1ZUm/oTlanJvA7rF9Wl41eHYe9o9XFqhiNbZmZmZmZmNfPDlpmZmZmZWQXcjHCEMpoR1qWrYeOu5rtMw45BkWZfbs5pVfE1alYtX2Nm3edmhGZmZmZmZjVzZGuELkW2mlLXMOhVrX8a9deSutZ0trTppZtl8nlsReQd3Mnn02iDjtO0HLsy/laOWqYtx8n/R9XDkS0zMzMzM7OaObI1QhcjW22pSTFru2l7CeYw01wzbdXzuVKeKo5l16MWbT6/un5srTx5zlNHtszMzMzMzGrmyNYIq81fO7Zab6FrNGw5g2q72lw7ZzZNXNvcXr4PmtmscmTLzMzMzMysZis1nYE2iyWPdqaGzjWKg/WOy5L11wBg3umXjkwHcx/DQfObPu6zWNvflXPeo0OWaxqPW9PnRFn3j2ksG7Np0pZ7TdvvFWXn05EtMzMzMzOzCrjP1ghdHI3QytNUDUxXan7MrBqzcA+YhX3sEpeHdVkbzl/32TIzMzMzM6uZ+2yNoQ1P0F3UtePWVD77+/Y0mZeq1HUutOmca1NeyjBt+9Mms3BMZ2Efq+Trr11mqTzauK9tyssgjmyZmZmZmZlVwH22Rhinz9Y4T/xtrCUo07TtX9kRpzzHZ9qO4SxyGZrVx6OAzjaX92xqstWM+2yZmZmZmZnVzA9bZmZmZmZmFXAzwhGmdeh3h9fH4+NWvTYNCuLyblZVx9/l2g4uh+HadB/sKp9fy7ThfGpDeVSdBzcjNDMzMzMzq5kjWyN0MbI16Mm9DTUKZmZWHd/n21GD33VdPY+qzvew9Xf1eI1jlvZ1HFMX2ZK0QNIhkm6U9LCkmyT9q6RnDEi7k6TLJT0i6RZJ+0ua10S+zczMzMxsdnQusiVJwHnAC4FPAL8ENgIOBK4Dtop0pyTtCJwEfB34FrAp8GngkIjYZ65tdSGy5ZoGMxtkWu8N2ejFrbusA8C6n72gqezUwhGb9pvW662NfKytjUZFtlaqOzMleAGwFfC+iDg8nXaWpCeAr5I8hF2bTv8McF5E7J7+fqak3wf2l/TFiLirzoybmZmZmdns6OLD1lPT79/2TX8w/X4KgKT1gJcBu/elOxo4AHgD8O8V5bFyrtmZfq7NtklM6zmT3a91Pzud+9hvWssyq+t/07qa7y7ysba6Zf8fg+LnYBcftq4CzgE+Lul64BqSZoSfAH4SEb9M022cfl+ZXTgibpK0OF3GzMzMzMysEp172IqIkPRGkgjVxZlZJwJvy/y+Rvr9wIDVPJCZvxxJu5NGw+az6sT5rYprdqafy9isfmVFlMeJ1HQ9utM0twYw64auXauT5rHx0QglbS8pcnzOyiz2NWALYA/gten35sB3JfX2Sen3oBFANGBakjji8IjYPCI2X5lVJt4/MzMzMzObTW2IbF0AvDhHusUAkt4EvBPYPiJOT+edI+lG4KfAnwE/BO5P5w2KYD09M99q1rUajWn10MItn/x5wfnJWDG98nANu82qss75cdYzy9dbGfs+y8evSZP2Z+lfh8tx+rWxjLP/E6121IWlrrvxh62IWEzS7yqvP0q/L+6b/vP0+8UkD1tXpb9vDDx51CStD6wKXF0wq2ZmZmZmZrk1/rA1ht5w7a8ETstMf1X6fTtARNwq6QrgXcARmXTvBpYCPykrQ66RKcbHqB2yNTeP9c3rLyNHuqrTXysMs32cfa7Vr6t/w7p+rnQh/4Py2D/NUUmbBmVHs7Ia77M1hu8DdwBHSdpT0raS9gSOAn4NHJ9Jux/wWkmHSdpG0oeA/Uleaux3bJmZmZmZWWUUMWj8iHZL36G1CNgOWAe4kyTKtSgibu9L++fAPwF/CPyGJMr1qYh4fK7t/IHWiFfpdeVm3syspbpQ025mZpMpcq9//HWbAbDKzcsPdZBn2Vn6m3JafPfSiNh80LwuNiMkIn4N/E3OtN8niYaZmZmZmZnVpovNCM3MzMzMzFqvk5EtMytuVDh/0CANw9La9HJ5V6+rg1FYO8xSsyyrTpHzZ97plybLVLydaebIlpmZmZmZWQUc2TKbEaNqmFz7tCLXIFsVmjyffE53n8vOrHsc2TIzMzMzM6uAI1tmZgO4Btmmjc9pM5tmbY3eO7JlZmZmZmZWAUe2rHRtrVmYVS4PMzMzm3Zt/T/HkS0zMzMzM7MKOLJlpWtrzcKscnmYLR7+dvQAABBxSURBVONIr5mZ1cmRLTMzMzMzswo4smVmZjPDES0zM6uTI1tmZmZmZmYVUEQ0nYfWknQPcEvT+WjYmsC9TWfCCnGZdYvLq3tcZt3i8uoel1m3uLzguRHxjEEz/LBlI0m6JCI2bzoflp/LrFtcXt3jMusWl1f3uMy6xeU1mpsRmpmZmZmZVcAPW2ZmZmZmZhXww5bN5fCmM2CFucy6xeXVPS6zbnF5dY/LrFtcXiO4z5aZmZmZmVkFHNkyMzMzMzOrgB+2zMzMzMzMKuCHLUPShyX9SNKdkkLSoiHpjkzn93++NCDtH0u6QNLDku6S9AVJT6t8Z2ZA3vJK0+4k6XJJj0i6RdL+kuYNSOfyqpmkm4dcTzsNSLubpGskLZF0raQ9msjzrJK0nqTvSnpI0m8lfV/Sc5rOl4GkbYZcRw/2pVtd0hGS7pX0O0mnSfqjpvI9KyQ9W9K/SLpQ0uK0bNYfkG6+pIPTv2sPp+lfMyDdUyTtm94/H5F0haS31rEvs6JAmQ267kLSy/rSzXyZ+WHLAHYDngn8IEfae4At+z5fzCaQ9FLgVOBu4E+B/YH3AkeWluPZlqu8JO0IfA+4GHgDcAhJWXy6L53LqzmnsOL1dHY2gaTdgMNIyvJPgO8Ah0ras96sziZJqwJnAH8I7Ar8FfAC4ExJv9dk3mw5f8/y19H2vRmSBJxAcv3sBbwVWJmkDJ9df1ZnyobAXwAPAOeOSPd1kr9tnyD5O3QncEr/P+7AgcAi4F9J/q5dBHxH0hvLzfZMy1tmkPyf0P837Lq+NC6ziPBnxj/AU9LvlYAAFg1JdyRwW471HQ/8Clg5M21huu6XN72/Xf8UKK/LgbP7pn0CeBRY2+XVeDneDBwzR5qVSB6Cv9k3/RvAvdky86eycvoH4HFgw8y05wGPAR9uOn+z/gG2Se9V249I85Y0zbaZaasB9wNfbnofpvnT+3uV/vy3aTms35dmk3T6ezPTVgKuBU7ITHsmsAQ4oG/504H/anpfp+WTp8zSeQF8co51ucwiHNkyiIgnylqXpJVJag+/HRFLM7O+TfJP/lvK2tasylNektYDXgYc0zfraJIa3Tek6Vxe7bYl8AwGl+MC4I9rz9HseTNwUURc35sQETcB5+ProyveDNwREWf2JkTEQ8CPcBlWKuf/F28GlgL/mVnuMeA4YEdJq6STdwSeyor3w2OAP5L0vMlzbGX+T4jLDHAzQivumWmb98ckXSdpn74+QBsA84ErswtFxCPADcBGNeZ1lm2cfveXw03AYpaVg8urWX+WtolfIumiAf21BpYjcFX67fKp3sasePwhKQMf//Y4VtLjku6T9B99fepGleFzJP1+PVm0ITYGboqIxX3TryL5R33DTLolwPUD0oGvxybsmf79WizpDEmv7pvvMiMJ05rl9QvgUpKLZD6wM3AQSf+Fv03TrJF+PzBg+fsz861ao8rhgcx8l1dzfkTSn+4mYC3gA8Dxkv4qInq1gMPK5/6++VadNRh+faxec15sRQ8Bnyfp6/hbYFNgP+BCSZtGxN0kZXjzgGV719HqwP9Wn1UbYtQ11pvf+34w0nZoI9JZPY4BfgzcATwX+ChwhqQdIuKsNI3LDD9sTR1J25MMdjCXsyNimyLrjoj+UQdPkvS/wAclfTYifgWol3xQ9opsbxZUWF55y8HlVYJxyjEi9upbx/EkHYcPYlmTi1HlY/Xx9dFSEXE5Sf/UnrMlnQP8nGTQjP1Jyspl2F55y8fl2CIR8VeZX8+V9EOSCPInWdbE3WWGH7am0QXAi3Ok6w/Xj+tbwAeBzUkGWRhVW7E6y0LHlqiqvEaVw9Mz811e5Zi4HCPicUnfAT4raZ2IuJPly+fOTPJeed2PVS0bCc5ancG18dawiLhM0nXAK9JJw6L0vciky7FZ9wODXqWwemZ+73t1SeqLlPSnswZExP9IOhH4m8xklxl+2Jo6aZvna2rcZH/N+w0k7XM3Xi6RNB94Psmw1ZaqsLx6D0kbAxf2JqbvylgVuDqd5PIqQYnl2H89Zcsx+7DVa+d+NVa1q+i7PlIb4ePfZtka9auA1w9IsxFwa0S4CWGzrgJ2lrRqX7+tjUgGaro+k24Vkr7G1/elA1+PbdAfyXKZ4QEybHJ/SXJhXQwQEY8CJwN/ISn7ML8LyQV3Qu05nEERcStwBfCuvlnvJhn16SdpOpdXS6TH/20k//zdlU6+kGSI90HleD/JiHhWrROALSQ9vzchrbTYGl8frSRpc+CFwM/SSScAz5L02kyaPwD+DJdhG5xAMkru23oT0vvh24GfRsSSdPLJJA9fg+6HV6YDQFlD0mvqTSy77sBlBjiyZTz5h2l9lj18byRpl/TnkyJisaTnkgw3fRxJ7cQqJANkvAc4LCJuyKxyEck/id+W9JV03QcD342ISyvdmRmQp7zSn/cDfizpMJLmnpuS9F84JPPPPLi8aifpnSRDTp8E/JpkgIz3A5sB7+yli4ilkj5O8hLj24HTgO2Avwb2Sh+WrVpfIxm85IeS9iepXDqQpNwOazJjBpKOJRlk5jLgQZL73L7A7cC/pMlOILnHHSPpoyTNBvclqYX/XN15njWZv0+bpd9vkHQPcE9EnB0Rv5D0n8CX0teR3ATsSfI+uyf/SY+IuyV9EdhX0v+QlPnbSe6JHsK/RHOVmaSPAC8CzmTZABkfAdbGZbaipl/05U/zH5KXFceQz/ppmjWAHwC3AI8AD5NcNB8g8wK8zDpfQ/LH7RHgN8CXgFWb3tdp+OQpr0zaPyeJcC0BbiV5qfE8l1fjZbgFcEZ6rJeSjKh2GrDjkPTvA65Ly/FXwN81vQ+z9CHpT/I9ktHu/ie9F67fdL78CUgemv4rvYaWkjwEHw6s05duDZKXgd9P0nfydGCTpvM/C58Rf6/OyqR5GvAF4K7079DPgG0GrGseSaXhLen98L+AXZrex2n7zFVmJFHh80laXiwF7iOp1Hily2zFj9IDYWZmZmZmZiVyny0zMzMzM7MK+GHLzMzMzMysAn7YMjMzMzMzq4AftszMzMzMzCrghy0zMzMzM7MK+GHLzMzMzMysAn7YMjMzMzMzq4AftszMzMzMzCrghy0zsxaRtKWk4yTdJulRSb+VdLGkAyWtM8b6Fkmq5O31knaS9OGqtjnJevLmbVi6Mkg6S9JZJa5vI0kh6fUlrKvUvJmZ2WB+2DIzawlJewPnA88A9ge2B94BnALsDnyjudwNtBMw6EHlCGDLmvPSL2/ehqVro5en35c1mgszM8ttpaYzYGZmIGlb4GDgkIj4UN/skyQdBLyt/pwVFxG3Abc1nY9B2py3HDYDbomIe5vOyLgkrRIRS5rOh5lZXRzZMjNrh32Ae9PvFUTE7yLiyN7vkjaUdLSkmyQ9LOlGSV+VtPqojRRZTtImko6XdF+a9lpJ+6bzjgR2BZ6VNm0LSTen81Zo/jdqXUXkyX/evI1K15uf/T0zfYUmeJLeIekaSUskXSVp5yH530TSCZIeSPN/vqRX59z9zRgjqpU3b3nzJ+md6foekfTfkt485JgsSo/pSySdIul/gW/n3U6B/LwwPbfuTvN0q6TvSHKFspk1zjciM7OGpf8Uvhb4fkQ8mnOxdUkiNB8EHgCeD+wHnMToJny5lpP0SuAs4HrgQ+kyLwBemiY5kKS54yuAN6fTBkYscqyriDz5z5u33PswiqTtgf8ATgT2Ttd5CLAycG0m3cuBc4HLgd2AxcAewGmStoqIS0dsQ8DLgM9Wkbe8+ZO0A3AscEK6vjWBLwHzgeuGZOOHwNfTvD+R9zgUOF4/Bh4E9iSpsHgW8EZcoWxmLeCHLTOz5i0g+Wf11v4Z/bXzEfFY+n0OcE4m3QUkDzPnSto0Ii4ftKECy/0zcB+wRUQsTqedkVnPDZLuAR6NiIvm2L+R6yoiT/7z5q3gPoxyAHAN8JaIeCLN1y+Bi1j+geZgkjLervdQLekU4Erg4yT9x4Z5IfB/gKEPZBPmLW/+DgCuBnaOiF6E8L/TfA172PpyRBzS+0XS6Tm2kys/ktYkeXB/S0SckNnmf+Q9QGZmVXKtj5lZ8zRworQ2sDT76T18SXqqpP3S5lwPp/PPTRd90dAN5VhO0qrA1sCxmYej8XasxHWl6xtrv6siaR5JZOy7vYcZgIj4GXBzJt3TSKKX3yGJ7qyUlqWA04DXzLGpzdLv3M0I8+Ytb/7S9W0OfK/3oJWu7zLgphFZOb7IdoqkI3mIvxH4jKTdJL0g7/ExM6uDH7bMzJp3L/AI8JwB01+Rfr7WN+8gYBFwDPAm4JXAn6fz5o/YVp7lVif5+1DGQBJlrgvG3++qrEnSJO83A+Zlp60BzCOJyCzt+3wAWF3SqL/JLwdui4i7K8hbrvxl1jcoD4O20XNnke2kxyFXuvShbwfgEpJz4zol/fj2HJEfM7PauBmhmVnDIuIxSecAO0h6aq/JVNpk8BIASX/at9g7gKMi4pO9CZJ+P8fm8iz3APAESd+XSZW5Lhh/v8fxCPDUAdMXkERUIHkgXgqsNSDdWsAt6c8PkhyHrwBHDdpYNvo0wGYUb0KYN2+58pdZ3zOHrG+FZrCp7GApuY6DpNzHKyJuBBam/do2IXkYO1TSzRHxkyF5MjOrhSNbZmbt8DmSyEHeARBWJfnHN+u9ZSyXNvc7D3h32pxrmCXAqPlF1pVX3v2eM2850t0CrJX2CwJA0gZkmitGxOPAxcAu2ciUpFcB62fS/Y6kueMmwGURcUn/Z1gG04eITSk4EmHevOXNX7q+S4C3pnnqrW8z4Hk585TrOIxzvCLxC5a9N+0lefJkZlYlR7bMzFogIk6X9I8kfU9eSlKbfxNJ07gXkkR0fseyKMHJwK7p4ATXkzSl2yrHpvIu9xHgbOBCSZ8naQb4fOBlEbFXmuZqYI20ydYlwCMR8d9jriuvvPnPm7dR6b5DMmLhsZK+QPIwvC9JhCfrn4CfAj+QdBjJiH8HAHf1pfswyeAep0j6OknzujVJmgjOi4h/HLLPGwCrkYlsSVqf5Pw4ICIWDVmuSN7y5q+3vuMlHZ7OX5Sub1Rkruh2cqVLr5VDgP8kOR/mAe8BHmPMQVjMzEoVEf74448//rTkQzKYxLeB24FHgd+SRCcOANbJpFsTOI6kmd4DJMNxv4LkYew9mXSLklt9seXStJsCPyJp+vUwyah2+2Tm/x7wrXQ9Adw8aJt51jXkWAxaT979zpW3Yeky83ciGf3uYeAK4PUkw9if1ZfunSSj+y0BrgJ2HpLuxWn+707T3kYyjPobRxyHt6d5y5b/xum0PXKcU7nyljd/wF8OWN/lwPGDyg9YaZzt5ElH0qTxmyQjIS4G7id5sN+x6WvZH3/88SciUMRy7500MzOzlpO0O/Ap4LlRwiiPE+bl2SRRpU9FxIFN5sXMrG3cjNDMzKx7Xgt8se4HrbTf3RdIhl+/l6Q56P8liSodUWdezMy6wJEtMzMzy0XSU0n6R21BMipjbyCL/SLiyibzZmbWRn7YMjMzMzMzq4CHfjczMzMzM6uAH7bMzMzMzMwq4IctMzMzMzOzCvhhy8zMzMzMrAJ+2DIzMzMzM6uAH7bMzMzMzMwq4IctMzMzMzOzCvx/K1f84WiVvu0AAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 1008x504 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure(figsize=(14,7))\n",
    "plt.hist2d(-l,b,bins=[np.linspace(-180,180,360),np.linspace(-90,90,180)],vmax=10)\n",
    "plt.tick_params(axis='both', which='major', labelsize=16)\n",
    "plt.xlabel('Galactic latitude $l$, degrees',fontsize=16)\n",
    "plt.ylabel('Galactic longitude $b$, degrees',fontsize=16)\n",
    "\n",
    "plt.savefig('Map.png',format='png',dpi=100)\n",
    "\n",
    "\n"
   ]
  },
  {
   "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.7.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}