This question originates from the definition of the Cox point process, but I suspect it might be a more general one.
If we define $$Q(\cdot) = \int_{\mathcal M} P_{\Lambda}(\cdot)Q_{\Psi}(d\Lambda)$$
Then $$\int_{\mathcal N} \mu(B) Q(d\mu) \stackrel{(\ast)}= \int_{\mathcal M} \int_{\mathcal N} \mu(B) P_{\Lambda}(d\mu) Q_{\Psi}(d\Lambda)$$
Where
$\mathcal M$ is a set of locally finite measures
$\mathcal N$ is a set of locally finite integer-valued measures
$P_\Lambda$ is the distribution of a Poisson process with intensity measure $\Lambda$
$\Psi$ is a random (diffusion) measure with distribution $Q_\Psi$
$B$ is a Borel set on the measurable space $X$ on which the measures in $\mathcal N$ are defined.
My question is: How to explain the equality $(\ast)$? Intuitivelly it makes sense. Possibly this could be contrasted with integration w.r.t. $\nu(E) = \int_E f\; d\mu$ which gives $\int_E g \; d\nu = \int_E fg \; d\mu$, if such contrast is helpful in answering the question.
Thank you.