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)