docs: 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'