If I'm following what C&T are doing, they are creating an orthonormal basis of one-forms
> [r/sqrt(r^2-N(r,z)^2),0,0,0];
> [0,exp(-v(r,z)/2),0,0];
> [0,0,exp(-v(r,z)/2),0,0];
> [ N(r,z)/sqrt( r^2-N(r,z)^2 ) ,0,0,sqrt(r^2-N(r,z)^2)];
which can be plugged directly into GRtII to give their metric. They are then finding the einstein tensor G in terms of the orthonormal basis, and equating all the pressure terms to zero.
This gives me the first two of their equations in (6). But I can't quite figure out how they are getting the third.
This is sort of a moot point, because I don't think C&T have really thought about the issue of the Komar mass being non-zero for a very short cylinder at z=0. C&T seem to actually admit that this happens, but don't seem to recognize properly the physical significance of this result.