So in Peskin and Schoreder, when computing the amplitude of $e^+e^-\rightarrow\mu^+\mu^-$, summing up over spin they write
\begin{align}\sum_{s,s'}\bar{v}^{s'}_a(p_2)\gamma^\mu_{ab}u^s_b(p_1)\bar{u}^{s}_c(p_1)\gamma^\nu_{cd}v^s_d(p_2) & =(\not{p}_2-m)_{da}\gamma^\mu_{ab}(\not{p}_1+m)_{bc}\gamma_{cd}^\nu\\ &= \textrm{tr}[(\not{p}_2-m)\gamma^\mu(\not{p}_1+m)\gamma^\nu] \end{align}
Why is that a trace?
As I understand that $a$, $b$, $c$, $d$ indexes are the matrix indexes and this is the combination to make a trace. Can someone please clarify, I can't find the answer anywhere.