Fermi-LAT via FermiPyLike¶
In this Example we show how to use the fermipy plugin in threeML. We perform a Binned likelihood analysis and a Bayesian analysis of the Crab, optimizing the parameters of the Crab Pulsar (PSR J0534+2200) keeping fixed the parameters of the Crab Nebula. In the model, the nebula is described by two sources, one representing the synchrotron spectrum, the othet the Inverse Compton emission. In this example we show how to download Fermi-LAT data, how to build a model starting from the 4FGL, how to free and fix parameters of the sources in the model, and how to perform a spectral analysis using the fermipy plugin.
[1]:
import shutil
from IPython.display import Image,display
import glob
import matplotlib as mpl
from matplotlib import pyplot as plt
from astropy.io import fits as pyfits
import numpy as np
from threeML import *
import scipy as sp
%matplotlib inline
/Users/omodei/miniconda/envs/threeML/lib/python2.7/site-packages/astromodels/core/parameter.py:555: UserWarning: We have set the min_value of K to 1e-99 because there was a postive transform
warnings.warn('We have set the min_value of %s to 1e-99 because there was a postive transform' % self.path)
/Users/omodei/miniconda/envs/threeML/lib/python2.7/site-packages/astromodels/core/parameter.py:555: UserWarning: We have set the min_value of xc to 1e-99 because there was a postive transform
warnings.warn('We have set the min_value of %s to 1e-99 because there was a postive transform' % self.path)
WARNING UserWarning: Using default configuration from /Users/omodei/miniconda/envs/threeML/lib/python2.7/site-packages/threeML/data/threeML_config.yml. You might want to copy it to /Users/omodei/.threeML/threeML_config.yml to customize it and avoid this warning.
WARNING RuntimeWarning: numpy.dtype size changed, may indicate binary incompatibility. Expected 96, got 88
WARNING RuntimeWarning: numpy.dtype size changed, may indicate binary incompatibility. Expected 96, got 88
WARNING RuntimeWarning: numpy.dtype size changed, may indicate binary incompatibility. Expected 96, got 88
WARNING RuntimeWarning: numpy.dtype size changed, may indicate binary incompatibility. Expected 96, got 88
WARNING RuntimeWarning: numpy.dtype size changed, may indicate binary incompatibility. Expected 96, got 88
WARNING RuntimeWarning: numpy.dtype size changed, may indicate binary incompatibility. Expected 96, got 88
WARNING RuntimeWarning: numpy.dtype size changed, may indicate binary incompatibility. Expected 96, got 88
WARNING RuntimeWarning: numpy.dtype size changed, may indicate binary incompatibility. Expected 96, got 88
WARNING RuntimeWarning: numpy.dtype size changed, may indicate binary incompatibility. Expected 96, got 88
WARNING RuntimeWarning: numpy.dtype size changed, may indicate binary incompatibility. Expected 96, got 88
WARNING RuntimeWarning: numpy.dtype size changed, may indicate binary incompatibility. Expected 96, got 88
WARNING RuntimeWarning: numpy.dtype size changed, may indicate binary incompatibility. Expected 96, got 88
WARNING RuntimeWarning: numpy.dtype size changed, may indicate binary incompatibility. Expected 96, got 88
INFO:keyring.backend:Loading KWallet
INFO:keyring.backend:Loading SecretService
INFO:keyring.backend:Loading Windows
INFO:keyring.backend:Loading chainer
INFO:keyring.backend:Loading macOS
WARNING RuntimeWarning: Env. variable OMP_NUM_THREADS is not set. Please set it to 1 for optimal performances in 3ML
WARNING RuntimeWarning: Env. variable MKL_NUM_THREADS is not set. Please set it to 1 for optimal performances in 3ML
WARNING RuntimeWarning: Env. variable NUMEXPR_NUM_THREADS is not set. Please set it to 1 for optimal performances in 3ML
The Fermi 4FGL catalog¶
Let’s interrogate the 4FGL to get the sources in a radius of 20.0 deg around the Crab
[2]:
lat_catalog = FermiLATSourceCatalog()
ra, dec, table = lat_catalog.search_around_source("Crab", radius=20.0)
table
Trying https://heasarc.gsfc.nasa.gov/cgi-bin/vo/cone/coneGet.pl?table=fermilpsc&
[2]:
<Table masked=True length=133>
| name | source_type | ra | dec | assoc_name | tevcat_assoc | Search_Offset |
|---|---|---|---|---|---|---|
| deg | deg | |||||
| object | str32 | float64 | float64 | object | object | float64 |
| 4FGL J0534.5+2200 | pulsar, identified by pulsations | 83.6367 | 22.0149 | PSR J0534+2200 | Crab pulsar | 0.2022 |
| 4FGL J0534.5+2201s | pulsar wind nebula | 83.6331 | 22.0199 | Crab Nebula | Crab | 0.3243 |
| 4FGL J0534.5+2201i | pulsar wind nebula | 83.6330 | 22.0200 | Crab Nebula | Crab | 0.3304 |
| 4FGL J0526.3+2246 | active galaxy of uncertain type | 81.5908 | 22.7778 | 2MASS J05262938+2247232 | 122.1994 | |
| 4FGL J0544.4+2238 | unknown | 86.1093 | 22.6418 | 142.4913 | ||
| 4FGL J0521.7+2112 | BL Lac type of blazar | 80.4445 | 21.2131 | TXS 0518+211 | VER J0521+211 | 184.2435 |
| 4FGL J0536.2+1733 | BL Lac type of blazar | 84.0719 | 17.5534 | TXS 0533+175 | 268.8090 | |
| 4FGL J0539.0+1644 | active galaxy of uncertain type | 84.7500 | 16.7432 | NVSS J053855+164612 | 322.5278 | |
| 4FGL J0550.9+2552 | unknown | 87.7457 | 25.8710 | 323.0606 | ||
| ... | ... | ... | ... | ... | ... | ... |
| 4FGL J0512.8+4041 | active galaxy of uncertain type | 78.2186 | 40.6945 | B3 0509+406 | 1153.9906 | |
| 4FGL J0639.6+3503 | active galaxy of uncertain type | 99.9024 | 35.0586 | B2 0635+35 | 1157.9315 | |
| 4FGL J0552.8+0313 | active galaxy of uncertain type | 88.2180 | 3.2322 | PKS 0550+032 | 1158.1471 | |
| 4FGL J0653.6+1636 | active galaxy of uncertain type | 103.4105 | 16.6106 | 2MASX J06533986+1636432 | 1164.7254 | |
| 4FGL J0658.7+2318 | unknown | 104.6808 | 23.3027 | 1166.9322 | ||
| 4FGL J0555.1+0304 | active galaxy of uncertain type | 88.7776 | 3.0710 | GB6 J0555+0304 | 1175.4617 | |
| 4FGL J0658.2+2709 | active galaxy of uncertain type | 104.5735 | 27.1501 | B2 0655+27A | 1181.6921 | |
| 4FGL J0642.4+1048 | unknown | 100.6081 | 10.8135 | 1183.9644 | ||
| 4FGL J0506.9+0323 | active galaxy of uncertain type | 76.7314 | 3.3917 | 2MASS J05065014+0323587 | 1187.4462 | |
| 4FGL J0409.2+2542 | unknown | 62.3144 | 25.7022 | 1189.1462 |
This gets a 3ML model (a Model instance) from the table above, where every source in the 4FGL becomes a Source instance. Note that by default all parameters of all sources are fixed.
[3]:
model = lat_catalog.get_model()
WARNING UserWarning: We have set the min_value of Crab_synch.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of Crab_IC.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _2MASS_J05262938p2247232.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0544d4p2238.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of TXS_0518p211.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of TXS_0533p175.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of NVSS_J053855p164612.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0550d9p2552.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0534d2p2751.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: Source 4FGL J0540.3+2756e is extended, support for extended source is not here yet. I will ignoreit
WARNING UserWarning: We have set the min_value of Sim_147.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _3C_138.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of MG2_J060149p2324.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of TXS_0506p196.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0533d9p2838.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4C_p27d15.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4C_p22d12.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of PKS_0507p17.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of NVSS_J052429p283918.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of TXS_0536p145.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of GB6_J0515p1527.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0609d0p2136.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _1RXS_J050107d1p242318.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of NVSS_J060900p203746.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0609d0p2006.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of PKS_0528p134.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _1RXS_J045931d5p192242.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0531d7p1241c.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of GB6_J0614p2426.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0616d5p2235.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: Source 4FGL J0617.2+2234e is extended, support for extended source is not here yet. I will ignoreit
WARNING UserWarning: We have set the min_value of IC_443.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of MG2_J061020p1701.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _2MASS_J05400502p1209550.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0613d1p1749.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of MG2_J045613p2702.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0618d9p2240c.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0527d1p1152.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of GB6_J0611p2803.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _2MASS_J05594100p3042282.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0620d1p2246.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0502d2p3016.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of NVSS_J060015p124344.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0447d2p2446.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0612d6p1520c.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of PKS_0459p135.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of GB6_J0549p3258.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of GB6_J0501p3048.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0453d3p2843.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of RX_J0620d6p2644.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of NVSS_J060922p140744.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of RX_J0621d7p1747.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0532d6p3358.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0537d5p0959.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0535d3p0934.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0512d3p3329.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of GB6_J0529p0934.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of SSTSL2_J043903d89p215240d5.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of GB6_J0629p2437.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _2MASX_J05054288p3321355.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of NVSS_J060835p114942.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0620d4p1445.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of MG2_J044950p3026.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0445d7p1535.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of PKS_0506p101.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of TXS_0628p203.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of NVSS_J045804p115142.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of B2_0437p27B.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of PMN_J0517p0858.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of TXS_0516p087.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0538d9p3549c.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of NVSS_J054218p354149.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of GB6_J0623p3045.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0622d5p3120.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0540d7p3611.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _1RXS_J044127d8p150455.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0618d7p1211.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of MG2_J050234p3436.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of OG_050.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of TXS_0437p145.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0436d9p2915.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0540d2p0655.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of PKS_0446p11.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of TXS_0429p174.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of RGB_J0640p245B.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of MG2_J043337p2905.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of B2_0619p33.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of RX_J0638d6p2751.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0534d0p3746c.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of MG1_J043022p1655.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of NVSS_J054529p061955.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of MG2_J042948p2843.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of B3_0548p378B.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of B2_0441p34.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of TXS_0420p210.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0435d5p3232.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of MS_0419d3p1943.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of NVSS_J043307p322840.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of GB6_J0649p2228.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of TXS_0506p056.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of B2_0557p38.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of PKS_0459p060.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4C_p06d21.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _1ES_0647p250.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of NVSS_J065035p205556.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0551d7p0446.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of GB6_J0648p1749.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0442d8p3609.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of B2_0552p39A.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _1RXS_J064814d1p160708.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of MG2_J065230p1934.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of PKS_0502p049.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of RX_J0648d7p1516.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of GB6_J0620p3806.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of GB6_J0653p2816.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of TXS_0442p071.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of MG1_J050533p0415.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of TXS_0431p092.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of B3_0509p406.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of B2_0635p35.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of PKS_0550p032.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _2MASX_J06533986p1636432.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0658d7p2318.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of GB6_J0555p0304.spectrum.main.Log_parabola.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of B2_0655p27A.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0642d4p1048.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _2MASS_J05065014p0323587.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
WARNING UserWarning: We have set the min_value of _4FGL_J0409d2p2542.spectrum.main.Powerlaw.K to 1e-99 because there was a postive transform
Let’s free all the normalizations within 3 deg from the center.
[4]:
model.free_point_sources_within_radius(3.0, normalization_only=True)
model.display()
Model summary:
Free parameters (5):
Fixed parameters (707):
(abridged. Use complete=True to see all fixed parameters)
Linked parameters (0):
(none)
Independent variables:
(none)
| N | |
|---|---|
| Point sources | 133 |
| Extended sources | 0 |
| Particle sources | 0 |
Free parameters (5):
| value | min_value | max_value | unit | |
|---|---|---|---|---|
| PSR_J0534p2200.spectrum.main.Super_cutoff_powerlaw.K | 1.61768e-13 | 1.61768e-16 | 1.61768e-10 | keV-1 s-1 cm-2 |
| Crab_synch.spectrum.main.Powerlaw.K | 2.3083e-10 | 2.3083e-13 | 2.3083e-07 | keV-1 s-1 cm-2 |
| Crab_IC.spectrum.main.Log_parabola.K | 5.5e-16 | 5.5e-19 | 5.5e-13 | keV-1 s-1 cm-2 |
| _2MASS_J05262938p2247232.spectrum.main.Powerlaw.K | 1.3747e-17 | 1.3747e-20 | 1.3747e-14 | keV-1 s-1 cm-2 |
| _4FGL_J0544d4p2238.spectrum.main.Powerlaw.K | 1.565e-16 | 1.565e-19 | 1.565e-13 | keV-1 s-1 cm-2 |
Fixed parameters (707):
(abridged. Use complete=True to see all fixed parameters)
Linked parameters (0):
(none)
Independent variables:
(none)
but then let’s fix the sync and the IC components of the Crab nebula (cannot fit them with just one month of data) (these two methods are equivalent)
[5]:
model['Crab_IC.spectrum.main.Log_parabola.K'].fix = True
model.Crab_synch.spectrum.main.Powerlaw.K.fix = True
However, let’s free the index of the Crab Pulsar
[6]:
model.PSR_J0534p2200.spectrum.main.Super_cutoff_powerlaw.index.free = True
model.display()
Model summary:
Free parameters (4):
Fixed parameters (708):
(abridged. Use complete=True to see all fixed parameters)
Linked parameters (0):
(none)
Independent variables:
(none)
| N | |
|---|---|
| Point sources | 133 |
| Extended sources | 0 |
| Particle sources | 0 |
Free parameters (4):
| value | min_value | max_value | unit | |
|---|---|---|---|---|
| PSR_J0534p2200.spectrum.main.Super_cutoff_powerlaw.K | 1.61768e-13 | 1.61768e-16 | 1.61768e-10 | keV-1 s-1 cm-2 |
| PSR_J0534p2200.spectrum.main.Super_cutoff_powerlaw.index | -1.8869 | -10 | 10 | |
| _2MASS_J05262938p2247232.spectrum.main.Powerlaw.K | 1.3747e-17 | 1.3747e-20 | 1.3747e-14 | keV-1 s-1 cm-2 |
| _4FGL_J0544d4p2238.spectrum.main.Powerlaw.K | 1.565e-16 | 1.565e-19 | 1.565e-13 | keV-1 s-1 cm-2 |
Fixed parameters (708):
(abridged. Use complete=True to see all fixed parameters)
Linked parameters (0):
(none)
Independent variables:
(none)
[7]:
# Download data from Jan 01 2010 to February 1 2010
tstart = "2010-01-01 00:00:00"
tstop = "2010-02-01 00:00:00"
# Note that this will understand if you already download these files, and will
# not do it twice unless you change your selection or the outdir
evfile, scfile = download_LAT_data(
ra,
dec,
20.0,
tstart,
tstop,
time_type="Gregorian",
destination_directory="Crab_data",
)
Existing event file /Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/Crab_data/L2004131748578D8991BF28_PH00.fits and Spacecraft file /Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/Crab_data/L2004131748578D8991BF28_SC00.fits correspond to the same selection. We assume you did not tamper with them, so we will return those instead of downloading them again. If you want to download them again, remove them from the outdir
Configuration for Fermipy¶
3ML provides and intreface into Fermipy via the FermipyLike plugin. We can use it to generate basic configuration files.
.. note:: Currently, the FermipyLike plugin does not provide an interface to handle extended sources. This will change[8]:
config = FermipyLike.get_basic_config(evfile=evfile, scfile=scfile, ra=ra, dec=dec)
# See what we just got
config.display()
binning:
roiwidth: 10.0
binsz: 0.1
binsperdec: 8
data:
evfile: /Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/Crab_data/L2004131748578D8991BF28_PH00.fits
scfile: /Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/Crab_data/L2004131748578D8991BF28_SC00.fits
selection:
emin: 100.0
emax: 100000.0
zmax: 100.0
evclass: 128
evtype: 3
filter: DATA_QUAL>0 && LAT_CONFIG==1
ra: 83.63309095468972
dec: 22.014494786634696
binning: roiwidth: 10.0 binsz: 0.1 binsperdec: 8 data: evfile: /Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/Crab_data/L2004131748578D8991BF28_PH00.fits scfile: /Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/Crab_data/L2004131748578D8991BF28_SC00.fits selection: emin: 100.0 emax: 100000.0 zmax: 100.0 evclass: 128 evtype: 3 filter: DATA_QUAL>0 && LAT_CONFIG==1 ra: 83.63309095468972 dec: 22.014494786634696
You can of course modify the configuration as a dictionary
[9]:
config["selection"]["emax"] = 300000.0
and even add sections
[10]:
config["gtlike"] = {"edisp": False}
config.display()
binning:
roiwidth: 10.0
binsz: 0.1
binsperdec: 8
data:
evfile: /Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/Crab_data/L2004131748578D8991BF28_PH00.fits
scfile: /Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/Crab_data/L2004131748578D8991BF28_SC00.fits
gtlike:
edisp: false
selection:
emin: 100.0
emax: 300000.0
zmax: 100.0
evclass: 128
evtype: 3
filter: DATA_QUAL>0 && LAT_CONFIG==1
ra: 83.63309095468972
dec: 22.014494786634696
binning: roiwidth: 10.0 binsz: 0.1 binsperdec: 8 data: evfile: /Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/Crab_data/L2004131748578D8991BF28_PH00.fits scfile: /Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/Crab_data/L2004131748578D8991BF28_SC00.fits gtlike: edisp: false selection: emin: 100.0 emax: 300000.0 zmax: 100.0 evclass: 128 evtype: 3 filter: DATA_QUAL>0 && LAT_CONFIG==1 ra: 83.63309095468972 dec: 22.014494786634696
FermipyLike¶
Let’s create an instance of the plugin/ Note that here no processing is made, because fermipy still doesn’t know about the model you want to use.
[11]:
LAT = FermipyLike("LAT", config)
WARNING RuntimeWarning: numpy.dtype size changed, may indicate binary incompatibility. Expected 96, got 88
WARNING RuntimeWarning: numpy.dtype size changed, may indicate binary incompatibility. Expected 96, got 88
WARNING RuntimeWarning: numpy.dtype size changed, may indicate binary incompatibility. Expected 96, got 88
The plugin modifies the configuration as needed to get the output files in a unique place, which will stay the same as long as your selection does not change.
[12]:
config.display()
binning:
roiwidth: 10.0
binsz: 0.1
binsperdec: 8
data:
evfile: /Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/Crab_data/L2004131748578D8991BF28_PH00.fits
scfile: /Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/Crab_data/L2004131748578D8991BF28_SC00.fits
fileio:
outdir: __6b76c18fe898c8f0dceae172dc081dd6
gtlike:
edisp: false
selection:
emin: 100.0
emax: 300000.0
zmax: 100.0
evclass: 128
evtype: 3
filter: DATA_QUAL>0 && LAT_CONFIG==1
ra: 83.63309095468972
dec: 22.014494786634696
binning: roiwidth: 10.0 binsz: 0.1 binsperdec: 8 data: evfile: /Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/Crab_data/L2004131748578D8991BF28_PH00.fits scfile: /Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/Crab_data/L2004131748578D8991BF28_SC00.fits fileio: outdir: __6b76c18fe898c8f0dceae172dc081dd6 gtlike: edisp: false selection: emin: 100.0 emax: 300000.0 zmax: 100.0 evclass: 128 evtype: 3 filter: DATA_QUAL>0 && LAT_CONFIG==1 ra: 83.63309095468972 dec: 22.014494786634696
Here is where the fermipy processing happens (the .setup method)
[14]:
fermipy_output_directory = config['fileio']['outdir']
print('Fermipy Output directoty: %s' % fermipy_output_directory)
#This remove the output directory, to start a fresh analysis...
if os.path.exists(fermipy_output_directory): shutil.rmtree(fermipy_output_directory)
# Here is where the fermipy processing happens (the .setup method)
data = DataList(LAT)
jl = JointLikelihood(model, data)
Fermipy Output directoty: __6b76c18fe898c8f0dceae172dc081dd6
Found Galactic template for IRF. P8R2_SOURCE_V6: /Users/omodei/miniconda/envs/threeML/share/fermitools/refdata/fermi/galdiffuse/gll_iem_v07.fits
Cutting the template around the ROI:
2020-04-17 14:54:16 INFO GTAnalysis.__init__():
--------------------------------------------------------------------------------
fermipy version 0.17.3
ScienceTools version unknown
Found Isotropic template for irf P8R2_SOURCE_V6: /Users/omodei/miniconda/envs/threeML/share/fermitools/refdata/fermi/galdiffuse/iso_P8R2_SOURCE_V6_v06.txt
2020-04-17 14:54:23 INFO GTAnalysis.setup(): Running setup.
2020-04-17 14:54:23 INFO GTBinnedAnalysis.setup(): Running setup for component 00
2020-04-17 14:54:23 INFO GTBinnedAnalysis.run_gtapp(): Running gtselect.
2020-04-17 14:54:23 INFO GTBinnedAnalysis.run_gtapp(): time -p gtselect infile=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/Crab_data/L2004131748578D8991BF28_PH00.fits outfile=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/ft1_00.fits ra=83.6330909547 dec=22.0144947866 rad=7.57106781187 tmin="INDEF" tmax="INDEF" emin=100.0 emax=300000.0 zmin=0.0 zmax=100.0 evclass=128 evtype=3 convtype=-1 phasemin=0.0 phasemax=1.0 evtable="EVENTS" chatter=3 clobber=yes debug=no gui=no mode="ql"
2020-04-17 14:54:24 INFO GTBinnedAnalysis.run_gtapp():
2020-04-17 14:54:24 INFO GTBinnedAnalysis.run_gtapp(): WARNING: version mismatch between CFITSIO header (v3.43) and linked library (v3.41).
2020-04-17 14:54:24 INFO GTBinnedAnalysis.run_gtapp():
2020-04-17 14:54:24 INFO GTBinnedAnalysis.run_gtapp(): This is gtselect version HEAD
2020-04-17 14:54:24 INFO GTBinnedAnalysis.run_gtapp(): Applying filter string: ((EVENT_CLASS&o200) != o0) && angsep(RA,DEC,83.63309095,22.01449479) < 7.571067812 && ((EVENT_TYPE&o3) != o0) && 100 < ENERGY && ENERGY <= 300000 && 0 < ZENITH_ANGLE && ZENITH_ANGLE <= 100 && gtifilter()
2020-04-17 14:54:24 INFO GTBinnedAnalysis.run_gtapp(): Done.
2020-04-17 14:54:24 INFO GTBinnedAnalysis.run_gtapp(): Finished gtselect. Execution time: 1.01 s
2020-04-17 14:54:24 INFO GTBinnedAnalysis.run_gtapp(): Running gtmktime.
2020-04-17 14:54:24 INFO GTBinnedAnalysis.run_gtapp(): time -p gtmktime scfile=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/Crab_data/L2004131748578D8991BF28_SC00.fits sctable="SC_DATA" filter="DATA_QUAL>0 && LAT_CONFIG==1" roicut=no evfile=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/ft1_00.fits evtable="EVENTS" outfile="/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/ft1_filtered_00.fits" apply_filter=yes overwrite=no header_obstimes=yes tstart=0.0 tstop=0.0 gtifile="default" chatter=2 clobber=yes debug=no gui=no mode="ql"
2020-04-17 14:54:24 INFO GTBinnedAnalysis.run_gtapp():
2020-04-17 14:54:24 INFO GTBinnedAnalysis.run_gtapp(): WARNING: version mismatch between CFITSIO header (v3.43) and linked library (v3.41).
2020-04-17 14:54:24 INFO GTBinnedAnalysis.run_gtapp():
2020-04-17 14:54:24 INFO GTBinnedAnalysis.run_gtapp(): Finished gtmktime. Execution time: 0.31 s
2020-04-17 14:54:24 INFO GTBinnedAnalysis.run_gtapp(): Running gtltcube.
2020-04-17 14:54:24 INFO GTBinnedAnalysis.run_gtapp(): time -p gtltcube evfile="/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/ft1_00.fits" evtable="EVENTS" scfile=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/Crab_data/L2004131748578D8991BF28_SC00.fits sctable="SC_DATA" outfile=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/ltcube_00.fits dcostheta=0.025 binsz=1.0 phibins=0 tmin=0.0 tmax=0.0 file_version="1" zmin=0.0 zmax=100.0 chatter=2 clobber=yes debug=no gui=no mode="ql"
2020-04-17 14:55:36 INFO GTBinnedAnalysis.run_gtapp():
2020-04-17 14:55:36 INFO GTBinnedAnalysis.run_gtapp(): WARNING: version mismatch between CFITSIO header (v3.43) and linked library (v3.41).
2020-04-17 14:55:36 INFO GTBinnedAnalysis.run_gtapp():
2020-04-17 14:55:36 INFO GTBinnedAnalysis.run_gtapp(): Working on file /Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/Crab_data/L2004131748578D8991BF28_SC00.fits
2020-04-17 14:55:36 INFO GTBinnedAnalysis.run_gtapp(): .....................!
2020-04-17 14:55:36 INFO GTBinnedAnalysis.run_gtapp(): Finished gtltcube. Execution time: 72.14 s
2020-04-17 14:55:37 INFO GTBinnedAnalysis.run_gtapp(): Running gtbin.
2020-04-17 14:55:37 INFO GTBinnedAnalysis.run_gtapp(): time -p gtbin evfile=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/ft1_00.fits scfile=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/Crab_data/L2004131748578D8991BF28_SC00.fits outfile=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/ccube_00.fits algorithm="ccube" ebinalg="LOG" emin=100.0 emax=300000.0 enumbins=28 ebinfile=NONE tbinalg="LIN" tbinfile=NONE nxpix=100 nypix=100 binsz=0.1 coordsys="CEL" xref=83.6330909547 yref=22.0144947866 axisrot=0.0 rafield="RA" decfield="DEC" proj="AIT" hpx_ordering_scheme="RING" hpx_order=3 hpx_ebin=yes hpx_region="" evtable="EVENTS" sctable="SC_DATA" efield="ENERGY" tfield="TIME" chatter=3 clobber=yes debug=no gui=no mode="ql"
2020-04-17 14:55:37 INFO GTBinnedAnalysis.run_gtapp():
2020-04-17 14:55:37 INFO GTBinnedAnalysis.run_gtapp(): WARNING: version mismatch between CFITSIO header (v3.43) and linked library (v3.41).
2020-04-17 14:55:37 INFO GTBinnedAnalysis.run_gtapp():
2020-04-17 14:55:37 INFO GTBinnedAnalysis.run_gtapp(): This is gtbin version HEAD
2020-04-17 14:55:37 INFO GTBinnedAnalysis.run_gtapp(): Finished gtbin. Execution time: 0.42 s
2020-04-17 14:55:37 INFO GTBinnedAnalysis.run_gtapp(): Running gtexpcube2.
2020-04-17 14:55:37 INFO GTBinnedAnalysis.run_gtapp(): time -p gtexpcube2 infile=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/ltcube_00.fits cmap=none outfile=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/bexpmap_00.fits irfs="P8R2_SOURCE_V6" evtype=3 edisp_bins=0 nxpix=360 nypix=180 binsz=1.0 coordsys="CEL" xref=0.0 yref=0.0 axisrot=0.0 proj="CAR" ebinalg="LOG" emin=100.0 emax=300000.0 enumbins=28 ebinfile="NONE" hpx_ordering_scheme="RING" hpx_order=6 bincalc="EDGE" ignorephi=no thmax=180.0 thmin=0.0 table="EXPOSURE" chatter=3 clobber=yes debug=no mode="ql"
2020-04-17 14:55:44 INFO GTBinnedAnalysis.run_gtapp():
2020-04-17 14:55:44 INFO GTBinnedAnalysis.run_gtapp(): WARNING: version mismatch between CFITSIO header (v3.43) and linked library (v3.41).
2020-04-17 14:55:44 INFO GTBinnedAnalysis.run_gtapp():
2020-04-17 14:55:44 INFO GTBinnedAnalysis.run_gtapp(): This is gtexpcube2 version HEAD
2020-04-17 14:55:44 INFO GTBinnedAnalysis.run_gtapp(): Using evtype=3 (i.e., FRONT/BACK irfs)
2020-04-17 14:55:44 INFO GTBinnedAnalysis.run_gtapp(): ResponseFunctions::load: IRF used: P8R2_SOURCE_V6
2020-04-17 14:55:44 INFO GTBinnedAnalysis.run_gtapp(): event_types: 0 1
2020-04-17 14:55:44 INFO GTBinnedAnalysis.run_gtapp(): Computing binned exposure map....................!
2020-04-17 14:55:44 INFO GTBinnedAnalysis.run_gtapp(): Finished gtexpcube2. Execution time: 7.09 s
2020-04-17 14:55:44 INFO GTBinnedAnalysis.run_gtapp(): Running gtexpcube2.
2020-04-17 14:55:44 INFO GTBinnedAnalysis.run_gtapp(): time -p gtexpcube2 infile=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/ltcube_00.fits cmap=none outfile=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/bexpmap_roi_00.fits irfs="P8R2_SOURCE_V6" evtype=3 edisp_bins=0 nxpix=100 nypix=100 binsz=0.1 coordsys="CEL" xref=83.6330909547 yref=22.0144947866 axisrot=0.0 proj="CAR" ebinalg="LOG" emin=100.0 emax=300000.0 enumbins=28 ebinfile="NONE" hpx_ordering_scheme="RING" hpx_order=6 bincalc="EDGE" ignorephi=no thmax=180.0 thmin=0.0 table="EXPOSURE" chatter=3 clobber=yes debug=no mode="ql"
2020-04-17 14:55:46 INFO GTBinnedAnalysis.run_gtapp():
2020-04-17 14:55:46 INFO GTBinnedAnalysis.run_gtapp(): WARNING: version mismatch between CFITSIO header (v3.43) and linked library (v3.41).
2020-04-17 14:55:46 INFO GTBinnedAnalysis.run_gtapp():
2020-04-17 14:55:46 INFO GTBinnedAnalysis.run_gtapp(): This is gtexpcube2 version HEAD
2020-04-17 14:55:46 INFO GTBinnedAnalysis.run_gtapp(): Using evtype=3 (i.e., FRONT/BACK irfs)
2020-04-17 14:55:46 INFO GTBinnedAnalysis.run_gtapp(): ResponseFunctions::load: IRF used: P8R2_SOURCE_V6
2020-04-17 14:55:46 INFO GTBinnedAnalysis.run_gtapp(): event_types: 0 1
2020-04-17 14:55:46 INFO GTBinnedAnalysis.run_gtapp(): Computing binned exposure map....................!
2020-04-17 14:55:46 INFO GTBinnedAnalysis.run_gtapp(): Finished gtexpcube2. Execution time: 1.50 s
WARNING: FITSFixedWarning: RADECSYS= 'FK5 '
the RADECSYS keyword is deprecated, use RADESYSa. [astropy.wcs.wcs]
WARNING:astropy:FITSFixedWarning: RADECSYS= 'FK5 '
the RADECSYS keyword is deprecated, use RADESYSa.
2020-04-17 14:55:46 INFO GTBinnedAnalysis.run_gtapp(): Running gtsrcmaps.
2020-04-17 14:55:46 INFO GTBinnedAnalysis.run_gtapp(): time -p gtsrcmaps scfile=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/Crab_data/L2004131748578D8991BF28_SC00.fits sctable="SC_DATA" expcube=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/ltcube_00.fits cmap=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/ccube_00.fits srcmdl=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/srcmdl_00.xml bexpmap=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/bexpmap_00.fits wmap=none outfile=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/srcmap_00.fits irfs="P8R2_SOURCE_V6" evtype=3 convol=yes resample=yes rfactor=2 minbinsz=0.05 ptsrc=yes psfcorr=yes emapbnds=no edisp_bins=0 copyall=no chatter=3 clobber=yes debug=no gui=no mode="ql"
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp():
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): WARNING: version mismatch between CFITSIO header (v3.43) and linked library (v3.41).
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp():
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): This is gtsrcmaps version HEAD
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Using evtype=3 (i.e., FRONT/BACK irfs)
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): ResponseFunctions::load: IRF used: P8R2_SOURCE_V6
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): event_types: 0 1
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named PSR_J0534p2200
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named Crab_synch
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named Crab_IC
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _2MASS_J05262938p2247232
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0544d4p2238
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named TXS_0518p211
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named TXS_0533p175
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named NVSS_J053855p164612
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0550d9p2552
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0534d2p2751
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named Sim_147
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _3C_138
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named MG2_J060149p2324
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named TXS_0506p196
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0533d9p2838
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4C_p27d15
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4C_p22d12
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named PKS_0507p17
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named NVSS_J052429p283918
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named TXS_0536p145
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named GB6_J0515p1527
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0609d0p2136
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _1RXS_J050107d1p242318
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named NVSS_J060900p203746
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0609d0p2006
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named PKS_0528p134
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _1RXS_J045931d5p192242
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0531d7p1241c
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named GB6_J0614p2426
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0616d5p2235
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named IC_443
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named MG2_J061020p1701
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _2MASS_J05400502p1209550
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0613d1p1749
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named MG2_J045613p2702
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named PSR_J0554p3107
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0618d9p2240c
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0527d1p1152
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named GB6_J0611p2803
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _2MASS_J05594100p3042282
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0620d1p2246
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0502d2p3016
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named NVSS_J060015p124344
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named PSR_J0621p2514
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0447d2p2446
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0612d6p1520c
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named PKS_0459p135
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named GB6_J0549p3258
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named GB6_J0501p3048
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0453d3p2843
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named RX_J0620d6p2644
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named NVSS_J060922p140744
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named RX_J0621d7p1747
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0532d6p3358
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0537d5p0959
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0535d3p0934
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0512d3p3329
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named GB6_J0529p0934
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named SSTSL2_J043903d89p215240d5
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named GB6_J0629p2437
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _2MASX_J05054288p3321355
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named NVSS_J060835p114942
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0620d4p1445
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named MG2_J044950p3026
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0445d7p1535
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named PKS_0506p101
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named TXS_0628p203
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named NVSS_J045804p115142
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named B2_0437p27B
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named PMN_J0517p0858
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named TXS_0516p087
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0538d9p3549c
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named NVSS_J054218p354149
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named GB6_J0623p3045
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0622d5p3120
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0540d7p3611
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _1RXS_J044127d8p150455
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0618d7p1211
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named MG2_J050234p3436
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named OG_050
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named PSR_J0633p1746
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named TXS_0437p145
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0436d9p2915
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0540d2p0655
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named PKS_0446p11
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named TXS_0429p174
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named RGB_J0640p245B
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named MG2_J043337p2905
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named B2_0619p33
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named RX_J0638d6p2751
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0534d0p3746c
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named MG1_J043022p1655
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named NVSS_J054529p061955
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named MG2_J042948p2843
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named B3_0548p378B
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named B2_0441p34
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named TXS_0420p210
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0435d5p3232
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named MS_0419d3p1943
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named NVSS_J043307p322840
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named PSR_J0605p3757
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named GB6_J0649p2228
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named TXS_0506p056
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named B2_0557p38
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named PKS_0459p060
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4C_p06d21
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _1ES_0647p250
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named NVSS_J065035p205556
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0551d7p0446
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named PSR_J0631p1036
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named GB6_J0648p1749
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0442d8p3609
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named B2_0552p39A
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _1RXS_J064814d1p160708
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named MG2_J065230p1934
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named PKS_0502p049
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named RX_J0648d7p1516
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named GB6_J0620p3806
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named PSR_J0622p3749
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named GB6_J0653p2816
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named TXS_0442p071
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named MG1_J050533p0415
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named TXS_0431p092
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named B3_0509p406
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named B2_0635p35
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named PKS_0550p032
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _2MASX_J06533986p1636432
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0658d7p2318
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named GB6_J0555p0304
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named B2_0655p27A
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0642d4p1048
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _2MASS_J05065014p0323587
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0409d2p2542
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named isodiff
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Creating source named galdiff
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for B2_0437p27B 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for B2_0441p34 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for B2_0552p39A 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for B2_0557p38 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for B2_0619p33 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for B2_0635p35 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for B2_0655p27A 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for B3_0509p406 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for B3_0548p378B 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for Crab_IC 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for Crab_synch 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for GB6_J0501p3048 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for GB6_J0515p1527 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for GB6_J0529p0934 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for GB6_J0549p3258 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for GB6_J0555p0304 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for GB6_J0611p2803 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for GB6_J0614p2426 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for GB6_J0620p3806 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for GB6_J0623p3045 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for GB6_J0629p2437 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for GB6_J0648p1749 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for GB6_J0649p2228 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for GB6_J0653p2816 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for IC_443 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for MG1_J043022p1655 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for MG1_J050533p0415 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for MG2_J042948p2843 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for MG2_J043337p2905 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for MG2_J044950p3026 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for MG2_J045613p2702 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for MG2_J050234p3436 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for MG2_J060149p2324 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for MG2_J061020p1701 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for MG2_J065230p1934 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for MS_0419d3p1943 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for NVSS_J043307p322840 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for NVSS_J045804p115142 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for NVSS_J052429p283918 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for NVSS_J053855p164612 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for NVSS_J054218p354149 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for NVSS_J054529p061955 29....................!
2020-04-17 14:56:49 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for NVSS_J060015p124344 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for NVSS_J060835p114942 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for NVSS_J060900p203746 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for NVSS_J060922p140744 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for NVSS_J065035p205556 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for OG_050 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for PKS_0446p11 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for PKS_0459p060 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for PKS_0459p135 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for PKS_0502p049 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for PKS_0506p101 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for PKS_0507p17 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for PKS_0528p134 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for PKS_0550p032 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for PMN_J0517p0858 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for PSR_J0534p2200 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for PSR_J0554p3107 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for PSR_J0605p3757 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for PSR_J0621p2514 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for PSR_J0622p3749 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for PSR_J0631p1036 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for PSR_J0633p1746 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for RGB_J0640p245B 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for RX_J0620d6p2644 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for RX_J0621d7p1747 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for RX_J0638d6p2751 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for RX_J0648d7p1516 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for SSTSL2_J043903d89p215240d5 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for Sim_147 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for TXS_0420p210 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for TXS_0429p174 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for TXS_0431p092 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for TXS_0437p145 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for TXS_0442p071 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for TXS_0506p056 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for TXS_0506p196 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for TXS_0516p087 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for TXS_0518p211 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for TXS_0533p175 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for TXS_0536p145 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for TXS_0628p203 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _1ES_0647p250 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _1RXS_J044127d8p150455 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _1RXS_J045931d5p192242 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _1RXS_J050107d1p242318 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _1RXS_J064814d1p160708 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _2MASS_J05065014p0323587 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _2MASS_J05262938p2247232 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _2MASS_J05400502p1209550 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _2MASS_J05594100p3042282 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _2MASX_J05054288p3321355 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _2MASX_J06533986p1636432 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _3C_138 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4C_p06d21 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4C_p22d12 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4C_p27d15 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0409d2p2542 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0435d5p3232 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0436d9p2915 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0442d8p3609 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0445d7p1535 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0447d2p2446 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0453d3p2843 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0502d2p3016 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0512d3p3329 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0527d1p1152 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0531d7p1241c 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0532d6p3358 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0533d9p2838 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0534d0p3746c 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0534d2p2751 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0535d3p0934 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0537d5p0959 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0538d9p3549c 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0540d2p0655 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0540d7p3611 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0544d4p2238 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0550d9p2552 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0551d7p0446 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0609d0p2006 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0609d0p2136 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0612d6p1520c 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0613d1p1749 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0616d5p2235 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0618d7p1211 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0618d9p2240c 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0620d1p2246 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0620d4p1445 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0622d5p3120 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0642d4p1048 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0658d7p2318 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for galdiff 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for isodiff 29....................!
2020-04-17 15:00:04 INFO GTBinnedAnalysis.run_gtapp(): Finished gtsrcmaps. Execution time: 257.87 s
2020-04-17 15:00:04 INFO GTBinnedAnalysis.setup(): Finished setup for component 00
2020-04-17 15:00:04 INFO GTBinnedAnalysis._create_binned_analysis(): Creating BinnedAnalysis for component 00.
WARNING: AstropyDeprecationWarning: "clobber" was deprecated in version 2.0 and will be removed in a future version. Use argument "overwrite" instead. [astropy.utils.decorators]
WARNING:astropy:AstropyDeprecationWarning: "clobber" was deprecated in version 2.0 and will be removed in a future version. Use argument "overwrite" instead.
2020-04-17 15:00:14 INFO GTAnalysis.setup(): Initializing source properties
2020-04-17 15:00:22 INFO GTAnalysis.setup(): Finished setup.
2020-04-17 15:00:22 INFO GTAnalysis.free_source(): Fixing parameters for PSR_J0534p2200 : ['Prefactor']
2020-04-17 15:00:22 INFO GTAnalysis.free_source(): Fixing parameters for Crab_synch : ['Prefactor']
2020-04-17 15:00:22 INFO GTAnalysis.free_source(): Fixing parameters for Crab_IC : ['Prefactor']
2020-04-17 15:00:22 INFO GTAnalysis.free_source(): Fixing parameters for _2MASS_J05262938p2247232: ['Prefactor']
2020-04-17 15:00:22 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0544d4p2238 : ['Prefactor']
2020-04-17 15:00:22 INFO GTAnalysis.free_source(): Fixing parameters for TXS_0518p211 : ['Prefactor']
2020-04-17 15:00:22 INFO GTAnalysis.free_source(): Fixing parameters for TXS_0533p175 : ['Prefactor']
2020-04-17 15:00:22 INFO GTAnalysis.free_source(): Fixing parameters for NVSS_J053855p164612 : ['Prefactor']
2020-04-17 15:00:22 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0550d9p2552 : ['Prefactor']
2020-04-17 15:00:22 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0534d2p2751 : ['Prefactor']
2020-04-17 15:00:22 INFO GTAnalysis.free_source(): Fixing parameters for Sim_147 : ['Prefactor']
2020-04-17 15:00:22 INFO GTAnalysis.free_source(): Fixing parameters for _3C_138 : ['Prefactor']
2020-04-17 15:00:22 INFO GTAnalysis.free_source(): Fixing parameters for MG2_J060149p2324 : ['Prefactor']
2020-04-17 15:00:22 INFO GTAnalysis.free_source(): Fixing parameters for TXS_0506p196 : ['Prefactor']
2020-04-17 15:00:22 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0533d9p2838 : ['Prefactor']
2020-04-17 15:00:22 INFO GTAnalysis.free_source(): Fixing parameters for _4C_p27d15 : ['Prefactor']
2020-04-17 15:00:22 INFO GTAnalysis.free_source(): Fixing parameters for _4C_p22d12 : ['Prefactor']
2020-04-17 15:00:22 INFO GTAnalysis.free_source(): Fixing parameters for PKS_0507p17 : ['Prefactor']
2020-04-17 15:00:22 INFO GTAnalysis.free_source(): Fixing parameters for NVSS_J052429p283918 : ['Prefactor']
2020-04-17 15:00:22 INFO GTAnalysis.free_source(): Fixing parameters for TXS_0536p145 : ['Prefactor']
2020-04-17 15:00:22 INFO GTAnalysis.free_source(): Fixing parameters for GB6_J0515p1527 : ['Prefactor']
2020-04-17 15:00:22 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0609d0p2136 : ['Prefactor']
2020-04-17 15:00:22 INFO GTAnalysis.free_source(): Fixing parameters for _1RXS_J050107d1p242318: ['Prefactor']
2020-04-17 15:00:22 INFO GTAnalysis.free_source(): Fixing parameters for NVSS_J060900p203746 : ['Prefactor']
2020-04-17 15:00:22 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0609d0p2006 : ['Prefactor']
2020-04-17 15:00:22 INFO GTAnalysis.free_source(): Fixing parameters for PKS_0528p134 : ['Prefactor']
2020-04-17 15:00:22 INFO GTAnalysis.free_source(): Fixing parameters for _1RXS_J045931d5p192242: ['Prefactor']
2020-04-17 15:00:22 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0531d7p1241c : ['Prefactor']
2020-04-17 15:00:22 INFO GTAnalysis.free_source(): Fixing parameters for GB6_J0614p2426 : ['Prefactor']
2020-04-17 15:00:22 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0616d5p2235 : ['Prefactor']
2020-04-17 15:00:22 INFO GTAnalysis.free_source(): Fixing parameters for IC_443 : ['Prefactor']
2020-04-17 15:00:22 INFO GTAnalysis.free_source(): Fixing parameters for MG2_J061020p1701 : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for _2MASS_J05400502p1209550: ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0613d1p1749 : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for MG2_J045613p2702 : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for PSR_J0554p3107 : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0618d9p2240c : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0527d1p1152 : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for GB6_J0611p2803 : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for _2MASS_J05594100p3042282: ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0620d1p2246 : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0502d2p3016 : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for NVSS_J060015p124344 : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for PSR_J0621p2514 : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0447d2p2446 : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0612d6p1520c : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for PKS_0459p135 : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for GB6_J0549p3258 : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for GB6_J0501p3048 : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0453d3p2843 : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for RX_J0620d6p2644 : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for NVSS_J060922p140744 : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for RX_J0621d7p1747 : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0532d6p3358 : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0537d5p0959 : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0535d3p0934 : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0512d3p3329 : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for GB6_J0529p0934 : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for SSTSL2_J043903d89p215240d5: ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for GB6_J0629p2437 : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for _2MASX_J05054288p3321355: ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for NVSS_J060835p114942 : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0620d4p1445 : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for MG2_J044950p3026 : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0445d7p1535 : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for PKS_0506p101 : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for TXS_0628p203 : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for NVSS_J045804p115142 : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for B2_0437p27B : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for PMN_J0517p0858 : ['Prefactor']
2020-04-17 15:00:23 INFO GTAnalysis.free_source(): Fixing parameters for TXS_0516p087 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0538d9p3549c : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for NVSS_J054218p354149 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for GB6_J0623p3045 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0622d5p3120 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0540d7p3611 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for _1RXS_J044127d8p150455: ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0618d7p1211 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for MG2_J050234p3436 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for OG_050 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for PSR_J0633p1746 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for TXS_0437p145 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0436d9p2915 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0540d2p0655 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for PKS_0446p11 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for TXS_0429p174 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for RGB_J0640p245B : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for MG2_J043337p2905 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for B2_0619p33 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for RX_J0638d6p2751 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0534d0p3746c : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for MG1_J043022p1655 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for NVSS_J054529p061955 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for MG2_J042948p2843 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for B3_0548p378B : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for B2_0441p34 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for TXS_0420p210 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0435d5p3232 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for MS_0419d3p1943 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for NVSS_J043307p322840 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for PSR_J0605p3757 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for GB6_J0649p2228 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for TXS_0506p056 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for B2_0557p38 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for PKS_0459p060 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for _4C_p06d21 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for _1ES_0647p250 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for NVSS_J065035p205556 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0551d7p0446 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for PSR_J0631p1036 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for GB6_J0648p1749 : ['Prefactor']
2020-04-17 15:00:24 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0442d8p3609 : ['Prefactor']
2020-04-17 15:00:25 INFO GTAnalysis.free_source(): Fixing parameters for B2_0552p39A : ['Prefactor']
2020-04-17 15:00:25 INFO GTAnalysis.free_source(): Fixing parameters for _1RXS_J064814d1p160708: ['Prefactor']
2020-04-17 15:00:25 INFO GTAnalysis.free_source(): Fixing parameters for MG2_J065230p1934 : ['Prefactor']
2020-04-17 15:00:25 INFO GTAnalysis.free_source(): Fixing parameters for PKS_0502p049 : ['Prefactor']
2020-04-17 15:00:25 INFO GTAnalysis.free_source(): Fixing parameters for RX_J0648d7p1516 : ['Prefactor']
2020-04-17 15:00:25 INFO GTAnalysis.free_source(): Fixing parameters for GB6_J0620p3806 : ['Prefactor']
2020-04-17 15:00:25 INFO GTAnalysis.free_source(): Fixing parameters for PSR_J0622p3749 : ['Prefactor']
2020-04-17 15:00:25 INFO GTAnalysis.free_source(): Fixing parameters for GB6_J0653p2816 : ['Prefactor']
2020-04-17 15:00:25 INFO GTAnalysis.free_source(): Fixing parameters for TXS_0442p071 : ['Prefactor']
2020-04-17 15:00:25 INFO GTAnalysis.free_source(): Fixing parameters for MG1_J050533p0415 : ['Prefactor']
2020-04-17 15:00:25 INFO GTAnalysis.free_source(): Fixing parameters for TXS_0431p092 : ['Prefactor']
2020-04-17 15:00:25 INFO GTAnalysis.free_source(): Fixing parameters for B3_0509p406 : ['Prefactor']
2020-04-17 15:00:25 INFO GTAnalysis.free_source(): Fixing parameters for B2_0635p35 : ['Prefactor']
2020-04-17 15:00:25 INFO GTAnalysis.free_source(): Fixing parameters for PKS_0550p032 : ['Prefactor']
2020-04-17 15:00:25 INFO GTAnalysis.free_source(): Fixing parameters for _2MASX_J06533986p1636432: ['Prefactor']
2020-04-17 15:00:25 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0658d7p2318 : ['Prefactor']
2020-04-17 15:00:25 INFO GTAnalysis.free_source(): Fixing parameters for GB6_J0555p0304 : ['Prefactor']
2020-04-17 15:00:25 INFO GTAnalysis.free_source(): Fixing parameters for B2_0655p27A : ['Prefactor']
2020-04-17 15:00:25 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0642d4p1048 : ['Prefactor']
2020-04-17 15:00:25 INFO GTAnalysis.free_source(): Fixing parameters for _2MASS_J05065014p0323587: ['Prefactor']
2020-04-17 15:00:25 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0409d2p2542 : ['Prefactor']
Performing the fit¶
[15]:
jl.set_minimizer("ROOT")
res = jl.fit()
WARNING UserWarning: get_number_of_data_points not implemented, values for statistical measurements such as AIC or BIC are unreliable
Best fit values:
| result | unit | |
|---|---|---|
| parameter | ||
| PSR_J0534p2200...K | (1.50 +/- 0.06) x 10^-13 | 1 / (cm2 keV s) |
| PSR_J0534p2200...index | -1.801 +/- 0.026 | |
| _2MASS_J05262938p2247232.spectrum.main.Powerlaw.K | (0.41 -0.32 +1.4) x 10^-16 | 1 / (cm2 keV s) |
| _4FGL_J0544d4p2238.spectrum.main.Powerlaw.K | (3.9 -2.3 +6) x 10^-16 | 1 / (cm2 keV s) |
Correlation matrix:
| 1.00 | 0.73 | -0.06 | -0.02 |
| 0.73 | 1.00 | 0.09 | 0.21 |
| -0.06 | 0.09 | 1.00 | 0.03 |
| -0.02 | 0.21 | 0.03 | 1.00 |
Values of -log(likelihood) at the minimum:
| -log(likelihood) | |
|---|---|
| LAT | 118116.776096 |
| total | 118116.776096 |
Values of statistical measures:
| statistical measures | |
|---|---|
| AIC | 236231.552191 |
| BIC | 236233.552191 |
Now let’s compute the errors on the best fit parameters
[16]:
res = jl.get_errors()
| result | unit | |
|---|---|---|
| parameter | ||
| PSR_J0534p2200...K | (1.50 +/- 0.06) x 10^-13 | 1 / (cm2 keV s) |
| PSR_J0534p2200...index | -1.801 +/- 0.026 | |
| _2MASS_J05262938p2247232.spectrum.main.Powerlaw.K | (4 -4 +7) x 10^-17 | 1 / (cm2 keV s) |
| _4FGL_J0544d4p2238.spectrum.main.Powerlaw.K | (3.9 -3.3 +4) x 10^-16 | 1 / (cm2 keV s) |
We might also want to look at the profile of the likelihood for each parameter.
[17]:
res = jl.get_contours(
model.PSR_J0534p2200.spectrum.main.Super_cutoff_powerlaw.index,-2.0,-1.6,30
)
[18]:
res[-1]
[18]:
Or we might want to produce a contour plot
[19]:
res = jl.get_contours(
'PSR_J0534p2200.spectrum.main.Super_cutoff_powerlaw.K',1.3e-13,1.7e-13, 20,
'PSR_J0534p2200.spectrum.main.Super_cutoff_powerlaw.index',-2.0,-1.6, 20
)
[20]:
res[-1]
[20]:
Pro-trick: We can also axcess the GTAnalysis object of fermipy:
[23]:
res = jl.fit()
LAT.gta.write_roi('test',make_plots=True)
plt.show()
Best fit values:
| result | unit | |
|---|---|---|
| parameter | ||
| PSR_J0534p2200...K | (1.50 +/- 0.06) x 10^-13 | 1 / (cm2 keV s) |
| PSR_J0534p2200...index | -1.802 +/- 0.026 | |
| _2MASS_J05262938p2247232.spectrum.main.Powerlaw.K | (0.41 -0.31 +1.4) x 10^-16 | 1 / (cm2 keV s) |
| _4FGL_J0544d4p2238.spectrum.main.Powerlaw.K | (3.9 -2.3 +6) x 10^-16 | 1 / (cm2 keV s) |
Correlation matrix:
| 1.00 | 0.73 | -0.06 | -0.02 |
| 0.73 | 1.00 | 0.09 | 0.21 |
| -0.06 | 0.09 | 1.00 | 0.02 |
| -0.02 | 0.21 | 0.02 | 1.00 |
Values of -log(likelihood) at the minimum:
| -log(likelihood) | |
|---|---|
| LAT | 118116.775933 |
| total | 118116.775933 |
Values of statistical measures:
| statistical measures | |
|---|---|
| AIC | 236231.551865 |
| BIC | 236233.551865 |
2020-04-17 15:19:36 INFO GTBinnedAnalysis.write_xml(): Writing /Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/test_00.xml...
2020-04-17 15:19:36 INFO GTAnalysis.write_fits(): Writing /Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/test.fits...
2020-04-17 15:19:41 INFO GTAnalysis.write_roi(): Writing /Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/test.npy...
All the plots are saved in the output directory as png files:
[22]:
pngs=glob.glob("%s/*png" % fermipy_output_directory)
for png in pngs:
print png
my_image=Image(png)
display(my_image)
__6b76c18fe898c8f0dceae172dc081dd6/test_model_map_2.000_5.477.png
__6b76c18fe898c8f0dceae172dc081dd6/test_counts_map_2.000_5.477.png
__6b76c18fe898c8f0dceae172dc081dd6/test_counts_spectrum.png
__6b76c18fe898c8f0dceae172dc081dd6/test_counts_map_yproj_2.000_5.477.png
__6b76c18fe898c8f0dceae172dc081dd6/test_counts_map_xproj_2.000_5.477.png
We can also plot the resulting model:
[24]:
energies=sp.logspace(1,8,100) *u.MeV
fig=plt.figure()
# we only want to visualize the relevant sources...
src_to_plot=['Crab','PSR_J0534p2200']
# Now loop over all point sources and plot them
for source_name, point_source in model.point_sources.iteritems():
for src in src_to_plot:
if src in source_name:
# Plot the sum of all components for this source
plt.loglog(energies, point_source(energies),label=source_name)
# If there is more than one component, plot them also separately
if len(point_source.components) > 1:
for component_name, component in point_source.components.iteritems():
plt.loglog(energies,component.shape(energies),'--',label="%s of %s" %(component_name, source_name))
# Add a legend
plt.legend(loc=0,frameon=False)
plt.xlabel("Energy (MeV)")
plt.ylabel(r"Flux (ph cm$^{-2}$ s$^{-1}$ keV$^{-1}$")
plt.ylim([1e-20,1e-3])
plt.show()
We can also do a bayesian analysis.
This will set priors based on the current defined min-max (log-uniform or uniform)
[25]:
model.PSR_J0534p2200.spectrum.main.Super_cutoff_powerlaw.K.set_uninformative_prior(
Log_uniform_prior
)
model.PSR_J0534p2200.spectrum.main.Super_cutoff_powerlaw.index.set_uninformative_prior(
Uniform_prior
)
model._2MASS_J05262938p2247232.spectrum.main.Powerlaw.K.set_uninformative_prior(
Log_uniform_prior
)
model._4FGL_J0544d4p2238.spectrum.main.Powerlaw.K.set_uninformative_prior(
Log_uniform_prior
)
[27]:
#It's better to remove the output rirectory,...
shutil.rmtree(fermipy_output_directory)
bayes = BayesianAnalysis(model, data)
Found Galactic template for IRF. P8R2_SOURCE_V6: /Users/omodei/miniconda/envs/threeML/share/fermitools/refdata/fermi/galdiffuse/gll_iem_v07.fits
Cutting the template around the ROI:
2020-04-17 15:29:17 INFO GTAnalysis.__init__():
--------------------------------------------------------------------------------
fermipy version 0.17.3
ScienceTools version unknown
Found Isotropic template for irf P8R2_SOURCE_V6: /Users/omodei/miniconda/envs/threeML/share/fermitools/refdata/fermi/galdiffuse/iso_P8R2_SOURCE_V6_v06.txt
2020-04-17 15:29:23 INFO GTAnalysis.setup(): Running setup.
2020-04-17 15:29:23 INFO GTBinnedAnalysis.setup(): Running setup for component 00
2020-04-17 15:29:23 INFO GTBinnedAnalysis.run_gtapp(): Running gtselect.
2020-04-17 15:29:23 INFO GTBinnedAnalysis.run_gtapp(): time -p gtselect infile=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/Crab_data/L2004131748578D8991BF28_PH00.fits outfile=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/ft1_00.fits ra=83.6330909547 dec=22.0144947866 rad=7.57106781187 tmin="INDEF" tmax="INDEF" emin=100.0 emax=300000.0 zmin=0.0 zmax=100.0 evclass=128 evtype=3 convtype=-1 phasemin=0.0 phasemax=1.0 evtable="EVENTS" chatter=3 clobber=yes debug=no gui=no mode="ql"
2020-04-17 15:29:24 INFO GTBinnedAnalysis.run_gtapp():
2020-04-17 15:29:24 INFO GTBinnedAnalysis.run_gtapp(): WARNING: version mismatch between CFITSIO header (v3.43) and linked library (v3.41).
2020-04-17 15:29:24 INFO GTBinnedAnalysis.run_gtapp():
2020-04-17 15:29:24 INFO GTBinnedAnalysis.run_gtapp(): This is gtselect version HEAD
2020-04-17 15:29:24 INFO GTBinnedAnalysis.run_gtapp(): Applying filter string: ((EVENT_CLASS&o200) != o0) && angsep(RA,DEC,83.63309095,22.01449479) < 7.571067812 && ((EVENT_TYPE&o3) != o0) && 100 < ENERGY && ENERGY <= 300000 && 0 < ZENITH_ANGLE && ZENITH_ANGLE <= 100 && gtifilter()
2020-04-17 15:29:24 INFO GTBinnedAnalysis.run_gtapp(): Done.
2020-04-17 15:29:24 INFO GTBinnedAnalysis.run_gtapp(): Finished gtselect. Execution time: 1.03 s
2020-04-17 15:29:24 INFO GTBinnedAnalysis.run_gtapp(): Running gtmktime.
2020-04-17 15:29:24 INFO GTBinnedAnalysis.run_gtapp(): time -p gtmktime scfile=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/Crab_data/L2004131748578D8991BF28_SC00.fits sctable="SC_DATA" filter="DATA_QUAL>0 && LAT_CONFIG==1" roicut=no evfile=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/ft1_00.fits evtable="EVENTS" outfile="/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/ft1_filtered_00.fits" apply_filter=yes overwrite=no header_obstimes=yes tstart=0.0 tstop=0.0 gtifile="default" chatter=2 clobber=yes debug=no gui=no mode="ql"
2020-04-17 15:29:24 INFO GTBinnedAnalysis.run_gtapp():
2020-04-17 15:29:24 INFO GTBinnedAnalysis.run_gtapp(): WARNING: version mismatch between CFITSIO header (v3.43) and linked library (v3.41).
2020-04-17 15:29:24 INFO GTBinnedAnalysis.run_gtapp():
2020-04-17 15:29:24 INFO GTBinnedAnalysis.run_gtapp(): Finished gtmktime. Execution time: 0.30 s
2020-04-17 15:29:24 INFO GTBinnedAnalysis.run_gtapp(): Running gtltcube.
2020-04-17 15:29:24 INFO GTBinnedAnalysis.run_gtapp(): time -p gtltcube evfile="/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/ft1_00.fits" evtable="EVENTS" scfile=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/Crab_data/L2004131748578D8991BF28_SC00.fits sctable="SC_DATA" outfile=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/ltcube_00.fits dcostheta=0.025 binsz=1.0 phibins=0 tmin=0.0 tmax=0.0 file_version="1" zmin=0.0 zmax=100.0 chatter=2 clobber=yes debug=no gui=no mode="ql"
2020-04-17 15:30:37 INFO GTBinnedAnalysis.run_gtapp():
2020-04-17 15:30:37 INFO GTBinnedAnalysis.run_gtapp(): WARNING: version mismatch between CFITSIO header (v3.43) and linked library (v3.41).
2020-04-17 15:30:37 INFO GTBinnedAnalysis.run_gtapp():
2020-04-17 15:30:37 INFO GTBinnedAnalysis.run_gtapp(): Working on file /Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/Crab_data/L2004131748578D8991BF28_SC00.fits
2020-04-17 15:30:37 INFO GTBinnedAnalysis.run_gtapp(): .....................!
2020-04-17 15:30:37 INFO GTBinnedAnalysis.run_gtapp(): Finished gtltcube. Execution time: 72.54 s
2020-04-17 15:30:38 INFO GTBinnedAnalysis.run_gtapp(): Running gtbin.
2020-04-17 15:30:38 INFO GTBinnedAnalysis.run_gtapp(): time -p gtbin evfile=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/ft1_00.fits scfile=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/Crab_data/L2004131748578D8991BF28_SC00.fits outfile=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/ccube_00.fits algorithm="ccube" ebinalg="LOG" emin=100.0 emax=300000.0 enumbins=28 ebinfile=NONE tbinalg="LIN" tbinfile=NONE nxpix=100 nypix=100 binsz=0.1 coordsys="CEL" xref=83.6330909547 yref=22.0144947866 axisrot=0.0 rafield="RA" decfield="DEC" proj="AIT" hpx_ordering_scheme="RING" hpx_order=3 hpx_ebin=yes hpx_region="" evtable="EVENTS" sctable="SC_DATA" efield="ENERGY" tfield="TIME" chatter=3 clobber=yes debug=no gui=no mode="ql"
2020-04-17 15:30:38 INFO GTBinnedAnalysis.run_gtapp():
2020-04-17 15:30:38 INFO GTBinnedAnalysis.run_gtapp(): WARNING: version mismatch between CFITSIO header (v3.43) and linked library (v3.41).
2020-04-17 15:30:38 INFO GTBinnedAnalysis.run_gtapp():
2020-04-17 15:30:38 INFO GTBinnedAnalysis.run_gtapp(): This is gtbin version HEAD
2020-04-17 15:30:38 INFO GTBinnedAnalysis.run_gtapp(): Finished gtbin. Execution time: 0.46 s
2020-04-17 15:30:38 INFO GTBinnedAnalysis.run_gtapp(): Running gtexpcube2.
2020-04-17 15:30:38 INFO GTBinnedAnalysis.run_gtapp(): time -p gtexpcube2 infile=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/ltcube_00.fits cmap=none outfile=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/bexpmap_00.fits irfs="P8R2_SOURCE_V6" evtype=3 edisp_bins=0 nxpix=360 nypix=180 binsz=1.0 coordsys="CEL" xref=0.0 yref=0.0 axisrot=0.0 proj="CAR" ebinalg="LOG" emin=100.0 emax=300000.0 enumbins=28 ebinfile="NONE" hpx_ordering_scheme="RING" hpx_order=6 bincalc="EDGE" ignorephi=no thmax=180.0 thmin=0.0 table="EXPOSURE" chatter=3 clobber=yes debug=no mode="ql"
2020-04-17 15:30:45 INFO GTBinnedAnalysis.run_gtapp():
2020-04-17 15:30:45 INFO GTBinnedAnalysis.run_gtapp(): WARNING: version mismatch between CFITSIO header (v3.43) and linked library (v3.41).
2020-04-17 15:30:45 INFO GTBinnedAnalysis.run_gtapp():
2020-04-17 15:30:45 INFO GTBinnedAnalysis.run_gtapp(): This is gtexpcube2 version HEAD
2020-04-17 15:30:45 INFO GTBinnedAnalysis.run_gtapp(): Using evtype=3 (i.e., FRONT/BACK irfs)
2020-04-17 15:30:45 INFO GTBinnedAnalysis.run_gtapp(): ResponseFunctions::load: IRF used: P8R2_SOURCE_V6
2020-04-17 15:30:45 INFO GTBinnedAnalysis.run_gtapp(): event_types: 0 1
2020-04-17 15:30:45 INFO GTBinnedAnalysis.run_gtapp(): Computing binned exposure map....................!
2020-04-17 15:30:45 INFO GTBinnedAnalysis.run_gtapp(): Finished gtexpcube2. Execution time: 7.04 s
2020-04-17 15:30:45 INFO GTBinnedAnalysis.run_gtapp(): Running gtexpcube2.
2020-04-17 15:30:45 INFO GTBinnedAnalysis.run_gtapp(): time -p gtexpcube2 infile=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/ltcube_00.fits cmap=none outfile=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/bexpmap_roi_00.fits irfs="P8R2_SOURCE_V6" evtype=3 edisp_bins=0 nxpix=100 nypix=100 binsz=0.1 coordsys="CEL" xref=83.6330909547 yref=22.0144947866 axisrot=0.0 proj="CAR" ebinalg="LOG" emin=100.0 emax=300000.0 enumbins=28 ebinfile="NONE" hpx_ordering_scheme="RING" hpx_order=6 bincalc="EDGE" ignorephi=no thmax=180.0 thmin=0.0 table="EXPOSURE" chatter=3 clobber=yes debug=no mode="ql"
2020-04-17 15:30:47 INFO GTBinnedAnalysis.run_gtapp():
2020-04-17 15:30:47 INFO GTBinnedAnalysis.run_gtapp(): WARNING: version mismatch between CFITSIO header (v3.43) and linked library (v3.41).
2020-04-17 15:30:47 INFO GTBinnedAnalysis.run_gtapp():
2020-04-17 15:30:47 INFO GTBinnedAnalysis.run_gtapp(): This is gtexpcube2 version HEAD
2020-04-17 15:30:47 INFO GTBinnedAnalysis.run_gtapp(): Using evtype=3 (i.e., FRONT/BACK irfs)
2020-04-17 15:30:47 INFO GTBinnedAnalysis.run_gtapp(): ResponseFunctions::load: IRF used: P8R2_SOURCE_V6
2020-04-17 15:30:47 INFO GTBinnedAnalysis.run_gtapp(): event_types: 0 1
2020-04-17 15:30:47 INFO GTBinnedAnalysis.run_gtapp(): Computing binned exposure map....................!
2020-04-17 15:30:47 INFO GTBinnedAnalysis.run_gtapp(): Finished gtexpcube2. Execution time: 1.50 s
2020-04-17 15:30:47 INFO GTBinnedAnalysis.run_gtapp(): Running gtsrcmaps.
2020-04-17 15:30:47 INFO GTBinnedAnalysis.run_gtapp(): time -p gtsrcmaps scfile=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/Crab_data/L2004131748578D8991BF28_SC00.fits sctable="SC_DATA" expcube=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/ltcube_00.fits cmap=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/ccube_00.fits srcmdl=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/srcmdl_00.xml bexpmap=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/bexpmap_00.fits wmap=none outfile=/Users/omodei/GRBWorkDir/MY_PYTHON_MODULES/gitrepository/threeML/docs/examples/__6b76c18fe898c8f0dceae172dc081dd6/srcmap_00.fits irfs="P8R2_SOURCE_V6" evtype=3 convol=yes resample=yes rfactor=2 minbinsz=0.05 ptsrc=yes psfcorr=yes emapbnds=no edisp_bins=0 copyall=no chatter=3 clobber=yes debug=no gui=no mode="ql"
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp():
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): WARNING: version mismatch between CFITSIO header (v3.43) and linked library (v3.41).
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp():
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): This is gtsrcmaps version HEAD
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Using evtype=3 (i.e., FRONT/BACK irfs)
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): ResponseFunctions::load: IRF used: P8R2_SOURCE_V6
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): event_types: 0 1
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named PSR_J0534p2200
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named Crab_synch
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named Crab_IC
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _2MASS_J05262938p2247232
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0544d4p2238
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named TXS_0518p211
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named TXS_0533p175
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named NVSS_J053855p164612
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0550d9p2552
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0534d2p2751
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named Sim_147
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _3C_138
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named MG2_J060149p2324
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named TXS_0506p196
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0533d9p2838
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4C_p27d15
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4C_p22d12
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named PKS_0507p17
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named NVSS_J052429p283918
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named TXS_0536p145
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named GB6_J0515p1527
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0609d0p2136
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _1RXS_J050107d1p242318
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named NVSS_J060900p203746
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0609d0p2006
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named PKS_0528p134
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _1RXS_J045931d5p192242
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0531d7p1241c
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named GB6_J0614p2426
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0616d5p2235
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named IC_443
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named MG2_J061020p1701
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _2MASS_J05400502p1209550
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0613d1p1749
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named MG2_J045613p2702
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named PSR_J0554p3107
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0618d9p2240c
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0527d1p1152
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named GB6_J0611p2803
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _2MASS_J05594100p3042282
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0620d1p2246
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0502d2p3016
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named NVSS_J060015p124344
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named PSR_J0621p2514
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0447d2p2446
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0612d6p1520c
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named PKS_0459p135
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named GB6_J0549p3258
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named GB6_J0501p3048
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0453d3p2843
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named RX_J0620d6p2644
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named NVSS_J060922p140744
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named RX_J0621d7p1747
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0532d6p3358
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0537d5p0959
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0535d3p0934
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0512d3p3329
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named GB6_J0529p0934
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named SSTSL2_J043903d89p215240d5
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named GB6_J0629p2437
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _2MASX_J05054288p3321355
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named NVSS_J060835p114942
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0620d4p1445
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named MG2_J044950p3026
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0445d7p1535
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named PKS_0506p101
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named TXS_0628p203
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named NVSS_J045804p115142
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named B2_0437p27B
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named PMN_J0517p0858
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named TXS_0516p087
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0538d9p3549c
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named NVSS_J054218p354149
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named GB6_J0623p3045
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0622d5p3120
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0540d7p3611
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _1RXS_J044127d8p150455
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0618d7p1211
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named MG2_J050234p3436
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named OG_050
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named PSR_J0633p1746
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named TXS_0437p145
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0436d9p2915
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0540d2p0655
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named PKS_0446p11
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named TXS_0429p174
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named RGB_J0640p245B
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named MG2_J043337p2905
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named B2_0619p33
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named RX_J0638d6p2751
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0534d0p3746c
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named MG1_J043022p1655
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named NVSS_J054529p061955
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named MG2_J042948p2843
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named B3_0548p378B
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named B2_0441p34
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named TXS_0420p210
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0435d5p3232
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named MS_0419d3p1943
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named NVSS_J043307p322840
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named PSR_J0605p3757
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named GB6_J0649p2228
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named TXS_0506p056
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named B2_0557p38
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named PKS_0459p060
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4C_p06d21
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _1ES_0647p250
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named NVSS_J065035p205556
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0551d7p0446
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named PSR_J0631p1036
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named GB6_J0648p1749
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0442d8p3609
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named B2_0552p39A
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _1RXS_J064814d1p160708
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named MG2_J065230p1934
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named PKS_0502p049
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named RX_J0648d7p1516
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named GB6_J0620p3806
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named PSR_J0622p3749
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named GB6_J0653p2816
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named TXS_0442p071
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named MG1_J050533p0415
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named TXS_0431p092
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named B3_0509p406
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named B2_0635p35
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named PKS_0550p032
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _2MASX_J06533986p1636432
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0658d7p2318
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named GB6_J0555p0304
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named B2_0655p27A
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0642d4p1048
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _2MASS_J05065014p0323587
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named _4FGL_J0409d2p2542
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named isodiff
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Creating source named galdiff
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for B2_0437p27B 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for B2_0441p34 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for B2_0552p39A 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for B2_0557p38 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for B2_0619p33 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for B2_0635p35 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for B2_0655p27A 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for B3_0509p406 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for B3_0548p378B 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for Crab_IC 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for Crab_synch 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for GB6_J0501p3048 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for GB6_J0515p1527 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for GB6_J0529p0934 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for GB6_J0549p3258 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for GB6_J0555p0304 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for GB6_J0611p2803 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for GB6_J0614p2426 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for GB6_J0620p3806 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for GB6_J0623p3045 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for GB6_J0629p2437 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for GB6_J0648p1749 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for GB6_J0649p2228 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for GB6_J0653p2816 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for IC_443 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for MG1_J043022p1655 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for MG1_J050533p0415 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for MG2_J042948p2843 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for MG2_J043337p2905 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for MG2_J044950p3026 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for MG2_J045613p2702 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for MG2_J050234p3436 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for MG2_J060149p2324 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for MG2_J061020p1701 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for MG2_J065230p1934 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for MS_0419d3p1943 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for NVSS_J043307p322840 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for NVSS_J045804p115142 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for NVSS_J052429p283918 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for NVSS_J053855p164612 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for NVSS_J054218p354149 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for NVSS_J054529p061955 29....................!
2020-04-17 15:31:51 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for NVSS_J060015p124344 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for NVSS_J060835p114942 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for NVSS_J060900p203746 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for NVSS_J060922p140744 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for NVSS_J065035p205556 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for OG_050 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for PKS_0446p11 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for PKS_0459p060 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for PKS_0459p135 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for PKS_0502p049 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for PKS_0506p101 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for PKS_0507p17 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for PKS_0528p134 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for PKS_0550p032 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for PMN_J0517p0858 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for PSR_J0534p2200 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for PSR_J0554p3107 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for PSR_J0605p3757 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for PSR_J0621p2514 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for PSR_J0622p3749 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for PSR_J0631p1036 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for PSR_J0633p1746 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for RGB_J0640p245B 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for RX_J0620d6p2644 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for RX_J0621d7p1747 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for RX_J0638d6p2751 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for RX_J0648d7p1516 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for SSTSL2_J043903d89p215240d5 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for Sim_147 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for TXS_0420p210 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for TXS_0429p174 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for TXS_0431p092 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for TXS_0437p145 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for TXS_0442p071 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for TXS_0506p056 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for TXS_0506p196 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for TXS_0516p087 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for TXS_0518p211 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for TXS_0533p175 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for TXS_0536p145 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for TXS_0628p203 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _1ES_0647p250 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _1RXS_J044127d8p150455 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _1RXS_J045931d5p192242 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _1RXS_J050107d1p242318 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _1RXS_J064814d1p160708 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _2MASS_J05065014p0323587 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _2MASS_J05262938p2247232 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _2MASS_J05400502p1209550 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _2MASS_J05594100p3042282 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _2MASX_J05054288p3321355 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _2MASX_J06533986p1636432 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _3C_138 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4C_p06d21 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4C_p22d12 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4C_p27d15 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0409d2p2542 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0435d5p3232 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0436d9p2915 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0442d8p3609 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0445d7p1535 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0447d2p2446 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0453d3p2843 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0502d2p3016 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0512d3p3329 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0527d1p1152 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0531d7p1241c 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0532d6p3358 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0533d9p2838 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0534d0p3746c 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0534d2p2751 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0535d3p0934 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0537d5p0959 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0538d9p3549c 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0540d2p0655 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0540d7p3611 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0544d4p2238 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0550d9p2552 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0551d7p0446 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0609d0p2006 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0609d0p2136 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0612d6p1520c 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0613d1p1749 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0616d5p2235 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0618d7p1211 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0618d9p2240c 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0620d1p2246 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0620d4p1445 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0622d5p3120 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0642d4p1048 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for _4FGL_J0658d7p2318 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for galdiff 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Generating SourceMap for isodiff 29....................!
2020-04-17 15:35:31 INFO GTBinnedAnalysis.run_gtapp(): Finished gtsrcmaps. Execution time: 284.42 s
2020-04-17 15:35:31 INFO GTBinnedAnalysis.setup(): Finished setup for component 00
2020-04-17 15:35:31 INFO GTBinnedAnalysis._create_binned_analysis(): Creating BinnedAnalysis for component 00.
2020-04-17 15:35:42 INFO GTAnalysis.setup(): Initializing source properties
2020-04-17 15:35:43 INFO GTAnalysis.setup(): Finished setup.
2020-04-17 15:35:43 INFO GTAnalysis.free_source(): Fixing parameters for PSR_J0534p2200 : ['Prefactor']
2020-04-17 15:35:43 INFO GTAnalysis.free_source(): Fixing parameters for Crab_synch : ['Prefactor']
2020-04-17 15:35:43 INFO GTAnalysis.free_source(): Fixing parameters for Crab_IC : ['Prefactor']
2020-04-17 15:35:43 INFO GTAnalysis.free_source(): Fixing parameters for _2MASS_J05262938p2247232: ['Prefactor']
2020-04-17 15:35:43 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0544d4p2238 : ['Prefactor']
2020-04-17 15:35:43 INFO GTAnalysis.free_source(): Fixing parameters for TXS_0518p211 : ['Prefactor']
2020-04-17 15:35:43 INFO GTAnalysis.free_source(): Fixing parameters for TXS_0533p175 : ['Prefactor']
2020-04-17 15:35:43 INFO GTAnalysis.free_source(): Fixing parameters for NVSS_J053855p164612 : ['Prefactor']
2020-04-17 15:35:43 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0550d9p2552 : ['Prefactor']
2020-04-17 15:35:43 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0534d2p2751 : ['Prefactor']
2020-04-17 15:35:43 INFO GTAnalysis.free_source(): Fixing parameters for Sim_147 : ['Prefactor']
2020-04-17 15:35:43 INFO GTAnalysis.free_source(): Fixing parameters for _3C_138 : ['Prefactor']
2020-04-17 15:35:43 INFO GTAnalysis.free_source(): Fixing parameters for MG2_J060149p2324 : ['Prefactor']
2020-04-17 15:35:43 INFO GTAnalysis.free_source(): Fixing parameters for TXS_0506p196 : ['Prefactor']
2020-04-17 15:35:43 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0533d9p2838 : ['Prefactor']
2020-04-17 15:35:43 INFO GTAnalysis.free_source(): Fixing parameters for _4C_p27d15 : ['Prefactor']
2020-04-17 15:35:43 INFO GTAnalysis.free_source(): Fixing parameters for _4C_p22d12 : ['Prefactor']
2020-04-17 15:35:43 INFO GTAnalysis.free_source(): Fixing parameters for PKS_0507p17 : ['Prefactor']
2020-04-17 15:35:43 INFO GTAnalysis.free_source(): Fixing parameters for NVSS_J052429p283918 : ['Prefactor']
2020-04-17 15:35:43 INFO GTAnalysis.free_source(): Fixing parameters for TXS_0536p145 : ['Prefactor']
2020-04-17 15:35:43 INFO GTAnalysis.free_source(): Fixing parameters for GB6_J0515p1527 : ['Prefactor']
2020-04-17 15:35:43 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0609d0p2136 : ['Prefactor']
2020-04-17 15:35:43 INFO GTAnalysis.free_source(): Fixing parameters for _1RXS_J050107d1p242318: ['Prefactor']
2020-04-17 15:35:43 INFO GTAnalysis.free_source(): Fixing parameters for NVSS_J060900p203746 : ['Prefactor']
2020-04-17 15:35:43 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0609d0p2006 : ['Prefactor']
2020-04-17 15:35:43 INFO GTAnalysis.free_source(): Fixing parameters for PKS_0528p134 : ['Prefactor']
2020-04-17 15:35:43 INFO GTAnalysis.free_source(): Fixing parameters for _1RXS_J045931d5p192242: ['Prefactor']
2020-04-17 15:35:43 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0531d7p1241c : ['Prefactor']
2020-04-17 15:35:43 INFO GTAnalysis.free_source(): Fixing parameters for GB6_J0614p2426 : ['Prefactor']
2020-04-17 15:35:43 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0616d5p2235 : ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for IC_443 : ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for MG2_J061020p1701 : ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for _2MASS_J05400502p1209550: ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0613d1p1749 : ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for MG2_J045613p2702 : ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for PSR_J0554p3107 : ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0618d9p2240c : ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0527d1p1152 : ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for GB6_J0611p2803 : ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for _2MASS_J05594100p3042282: ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0620d1p2246 : ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0502d2p3016 : ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for NVSS_J060015p124344 : ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for PSR_J0621p2514 : ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0447d2p2446 : ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0612d6p1520c : ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for PKS_0459p135 : ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for GB6_J0549p3258 : ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for GB6_J0501p3048 : ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0453d3p2843 : ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for RX_J0620d6p2644 : ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for NVSS_J060922p140744 : ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for RX_J0621d7p1747 : ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0532d6p3358 : ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0537d5p0959 : ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0535d3p0934 : ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0512d3p3329 : ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for GB6_J0529p0934 : ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for SSTSL2_J043903d89p215240d5: ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for GB6_J0629p2437 : ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for _2MASX_J05054288p3321355: ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for NVSS_J060835p114942 : ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0620d4p1445 : ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for MG2_J044950p3026 : ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0445d7p1535 : ['Prefactor']
2020-04-17 15:35:44 INFO GTAnalysis.free_source(): Fixing parameters for PKS_0506p101 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for TXS_0628p203 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for NVSS_J045804p115142 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for B2_0437p27B : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for PMN_J0517p0858 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for TXS_0516p087 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0538d9p3549c : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for NVSS_J054218p354149 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for GB6_J0623p3045 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0622d5p3120 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0540d7p3611 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for _1RXS_J044127d8p150455: ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0618d7p1211 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for MG2_J050234p3436 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for OG_050 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for PSR_J0633p1746 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for TXS_0437p145 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0436d9p2915 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0540d2p0655 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for PKS_0446p11 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for TXS_0429p174 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for RGB_J0640p245B : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for MG2_J043337p2905 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for B2_0619p33 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for RX_J0638d6p2751 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0534d0p3746c : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for MG1_J043022p1655 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for NVSS_J054529p061955 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for MG2_J042948p2843 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for B3_0548p378B : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for B2_0441p34 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for TXS_0420p210 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0435d5p3232 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for MS_0419d3p1943 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for NVSS_J043307p322840 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for PSR_J0605p3757 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for GB6_J0649p2228 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for TXS_0506p056 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for B2_0557p38 : ['Prefactor']
2020-04-17 15:35:45 INFO GTAnalysis.free_source(): Fixing parameters for PKS_0459p060 : ['Prefactor']
2020-04-17 15:35:46 INFO GTAnalysis.free_source(): Fixing parameters for _4C_p06d21 : ['Prefactor']
2020-04-17 15:35:46 INFO GTAnalysis.free_source(): Fixing parameters for _1ES_0647p250 : ['Prefactor']
2020-04-17 15:35:46 INFO GTAnalysis.free_source(): Fixing parameters for NVSS_J065035p205556 : ['Prefactor']
2020-04-17 15:35:46 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0551d7p0446 : ['Prefactor']
2020-04-17 15:35:46 INFO GTAnalysis.free_source(): Fixing parameters for PSR_J0631p1036 : ['Prefactor']
2020-04-17 15:35:46 INFO GTAnalysis.free_source(): Fixing parameters for GB6_J0648p1749 : ['Prefactor']
2020-04-17 15:35:46 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0442d8p3609 : ['Prefactor']
2020-04-17 15:35:46 INFO GTAnalysis.free_source(): Fixing parameters for B2_0552p39A : ['Prefactor']
2020-04-17 15:35:46 INFO GTAnalysis.free_source(): Fixing parameters for _1RXS_J064814d1p160708: ['Prefactor']
2020-04-17 15:35:46 INFO GTAnalysis.free_source(): Fixing parameters for MG2_J065230p1934 : ['Prefactor']
2020-04-17 15:35:46 INFO GTAnalysis.free_source(): Fixing parameters for PKS_0502p049 : ['Prefactor']
2020-04-17 15:35:46 INFO GTAnalysis.free_source(): Fixing parameters for RX_J0648d7p1516 : ['Prefactor']
2020-04-17 15:35:46 INFO GTAnalysis.free_source(): Fixing parameters for GB6_J0620p3806 : ['Prefactor']
2020-04-17 15:35:46 INFO GTAnalysis.free_source(): Fixing parameters for PSR_J0622p3749 : ['Prefactor']
2020-04-17 15:35:46 INFO GTAnalysis.free_source(): Fixing parameters for GB6_J0653p2816 : ['Prefactor']
2020-04-17 15:35:46 INFO GTAnalysis.free_source(): Fixing parameters for TXS_0442p071 : ['Prefactor']
2020-04-17 15:35:46 INFO GTAnalysis.free_source(): Fixing parameters for MG1_J050533p0415 : ['Prefactor']
2020-04-17 15:35:46 INFO GTAnalysis.free_source(): Fixing parameters for TXS_0431p092 : ['Prefactor']
2020-04-17 15:35:46 INFO GTAnalysis.free_source(): Fixing parameters for B3_0509p406 : ['Prefactor']
2020-04-17 15:35:46 INFO GTAnalysis.free_source(): Fixing parameters for B2_0635p35 : ['Prefactor']
2020-04-17 15:35:46 INFO GTAnalysis.free_source(): Fixing parameters for PKS_0550p032 : ['Prefactor']
2020-04-17 15:35:46 INFO GTAnalysis.free_source(): Fixing parameters for _2MASX_J06533986p1636432: ['Prefactor']
2020-04-17 15:35:46 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0658d7p2318 : ['Prefactor']
2020-04-17 15:35:46 INFO GTAnalysis.free_source(): Fixing parameters for GB6_J0555p0304 : ['Prefactor']
2020-04-17 15:35:46 INFO GTAnalysis.free_source(): Fixing parameters for B2_0655p27A : ['Prefactor']
2020-04-17 15:35:46 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0642d4p1048 : ['Prefactor']
2020-04-17 15:35:46 INFO GTAnalysis.free_source(): Fixing parameters for _2MASS_J05065014p0323587: ['Prefactor']
2020-04-17 15:35:46 INFO GTAnalysis.free_source(): Fixing parameters for _4FGL_J0409d2p2542 : ['Prefactor']
[30]:
bayes.set_sampler('emcee')
n_walkers = 10
burn_in = 10
n_samples = 500
bayes.sampler.setup(n_iterations=n_samples,n_burn_in=burn_in,n_walkers=n_walkers )
res = bayes.sample()
100%|██████████| 10/10 [00:02<00:00, 4.89it/s]
100%|██████████| 500/500 [01:27<00:00, 6.22it/s]
Mean acceptance fraction: 0.4832
Maximum a posteriori probability (MAP) point:
| result | unit | |
|---|---|---|
| parameter | ||
| PSR_J0534p2200...K | (1.49 +/- 0.05) x 10^-13 | 1 / (cm2 keV s) |
| PSR_J0534p2200...index | -1.806 +/- 0.024 | |
| _2MASS_J05262938p2247232.spectrum.main.Powerlaw.K | (5 +/- 4) x 10^-17 | 1 / (cm2 keV s) |
| _4FGL_J0544d4p2238.spectrum.main.Powerlaw.K | (3.9 +/- 3.0) x 10^-16 | 1 / (cm2 keV s) |
Values of -log(posterior) at the minimum:
| -log(posterior) | |
|---|---|
| LAT | -118069.424919 |
| total | -118069.424919 |
Values of statistical measures:
| statistical measures | |
|---|---|
| AIC | 236136.849838 |
| BIC | 236138.849838 |
| DIC | 236148.439706 |
| PDIC | 2.004265 |
[32]:
#credible_intervals = bayes.get_credible_intervals()
You can access to the parameter range like this (HPD):
[46]:
this_K = bayes.results.get_variates(
'PSR_J0534p2200.spectrum.main.Super_cutoff_powerlaw.K'
)
this_idx = bayes.results.get_variates(
'PSR_J0534p2200.spectrum.main.Super_cutoff_powerlaw.index'
)
print('Highest_posterior_density_intervals :')
print('K (68%%): %10.2e,%10.2e' % this_K.highest_posterior_density_interval(cl=.68))
print('K (95%%): %10.2e,%10.2e' % this_K.highest_posterior_density_interval(cl=.95))
print('Index (68%%): %10.2e,%10.2e' % this_idx.highest_posterior_density_interval(cl=.68))
print('Index (95%%): %10.2e,%10.2e' % this_idx.highest_posterior_density_interval(cl=.95))
Highest_posterior_density_intervals :
K (68%): 1.44e-13, 1.54e-13
K (95%): 1.40e-13, 1.60e-13
Index (68%): -1.83e+00, -1.78e+00
Index (95%): -1.85e+00, -1.76e+00
[47]:
corner_figure = bayes.results.corner_plot()
plt.show()
[ ]: