We present a robust and highly efficient Shannon wavelet pricing method for plain-vanilla foreign exchange European options under the jump-extended Heston model with multi-factor CIR interest rate dynamics. Under a Monte Carlo and partial differential equation hybrid computational framework, the option price can be expressed as an expectation, conditional on the variance factor, of a convolution product that involves the densities of the time-integrated domestic and foreign multi-factor CIR interest rate processes. We propose an efficient treatment to this convolution product that effectively results in a significant dimension reduction, from two multi-factor interest rate processes to only a single-factor process. By means of a state-of-the-art Shannon wavelet inverse Fourier technique, the resulting convolution product is approximated analytically and the conditional expectation can be computed very efficiently. We develop sharp approximation error bounds for the option price and hedging parameters. Numerical experiments confirm the robustness and efficiency of the method
