Abstract: The use of a detailed Monte Carlo simulation code for X-ray fluorescence spectrometers is demonstrated for calculating the outcome of X-ray scattering experiments in the incident energy range 40-80 keV. The code was validated by comparisons of experimental and simulated spectral distributions in the case of thick, homogeneous samples in which multiple photon scattering occurs with high probability. The experimental spectral distributions were collected at beamline BW5 of HASYLAB, Germany, where a highly energetic, monochromatic synchrotron beam is available. With respect to heterogeneous samples, the code was employed to evaluate the use of Rayleigh and Compton scatter signals for obtaining three dimensional information on the sample dark matrix composition.