Recent Fermi-LAT observations of dwarf spheroidal galaxies in the Milky Way have placed strong limits on the gamma-ray flux from dark matter annihilation. In order to produce the strongest limit on the dark matter annihilation cross-section, the observations of each dwarf galaxy have typically been "stacked" in a joint-likelihood analysis, utilizing optical observations to constrain the dark matter density profile in each dwarf. These limits have typically been computed only for singular annihilation final states, such as $b\bar{b}$ or $\tau^+\tau^-$. In this paper, we generalize this approach by producing an independent joint-likelihood analysis to set constraints on models where the dark matter particle annihilates to multiple final state fermions. We interpret these results in the context of the most popular simplified models, including those with s- and t-channel dark matter annihilation through scalar and vector mediators. We present our results as constraints on the minimum dark matter mass and the mediator sector parameters. Additionally, we compare our simplified model results to those of Effective Field Theory contact interactions in the high-mass limit.

Using limits on photon flux from Dwarf Spheroidal galaxies, we place bounds on the parameter space of models in which Dark Matter annihilates into multiple final state particle pair channels. We derive constraints on effective operator models with Dark Matter couplings to third generation fermions and to pairs of Standard Model vector bosons. We present limits in various slices of model parameter space along with estimations of the region of maximal validity of the effective operator approach for indirect detection. We visualize our bounds for models with multiple final state annihilations by projecting parameter space constraints onto triangles, a technique familiar from collider physics; and we compare our bounds to collider limits on equivalent models.