With[{R = 106, r = 50}, ContourPlot3D[(Sqrt[x^2 + y^2] - R)^2 + z^2 == r^2, {x, -(R + r), R + r}, {y, -(R + r), R + r}, {z, -r, r}]]

BoxRatios->{1,1,1/3}