Observations made in a JET experiment aimed at accelerating deuterons to the MeV range by third harmonic radio-frequency (RF) heating coupled into a deuterium beam are reported. Measurements are based on a set of advanced neutron and gamma-ray spectrometers that, for the first time, observe the plasma simultaneously along vertical and oblique lines of sight. Parameters of the fast ion energy distribution, such as the high energy cut-off of the deuteron distribution function and the RF coupling constant, are determined from data within a uniform analysis framework for neutron and gamma-ray spectroscopy based on a one-dimensional model and by a consistency check among the individual measurement techniques. A systematic difference is seen between the two lines of sight and is interpreted to originate from the sensitivity of the oblique detectors to the pitch-angle structure of the distribution around the resonance, which is not correctly portrayed within the adopted one dimensional model. A framework to calculate neutron and gamma-ray emission from a spatially resolved, two-dimensional deuteron distribution specified by energy/pitch is thus developed and used for a first comparison with predictions from ab initio models of RF heating at multiple harmonics. The results presented in this paper are of relevance for the development of advanced diagnostic techniques for MeV range ions in high performance fusion plasmas, with applications to the experimental validation of RF heating codes and, more generally, to studies of the energy distribution of ions in the MeV range in high performance deuterium and deuterium-tritium plasmas.