Charged particle sputtering, micrometeoroid impact vaporization, and photon-stimulated desorption are fundamental processes operating at airless surfaces throughout the solar system. At larger bodies, such as Earth's Moon and several of the outer planet moons, these processes generate tenuous surface-bound exospheres that have been observed by a variety of methods. Phobos and Deimos, in contrast, are too gravitationally weak to keep ejected neutrals bound and, thus, are suspected to generate neutral tori in orbit around Mars. While these tori have not yet been detected, the distribution and density of both the neutral and ionized components are of fundamental interest. We combine a neutral Monte Carlo model and a hybrid plasma model to investigate both the neutral and ionized components of the Phobos torus. We show that the spatial distribution of the neutral torus is highly dependent on each individual species (due to ionization rates that span nearly 4 orders of magnitude) and on the location of Phobos with respect to Mars. Additionally, we present the flux distribution of torus pickup ions throughout the Martian system and estimate typical pickup ion fluxes. We find that the predicted pickup ion fluxes are too low to perturb the ambient plasma, consistent with previous null detections by spacecraft around Mars.