docs: def unitResponse(self,band):
          """This is used internally for :ref:`pysynphot-formula-effstim`
          calculations."""
          #sum = asumr(band,nwave)
          total = band.throughput.sum()
          return 2.5*math.log10(total)