From 1ab07d335f6bffd40be881e441199be91b1f1b2d Mon Sep 17 00:00:00 2001
From: Fayorg <ebr@aketechnology.com>
Date: Sun, 25 Feb 2024 22:12:45 +0100
Subject: [PATCH] add: phase shift corrected with average

---
 main.cpp | 9 ++++++---
 1 file changed, 6 insertions(+), 3 deletions(-)

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");