docs: def unitResponse(self,band):
          """This is used internally for :ref:`pysynphot-formula-effstim`
          calculations."""
          wave=band.wave
          total = band.trapezoidIntegration(wave,band.throughput/wave)
          modtot = total * (1.0e-29/H)
          return 1.0/modtot