Loads a top hat filter given wavelength min and max values

        Parameters
        ----------
        wave_min: astropy.units.quantity (optional)
            The minimum wavelength to use
        wave_max: astropy.units.quantity (optional)
            The maximum wavelength to use
        n_pixels: int
            The number of pixels for the filter

def load_TopHat(self, wave_min, wave_max, pixels_per_bin=100):
        """
        Loads a top hat filter given wavelength min and max values

        Parameters
        ----------
        wave_min: astropy.units.quantity (optional)
            The minimum wavelength to use
        wave_max: astropy.units.quantity (optional)
            The maximum wavelength to use
        n_pixels: int
            The number of pixels for the filter
        """
        # Get min, max, effective wavelengths and width
        self.pixels_per_bin = pixels_per_bin
        self.n_bins = 1
        self._wave_units = q.AA
        wave_min = wave_min.to(self.wave_units)
        wave_max = wave_max.to(self.wave_units)

        # Create the RSR curve
        self._wave = np.linspace(wave_min, wave_max, pixels_per_bin)
        self._throughput = np.ones_like(self.wave)
        self.raw = np.array([self.wave.value, self.throughput])

        # Calculate the effective wavelength
        wave_eff = ((wave_min + wave_max) / 2.).value
        width = (wave_max - wave_min).value

        # Add the attributes
        self.path = ''
        self.refs = ''
        self.Band = 'Top Hat'
        self.CalibrationReference = ''
        self.FWHM = width
        self.Facility = '-'
        self.FilterProfileService = '-'
        self.MagSys = '-'
        self.PhotCalID = ''
        self.PhotSystem = ''
        self.ProfileReference = ''
        self.WavelengthMin = wave_min.value
        self.WavelengthMax = wave_max.value
        self.WavelengthCen = wave_eff
        self.WavelengthEff = wave_eff
        self.WavelengthMean = wave_eff
        self.WavelengthPeak = wave_eff
        self.WavelengthPhot = wave_eff
        self.WavelengthPivot = wave_eff
        self.WavelengthUCD = ''
        self.WidthEff = width
        self.ZeroPoint = 0
        self.ZeroPointType = ''
        self.ZeroPointUnit = 'Jy'
        self.filterID = 'Top Hat'