code: def _local_pauli_eig_meas(op, idx):
      
      if op == 'X':
          return Program(RY(-pi / 2, idx))
      elif op == 'Y':
          return Program(RX(pi / 2, idx))
      elif op == 'Z':
          return Program()
      raise ValueError(f'Unknown operation {op}')
docs: Generate gate sequence to measure in the eigenbasis of a Pauli operator, assuming
      we are only able to measure in the Z eigenbasis. (Note: The unitary operations of this
      Program are essentially the Hermitian conjugates of those in :py:func:`_one_q_pauli_prep`)