def setUp(self):
self.sp = S.FlatSpectrum(10, fluxunits='flam')
self.wave = S.Waveset(1000, 11000, 1000)
self.ref = S.ArraySpectrum(wave=self.wave,
flux=self.sp.flux[0] *
N.ones(self.wave.shape),
fluxunits=self.sp.fluxunits)
def setUp(self):
self.inwave = S.Waveset(1300, 1800)
self.influx = N.random.lognormal(size=len(self.inwave)) * 1e-15
self.sp = S.ArraySpectrum(wave=self.inwave,
flux=self.influx,
waveunits='nm',
fluxunits='fnu',
name='Tab2 spectrum')
def test_add_spectrum(self):
# is not allowed
wv = S.Waveset(1000, 10000, 1)
other = S.ArraySpectrum(wave=wv,
flux=np.zeros(wv.shape) + 10,
waveunits='angstroms',
fluxunits='counts')
self.assertRaises(NotImplementedError, self.obs.__add__, other)
def setUp(self):
#This case is deliberately constructed to have a waveset
#that partially overlaps, but for which all the flux is
#fully contained.
self.bp = S.ObsBandpass('acs,hrc,f555w')
w = wave = S.Waveset(1000, 10000)
f = N.zeros(w.shape)
f[4000:4010] = 1.0
self.sp = S.ArraySpectrum(wave=w, flux=f, fluxunits='flam')
def setUp(self):
self.inwave = S.Waveset(1300, 1800)
self.influx = -2.5 * N.log10(self.inwave**2)
self.sp = S.ArraySpectrum(wave=self.inwave,
flux=self.influx,
fluxunits='abmag')
def setUp(self):
self.wave = S.Waveset(3000, 5000, 10)
def setUp(self):
sp = S.FlatSpectrum(1, fluxunits='counts')
bp = S.Box(5000, 100)
binset = S.Waveset(1000, 10000, 1)
self.obs = S.Observation(sp, bp, binset=binset)
self.obs.convert('counts')
