code: def _bindex_blackbody(self, Teff, photon_weighted=False):
          
  
          if photon_weighted:
              num   = lambda w: w*self._planck(w, Teff)*self.ptf(w)*self._planck_spi(w, Teff)
              denom = lambda w: w*self._planck(w, Teff)*self.ptf(w)
              return integrate.quad(num, self.wl[0], self.wl[-1], epsabs=1e10, epsrel=1e-8)[0]/integrate.quad(denom, self.wl[0], self.wl[-1], epsabs=1e10, epsrel=1e-6)[0]
          else:
              num   = lambda w: self._planck(w, Teff)*self.ptf(w)*self._planck_spi(w, Teff)
              denom = lambda w: self._planck(w, Teff)*self.ptf(w)
              return integrate.quad(num, self.wl[0], self.wl[-1], epsabs=1e10, epsrel=1e-8)[0]/integrate.quad(denom, self.wl[0], self.wl[-1], epsabs=1e10, epsrel=1e-6)[0]
docs: Computes the mean boosting index using blackbody atmosphere:
  
          B_pb^E = \int_\lambda I(\lambda) P(\lambda) B(\lambda) d\lambda / \int_\lambda I(\lambda) P(\lambda) d\lambda
          B_pb^P = \int_\lambda \lambda I(\lambda) P(\lambda) B(\lambda) d\lambda / \int_\lambda \lambda I(\lambda) P(\lambda) d\lambda
  
          Superscripts E and P stand for energy and photon, respectively.
  
          @Teff: effective temperature in K
          @photon_weighted: photon/energy switch
  
          Returns: mean boosting index using blackbody atmosphere.