I suggest that you do it algebraically, not with a calculator. Then you should be able to do it according to your original idea. The other benefit is that you see immediately from the formula how it would be for other values of the energy E ( = 2TeV in your case )
Otherwise, I think it can be done without using a boost, like this:
Assume that the one moving proton has a 4-momentum P1. The stationary proton's 4-momentum is denoted P2.
The total 4-momentum is therefore P = P1 + P2.
Now, the energy of two particle system in its rest frame (COM frame) is the same as the system's rest mass. This can be found by taking minus the square root of the square of P. This is analogous to the rest mass of a single particle, which satisfies P^2 = -M^2 where M is the proton rest mass.
You would demand that this is equal to E:
(P1 + P2)^2 = -E^2
thus
P1^2 + P2^2 + 2P1.P2 = -E^2
Then I use:
P1^2 = P2^2 = -M^2
and
P1.P2 = -M*sqrt(M^2+p^2)
where I used the components of the vectors P1 and P2 in the lab frame:
P1 = (sqrt(M^2+p^2), p)
P2 = (M, 0)
I have restricted to one spatial dimension, and written a lowercase p for the spatial component of the momentum.
The equation can then be reduces to:
-2M^2 - 2Msqrt(M^2+p^2) = -E^2
or
sqrt(M^2+p^2) = (E^2 - 2M^2) / 2M
But the left hand side is actually the energy of the incoming proton that has 3-momentum p, so that would be the answer to the question.
Please note that I may have made some errors here, so please check that you agree with everything yourself. Please also double-check signs. I have used the (-,+,+,+) signature. Maybe you use the opposite one. If so, you could correct for that.
Torquil