diff --git a/Documentation/Examples/vertical_EAS.cc b/Documentation/Examples/vertical_EAS.cc
index 5cf491d29abc901763189aeab63c1adfe9875734..f9b44e638291abd20119ee1f5be11623d695d65d 100644
--- a/Documentation/Examples/vertical_EAS.cc
+++ b/Documentation/Examples/vertical_EAS.cc
@@ -189,7 +189,7 @@ int main(int argc, char** argv) {
 
   // setup processes, decays and interactions
 
-  process::particle_cut::ParticleCut cut{60_GeV, true, true};
+  process::particle_cut::ParticleCut cut{60_GeV, false, true};
   process::proposal::Interaction proposal(env, cut.GetECut());
   process::proposal::ContinuousProcess em_continuous(env, cut.GetECut());
   process::interaction_counter::InteractionCounter proposalCounted(proposal);
diff --git a/Tools/plot_tracks.sh b/Tools/plot_tracks.sh
index c8e7fbcd4263146968793b485ea2a411f1aa096b..a600b35573ccdf6654d32f531f2e8a94a4f4fa5a 100755
--- a/Tools/plot_tracks.sh
+++ b/Tools/plot_tracks.sh
@@ -2,8 +2,6 @@
 
 # (c) Copyright 2018 CORSIKA Project, corsika-project@lists.kit.edu
 #
-# See file AUTHORS for a list of contributors.
-#
 # This software is distributed under the terms of the GNU General Public
 # Licence version 3 (GPL Version 3). See file LICENSE for a full version of
 # the license.
@@ -19,6 +17,8 @@ fi
 directory=$(dirname "$tracks_dat")
 hadrons_dat="$directory/hadrons.dat"
 muons_dat="$directory/muons.dat"
+electrons_dat="$directory/electrons.dat"
+gammas_dat="$directory/gammas.dat"
 
 if [ ! -e "$tracks_dat" ]; then
   echo "$tracks_dat does not exist" >&2
@@ -28,8 +28,10 @@ fi
 if [ "$muons_dat" -ot "$tracks_dat" ] || [ "$hadrons_dat" -ot "$tracks_dat" ]; then
   echo "splitting $tracks_dat into $muons_dat and $hadrons_dat."
   cat "$tracks_dat" | egrep '^\s+-*13\s' > "$muons_dat"
-  cat "$tracks_dat" | egrep -v '^\s+-*13\s' > "$hadrons_dat"
-fi
+  cat "$tracks_dat" | egrep '^\s+-*11\s' > "$electrons_dat"
+  cat "$tracks_dat" | egrep '^\s+-*22\s' > "$gammas_dat"
+  cat "$tracks_dat" | egrep -v '^\s+-*13\s' | egrep -v '^\s+-*11\s' | egrep -v '^\s+-*22\s' > "$hadrons_dat"
+#fi
 
 output="$2"
 if [ -z "$output" ]; then
@@ -55,7 +57,7 @@ do for [t=0:359:1] {
 #	set output sprintf("%03d_$output", t)
 
 	set view 80, t
-        splot "$muons_dat" u 3:4:5:6:7:8 w vectors nohead lt rgb "red" t "", "$hadrons_dat" u 3:4:5:6:7:8 w vectors nohead  lc rgb "black" t ""
+        splot "$gammas_dat" u 3:4:5:6:7:8 w vectors nohead lt rgb "orange" t "", "$electrons_dat" u 3:4:5:6:7:8 w vectors nohead lt rgb "blue" t "", "$muons_dat" u 3:4:5:6:7:8 w vectors nohead lt rgb "red" t "", "$hadrons_dat" u 3:4:5:6:7:8 w vectors nohead  lc rgb "black" t ""
 }
 EOF