err[m] = np.sqrt(shift[0]*shift[0] + shift[1]*shift[1])