How do the eigenfunctions of the total angular momentum operator
analytically look like?
I mean the operator is given by $J = L+S$ so the eigenfunctions have
to be tensor-product states, right? Can we explicitely say what they
are?
The eigenfunctions of $J$ are going to be made up of linear combinations of tensor-product states of the eigenfunctions of $L$ and the eigenfunctions of $S$. In general, the linear combinations will involved more than just one tensor-product of the eigenfunctions of $L$ and $S$, however the "stretch-state" is the exception and is the starting point for constructing the others.
For example, if $L$ and $S$ are both spin 1/2 (yes, I know that $L$ usually stands for "orbital" angular momentum, but in this example $L$ and $S$ are both spin 1/2) then total $J$ can be either spin 1=|1/2+1/2| or spin 0=|1/2-1/2|. One of the eigenfunctions of $J=1$ is given by a tensor-product state
$$
|1,1\rangle=|1/2,1/2\rangle\otimes|1/2,1/2\rangle\;,
$$
which is the "stretch-state". The other states can be obtained by applying the lowering operator
$$
J_- = L_{-}\otimes 1 +1\otimes S_{-}\;,
$$
to the stretch-state and normalizing. E.g., we find (I'm putting in the square root of two on the LHS by hand to indicate that the RHS is not generated as a normalized state)
$$
\sqrt{2}|1,0\rangle=|1/2,-1/2\rangle\otimes|1/2,1/2\rangle+|1/2,1/2\rangle\otimes|1/2,-1/2\rangle
$$
and
$$
|1,-1\rangle=|1/2,-1/2\rangle\otimes|1/2,-1/2\rangle\;.
$$
And the $|0,0\rangle$ state is generated by creating a state with $J_z=0$ that is orthogonal to the $|1,0\rangle$ state. It is
$$
\sqrt{2}|0,0\rangle=|1/2,-1/2\rangle\otimes|1/2,1/2\rangle-|1/2,1/2\rangle\otimes|1/2,-1/2\rangle\;.
$$
So, you see that two $(|1,1\rangle$ and $|1,-1\rangle)$ of the eigenfunctions of J in this case are simple tensor-products of the eigenfunctions of L and S, and the other two $(|1,0\rangle$ and $|0,0\rangle)$ are linear combinations of more than one tensor-product of the eigenfunctions of L and S.