code: def throughput(self, points):
          
          # Test shape
          if not points.shape == self.wave.shape:
              raise ValueError("Throughput and wavelength must be same shape.")
  
          self._throughput = points
docs: A setter for the throughput
  
          Parameters
          ----------
          throughput: sequence
              The array of throughput points