diff --git a/tools/RegressionSceneData.py b/tools/RegressionSceneData.py
index ea8fc56..d43db5f 100644
--- a/tools/RegressionSceneData.py
+++ b/tools/RegressionSceneData.py
@@ -6,9 +6,6 @@
 import pathlib
 
 import Sofa
-import Sofa.Gui
-
-debug_info = False
 
 def is_simulated(node):
     if node.hasODESolver():
@@ -23,6 +20,37 @@ def is_simulated(node):
     return False
 
 
+class ReplayState(Sofa.Core.Controller):
+    def __init__(self, node, slave_mo, state_filename, **kwargs):
+        super().__init__(**kwargs)
+        self.node = node
+        self.slave_mo = slave_mo
+        self.keyframes = []
+        self.frame_step = 0
+        self.t_sim = 0.0
+
+        with gzip.open(state_filename, 'r') as zipfile:
+            self.ref_data = json.loads(zipfile.read().decode('utf-8'))
+            for key in self.ref_data:
+                self.keyframes.append(float(key))
+        
+        if (self.keyframes[0] == 0.0): # frame 0.0
+            tmp_position = np.asarray(self.ref_data[str(self.keyframes[0])])
+            self.slave_mo.position = tmp_position.tolist()
+            self.frame_step = 1
+        #print('nbr Frames', len(self.keyframes))
+           
+    def onAnimateEndEvent(self, event):
+       dt = float(self.node.getRootContext().dt.value)
+       self.t_sim += dt
+
+       if abs(self.t_sim - self.keyframes[self.frame_step]) < 0.000001:
+           tmp_position = np.asarray(self.ref_data[str(self.keyframes[self.frame_step])])
+           self.slave_mo.position = tmp_position.tolist()
+           #print('tmp_position: ', tmp_position.shape)
+           self.frame_step += 1
+
+
 class NumpyArrayEncoder(JSONEncoder):
     def default(self, obj):
         if isinstance(obj, np.ndarray):
@@ -32,7 +60,7 @@ def default(self, obj):
 
 class RegressionSceneData:
     def __init__(self, file_scene_path: str = None, file_ref_path: str = None, steps = 1000,
-                 epsilon = 0.0001, meca_in_mapping = True, dump_number_step = 1, disable_progress_bar = False):
+                 epsilon = 0.0001, meca_in_mapping = True, dump_number_step = 1, disable_progress_bar = False, verbose = False):
         """
         /// Path to the file scene to test
         std::string m_fileScenePath;
@@ -63,6 +91,7 @@ def __init__(self, file_scene_path: str = None, file_ref_path: str = None, steps
         self.regression_failed = False
         self.root_node = None
         self.disable_progress_bar = disable_progress_bar
+        self.verbose = verbose
 
     def print_info(self):
         print("Test scene: " + self.file_scene_path + " vs " + self.file_ref_path + " using: " + str(self.steps)
@@ -99,9 +128,12 @@ def parse_node(self, node, level = 0):
     def add_compare_state(self):
         counter = 0
         for meca_obj in self.meca_objs:
-            _filename = self.file_ref_path + ".reference_" + str(counter) + "_" + meca_obj.name.value + "_mstate" + ".txt.gz"
+            #_filename = self.file_ref_path + ".reference_" + str(counter) + "_" + meca_obj.name.value + "_mstate" + ".txt.gz"
+            _filename = self.file_ref_path + ".reference_mstate_" + str(counter) + "_" + meca_obj.name.value + ".json.gz"
             
-            meca_obj.getContext().addObject('CompareState', filename=_filename)
+            compareNode = meca_obj.getContext().addChild("CompareStateNode_"+str(counter))
+            cloudPoint = compareNode.addObject('VisualPointCloud', pointSize=10, drawMode="Point", color="green")
+            compareNode.addObject(ReplayState(node=compareNode, slave_mo=cloudPoint, state_filename=_filename))
             counter = counter+1
 
 
@@ -217,7 +249,7 @@ def compare_references(self):
                     full_dist = np.linalg.norm(data_diff)
                     error_by_dof = full_dist / float(data_diff.size)
                     
-                    if debug_info:
+                    if self.verbose:
                         print (str(step) + "| " + self.meca_objs[meca_id].name.value + " | full_dist: " + str(full_dist) + " | error_by_dof: " + str(error_by_dof) + " | nbrDofs: " + str(data_ref.size))
 
                     self.total_error[meca_id] = self.total_error[meca_id] + full_dist
@@ -243,10 +275,17 @@ def compare_references(self):
         return True
     
 
-    def replayReferences(self):
-        Sofa.Gui.GUIManager.Init("myscene", "qglviewer")
+    def replay_references(self):
+        
+        # Import the GUI package
+        import SofaImGui
+        import Sofa.Gui
+        #supported_gui = Sofa.Gui.GUIManager.ListSupportedGUI(",")
+        self.root_node.bbox = "-10 -10 -10 10 10 10"
+        #print ("Supported GUIs are " + supported_gui)
+        Sofa.Gui.GUIManager.Init("myscene", "imgui")
         Sofa.Gui.GUIManager.createGUI(self.root_node, __file__)
-        Sofa.Gui.GUIManager.SetDimension(1080, 1080)
+        Sofa.Gui.GUIManager.SetDimension(1920, 1080)
         Sofa.Gui.GUIManager.MainLoop(self.root_node)
         Sofa.Gui.GUIManager.closeGUI()
 
diff --git a/tools/RegressionSceneList.py b/tools/RegressionSceneList.py
index 3ae6a5f..867cd81 100644
--- a/tools/RegressionSceneList.py
+++ b/tools/RegressionSceneList.py
@@ -61,7 +61,7 @@ def process_file(self):
             if len(values) == 5:
                 scene_data = RegressionSceneData.RegressionSceneData(full_file_path, full_ref_file_path,
                                                                      values[1], values[2], values[3], values[4],
-                                                                     self.disable_progress_bar)
+                                                                     self.disable_progress_bar, self.verbose)
             
                 #scene_data.printInfo()
                 self.scenes_data_sets.append(scene_data)
@@ -119,6 +119,7 @@ def compare_all_references(self):
 
     def replay_references(self, id_scene):
         self.scenes_data_sets[id_scene].load_scene()
+        self.scenes_data_sets[id_scene].add_compare_state()
         self.scenes_data_sets[id_scene].replay_references()