y = np.real(np.trace(np.dot(rho,np.array([[0,-1j],[1j,0]]))))