diff --git a/main.cpp b/main.cpp
index 83ff7d3..4389d5d 100644
--- a/main.cpp
+++ b/main.cpp
@@ -203,12 +203,15 @@ int main() {
     }
 
     cout << "[INFO] Starting the calculate the phase shift" << endl;
-    double delta = 0;
+    double deltas;
+    int ndeltas = 0;
     for(int i = 0; i < channelsZeros[0].size(); i++) {
         if(channelsZeros[1].size() <= i) continue;
-        cout << "[INFO] Phase shift for zero " << i << " is " << ((channelsXData[0][channelsZeros[0][i]] - channelsXData[1][channelsZeros[1][i]]) / channelOneT) * 2 * M_PI << endl;
-        delta += channelsZeros[1][i] - channelsZeros[0][i];
+        // cout << "[INFO] Phase shift for zero " << i << " is " << ((channelsXData[0][channelsZeros[0][i]] - channelsXData[1][channelsZeros[1][i]]) / channelOneT) * 360 << endl;
+        deltas += channelsXData[1][channelsZeros[1][i]] - channelsXData[0][channelsZeros[0][i]];
+        ndeltas++;
     }
+    cout << "[INFO] Average phase shift (CHA0 to CHA1) " << ((deltas/ndeltas) / channelOneT) * 360 << endl;
 
     plt::plot(channelsXData[0], channelsYData[0], channelsXData[1], channelsYData[1]);
     plt::xlabel("X Label");