Improving the modified XFEM for optimal high‐order approximation