In this contribution calculations of leading and next-to-leading order (hard photon radiation) matrix element for Drell--Yan process (qq \rightarrow Z \rightarrow ll) are presented. Two different spin amplitude formalisms are used in the calculations. Results are compared numerically and cross-checked against the published calculations based on the Dirac-trace methods.

PACS numbers: 11.80.Cr