This article details the derivation of marginal likelihood for stochastic field amplitude, using Bessel and Gamma functions to analytically integrate random variables. The resulting normalized likelihood is split into three components: two for sum/difference peaks and one for the Compton peak, with the forms of the first two being equivalent.

(1) Dorian W. P. Amaral, Department of Physics and Astronomy, Rice University and These authors contributed approximately equally to this work;

(2) Mudit Jain, Department of Physics and Astronomy, Rice University, Theoretical Particle Physics and Cosmology, King’s College London and These authors contributed approximately equally to this work;

(3) Mustafa A. Amin, Department of Physics and Astronomy, Rice University;

(4) Christopher Tunnell, Department of Physics and Astronomy, Rice University.

B Derivation of Marginal Likelihood with Stochastic Field Amplitude

The full signal in time space is given by

Using the series representation of the Bessel function, together with Gamma function identities, the 5 random variables can be integrated out analytically. We arrive at the following marginalized (and normalized) likelihood:

where

This likelihood can be split into three individual likelihoods for the sum/difference peaks and the Compton peak, as given in Eq. (3.2). The form of the likelihoods for the sum and difference peaks is equivalent.