Monte Carlo simulations for particle and -ray emissions from a compound nucleus based on the Hauser-Feshbach statistical theory with pre-equilibrium emission are performed. The simulation yields reliable nuclear-reaction-wise energy spectra, or so-called exclusive spectra, for emitted neutrons and -rays, which are required in particle transport calculations for nuclear applications. The Monte Carlo method is applied to neutron-induced nuclear reactions on 56 Fe, and the results are compared with a traditional deterministic method. The neutron and -ray emission correlation is examined by gating on an 847 keV -ray that is produced by an inelastic scattering process. The partial -ray energy spectra for different -ray multiplicities are inferred using this Monte Carlo method. In addition, we investigate a correlation between two neutrons in the (n,2n) reaction.