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