Construction of high order maps for large proton accelerators