  In the last two years substantial progress for the simulation of the process:
$\tau^- \to \pi^-\pi^+\pi^- \nu_\tau$ by the Monte Carlo generator TAUOLA was
achieved. It is related to a new parametrization of the corresponding hadronic
current based on the Resonance Chiral Lagrangian and the recent availability of
the unfolded distributions from the BaBar Collaboration analysis for all
invariant hadronic masses. The theoretical model parameters were fitted to the
one-dimensional distributions provided by the BaBar Collaboration and results
of the fit are discussed. A set of the hadronic currents for other final states
with two and three pseudoscalars is also installed in TAUOLA and the
preliminary results for fitting $K^+K^-\pi^-$ and $\pi^0\pi^-$ to BaBar and
Belle data are presented.
