diff --git a/pes_to_spec/model.py b/pes_to_spec/model.py
index 644f4696169d26006964f2e65bd38fb158d1c7e9..b696dba65c23de698e0740f721bb639f4062597e 100644
--- a/pes_to_spec/model.py
+++ b/pes_to_spec/model.py
@@ -530,6 +530,9 @@ class SelectRelevantLowResolution(TransformerMixin, BaseEstimator):
                 chi2[ch] = np.sum(((meanX[ch] - self.mean[ch][:,0,:])/self.std[ch][:,0,:])**2)
         shift = {ch: shifts[np.argmin(chi2[ch])]
                 for ch in y.keys()}
+        print("Shifts", shift)
+        y = {ch: np.roll(y[ch], shift[ch], axis=-1)
+             for ch in y.keys()}
 
         if not keep_dictionary_structure:
             selected = list(y.values())