Let $W_1, W_2 \in \mathbb{R}^{n \times n}$, $x, b_1, b_2 \in \mathbb{R}^{n}$ and $\sigma \in \mathbb{R}^{+}$, I am trying to integrate a Gaussian over the intersection of two half-space: $$ f(x) = \frac{1}{\sqrt{2 \pi} \sigma} \int_{\Omega_1 \cap \Omega_2} \exp\left(-\frac{\left\Vert x - t \right\Vert_2^2}{2 \sigma^2}\right) d t $$ where $$ \Omega_1 = \left\{ z \in \mathbb{R}^{n}: W_1 z + b_1 \geq 0 \right\} $$ and $$ \Omega_2 = \left\{ z \in \mathbb{R}^{n}: W_2 z + b_2 \geq 0 \right\} $$
Related questions have been answered here and here but I still haven't figure out a proper way to compute this integral. I also wondering if the result from this paper can be of help.
Thank you for your help.