In this paper (cf. Lychagin and Yumaguzhin, in Anal Math Phys, 2016) a class of totally geodesics solutions for the vacuum Einstein equations is introduced. It consists of Einstein metrics of signature (1,3) such that 2-dimensional distributions, defined by the Weyl tensor, are completely integrable and totally geodesic. The complete and explicit description of metrics from these class is given. It is shown that these metrics depend on two functions in one variable and one harmonic function.