def setparms(self):
self.sp = parse_spec(
'rn(pl(4000.0,-1.5,flam),box(1850.0,1.0),1.00E-14,flam)')
self.bp = ObsBandpass('cos,nuv,g185m,c2010')
self.ref_rate = 71.8377
self.cmd = 'SpecSourcerateSpec'
self.fname = 'specAV95.fits'
def setparms(self):
self.sp = parse_spec(
'rn(icat(k93models,25400,0.0,3.9),box(2000.0,1.0),2.00E-13,flam)')
self.bp = ObsBandpass('cos,nuv,g225m,c2186')
self.ref_rate = 1460.12
self.cmd = 'SpecSourcerateSpec'
self.fname = 'specAV105.fits'
def setparms(self):
self.sp = parse_spec(
'rn(icat(k93models,44500,0.0,5.0),box(2000.0,1.0),2.00E-13,flam)')
self.bp = ObsBandpass('cos,nuv,g185m,c1835')
self.ref_rate = 1559.37
self.cmd = 'SpecSourcerateSpec'
self.fname = 'specAV6.fits'
def setparms(self):
self.sp = parse_spec(
'rn(z(spec(qso_template.fits),0.2),box(1850.0,1.0),1.00E-14,flam)')
self.bp = ObsBandpass('cos,nuv,g185m,c2010')
self.ref_rate = 47.2336
self.cmd = 'SpecSourcerateSpec'
self.fname = 'specAV93.fits'
def setparms(self):
self.sp = parse_spec(
'rn(icat(k93models,44500,0.0,5.0)*ebmvx(0.1,gal1),box(1850.0,1.0),2.00E-13,flam)'
)
self.bp = ObsBandpass('cos,nuv,g185m,c1786')
self.ref_rate = 1180.17
self.cmd = 'SpecSourcerateSpec'
self.fname = 'specAV98.fits'
def setparms(self):
self.sp = parse_spec(
'spec(earthshine.fits)*0.5+rn(spec(Zodi.fits),band(johnson,v),22.7,vegamag)+(spec(el1215a.fits)+spec(el1302a.fits)+spec(el1356a.fits)+spec(el2471a.fits))*0.5'
)
self.bp = ObsBandpass('cos,nuv,g185m,c1850')
self.ref_rate = 6.7e-05
self.cmd = 'SpecSourcerateSpec'
self.fname = 'specAV7.fits'
self.accuracy = 1e-2
def setUp(self):
tmg = locations._refTable(os.path.join('mtab', 'u921351jm_tmg.fits'))
tmc = locations._refTable(os.path.join('mtab', 'ub31649mm_tmc.fits'))
spec = FlatSpectrum(1, fluxunits='flam')
bandpass = ObsBandpass('acs,hrc,f555w', graphtable=tmg, comptable=tmc)
self.obs = Observation(spec, bandpass)
self.obs.initbinflux()
def setparms(self):
self.sp = parse_spec(
'spec(earthshine.fits)*2.0+rn(spec(Zodi.fits),band(johnson,v),22.7,vegamag)+(spec(el1215a.fits)+spec(el1302a.fits)+spec(el1356a.fits)+spec(el2471a.fits))*0.5'
)
self.bp = ObsBandpass('cos,nuv,g225m,c2186')
self.ref_rate = 0.000906
self.cmd = 'SpecSourcerateSpec'
self.fname = 'specAV109.fits'
self.accuracy = 1e-2
def test_photbw_stis_02x29_mirror_fuvmama(): mode = 'stis,0.2x29,mirror,fuvmama' # from Synphot ref_photbw = 134.79 band = ObsBandpass(mode) test_photbw = band.photbw() np.testing.assert_allclose(ref_photbw, test_photbw, rtol=0.001)
def get_filter_info(name, vega=vega):
filter = ObsBandpass(name)
vega_obs = obs.Observation(vega, filter, binset=filter.wave, force='taper')
vega_flux = vega_obs.binflux.sum()
vega_mag = 0.03
filter.flux0 = vega_flux
filter.mag0 = vega_mag
return filter
def test_photbw_acs_sbc_f122m(): mode = 'acs,sbc,f122m' # from Synphot ref_photbw = 91.063 band = ObsBandpass(mode) test_photbw = band.photbw() np.testing.assert_allclose(ref_photbw, test_photbw, rtol=0.001)
def test_photbw_wfc3_uvis2_f336w(): mode = 'wfc3,uvis2,f336w' # from Synphot ref_photbw = 158.36 band = ObsBandpass(mode) test_photbw = band.photbw() np.testing.assert_allclose(ref_photbw, test_photbw, rtol=0.001)
def test_photbw_acs_hrc_f555w(): mode = 'acs,hrc,f555w' # from Synphot ref_photbw = 357.17 band = ObsBandpass(mode) test_photbw = band.photbw() np.testing.assert_allclose(ref_photbw, test_photbw, rtol=0.001)
def test_photbw_wfc3_ir_f164n(): mode = 'wfc3,ir,f164n' # from Synphot ref_photbw = 700.05 band = ObsBandpass(mode) test_photbw = band.photbw() np.testing.assert_allclose(ref_photbw, test_photbw, rtol=0.001)
def test_photbw_nicmos_1_f090m_dn(): mode = 'nicmos,1,f090m,dn' # from Synphot ref_photbw = 559.59 band = ObsBandpass(mode) test_photbw = band.photbw() np.testing.assert_allclose(ref_photbw, test_photbw, rtol=0.001)
def test_photbw_cos_boa_nuv_mirrora(): mode = 'cos,boa,nuv,mirrora' # from Synphot ref_photbw = 370.65 band = ObsBandpass(mode) test_photbw = band.photbw() np.testing.assert_allclose(ref_photbw, test_photbw, rtol=0.001)
def test_photbw_acs_wfc1_f775w_pol_v(): mode = 'acs,wfc1,f775w,pol_v' # from Synphot ref_photbw = 444.05 band = ObsBandpass(mode) test_photbw = band.photbw() np.testing.assert_allclose(ref_photbw, test_photbw, rtol=0.001)
def get_filter_info(name, vega=vega):
if name.startswith('nirc2'):
tmp = name.split(',')
filterName = tmp[-1]
filter = nirc2syn.filters[filterName]
flux0 = nirc2syn.filter_flux0[filterName]
mag0 = nirc2syn.filter_mag0[filterName]
else:
filter = ObsBandpass(name)
vega_obs = obs.Observation(vega, filter, binset=filter.wave, force='taper')
vega_flux = vega_obs.binflux.sum()
vega_mag = 0.03
filter.flux0 = vega_flux
filter.mag0 = vega_mag
return filter
def ABtoVega(instrument, bandpass):
bp = ObsBandpass(str(instrument) + ',wfc1,'\
+ str(bandpass) + ',mjd#57754')
spec_bb = BlackBody(10000)
spec_bb_norm = spec_bb.renorm(1, 'counts', bp)
obs = observation.Observation(spec_bb_norm, bp)
# Get photometric calibration information.
photflam = obs.effstim('flam')
photplam = bp.pivot()
zp_vega = obs.effstim('vegamag')
zp_st = obs.effstim('stmag')
zp_ab = obs.effstim('abmag')
difference = zp_vega - zp_ab
return difference
def setparms(self):
self.sp = parse_spec('rn(bb(40000.0),box(1850.0,1.0),1.00E-14,flam)')
self.bp = ObsBandpass('cos,nuv,g185m,c2010')
self.ref_rate = 64.6011
self.cmd = 'SpecSourcerateSpec'
self.fname = 'specAV97.fits'
