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)
