Birk Emil Karlsen-Baeck
--- a/blond/interfaces/xsuite.py 0 → 100644
+++ b/blond/interfaces/xsuite.py 0 → 100644
+import numpy as np
+from scipy.constants import c as clight
+
+from xtrack import ReferenceEnergyIncrease, ZetaShift
+from blond.trackers.tracker import RingAndRFTracker
+from blond.impedances.impedance import TotalInducedVoltage
+
+if TYPE_CHECKING:
+    from xtrack import Particles
+    from xtrack import Line
+    from blond.beam.beam import Beam
+    from blond.beam.profile import Profile
+    from typing import Any, Union
+
+
+def blond_to_xsuite_transform(dt: Union[float, np.ndarray], de: Union[float, np.ndarray], beta0: float,
--- a/blond/interfaces/xsuite.py 0 → 100644
+++ b/blond/interfaces/xsuite.py 0 → 100644
+        Synchronous energy.
+    :param omega_rf: float.
+        The rf angular frequency.
+    :param phi_s: float.
+        Synchronous phase in radians equivalent to Xsuite's phi_rf (below transition energy input should be phi_s-phi_rf).
+        The default value is 0.
+
+    :return zeta, ptau: as numpy-arrays (or single variable)
+    '''
+
+    ptau = de / (beta0 * energy0)
+    zeta = -(dt - phi_s / omega_rf) * beta0 * clight
+    return zeta, ptau
+
+
+def xsuite_to_blond_transform(zeta: Union[float, np.ndarray], ptau: Union[float, np.ndarray],