I am recently puzzled by one question in CFT. I want to compute the correlator of $$\langle T(z)J(w)O_1(x_1) O_2(x_2)\rangle$$ where $T$ is the stress tensor, $J$ is the $U(1)$ Kac-Moody current, $O_1,O_2$ have the same weights $h$ but opposite charges $q$. How can I compute the correlator?
I think I can use the Ward identity corresponding to stress tensor and Kac-Moody current. Starting with $$\langle O_1(x_1) O_2(x_2)\rangle=\frac{1}{(x_1 - x_2)^{2 h}} $$ I can compute $$\langle TO_1O_2\rangle=\frac{h \left(x_1-x_2\right){}^{2-2 h}}{\left(x_1-z\right){}^2 \left(x_2-z\right){}^2} $$ using stress Ward identity. Then I use the Ward identity for $J$ to compute
$$F_1=\langle JTO_1O_2\rangle=\frac{h q \left(x_1-x_2\right){}^{3-2 h}}{\left(w-x_1\right) \left(w-x_2\right) \left(x_1-z\right){}^2 \left(x_2-z\right){}^2}$$
On the other hand, I can change the order: I can first compute $$\langle JO_1O_2\rangle=\frac{q \left(x_1-x_2\right){}^{1-2 h}}{\left(w-x_1\right) \left(w-x_2\right)} $$
and then compute
$$F_2=\langle TJO_1O_2\rangle=\frac{q \left(x_1-x_2\right){}^{1-2 h} \left(\frac{h \left(x_1-x_2\right){}^2}{\left(w-x_1\right) \left(w-x_2\right)}+\frac{\left(z-x_1\right) \left(z-x_2\right)}{(w-z)^2}\right)}{\left(x_1-z\right){}^2 \left(x_2-z\right){}^2}$$
The problem is that I find the two results do not agree $F_1\neq F_2$. I am not sure what is wrong? Is the Ward identity correct?
In general I am wondering how to compute the correlators of involving different types of currents, like $J,T$ here...Intuitively I am not sure why taking different orders for Ward identity should give the same result, although I believe this is should be guaranteed by associativity of the algebra...