Likelihood preservation and statistical reproduction of searches for new physics

Likelihoods associated with statistical fits in searches for new physics are beginning to be published by LHC experiments on HEPData. The first of these is the search for bottom-squark pair production by ATLAS. These likelihoods adhere to a specification first defined by the $\texttt{HistFactory}$ p.d.f. template. This is per-se independent of its implementation in $\texttt{ROOT}$ and it is useful to be able to run statistical analysis outside of the $\texttt{ROOT}$ and $\texttt{RooStats}$/$\texttt{RooFit}$ framework. We introduce a JSON schema that fully describes the $\texttt{HistFactory}$ statistical model and is sufficient to reproduce key results from published ATLAS analyses. Using two independent implementations of the model, one in $\texttt{ROOT}$ and one in pure Python, we reproduce the sbottom multi-$b$ limits using the published likelihoods on HEPData underscoring the implementation independence and long-term viability of the archived data.