I am looking for proof of the "well-known" result that for a $4$-dimensional Riemannian manifold $(M, g)$, we have an isomorphism $$ \operatorname{Sym}^2_0(T^*M) \simeq \Lambda_- \otimes \Lambda_+ $$ between the traceless, symmetric, covariant $2$-tensors on $M$ and the tensor product of the spaces of anti-self-dual and self-dual $2$-forms on $M$.
I saw in some papers on Einstein manifolds that the isomorphism is the map $$ \Lambda_- \otimes \Lambda_+ \ni \omega_- \otimes \omega_+ \mapsto \omega_- \omega_+ \in \operatorname{Sym}^2_0(T^*M), $$ where $\omega_- \omega_+$ is given, in coordinates, as $$ (\omega_- \omega_+)_{ij} = (\omega_-)_{ik} g^{kl} (\omega_+)_{lj} $$ (this is the notation from, for instance, the proof of Proposition 4.5.3 in "Wormholes in ACH Einstein manifolds" by Biquard and Rollin). However, none of these references provides a proof of this result. While I understand the definition of $\omega_- \omega_+$, I do not know how to prove that it is the desired isomorphism; I am not sure how to even check that $\omega_- \omega_+$ is a symmetric tensor.
Remark. I also posted this question on Stack Exchange here a month ago, but I have received no answers.