This article is devoted to establishing a least square based weak Galerkin method for second order elliptic equations in non-divergence form using a discrete weak Hessian operator.Naturally,the resulting linear system is symmetric and positive definite,and thus the algorithm is easy to implement and analyze.Convergence analysis in the H2 equivalent norm is established on an arbitrary shape regular polygonal mesh.A superconvergence result is proved when the coefficient matrix is constant or piecewise constant.Numerical examples are performed which not only verify the theoretical results but also reveal some unexpected superconvergence phenomena.