diff --git a/main.cpp b/main.cpp
index b659fa7..83ff7d3 100644
--- a/main.cpp
+++ b/main.cpp
@@ -110,6 +110,7 @@ int main() {
 
     vector<vector<double>> channelsXData, channelsYData, channelsZeros;
     vector<double> channelsP2P;
+    double channelOneT = 0;
 
     for(int channel = 0; channel < channels.size(); channel++) {
         string line;
@@ -170,6 +171,10 @@ int main() {
         avgFrequency *= 2;
         cout << "[INFO] Average peak time is " << avgFrequency << "s" << endl;
         cout << "[INFO] Average frequency is " << 1/avgFrequency << "Hz" << endl;
+        if(channel == 0) {
+            channelOneT = avgFrequency;
+            cout << "Set chan 1 freq: " << channelOneT << endl;
+        }
 
         vector<double> mins, maxs;
         for(int i = 1; i < channelsZeros[channel].size(); i++) {
@@ -197,6 +202,14 @@ int main() {
         channelsP2P.push_back(abs(avgMin) + abs(avgMax));
     }
 
+    cout << "[INFO] Starting the calculate the phase shift" << endl;
+    double delta = 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];
+    }
+
     plt::plot(channelsXData[0], channelsYData[0], channelsXData[1], channelsYData[1]);
     plt::xlabel("X Label");
     plt::ylabel("Y Label");