Bandpasses¶
Constructing a Bandpass¶
Bandpass objects represent the transmission fraction of an astronomical filter as a function of dispersion (photon wavelength, frequency or energy). They are basically simple containers for arrays of these values, with a couple special features. To get a bandpass that is in the registry (built-in):
>>> import sncosmo
>>> band = sncosmo.get_bandpass('sdssi')
>>> band
<Bandpass 'sdssi' at 0x...>
To create a Bandpass directly, you can supply arrays of wavelength and transmission values:
>>> wavelength = [4000., 5000.]
>>> transmission = [1., 1.]
>>> sncosmo.Bandpass(wavelength, transmission, name='tophatg')
<Bandpass 'tophatg' at 0x...>
By default, the first argument is assumed to be wavelength in Angstroms.
To specify a different dispersion unit, use a unit from the
astropy.units
package:
>>> import astropy.units as u
>>> wavelength = [400., 500.]
>>> transmission = [1., 1.]
>>> Bandpass(wavelength, transmission, wave_unit=u.nm)
<Bandpass 'tophatg' at 0x...>
Using a Bandpass¶
A Bandpass acts like a continuous 1-d function, returning the transmission at supplied wavelengths (always in Angstroms):
>>> band([4100., 4250., 4300.])
array([ 0., 1., 1.])
Note that the transmission is zero outside the defined wavelength range. Linear interpolation is used between the defined wavelengths.
Bnadpasses have a few other useful properties. You can get the range of wavelengths where the transmission is non-zero:
>>> band.minwave(), band.maxwave()
(4000.0, 5000.0)
Or the transmission-weighted effective wavelength:
>>> band.wave_eff
4500.0
Or the name:
>>> band.name
'tophatg'
Adding Bandpasses to the Registry¶
You can create your own bandpasses and use them like built-ins by adding them to the registry. Suppose we want to register the ‘tophatg’ bandpass we created:
>>> sncosmo.register(band, 'tophatg')
Or if band.name
has been set:
>>> sncosmo.register(band) # registers band under band.name
After doing this, we can get the bandpass object by doing
>>> band = sncosmo.get_bandpass('tophatg')
Also, we can pass the string 'tophatg'
to any function that
takes a Bandpass
object. This means that you can create
and register bandpasses at the top of a script, then just keep track
of string identifiers throughout the rest of the script.