Bayesian analysis methods are increasingly being used in theoretical nuclear physics [1]}. A specific example is the determination of parameters in a chiral effective field theory description of the low-energy, strong interaction between neutrons and protons (see e.g. [2]}, [3]}, [4]}, and references therein). In short, we seek to find the parameter vector \(\theta \)  that minimizes the deviation between the model and experimental data. In this specific case, the calibration data corresponds to the 3S1 phase shifts from an analysis of neutron-proton scattering cross sections [5]}). See, e.g. [1]} for more details on the definition of the likelihood function.
