I'm studying G2++ Model in Brigo(2007)'s book.
The model constructed as follows,
$$ r(t) = x(t) + y(t) + φ(t), \quad r(0) = r_0\\ $$ with the dynamics of $dx(t)$ and $dy(t)$ described by: \begin{align} dx(t) &= -ax(t)dt + σdW_1(t), \quad x(0) = 0,\\ dy(t) &= -by(t)dt + ηdW_2(t),\quad y(0) = 0,\\ \end{align} and $dW_1(t)\cdot dW_2(t) = ρdt$.
Problem: When we calculate the variance, There is something that I cannot derive, which is described below:
$$∫^T_t(T-u)dx(u) = -a∫^T_t(T-u)x(u)du + σ∫^T_t(T-u)dW_1(u)$$ Then, $$∫^T_t(T-u)x(u)du = x(t)∫^T_t(T-u)e^{-a(u-t)}du + σ∫^T_t(T-u)∫^T_te^{-a(u-s)}dW_1(s)du.$$
I tried to derive the above integral, $\int_t^T (T-u) x(u) \: du$, but failed.
I want to know how to derive this equation. Especially, where in the world $e^{-a(u-t)}$ and $σ∫^T_t(T-u)∫^T_te^{-a(u-s)}dW_1(s)du$ came from, from the above equation?