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')