We study the origin and properties of unbound stars in the kinematic samples of dwarf spheroidal galaxies. For this purpose we have run a high resolution N-body simulation of a two-component dwarf galaxy orbiting in a Milky Way potential. We create mock kinematic data sets by observing the dwarf in different directions. When the dwarf is observed along the tidal tails the kinematic samples are strongly contaminated by unbound stars from the tails. However, most of the unbound stars can be removed by the method of interloper rejection proposed by den Hartog & Katgert. We model the velocity dispersion profiles of the cleaned-up kinematic samples using solutions of the Jeans equation. We show that even for such a strongly stripped dwarf the Jeans analysis, when applied to cleaned samples, allows us to reproduce the mass and mass-to-light ratio of the dwarf with accuracy typically better than 25%.