Obtaining the propagator for the Faddeev-Popov (FP) ghosts from the path integral language is straightforward. It is simply
$$\langle T(c(x) \bar c(y))\rangle~=~\int\frac{d^4 p}{(2\pi)^4}\frac{i e^{-ip.(x-y)}}{p^2-\xi m^2+i\epsilon}.$$
But I am unable to properly derive it following canonical quantization route.
The problem is that for anti-commuting fields, the time-ordered product is defined as
$$\langle T(c(x) \bar c(y))\rangle=\theta(x^0-y^0)\langle c(x) \bar c(y)\rangle-\theta(y^0-x^0)\langle \bar c(y) c(x)\rangle$$
with a minus sign between the two terms. This minus sign is preventing me from closing the contours in the correct way to obtain the expression in my first equation.
The only way I can save this is to say that FP ghosts are special, and their time-ordered product is defined with a plus sign instead of the minus sign. Is this legitimate? What is the right way to get the Ghost propagator following canonical quantization route?