  We study parton energy-momentum exchange with the quark gluon plasma (QGP)
within a multistage approach composed of in-medium DGLAP evolution at high
virtuality, and (linearized) Boltzmann Transport formalism at lower virtuality.
This multistage simulation is then calibrated in comparison with high $p_T$
charged hadrons, D-mesons, and the inclusive jet nuclear modification factors,
using Bayesian model-to-data comparison, to extract the virtuality-dependent
transverse momentum broadening transport coefficient $\hat{q}$. To facilitate
this undertaking, we develop a quantitative metric for validating the Bayesian
workflow, which is used to analyze the sensitivity of various model parameters
to individual observables. The usefulness of this new metric in improving
Bayesian model emulation is shown to be highly beneficial for future such
analyses.
