Optimality systems for optimal control problems governed by elliptic variational inequalities are derived. Existence of appropriately defined Lagrange multipliers is proved. A primal-dual active set method is proposed to solve the optimality systems numerically. Examples with and without lack of strict complementarity are included.