  Breakup reactions are one of the favored probes to study loosely bound
nuclei, particularly those in the limit of stability forming a halo. In order
to interpret such breakup experiments, the continuum discretized coupled
channel method is typically used. In this study, the first Bayesian analysis of
a breakup reaction model is performed. We use a combination of statistical
methods together with a three-body reaction model (the continuum discretized
coupled channel method) to quantify the uncertainties on the breakup
observables due to the parameters in the effective potential describing the
loosely bound projectile of interest. The combination of tools we develop opens
the path for a Bayesian analysis of not only breakup processes, but also a wide
array of complex processes that require computationally intensive reaction
models.
