code: def rotoreflection(axis, angle, origin=(0, 0, 0)):
          
          rot = SymmOp.from_origin_axis_angle(origin, axis, angle)
          refl = SymmOp.reflection(axis, origin)
          m = np.dot(rot.affine_matrix, refl.affine_matrix)
          return SymmOp(m)
docs: Returns a roto-reflection symmetry operation
  
          Args:
              axis (3x1 array): Axis of rotation / mirror normal
              angle (float): Angle in degrees
              origin (3x1 array): Point left invariant by roto-reflection.
                  Defaults to (0, 0, 0).
  
          Return:
              Roto-reflection operation