I am reading Weinberg's QFT book and in 10.4 he introduced a derivation of Ward-Takahashi identity (where $T$ is the time ordering): $$\begin{align} \frac{\partial}{\partial x^\mu}T{\{J^\mu(x)\Psi_n(y)\bar\Psi_m(z)\}} &= T{ \{ \partial_\mu J^\mu(x)\Psi_n(y)\bar\Psi_m(z)\}} \\ &+ \delta(x^0 - y^0) T\{ [J^0(x),\Psi_n(y)]\bar\Psi_m(z)\} \\ &+ \delta(x^0 - z^0) T\{ \Psi_n(y),[J^0(x),\bar\Psi_m(z)]\}. \end{align} \tag{10.4.21}$$ I understand the deltas come from time derivative of step functions in the time ordering operators,but after taking the derivatives, isn't the step function disappear? how do we still have time operators. how does $J^0$ appears in the second and third term on the right hand side?
I actually natively tried to express time ordering as a function of step function, $J^\mu(x),\Psi_n(y)$, and $\bar\Psi_m(z)$ and take the multivariable chain rule of it without luck.