An alternative is to use this identity for quadratic forms:
$$ Q(\lambda x+(1-\lambda)y) = \lambda Q(x) + (1-\lambda) Q(y) - \lambda (1-\lambda) Q(x-y). $$
It's then pretty obvious that $\lambda (1-\lambda) Q(x-y)>0$ if $0<\lambda<1$ and $Q>0$, so the convexity is obvious.
Now, you may be concerned about this identity applying here. $\Phi(x)=(Tx,x)$ is a quadratic form for real scalars:
$$ \Phi(\lambda x) = (T(\lambda x,\lambda x) = \lvert \lambda \rvert^2 (Tx,x) = \lambda^2 (Tx,x) = \lambda^2 \Phi(x), $$
and
$$ \Phi(x+y)-\Phi(x)-\Phi(y) = (T(x+y),x+y)-(Tx,x)-(Ty,y) = (Tx,y)+(Ty,x), $$
and $(T(\lambda x),y)+(Ty,\lambda x)=\bar{\lambda}(Tx,y)+\lambda(Ty,x)= \lambda((Tx,y)+(Ty,x))$. We are considering real $\lambda$, so this is all fine.
The identity itself is quite straightforward to prove: if $Q(x+y)-Q(x)-Q(y)=2B(x,y)$ is the bilinear,
$$ Q(\lambda x + (1-\lambda)y) = \lambda^2 Q(x)+(1-\lambda)^2 Q(y) + \lambda (1-\lambda) 2B(x,y) \\
= \lambda^2 Q(x)+(1-\lambda)^2 Q(y) + \lambda (1-\lambda) \left( Q(x)+Q(y)-Q(x-y) \right) \\
= \left( \lambda^2+\lambda(1-\lambda) \right)Q(x) + \left( (1-\lambda)^2 +\lambda (1-\lambda) \right)Q(y) -\lambda(1-\lambda)Q(x-y) \\
= \lambda Q(x) + (1-\lambda) Q(y) -\lambda(1-\lambda)Q(x-y). $$