Given a $2\rightarrow2$ body scattering problem in CM frame, i.e., $\gamma_1\gamma_2\rightarrow\gamma_3\gamma_4$ process, and the coordinate is chosen such that $\gamma_1$ is moving along $+z$-axis, with scattering angle $\theta$ of $\gamma_3$ in the $xz$-plane.
Then the polarization vectors for for $\gamma_1$ is by convention:
\begin{align} \hat{e}_+&=-\frac{1}{\sqrt{2}}(\hat{x}+i\hat{y})\\ \hat{e}_-&=\frac{1}{\sqrt{2}}(\hat{x}-i\hat{y}) \end{align}
and for $\gamma_3$, one just needs to apply $R_y(\theta)$ on the vectors above.
However, how do I find the polarization vectors for $\gamma_2$ and $\gamma_4$, where they have momentum $\vec{k}_2=-\vec{k}_1$ and $\vec{k}_3=-\vec{k}_4$?
My idea is to apply a rotation in $y$-direction with angle $\pi$ and $\pi-\theta$ on polarization vectors of $\gamma_1$ and $\gamma_3$ correspondingly.
In short, the question would be how to find the polarization vectors with inverted momentum $\vec{k}$?