Feynman Green Function in 1+1D for d'Alembertian I'm trying to obtain the Feynman Green Function (i.e. I'm using the Feynman Causal prescription to compute the green function) for the D'Alembertian in 1+1D, I'm finding
$$G^{(2)}_F (t; \vec x) = \frac{1}{4\pi} \frac{\Theta(t^2-\vec x^2)}{\sqrt{t^2-\vec x^2}} - \frac{i}{4\pi^2} \int_{-\infty}^\infty dz' \text{P.V.}\left(\frac{1}{t^2-\vec x^2 - z'^2}\right).$$
This result is right ou wrong? I don't have any reference with the answer. Supposing the result is right, there is some way to simplify it more?

The way I find the result.
I've used the usual procedure (as in Eleftherios Economou  and in Morse Feshbach) to find the Green function in 1+1D as a potential generated by an infinite line of charge in 2+1D.
$$G^{(2)}_F (t; \vec x) = \int dt' d^3 r' G^{(3)}_C (t-t'; \vec r-\vec r') J(t',\vec r').$$
The Feynman green function in 2+1D is,
$$G^{(3)}_F (t-t'; \vec r-\vec r') = \frac{1}{4\pi} \left[\delta((t-t')^2-||\vec r - \vec r'||) - \frac{i}{\pi} \text{P.V.}\left(\frac{1}{(t-t')^2-||\vec r-\vec r'||}\right)\right]$$
as can be checked on Bogoliubov-Shirkov (Appendix II, pag 605, A2b.6)
And the source is,
$$J(t',\vec r') = \delta(x') \delta(y') \delta(t')$$
So that,
$$G^{(2)}_F (t; \vec x) = \int_{-\infty}^\infty dz'\frac{1}{4\pi} \left[\delta(t^2-\vec x^2 - z'^2) - \frac{i}{\pi} \text{P.V.}\left(\frac{1}{t^2-\vec x^2 - z'^2}\right)\right]$$
Using a basic Dirac Delta property,
$$\delta(x^2-a^2) = \frac{1}{2|a|} \left(\delta(x+a) + \delta(x-a) \right)$$ we get for the first integral,
$$\frac{1}{8\pi} \int_{-\infty}^\infty dz'\frac{1}{||\sqrt{t^2-\vec x^2}||} \left(\delta(z-\sqrt{t^2-\vec x^2}) + \delta(z+\sqrt{t^2-\vec x^2})\right)$$
For $T^2>\vec x^2$ (time-like interval) the points $\pm \sqrt{t^2-\vec x^2}$ are real and belong to interval $(-\infty,\infty)$. So we have (for the first integral),
$$\frac{1}{4\pi} \frac{\Theta(t^2-\vec x^2)}{\sqrt{t^2-\vec x^2}}$$
And, finally,
$$G^{(2)}_F (t; \vec x) = \frac{1}{4\pi} \frac{\Theta(t^2-\vec x^2)}{\sqrt{t^2-\vec x^2}} + \int_{-\infty}^\infty dz' \text{P.V.}\left(\frac{1}{t^2-\vec x^2 - z'^2}\right)$$

Just putting here what happens with my solution if we are put of the light cone singularity ($t=\pm \vec x$). I think that we can forget about the principal value at this case.
If I take the integral and solve it, I get
$$\int_{-\infty}^\infty dz' \left(\frac{1}{t^2-\vec x^2 - z'^2}\right) = \frac{i\pi}{\sqrt{t^2-\vec x^2}}.$$
So I get,
$$G^{(2)}_F (t; \vec x) = \frac{1}{4\pi} \frac{1}{\sqrt{t^2-\vec x^2}} \left(\Theta(t^2-\vec x^2) + 1\right)$$
 A: First of all, the retarded (causal) Green's function for the dalembertian is ($\vec{R} \equiv \vec{x} - \vec{x}'$, $T \equiv t - t'$, $R\equiv |\vec{R}|$), $$G^{(3)}_R(R,T) = \frac{\Theta(T)\delta(T-R/c)}{4\pi R}.$$ This leads to the retarded potential in electromagnetism. It depends only on the difference of position vectors because of the symmetry of the boundary condition in unbounded space, and on the time difference, because the defining equations and boundary condition are invariant if we substitute $t \to t-t'$. 
Using the method of embedding or descent with appropriate boundary conditions, Green's function in dimension 2 is given by ($r \equiv \sqrt{ (x-x')^2 + (y-y')^2}$), $$G^{(2)}_R(r,T) = \frac{\Theta(T)}{4\pi} \int_{-\infty}^\infty \frac{\delta ( T - R/c)}{R} d z = \frac{\Theta(T)}{4\pi} \int_{-\infty}^\infty \frac{\delta\left(  T - \sqrt{ r^2 + (z-z')^2}/c \right)}{\sqrt{ r^2 + (z-z')^2}} d (z-z'),$$ which equals $$G^{(2)}_R(r,T) = \frac{\Theta(T - r/c)}{2\pi \sqrt{T^2- r^2/c^2}},$$ for $T>0$, otherwise $0$.
The principal value solution can be found as the linear combination of the retarded and the advanced functions, according to the expansion $$-\text{PV} \int_{-\infty}^\infty \frac{\phi_p(x) \phi^*_p(x')}{\lambda_p} d p = G^{\pm}  \mp \Lambda$$ in terms of eigenfunctions, where $\Lambda$ is the integral comprising the eigenfunctions with zero eigenvalue. You could also see this with the Sokhotski-Plemelj theorem.
