diff --git a/benchmarl/eval_results.py b/benchmarl/eval_results.py
new file mode 100644
index 00000000..59503b40
--- /dev/null
+++ b/benchmarl/eval_results.py
@@ -0,0 +1,161 @@
+import collections
+import importlib
+import json
+from os import walk
+from pathlib import Path
+from typing import Dict, List, Optional
+
+_has_marl_eval = importlib.util.find_spec("marl_eval") is not None
+if _has_marl_eval:
+    from marl_eval.plotting_tools.plotting import (
+        aggregate_scores,
+        performance_profiles,
+        plot_single_task,
+        probability_of_improvement,
+        sample_efficiency_curves,
+    )
+    from marl_eval.utils.data_processing_utils import (
+        create_matrices_for_rliable,
+        data_process_pipeline,
+    )
+    from matplotlib import pyplot as plt
+
+
+def get_raw_dict_from_multirun_folder(multirun_folder: str) -> Dict:
+    return _load_and_merge_json_dicts(_get_json_files_from_multirun(multirun_folder))
+
+
+def _get_json_files_from_multirun(multirun_folder: str) -> List[str]:
+    files = []
+    for dirpath, _, filenames in walk(multirun_folder):
+        for file_name in filenames:
+            if file_name.endswith(".json") and "wandb" not in file_name:
+                files.append(str(Path(dirpath) / Path(file_name)))
+    return files
+
+
+def _load_and_merge_json_dicts(
+    json_input_files: List[str], json_output_file: Optional[str] = None
+) -> Dict:
+    def update(d, u):
+        for k, v in u.items():
+            if isinstance(v, collections.abc.Mapping):
+                d[k] = update(d.get(k, {}), v)
+            else:
+                d[k] = v
+        return d
+
+    dicts = []
+    for file in json_input_files:
+        with open(file, "r") as f:
+            dicts.append(json.load(f))
+    full_dict = {}
+    for single_dict in dicts:
+        update(full_dict, single_dict)
+
+    if json_output_file is not None:
+        with open(json_output_file, "w+") as f:
+            json.dump(full_dict, f, indent=4)
+
+    return full_dict
+
+
+class Plotting:
+
+    METRICS_TO_NORMALIZE = ["return"]
+    METRIC_TO_PLOT = "return"
+
+    @staticmethod
+    def process_data(raw_data: Dict):
+        # Call data_process_pipeline to normalize the choosen metrics and to clean the data
+        return data_process_pipeline(
+            raw_data=raw_data, metrics_to_normalize=Plotting.METRICS_TO_NORMALIZE
+        )
+
+    @staticmethod
+    def create_matrices(processed_data, env_name: str):
+        return create_matrices_for_rliable(
+            data_dictionary=processed_data,
+            environment_name=env_name,
+            metrics_to_normalize=Plotting.METRICS_TO_NORMALIZE,
+        )
+
+    ############################
+    # Environment level plotting
+    ############################
+
+    @staticmethod
+    def performance_profile_figure(environment_comparison_matrix):
+        return performance_profiles(
+            environment_comparison_matrix,
+            metric_name=Plotting.METRIC_TO_PLOT,
+            metrics_to_normalize=Plotting.METRICS_TO_NORMALIZE,
+        )
+
+    @staticmethod
+    def aggregate_scores(environment_comparison_matrix):
+        return aggregate_scores(
+            dictionary=environment_comparison_matrix,
+            metric_name=Plotting.METRIC_TO_PLOT,
+            metrics_to_normalize=Plotting.METRICS_TO_NORMALIZE,
+            save_tabular_as_latex=True,
+        )
+
+    @staticmethod
+    def probability_of_improvement(
+        environment_comparison_matrix, algorithms_to_compare: List[List[str]]
+    ):
+        return probability_of_improvement(
+            environment_comparison_matrix,
+            metric_name=Plotting.METRIC_TO_PLOT,
+            metrics_to_normalize=Plotting.METRICS_TO_NORMALIZE,
+            algorithms_to_compare=algorithms_to_compare,
+        )
+
+    @staticmethod
+    def environemnt_sample_efficiency_curves(sample_effeciency_matrix):
+        return sample_efficiency_curves(
+            dictionary=sample_effeciency_matrix,
+            metric_name=Plotting.METRIC_TO_PLOT,
+            metrics_to_normalize=Plotting.METRICS_TO_NORMALIZE,
+        )
+
+    ############################
+    # Task level plotting
+    ############################
+
+    @staticmethod
+    def task_sample_efficiency_curves(processed_data, task, env):
+        return plot_single_task(
+            processed_data=processed_data,
+            environment_name=env,
+            task_name=task,
+            metric_name="return",
+            metrics_to_normalize=Plotting.METRICS_TO_NORMALIZE,
+        )
+
+
+if __name__ == "__main__":
+    raw_dict = get_raw_dict_from_multirun_folder(
+        multirun_folder="/Users/matbet/PycharmProjects/BenchMARL/benchmarl/multirun/2023-09-22/17-21-34"
+    )
+    processed_data = Plotting.process_data(raw_dict)
+    (
+        environment_comparison_matrix,
+        sample_efficiency_matrix,
+    ) = Plotting.create_matrices(processed_data, env_name="vmas")
+
+    Plotting.performance_profile_figure(
+        environment_comparison_matrix=environment_comparison_matrix
+    )
+    Plotting.aggregate_scores(
+        environment_comparison_matrix=environment_comparison_matrix
+    )
+    Plotting.environemnt_sample_efficiency_curves(
+        sample_effeciency_matrix=sample_efficiency_matrix
+    )
+
+    Plotting.task_sample_efficiency_curves(
+        processed_data=processed_data, env="vmas", task="navigation"
+    )
+    plt.show()
diff --git a/benchmarl/experiment/logger.py b/benchmarl/experiment/logger.py
index 47267549..08a1194e 100644
--- a/benchmarl/experiment/logger.py
+++ b/benchmarl/experiment/logger.py
@@ -264,7 +264,7 @@ def __init__(
         seed: int,
     ):
         self.path = Path(folder) / Path(name)
-        self.run_data = {}
+        self.run_data = {"absolute_metrics": {}}
         self.data = {
             environment_name: {
                 task_name: {algorithm_name: {f"seed_{seed}": self.run_data}}
@@ -280,5 +280,13 @@ def write(self, total_frames: int, metrics: Dict[str, Any], step: int):
         else:
             self.run_data[step_str] = metrics
 
+        # Store the maximum of each metric
+        for metric_name in metrics.keys():
+            max_metric = max(metrics[metric_name])
+            if metric_name in self.run_data["absolute_metrics"]:
+                prev_max_metric = self.run_data["absolute_metrics"][metric_name][0]
+                max_metric = max(max_metric, prev_max_metric)
+            self.run_data["absolute_metrics"][metric_name] = [max_metric]
+
         with open(self.path, "w+") as f:
             json.dump(self.data, f, indent=4)