diff --git a/mbtrack2/utilities/beamloading.py b/mbtrack2/utilities/beamloading.py
index d8d64b1e842ddf6c5d7cc26016f6aaf17224aedf..041741b61b2975426f819aa697a90ea833125843 100644
--- a/mbtrack2/utilities/beamloading.py
+++ b/mbtrack2/utilities/beamloading.py
@@ -107,7 +107,6 @@ class BeamLoadingEquilibrium():
 
     def center_of_mass(self):
         """Return center of mass position in [s]"""
-        self.update_rho()
         CM = np.average(self.z0, weights=self.rho0)
         return CM / c
 
@@ -166,6 +165,7 @@ class BeamLoadingEquilibrium():
         else:
             self.F = x[::2]
         self.PHI = x[1::2]
+        self.update_rho()
 
         # Compute system
         if self.auto_set_MC_theta: