In one dimension, given a particle in a quantum state $| \psi\rangle$, the probability density of position is given as $| \psi(x) |^2 = \psi^*(x) \psi(x) =\langle x | \psi \rangle\langle \psi | x \rangle $ where $| x \rangle$ is the base ket for the position operator. I have just learnt about tensor products role in describing states of multi-particles or singular particles where we consider different observable spaces simultaneously.
Consider the three dimensional problem: An electron in a hydrogen atom occupies the combined spin and position state given as $$R_{21} \bigg(\sqrt{\frac{1}{3}} Y_{1}^{0} \otimes \chi_{+} + \sqrt{\frac{2}{3}}Y_{1}^{1} \otimes \chi_{-} \bigg)$$
where $\chi_{+}$ and $\chi_{-}$ are the eigenstates of the z-component of spin (spin up and spin down). If you measure the position of the particle, what is the probability density of finding it at $r, \theta, \phi$?
The given answer suggests you use the scalar product defined for tensor products of Hilbert spaces (i.e. $\langle \phi_1 \otimes \phi_2, \psi_1 \otimes \psi_2 \rangle = \langle \phi_1, \psi_1 \rangle_1 \langle \phi_2, \psi_2 \rangle_2 $ ) so we get $$\begin{align} |\psi|^2 &= |R_{21}|^2\bigg\langle \sqrt{\frac{1}{3} }Y_{1}^{0*} \otimes \chi_{+}^{\dagger} + \sqrt{\frac{2}{3}}Y_{1}^{1*} \otimes \chi_{-}^{\dagger}, \sqrt{\frac{1}{3} }Y_{1}^{0} \otimes \chi_{+} + \sqrt{\frac{2}{3}}Y_{1}^{1} \otimes \chi_{-} \bigg\rangle \end{align}\\ =|R_{21}|^2\bigg( \frac{1}{3}|Y_{1}^{0}|^2 \chi_{+}^{\dagger}\chi_{+} + \frac{\sqrt{2}}{3}Y_{1}^{1*}Y_{1}^{0} \chi_{-}^{\dagger}\chi_{+} + \frac{\sqrt{2}}{3}Y_{1}^{0*}Y_{1}^{1} \chi_{+}^{\dagger}\chi_{-} + \frac{2}{3}Y_{1}^{1}\chi_{-}^{\dagger}\chi_{-} \bigg) \\ = \frac{1}{3}|R_{21}|^{2} \big( |Y_{1}^{0}|^2 + 2|Y_{1}^{1}|^2 \big).$$
Question: But if this is correct then what we have actually done is taken the scalar product $\langle \psi| \psi \rangle$ (which we usually want to be normalized to equal $1$) of a state $| \psi \rangle = R_{21} \bigg(\sqrt{\frac{1}{2}} Y_{1}^{0} \otimes \chi_{+} + \sqrt{\frac{2}{3}}Y_{1}^{1} \otimes \chi_{-} \bigg)$, which would usually give the summation of the probability coefficients rather than the probability density? What am I missing here? Why is this the way to get the probability density or is there a way which is more descriptive?
See my proposed answer below.
Thanks for any assistance.