We present MCSED, a new spectral energy distribution (SED)-fitting code, which mates flexible stellar evolution calculations with the Markov Chain Monte Carlo algorithms of the software package emcee. MCSED takes broad, intermediate, and narrow-band photometry, emission-line fluxes, and/or absorption line spectral indices, and returns probability distributions and co-variance plots for all model parameters. MCSED includes a variety of dust attenuation curves with parameters for varying the UV slopes and bump strengths, a prescription for continuum and PAH emission from dust, models for continuum and line emission from ionized gas, options for fixed and variable stellar metallicity, and a selection of star formation rate (SFR) histories. The code is well-suited for exploring parameter inter-dependencies in sets of galaxies with known redshifts, for which there is multi-band photometry and/or spectroscopy. We apply MCSED to a sample of $\sim2000$ $1.90<z<2.35$ galaxies in the five CANDELS fields, which were selected via their strong [O III] $\lambda5007$ emission, and explore the systematic behavior of their SEDs. We find the galaxies become redder with stellar mass, due to both increasing internal attenuation and a greater population of older stars. The slope of the UV extinction curve also changes with stellar mass, and at least some galaxies exhibit an extinction excess at 2175 Angstroms. Finally, we demonstrate that below $M\lesssim10^9\,M_{\odot}$), the shape of the star-forming galaxy main sequence is highly dependent on the galaxies' assumed SFR history, as calculations which assume a constant SFR produce stellar masses that are $\sim1$ dex smaller than those found using more realistic SFR histories.

Euclid, WFIRST, and HETDEX will make emission-line selected galaxies the largest observed constituent in the $z > 1$ universe. However, we only have a limited understanding of the physical properties of galaxies selected via their Ly$\alpha$ or rest-frame optical emission lines. To begin addressing this problem, we present the basic properties of $\sim 2,000$ AEGIS, COSMOS, GOODS-N, GOODS-S, and UDS galaxies identified in the redshift range $1.90 < z < 2.35$ via their [O II], H$\beta$, and [O III] emission lines. For these $z \sim 2$ galaxies, [O III] is generally much brighter than [O II] and H$\beta$, with typical rest-frame equivalent widths of several hundred Angstroms. Moreover, these strong emission-line systems span an extremely wide range of stellar mass ($\sim 3$ dex), star-formation rate ($\sim 2$ dex), and [O III] luminosity ($\sim 2$ dex). Comparing the distributions of these properties to those of continuum selected galaxies, we find that emission-line galaxies have systematically lower stellar masses and lower optical/UV dust attenuations. These measurements lay the groundwork for an extensive comparison between these rest-frame optical emission-line galaxies and Ly$\alpha$ emitters identified in the HETDEX survey.