Atomic and maximal function characterizations of Musielak–Orlicz–Hardy spaces associated to non-negative self-adjoint operators on spaces of homogeneous type