diff --git a/mixnet/ordering/src/bin/aggregate.rs b/mixnet/ordering/src/bin/aggregate.rs
index d227de8..7c7211a 100644
--- a/mixnet/ordering/src/bin/aggregate.rs
+++ b/mixnet/ordering/src/bin/aggregate.rs
@@ -49,6 +49,7 @@ fn aggregate(path: &str) {
             //     "data_msg_count_",
             // );
             add_stats_columns(&mut df, entry.path().join("latency_stats.csv"), "latency_");
+            add_stats_columns(&mut df, entry.path().join("hops_stats.csv"), "hops_");
 
             for prefix in ["strong", "casual", "weak"] {
                 let coeff_stats_path = entry.path().join(format!("{}_coeff_stats.csv", prefix));
diff --git a/mixnet/ordering/src/bin/hops.rs b/mixnet/ordering/src/bin/hops.rs
new file mode 100644
index 0000000..e9bbd8a
--- /dev/null
+++ b/mixnet/ordering/src/bin/hops.rs
@@ -0,0 +1,74 @@
+use glob::glob;
+use polars::prelude::*;
+use std::env;
+use std::fs::File;
+use walkdir::WalkDir;
+
+fn aggregate(path: &str) {
+    for entry in WalkDir::new(path)
+        .into_iter()
+        .filter_map(|e| e.ok())
+        .filter(|e| e.file_type().is_dir())
+    {
+        let dir_name = entry.path().file_name().unwrap().to_string_lossy();
+        if dir_name.starts_with("paramset_") {
+            let mut aggregated_series = Series::new_empty("", &DataType::Int64);
+            let pattern = format!("{}/**/hops.csv", entry.path().display());
+
+            for file in glob(&pattern).unwrap().filter_map(Result::ok) {
+                let df = CsvReadOptions::default()
+                    .with_has_header(true)
+                    .try_into_reader_with_file_path(Some(file.clone()))
+                    .unwrap()
+                    .finish()
+                    .unwrap();
+                aggregated_series
+                    .extend(
+                        &df.column("hops")
+                            .unwrap()
+                            .i64()
+                            .unwrap()
+                            .clone()
+                            .into_series(),
+                    )
+                    .unwrap();
+
+                println!("Processed {}", file.display());
+            }
+
+            let output_file = format!("{}/hops_stats.csv", entry.path().display());
+            save_stats(&aggregated_series, &output_file);
+        }
+    }
+}
+
+fn save_stats(aggregated: &Series, outpath: &str) {
+    let min = aggregated.min::<f64>().unwrap();
+    let max = aggregated.max::<f64>().unwrap();
+    let mean = aggregated.mean().unwrap();
+    let median = aggregated.median().unwrap();
+    let std = aggregated.std(1).unwrap();
+
+    let mut df = DataFrame::new(vec![
+        Series::new("min", &[min]),
+        Series::new("median", &[median]),
+        Series::new("mean", &[mean]),
+        Series::new("std", &[std]),
+        Series::new("max", &[max]),
+    ])
+    .unwrap();
+
+    let mut file = File::create(outpath).unwrap();
+    CsvWriter::new(&mut file).finish(&mut df).unwrap();
+    println!("Saved {}", outpath);
+}
+
+fn main() {
+    let args: Vec<String> = env::args().collect();
+    if args.len() < 2 {
+        eprintln!("Usage: {} <path>", args[0]);
+        std::process::exit(1);
+    }
+    let path = &args[1];
+    aggregate(path);
+}