diff --git a/Framework/Utilities/SaveBoostHistogram.hpp b/Framework/Utilities/SaveBoostHistogram.hpp
index 813099c4c3a13f9ba3623dedf6e0590e8b12597e..1a30b491c2f2841b5227cc3d9ba6295795db874e 100644
--- a/Framework/Utilities/SaveBoostHistogram.hpp
+++ b/Framework/Utilities/SaveBoostHistogram.hpp
@@ -72,16 +72,12 @@ namespace corsika::utl {
         cnpy::npz_save(filename, std::string{"bins_"} + std::to_string(i), bins.data(),
                        {bins.size()}, "a");
       }
-
-      cnpy::npz_save(filename, std::string{"axistypes"}, axis_types.data(),
-                     {axis_types.size()}, "a");
-
-      cnpy::npz_save(filename, std::string{"overflow"}, overflow.get(),
-                     {axis_types.size()}, "a");
-      cnpy::npz_save(filename, std::string{"underflow"}, underflow.get(),
-                     {axis_types.size()}, "a");
     }
 
+    cnpy::npz_save(filename, std::string{"axistypes"}, axis_types.data(), {rank}, "a");
+    cnpy::npz_save(filename, std::string{"overflow"}, overflow.get(), {rank}, "a");
+    cnpy::npz_save(filename, std::string{"underflow"}, underflow.get(), {rank}, "a");
+
     auto const prod_axis_size = std::accumulate(axes_dims.cbegin(), axes_dims.cend(),
                                                 unsigned{1}, std::multiplies<>());
     auto temp = std::make_unique<float[]>(prod_axis_size);