Tutorial1Part2
Contents
|
Working with Bandpasses
Create a bandpass from an obsmode string
ÂÂ >>> bp1=S.ObsBandpass('acs,hrc,f555w')ÂÂ >>> bp1.showfiles()
Read bandpasses
from a file
ÂÂ >>> ota=S.FileBandpass('/grp/hst/cdbs/comp/ota/hst_ota_007_syn.fits')ÂÂ >>> ccd=S.FileBandpass('/grp/hst/cdbs/comp/acs/acs_hrc_ccd_013_syn.fits')
Plot and compare several bandpasses
ÂÂ >>> P.plot(bp1.wave,bp1.throughput)ÂÂ >>> P.plot(ota.wave,ota.throughput) ÂÂ >>> P.plot(ccd.wave,ccd.throughput)ÂÂ >>> P.legend(['acs,hrc,f555w','ota','ccd'])ÂÂ >>> P.savefig('tut1_bp.ps')ÂÂ >>> bp2=ota*ccdÂÂ >>> P.plot(bp2.wave,bp2.throughput)
Exercise
Create a bandpass from another
obsmode, and plot it.
Back to Tutorial1Part1 * On to Tutorial1Part3