From 0d53d8d06f7a8e0793edfa6b4bdfd7b0081fc162 Mon Sep 17 00:00:00 2001
From: ralfulrich <ralf.ulrich@kit.edu>
Date: Fri, 9 Oct 2020 11:18:12 +0200
Subject: [PATCH] make track plot also with gammas and electrons
---
Documentation/Examples/vertical_EAS.cc | 2 +-
Tools/plot_tracks.sh | 12 +++++++-----
2 files changed, 8 insertions(+), 6 deletions(-)
diff --git a/Documentation/Examples/vertical_EAS.cc b/Documentation/Examples/vertical_EAS.cc
index 5cf491d29..f9b44e638 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 c8e7fbcd4..a600b3557 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
--
GitLab