The matrix $M$ below is the matrix of a linear operator $T$ over $\mathbb{R}^4$ with respect to the canonical basis $\{e_1, e_2, e_3, e_4\}$. $$ M = \sqrt{2} \begin{bmatrix} 1/2 & -1/3 & -1/3 & 1/6 \\ 1/3 & 1/2 & -1/6 & -1/3 \\ 1/3 & 1/6 & 1/2 & 1/3 \\ -1/6& 1/3 & -1/3 &1/2 \end{bmatrix}.$$ I proved that $T$ is an isometry showing that the set $\{T(e_1), T(e_2), T(e_3), T(e_4)\}=\{Me_1, Me_2, Me_3, Me_4\}$ is orthonormal. So, by Theorem 9.36 of Axler [1], there is an orthonormal basis of $\mathbb{R}^4$ with respect to which $T$ has a block diagonal matrix such that each block on the diagonal is a 1-by-1 matrix containing $1$ or $-1$ or is a 2-by-2 matrix of the form $$ \begin{bmatrix} \cos \theta & -\sin \theta \\ \sin \theta & \cos \theta \end{bmatrix},$$ with $\theta \in (0, \pi)$. How can I identify this basis?
[1] Axler, S. Linear Algebra Done Right. $3^{rd}$ edition. Springer: New York, 2015.