Let X , d , μ be a metric measure space endowed with a metric d and a non-negative Borel doubling measure μ . Let L be a non-negative self-adjoint operator on L 2 X . Assume that the (heat) kernel associated to the semigroup e − t L satisfies a Gaussian upper bound. In this paper, we prove that the Musielak–Orlicz–Hardy space H φ , L X associated with L in terms of the Lusin-area function and the Musielak–Orlicz–Hardy space H L , G , φ X associated with L in terms of the Littlewood–Paley function coincide and their norms are equivalent. To do this, we first establish the discrete characterization of these two spaces. It improves the known results in the literature.