def test_clarsach_rmf(self): arf_c = ARF(self.arffile) rmf_c = RMF(self.rmffile) m_arf_c = arf_c.apply_arf(self.m, self.exposure) m_rmf_c = rmf_c.apply_rmf(m_arf_c) assert np.allclose(self.sherpa_rmf, m_rmf_c)
def test_clarsach_rmf_with_matrix(self): arf_c = ARF(self.arffile) rmf_c = RMF(self.rmffile) m_arf_c = arf_c.specresp * self.m * self.exposure m_rmf_c = numpy.dot(m_arf_c, rmf_c.get_dense_matrix()) assert np.allclose(self.sherpa_rmf, m_rmf_c)
def setup_class(cls): cls.filename = "data/PCU2.rsp" cls.rmf = RMF(cls.filename) cls.h = fits.open(cls.filename) return
def _read_chandra(self, filename): this_dir = os.path.dirname(os.path.abspath(filename)) ff = fits.open(filename) data = ff[1].data hdr = ff[1].header self.bin_lo = data['BIN_LO'] self.bin_hi = data['BIN_HI'] self.bin_unit = data.columns['BIN_LO'].unit self.counts = data['COUNTS'] self.rmf_file = this_dir + "/" + hdr['RESPFILE'] self.arf_file = this_dir + "/" + hdr['ANCRFILE'] self.rmf = RMF(self.rmf_file) self.arf = ARF(self.arf_file) if "EXPOSURE" in list(hdr.keys()): self.exposure = hdr['EXPOSURE'] # seconds else: self.exposure = 1.0 ff.close() return
def test_clarsach_rmf(self): rmf_c = RMF(self.rmffile) m_rmf_c = rmf_c.apply_rmf(self.m) assert np.allclose(self.sherpa_rmf, m_rmf_c)
def test_clarsach_rmf(self): rmf_c = RMF("data/NICER_May2014_rbn.rsp") m_rmf_c = rmf_c.apply_rmf(self.m) assert np.allclose(self.sherpa_rmf, m_rmf_c)
def test_failure_when_arrays_are_not_the_same_length(self): rmf = RMF(self.filename) with pytest.raises(ValueError): _, _, _, _ = rmf._flatten_arrays(rmf.n_grp, rmf.f_chan, rmf.n_chan[:-1], rmf.matrix)
def test_rmf_exists(self): rmf = RMF(self.filename)