Simultaneous reconstruction of absorption, scattering and anisotropy factor distributions in quantitative photoacoustic tomography