Suppose that $X$ and $Y$ are normally distributed with mean zero and nonzero covariance. I want to know the distributions of $X | X - Y > c$ and $Y | X - Y > c$, which I believe should be jointly distributed as a multivariate normal distribution. I can compute the conditional expectations, using this answer and this answer, but I also want to know the variances and covariance.
That is, I want to compute $$\text E[X | X - Y > c], \quad \text E[Y | X - Y > c],$$ $$\text {Var}(X | X - Y > c), \quad \text {Var}(Y | X - Y > c),$$ and $$\text {Cov}(X, Y | X - Y > c).$$ How do I go about this?