diff --git a/src/qibocal/protocols/characterization/couplers/coupler_chevron.py b/src/qibocal/protocols/characterization/couplers/coupler_chevron.py
index 2c63e4361..c091ea077 100644
--- a/src/qibocal/protocols/characterization/couplers/coupler_chevron.py
+++ b/src/qibocal/protocols/characterization/couplers/coupler_chevron.py
@@ -66,10 +66,11 @@ def _aquisition(
 
     # create a DataUnits object to store the results,
     data = ChevronData()
-    # sort high and low frequency qubit
-    for pair in qubits:
-        sequence = PulseSequence()
+    sequence = PulseSequence()
 
+    for pair in qubits:
+        print(pair)
+        # sort high and low frequency qubit
         ordered_pair = order_pair(pair, platform.qubits)
         q_lf = ordered_pair[0] # low frequency qubit
         q_hf = ordered_pair[1] # high frequency qubit
@@ -82,6 +83,7 @@ def _aquisition(
         sequence.add(RX_q_lf)
 
         # Coupler Flux pulse applied to the coupler between qubits: q_lf, q_hf
+        print(type(q_lf))
         flux_coupler_pulse = platform.create_coupler_pulse(q_lf, start=sequence.finish)
         sequence.add(flux_coupler_pulse)
 
@@ -137,12 +139,12 @@ def _aquisition(
 
         # TODO: Explore probabilities instead of magnitude
         data.register_qubit(
-            ordered_pair[0],
-            ordered_pair[1],
+            q_lf,
+            q_hf,
             delta_duration_range,
             delta_amplitude_range,
-            results[ordered_pair[0]].magnitude,
-            results[ordered_pair[1]].magnitude,
+            results[q_lf].magnitude,
+            results[q_hf].magnitude,
         )
 
     return data
diff --git a/src/qibocal/protocols/characterization/test_state_1_flux.py b/src/qibocal/protocols/characterization/test_state_1_flux.py
index 405fb1a29..82ee42a8b 100644
--- a/src/qibocal/protocols/characterization/test_state_1_flux.py
+++ b/src/qibocal/protocols/characterization/test_state_1_flux.py
@@ -187,6 +187,7 @@ def _acquisition(
         state1_sequence.add(RX_pulses[qubit])
         
         # Coupler Flux pulse applied to the coupler
+        print(type(qubit))
         flux_coupler_pulse = platform.create_coupler_pulse(qubit, start=state1_sequence.finish)
         state1_sequence.add(flux_coupler_pulse)
 
@@ -228,6 +229,10 @@ def _acquisition(
             acquisition_type=AcquisitionType.INTEGRATION,
         ),
     )
+
+    state1_sequence.plot("./state1_sequence.png")
+    print(state1_sequence)
+    
     # retrieve and store the results for every qubit
     for qubit in qubits:
         result = state1_results[ro_pulses[qubit].serial]