def test_from_events_works_aps(self):
pds_ev = AveragedPowerspectrum.from_events(
self.events,
segment_size=self.segment_size,
dt=self.dt,
norm="leahy",
silent=True)
assert np.allclose(self.leahy_pds.power, pds_ev.power)
def test_from_lc_with_err_works(self, norm):
lc = self.events.to_lc(self.dt)
lc._counts_err = np.sqrt(lc.counts.mean()) + np.zeros_like(lc.counts)
pds = AveragedPowerspectrum.from_lightcurve(
lc, segment_size=self.segment_size, norm=norm, silent=True,
gti=lc.gti)
pds_ev = AveragedPowerspectrum.from_events(
self.events, segment_size=self.segment_size, dt=self.dt, norm=norm, silent=True, gti=self.events.gti)
for attr in ["power", "freq", "m", "n", "nphots", "segment_size"]:
assert np.allclose(getattr(pds, attr), getattr(pds_ev, attr))