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