Evaluating massive planar two-loop tensor vertex integrals

Using the parallel/orthogonal space method, we calculate the planar two-loop three-point diagram and two rotated reduced planar two-loop three-point diagrams. Together with the crossed topology, these diagrams are the most complicated ones in the two-loop corrections necessary, for instance, for the decay of the Z0 boson. Instead of calculating particular decay processes, we present a new algorithm which allows us to perform arbitrary next-to-next-to-leading order (NNLO) calculations for massive planar two-loop vertex functions in the general mass case. All integration steps up to the last two are performed analytically and will be implemented under XLOOPS as part of the Mainz XLOOPS-GiNaC project. The last two integrations are done numerically using methods like VEGAS and Divonne. Thresholds originating from Landau singularities are found and discussed in detail. In order to demonstrate the numerical stability of our methods we consider particular Feynman integrals which contribute to different physical processes. Our results can be generalized to the case of the crossed topology. (orig.)