def pos_in_box(pos, lbox):
  """Positions in [-lbox/2, lbox/2)

  Args:
    pos (np.array): positions in open BC
    lbox (float): cubic box side length
  Return:
    np.array: positions in box centered at 0
  """
  return (pos+lbox/2.) % lbox - lbox/2.