I am a bit confused by some point in the calculation of the electron-electron scattering amplitude in Zee's QFT in a Nutshell, Section II.6. He extracts this quantify to work it out separately:
$$\tau^{\mu\nu}(P,p) = \frac{1}{2} \sum_s \sum_S \bar{u}(P,S) \gamma^\mu u(p,s) \bar{u}(p,s)\gamma^\nu u(P,S)$$ where the $u$ function is the electron function from the planewave solution for the Dirac equation. All Zee tells me about $u$ are its properties:
$$\bar{u}(p,s)u(p,s) = 1$$ $$\sum_s u(p,s)\bar{u}(p,s) = \frac{\gamma^\mu p_\mu +m}{2m}$$
I can perform an intermediate step (that Zee skips) using the second property of $u$:
$$\tau^{\mu\nu}(P,p) = \frac{1}{2(2m)} \sum_S \bar{u}(P,S) \gamma^\mu (\gamma^\rho p_\rho + m)\gamma^\nu u(P,S)$$
But Zee goes one step further and introduces a trace without any explanation, in order to obtain:
$$\tau^{\mu\nu}(P,p) = \frac{1}{2(2m)^2} \mathrm{Tr}\left[(\gamma^\lambda P_\lambda + m) \gamma^\mu (\gamma^\rho p_\rho + m)\gamma^\nu \right]$$
This got me completely confused. Where did the trace come from? How did the $(\gamma^\lambda P_\lambda + m)/(2m)$ appear, considering that the $\bar{u}$ is on the left in the intermediate step and there are also gamma matrices in between them?