I'm reading this book The supersymmetric method in random matrix theory and applications to QCD. In page 302-303, the author calculate the following integral $$ Z=\int d\psi \, dH \, P(H)\exp\left(-\sum_{kl} i\phi^*_k(z-H)_{kl} \phi_l+i\chi^*_k (z+J-H)_{kl}\chi_l\right) $$ where $H$ is a $N\times N$ Hermitian matrix, $\phi_i$ and $\chi_i$ are respectively the commuting and anti-commuting (Grassmann) variables, the measurement $d\psi={-2\pi i}^{-1}\prod_{j=1}^N d\phi_j \, d\phi^*_j \, d\chi_j \, d\chi^*_j$. The distribution is $$ P(H)=e^{-\frac{N}{2} \operatorname{Tr}H^\dagger H}. $$ Integrating over the matrix $H$ results in $$ Z=\int d\psi \, \exp\left(-\frac{1}{2N} \operatorname{Str} \begin{pmatrix}\sum_j \phi^*_j\phi_j & \sum_j \chi^*_j \phi_j \\ \sum_j \chi_j \phi^*_j & \sum_j \chi^*_j \chi_j \end{pmatrix}^2-i\sum_j (\phi^*_jz \phi_j + \chi^*_j (z+J) \phi_j ) \right). $$ where the superetrace $\operatorname{Str} \begin{pmatrix}a & b \\ c & d\end{pmatrix}=\operatorname{Tr}a-\operatorname{Tr}d$. I'm not familiar with matrix integral and functional integral. How can I reach to the last equation?
I didn't realize that there's a very useful and introduction-level paper by Alexander D. Mirlin before I post my question, and this question will be closed.
