In order to determine the observable signatures of modified gravity theories, it is important to consider the effect of baryonic physics.We used a modified version of the ISIS code to run cosmological hydrodynamic simulations in order to study degeneracies between modified gravity and radiative hydrodynamic processes. One of the simulations was the standard ΛCDM model and four were variations of the Symmetron model. For each model we ran three variations of baryonic processes: nonradiative hydrodynamics; cooling and star formation; and cooling, star formation, and supernova feedback. We constructed stacked gas density, temperature, and dark matter density profiles of the halos in the simulations, and studied 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.