sys.cx(pair1[1], pair2[1]) sys.cx(pair1[1], pair3[1]) sys.barrier() #distillation(1) sys.measure(pair2, check1) sys.measure(pair3, check2) sys.barrier() #checking fidelity through bell-measurement sys.cx(pair1[0], pair1[1]) sys.h(pair1[0]) sys.barrier(pair1) sys.measure(pair1, result) print(sys.draw()) style = {'fold': 24} sys.draw(output='mpl', style=style, plot_barriers=True) from qiskit import BasicAer backend = BasicAer.get_backend('qasm_simulator') job = execute(sys, backend=backend, shots=10000) result = job.result() from qiskit.tools.visualization import plot_histogram counts = result.get_counts(sys) plot_histogram(counts)