Improving NLO-parton shower matched simulations with higher order matrix elements