In Peskin and Schroeder the Dirac equation is solved in the rest frame for solutions with positive frequency: $$\psi(x) = u(p) e^{-ip\cdot x}$$ $$u(p_0) = \sqrt{m} \begin{pmatrix} \xi \\ \xi \end{pmatrix},$$ for any numerical two-component spinor $\xi.$ Boosting to any other frame yields the solution: $$u(p) = \begin{pmatrix} \sqrt{p \cdot \sigma}\ \xi \\ \sqrt{p\cdot \bar{\sigma}}\ \xi \end{pmatrix},$$ where in taking the square root of a matrix, we take the positive root of each eigenvalue.
Then, they summarize: The general solution of the Dirac equation can be written as a linear combination of plane waves. The positive frequency waves are of the form $$\psi(x) = u(p)e^{-ip\cdot x}, \ \ \ p^2 = m^2, \ \ \ p^0 >0.$$ And there are two linearly independent solutions for $u(p),$ $$u^s (p) = \begin{pmatrix} \sqrt{p \cdot \sigma}\ \xi^s \\ \sqrt{p\cdot \bar{\sigma}}\ \xi^s \end{pmatrix}, \ \ \ s=1, 2, $$ which are normalized: $$\bar{u}^r (p) u^s (p) = 2m \delta^{rs}.$$
Next, we can consider the unpolarized cross section for $e^+e^{-} \to \mu^+ \mu^-$ to lowest order. The amplitude is given by: $$\bar{v}^{s'} \left(p'\right) \left(-ie\gamma^{\mu}\right)u^s\left(p\right)\left(\frac{-ig_{\mu \nu}}{q^2}\right)\bar{u}^{r} \left(k\right) \left(-ie\gamma^{\nu}\right)v^{r'}\left(k'\right)$$
Then, I quote
In most experiments the electron and positron beams are unpolarized, so the measured cross section is an average over the electron and positron spins $s$ and $s'$. Muon detectors are normally blind to polarization, so the measured cross section is a sum over the muon spins $r$ and $r'.$ ... We want to compute $$\frac{1}{2}\sum_s \frac{1}{2} \sum_{s'} \sum_r \sum_{r'}|M(s, s' \to r, r')|^2.$$
Why, in order to take the average and the sum, do we only need to sum, rather than integrate, over the spin phase space? Doesn't each incoming particle have an infinite number of spinors? Assuming it is unpolarized, the probability distribution will be uniform, but in principle it still seems like we should integrate over some $\theta$ for spinors $\xi = \begin{pmatrix} \cos \theta \\ \sin \theta \end{pmatrix}$. Instead, what has been done is to assume each incoming particle is in a definite state of "spin-up" or "spin-down" and assign a $50/50$ probability to each.