docs: def ElectronIC(pars, data):
      """
      Define particle distribution model, radiative model, and return model flux
      at data energy values
      """
  
      ECPL = ExponentialCutoffPowerLaw(
          pars[0] / u.eV, 10.0 * u.TeV, pars[1], 10 ** pars[2] * u.TeV
      )
      IC = InverseCompton(ECPL, seed_photon_fields=["CMB"])
  
      return IC.flux(data, distance=1.0 * u.kpc)