def __call__(self, wave):
        """Input wavelengths assumed to be in Angstrom."""

        if N.isscalar(wave):
            if (wave >= self.lower) & (wave <= self.upper):
                thru = 1.0
            else:
                thru = 0.0
        else:
            wave = N.asarray(wave)
            thru = N.zeros(wave.shape, dtype=N.float64)
            thru[(wave >= self.lower) & (wave <= self.upper)] = 1.0

        return thru
def __call__ ( self , wave ) : if N . isscalar ( wave ) : if ( wave >= self . lower ) & ( wave <= self . upper ) : thru = 1.0 else : thru = 0.0 else : wave = N . asarray ( wave ) thru = N . zeros ( wave . shape , dtype = N . float64 ) thru [ ( wave >= self . lower ) & ( wave <= self . upper ) ] = 1.0 return thru
Input wavelengths assumed to be in Angstrom.