code: def compute_pblum_scale(self, dataset, pblum, **kwargs):
          
          logger.debug("{}.compute_pblum_scale(dataset={}, pblum={})".format(self.component, dataset, pblum))
  
          abs_luminosity = self.compute_luminosity(dataset, **kwargs)
  
          # We now want to remember the scale for all intensities such that the
          # luminosity in relative units gives the provided pblum
          pblum_scale = pblum / abs_luminosity
  
          self.set_pblum_scale(dataset, pblum_scale)
docs: intensities should already be computed for this dataset at the time for which pblum is being provided
  
          TODO: add documentation