diff --git a/Processes/InteractionCounter/InteractionHistogram.cc b/Processes/InteractionCounter/InteractionHistogram.cc
index 556b58f2788f1063a594bb41fb06a552ed5c563d..61ed9d5c8f6ce972146bb656a92ba0d6e533f3bd 100644
--- a/Processes/InteractionCounter/InteractionHistogram.cc
+++ b/Processes/InteractionCounter/InteractionHistogram.cc
@@ -19,7 +19,6 @@ void InteractionHistogram::fill(particles::Code projectile_id,
                                 units::si::HEPEnergyType lab_energy,
                                 units::si::HEPEnergyType mass_target, int A, int Z) {
   using namespace units::si;
-
   if (projectile_id == particles::Code::Nucleus) {
     auto const sqrtS =
         sqrt(A * A * (units::constants::nucleonMass * units::constants::nucleonMass) +
@@ -52,7 +51,6 @@ void InteractionHistogram::fill(particles::Code projectile_id,
 
     interaction_histogram_cms_(
         static_cast<corsika::particles::CodeIntType>(projectile_id), sqrtS / 1_GeV);
-
     interaction_histogram_lab_(
         static_cast<corsika::particles::CodeIntType>(projectile_id), lab_energy / 1_GeV);
   }
@@ -82,6 +80,9 @@ InteractionHistogram InteractionHistogram::operator+(InteractionHistogram other)
   other.nuclear_inthist_lab_ += nuclear_inthist_lab_;
   other.nuclear_inthist_cms_ += nuclear_inthist_cms_;
 
+  other.interaction_histogram_lab_ += interaction_histogram_lab_;
+  other.interaction_histogram_cms_ += interaction_histogram_cms_;
+
   return other;
 }