diff --git a/examples/boundary_example.cpp b/examples/boundary_example.cpp
index d9522ce15535b4edf56c3f99c4f1857862d1ceff..e5f707bdef84d99c6a553ad32b58ff0e949f8b86 100644
--- a/examples/boundary_example.cpp
+++ b/examples/boundary_example.cpp
@@ -164,7 +164,9 @@ int main() {
   // define air shower object, run simulation
   Cascade EAS(env, tracking, sequence, output, stack);
 
+  output.startOfShower();
   EAS.run();
+  output.endOfShower();
 
   CORSIKA_LOG_INFO("Result: E0={}GeV", E0 / 1_GeV);
   cut.showResults();
diff --git a/examples/cascade_example.cpp b/examples/cascade_example.cpp
index 2852428d2c074b0afaa22fcdc58ec5eeeaa0825c..623f443120ec69c52e5bacbecf6dfc02039e87e9 100644
--- a/examples/cascade_example.cpp
+++ b/examples/cascade_example.cpp
@@ -153,7 +153,9 @@ int main() {
   // define air shower object, run simulation
   Cascade EAS(env, tracking, sequence, output, stack);
 
+  output.startOfShower();
   EAS.run();
+  output.endOfShower();
 
   eLoss.printProfile(); // print longitudinal profile
 
diff --git a/examples/cascade_proton_example.cpp b/examples/cascade_proton_example.cpp
index e6ab7d76fff01d3b0081e7e5a1bb784a01fb59a3..5491258ac68006e5e7d483878543b85a3a1d2de8 100644
--- a/examples/cascade_proton_example.cpp
+++ b/examples/cascade_proton_example.cpp
@@ -146,7 +146,9 @@ int main() {
 
   // define air shower object, run simulation
   Cascade EAS(env, tracking, sequence, output, stack);
+  output.startOfShower();
   EAS.run();
+  output.endOfShower();
 
   cout << "Result: E0=" << E0 / 1_GeV << endl;
   cut.showResults();
diff --git a/examples/em_shower.cpp b/examples/em_shower.cpp
index 5dc3e53101a904556c4fa57b1646fd8d10e630e6..9e8d2a39e6740bf8e6403f75be20b9ff80b9741b 100644
--- a/examples/em_shower.cpp
+++ b/examples/em_shower.cpp
@@ -170,7 +170,9 @@ int main(int argc, char** argv) {
   //  EAS.setNodes();
   //  EAS.forceInteraction();
 
+  output.startOfShower();
   EAS.run();
+  output.endOfShower();
 
   cut.showResults();
   emContinuous.showResults();
diff --git a/examples/hybrid_MC.cpp b/examples/hybrid_MC.cpp
index 3a4f6757708a4a2ff081d8240b2bdf537a90a6f6..dc4f7e8ee784fa1e5188785a38fa3b6dbc3d4be1 100644
--- a/examples/hybrid_MC.cpp
+++ b/examples/hybrid_MC.cpp
@@ -250,7 +250,9 @@ int main(int argc, char** argv) {
   //  EAS.SetNodes();
   //  EAS.forceInteraction();
 
+  output.startOfShower();
   EAS.run();
+  output.endOfShower();
 
   cut.showResults();
   eLoss.showResults();
diff --git a/examples/vertical_EAS.cpp b/examples/vertical_EAS.cpp
index 2fc5ad8f0d855c5ed8e31c9fb38006d801c7bb21..b8e16ab432e45620a2007888c20a4a633848e95f 100644
--- a/examples/vertical_EAS.cpp
+++ b/examples/vertical_EAS.cpp
@@ -295,7 +295,9 @@ int main(int argc, char** argv) {
   // define air shower object, run simulation
   setup::Tracking tracking;
   Cascade EAS(env, tracking, sequence, output, stack);
+  output.startOfShower();
   EAS.run();
+  output.endOfShower();
 
   cut.showResults();
   // emContinuous.showResults();