We perform in-depth dynamical modelling of the luminous and dark matter (DM) content of the elliptical galaxy NGC 1407. Our strategy consists of solving the spherical Jeans equations for three independent dynamical tracers: stars, blue GCs and red GCs in a self-consistent manner. We adopt a maximum-likelihood Markov-Chain Monte Carlo fitting technique in the attempt to constrain the inner slope of the DM density profile (the cusp/core problem), and the stellar initial mass function (IMF) of the galaxy. We find the inner logarithmic slope of the DM density profiles to be γ = 0.6±0.4, which is consistent with either a DM cusp (γ = 1) or with a DM core (γ = 0). Our findings are consistent with a Salpeter IMF, and marginally consistent with a Kroupa IMF. We infer tangential orbits for the blue GCs, and radial anisotropy for red GCs and stars. The modelling results are consistent with the virial mass-concentration relation predicted by ΛCDM simulations. The virial mass of NGC 1407 is log M vir = 13.3 ± 0.2M ⊙ , whereas the stellar mass is log M * = 11.8 ± 0.1M ⊙ . The overall uncertainties on the mass of NGC 1407 are only 5 per cent at the projected stellar effective radius. We attribute the disagreement between our results and previous X-ray results to the gas not being in hydrostatic equilibrium in the central regions of the galaxy. The halo of NGC 1407 is found be DM dominated, with a dynamical mass-to-light ratio of M/L = 260 +174 −100 M ⊙ /L ⊙,B . However, this value can be larger up to a factor of 3 depending on the assumed prior on the DM scale radius.