docs: def unitResponse(self,band):
          """This is used internally for :ref:`pysynphot-formula-effstim`
          calculations."""
          #sumfilt(wave,1,band)
          # SUMFILT = Sum [ FILT(I) * WAVE(I) ** NPOW * DWAVE(I) ]
          wave=band.wave
          total = band.trapezoidIntegration(wave,band.throughput*wave)
          modtot = total/(H*C)
          return 2.5*math.log10(modtot) + STZERO