Monte Carlo ICRH simulations in fully shaped anisotropic plasmas