  Some theories of star formation suggest massive stars may only form in
clustered environments, which would create a deficit of massive stars in low
density environments. Observationally, Massey (2002) finds such a deficit in
samples of the field population in the Small and Large Magellanic Clouds, with
an IMF slope of {\Gamma} ~ 4. These IMF measurements represent some of the
largest known deviations from the standard Salpeter IMF slope of {\Gamma}=1.35.
Here, we carry out a comprehensive investigation of the mass function above 20
solar masses for the entire field population of the Small Magellanic Cloud,
based on data from the Runaways and Isolated O Type Star Spectroscopic Survey
of the SMC (RIOTS4). This is a spatially complete census of the entire field OB
star population of the SMC obtained with the IMACS multi-object spectrograph
and MIKE echelle spectrograph on the Magellan telescopes. Based on Monte-Carlo
simulations of the evolved present-day mass function, we find the slope of the
field IMF above 20 solar masses is {\Gamma}=2.3+/-0.4. We extend our IMF
measurement to lower masses using BV photometry from the OGLE II survey. We use
a statistical approach to generate a probability distribution for the mass of
each star from the OGLE photometry, and we again find {\Gamma}=2.3+/-0.6 for
stellar masses from 7-20 solar masses. The discovery and removal of ten
runaways in our RIOTS4 sample steepens the field IMF slope to
{\Gamma}=2.8+/-0.5. We discuss the possible effects of binarity and
star-formation history on our results, and conclude that the steep field
massive star IMF is most likely a real effect.
