from __future__ import absolute_import, division, print_function
import os
import numpy as np
from numpy.testing import assert_allclose
from astropy.tests.helper import pytest
from astropy.table import Table
from fermipy.tests.utils import requires_dependency, requires_st_version
from fermipy import spectrum
try:
from fermipy import gtanalysis
except ImportError:
pass
# Skip tests in this file if Fermi ST aren't available
pytestmark = requires_dependency('Fermi ST')
@pytest.fixture(scope='module')
def create_draco_analysis(request, tmpdir_factory):
path = tmpdir_factory.mktemp('draco')
url = 'https://raw.githubusercontent.com/fermiPy/fermipy-extras/master/data/fermipy_test_draco.tar.gz'
outfile = path.join('fermipy_test_draco.tar.gz')
dirname = path.join()
os.system('curl -o %s -OL %s' % (outfile, url))
os.system('cd %s;tar xzf %s' % (dirname, outfile))
request.addfinalizer(lambda: path.remove(rec=1))
cfgfile = path.join('fermipy_test_draco', 'config.yaml')
gta = gtanalysis.GTAnalysis(str(cfgfile))
gta.setup()
from __future__ import absolute_import, division, print_function
import os
import numpy as np
from numpy.testing import assert_allclose
from astropy.tests.helper import pytest
from astropy.coordinates import SkyCoord
from fermipy.tests.utils import requires_dependency
from fermipy import spectrum
try:
from fermipy import irfs
except ImportError:
pass
# Skip tests in this file if Fermi ST aren't available
pytestmark = requires_dependency('Fermi ST')
def test_calc_exp():
ltc = irfs.LTCube.create_from_obs_time(3.1536E8)
c = SkyCoord(10.0, 10.0, unit='deg')
egy = 10**np.linspace(1.0, 6.0, 6)
cth_bins = np.array([0.2, 0.6, 1.0])
exp = irfs.calc_exp(c, ltc, 'P8R2_SOURCE_V6', ['FRONT', 'BACK'],
egy, cth_bins)
assert_allclose(exp,
np.array([[1.99488297e+07, 1.19873696e+08],
[2.62408989e+10, 1.45160418e+11],