So for the first part of C , all I'm doing is following the computation of Gordon in the link and replacing :
$(t-1) \rightarrow{t}, t \rightarrow{T}, 0 \rightarrow{t} , s\rightarrow{s}$
So I end up with:
\begin{align*}
\ln{C} &= \int_{t}^{T} \int_{t}^{s}e^{-\alpha(T-s)-\beta(s-r)}dW(r) \,ds +... \\
&= \int_{t}^{t} \int_{t}^{T}e^{-\alpha(T-s)-\beta(s-r)}ds\,dW(r) + \int_{t}^{T} \int_{r}^{T}e^{-\alpha(T-s)-\beta(s-r)}ds\,dW(r)+...\\
&= 0 + \frac{e^{-\alpha T}}{\alpha-\beta}\int_{t}^T e^{\beta r} \left[e^{(\alpha -\beta)T}-e^{(\alpha -\beta)r} \right] dW(r)+...\\
&= 0 + \frac{e^{-\beta T}}{\alpha -\beta}\int_{t}^T e^{\beta r} dW(r) - \frac{e^{-\alpha T}}{\alpha -\beta}\int_{t}^T e^{\alpha r} dW(r)+...
\end{align*}
From there its not too complicated to compute variance of the process, but I have no idea if this work like that as I dont fully understand the decomposition in the first place.
EDIT: This is working, tested with some closed form solutions on special cases, limit cases and Monte Carlo simulation.