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)