This paper presents a comparative study using the first-order formula (also called “sandwich”) and the multivariate sampling methodologies for propagating uncertainty in the resolved resonance region. A distinctive aspect of this work is the generation of random cross sections by sampling Resonance Parameters (RPs) taking in consideration their correlations provided in nuclear data libraries via the Resonance Parameter Covariance Matrix (RPCM). SCOOBY (Sampling COvariance OBservatorY), a newly developed sampling tool, is presented in this paper. This tool relies on the GAIA-2 nuclear data processing code to read and correct the RPCM and generate random cross sections. The study compares the sandwich method that relies on sensitivity coefficients and covariance matrices, with the sampling method, which involves numerous Monte Carlo simulations with random cross sections. The comparison is tested on an ICSBEP benchmark, PU-MET-MIXED-002, chosen for its sensitivity to 239Pu cross sections. The results indicate that both methods quantify similar uncertainties, confirming the reliability of both the SCOOBY module and the GAIA-2 code. By using the PU-MET-MIXED-002 benchmark and sampling 239Pu resonance parameters, this work demonstrates the effectiveness of these methodologies in estimating uncertainty in criticality safety calculations. The findings highlight the robustness of these approaches in uncertainty propagation, suggesting the need for further research on additional benchmarks and expanded uncertainty propagation studies.