def test_observation_cta_1dc():
ontime = 5.0 * u.hr
pointing = SkyCoord(0, 0, unit="deg", frame="galactic")
irfs = load_cta_irfs(
"$GAMMAPY_DATA/cta-1dc/caldb/data/cta/1dc/bcf/South_z20_50h/irf_file.fits"
)
t_ref = Time('2020-01-01T00:00:00')
tstart = 20 * u.hour
location = EarthLocation(lon="-70d18m58.84s",
lat="-24d41m0.34s",
height="2000m")
obs = Observation.create(
pointing,
irfs=irfs,
deadtime_fraction=0.1,
tstart=tstart,
tstop=tstart + ontime,
reference_time=t_ref,
location=location,
)
assert_skycoord_allclose(obs.pointing_radec, pointing.icrs)
assert_allclose(obs.observation_live_time_duration, 0.9 * ontime)
assert_allclose(obs.target_radec.ra, np.nan)
assert not np.isnan(obs.pointing_zen)
assert_allclose(obs.muoneff, 1)
def test_observation():
livetime = 5.0 * u.hr
pointing = SkyCoord(0, 0, unit="deg", frame="galactic")
irfs = load_cta_irfs(
"$GAMMAPY_DATA/cta-1dc/caldb/data/cta/1dc/bcf/South_z20_50h/irf_file.fits"
)
obs = Observation.create(pointing,
livetime=livetime,
irfs=irfs,
deadtime_fraction=0.1)
assert_skycoord_allclose(obs.pointing_radec, pointing.icrs)
assert_allclose(obs.observation_live_time_duration, 0.9 * livetime)
def test_data_store_observation(data_store):
"""Test DataStoreObservation class"""
obs = data_store.obs(23523)
assert_time_allclose(obs.tstart, Time(53343.92234009259, scale="tt", format="mjd"))
assert_time_allclose(obs.tstop, Time(53343.94186555556, scale="tt", format="mjd"))
c = SkyCoord(83.63333129882812, 21.51444435119629, unit="deg")
assert_skycoord_allclose(obs.pointing_radec, c)
c = SkyCoord(22.481705, 41.38979, unit="deg")
assert_skycoord_allclose(obs.pointing_altaz, c)
c = SkyCoord(83.63333129882812, 22.01444435119629, unit="deg")
assert_skycoord_allclose(obs.target_radec, c)