As successful as the Standard Model of particle physics has been it still has several major shortcomings which range from unanswered theoretical questions to a lack of any explanation for observed phenomena such as dark matter. One proposed theory for physics beyond the Standard Model which provides solutions for some of these issues is supersymmetry. This dissertation presents a search for supersymmetry using 2.3 fb^(-1) of proton-proton collision data. This data was collected at a center-of-mass energy of 13 TeV by the CMS detector at the LHC during 2015. This search focuses on top squark pair production where the produced stops both decay to an all hadronic final state. These decays are characterized by multiple jets and missing transverse momentum. A baseline search region is defined to be sensitive to signal processes which occur at rates many orders of magnitude lower than Standard Model processes. The sensitivity to various signal models is improved by dividing this baseline region into distinct categories. Events with an unreconstructed lepton from leptonic W boson decays constitute the primary background. There are also significant contributions from events where a Z boson decays invisibly to neutrinos especially in bins with higher missing transverse momentum. Events with multijet production where one jet has been severely mis-measured as well as those with a pair of top quarks and an invisibly decaying Z boson also have a small presence in the search region. The contributions from these processes to the search region is estimated using data control regions. No statistically significant deviations from the predicted background yields are observed. The results are interpreted in terms of exclusion limits using the Simplified Model Spectrum framework. Stop and neutralino masses are probed up to 780 GeV and 260 GeV respectively.