diff --git a/lib/damask/orientation.py b/lib/damask/orientation.py
index e8bb49e13..9f6c80bfd 100644
--- a/lib/damask/orientation.py
+++ b/lib/damask/orientation.py
@@ -966,7 +966,8 @@ class Orientation:
 
     eig, vec = np.linalg.eig(M/N)
 
-    return Orientation(quaternion = Quaternion(quatArray = np.real(vec.T[eig.argmax()])))
+    return Orientation(quaternion = Quaternion(quatArray = np.real(vec.T[eig.argmax()])),
+                       symmetry = reference.symmetry.lattice)
 
 
   def related(self,