In some references (see for example Ballentine ch. 18.5) the ground state of the BCS theory is assumed to be
\begin{equation} |BCS\rangle = \prod_{\bf k} (u_{\bf k}+v_{\bf k}\hat{c}^{\dagger}_{\bf k,\uparrow}\hat{c}^{\dagger}_{-\bf k,\downarrow})|0\rangle. \end{equation}
with the normalization constraint $u_{\bf k}^2 + v_{\bf k}^2 = 1$.
The same can be found in the original paper by Bardeen, Cooper and Schrieffer where you can read
\begin{equation} |BCS\rangle = \prod_{\bf k}( (1-h_{\bf{k}})^{1/2} + h_{\bf{k}}^{1/2} b^{\dagger}_k |0\rangle. \end{equation}
But in other books like the one by Annett those parameters are considered complex.
Why the parameter $u_{\bf k}$ and $v_{\bf k}$ are taken real? Shouldn't they be taken complex amplitudes? Is this state a classical mixture?
My guess is that this choice is justified a posteriori, but I cannot find any detail about it in any reference I have read.