Finding the Higgs boson in decays to $Z \gamma$ using the matrix element method at Next-to-Leading Order