In order to determine the observable signatures of modified gravity theories, it is important to consider the effect of baryonic physics. We use a modified version of the ISIS code to run cosmological hydrodynamic simulations to study degeneracies between modified gravity and radiative hydrodynamical processes. Of these, one was the standard $\Lambda$CDM model and four were variations of the Symmetron model. For each model we ran three variations of baryonic processes: non-radiative hydrodynamics; cooling and star formation; and cooling, star formation, and supernova feedback. We construct stacked gas density, temperature, and dark matter density profiles of the halos in the simulations, and study the differences between them. We find that both radiative variations of the models show degeneracies between their processes and at least two of the three parameters defining the Symmetron model.