IAP GITLAB

Skip to content
Snippets Groups Projects
Commit 0d53d8d0 authored by ralfulrich's avatar ralfulrich
Browse files

make track plot also with gammas and electrons

parent aa62eb17
No related branches found
No related tags found
No related merge requests found
...@@ -189,7 +189,7 @@ int main(int argc, char** argv) { ...@@ -189,7 +189,7 @@ int main(int argc, char** argv) {
// setup processes, decays and interactions // 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::Interaction proposal(env, cut.GetECut());
process::proposal::ContinuousProcess em_continuous(env, cut.GetECut()); process::proposal::ContinuousProcess em_continuous(env, cut.GetECut());
process::interaction_counter::InteractionCounter proposalCounted(proposal); process::interaction_counter::InteractionCounter proposalCounted(proposal);
......
...@@ -2,8 +2,6 @@ ...@@ -2,8 +2,6 @@
# (c) Copyright 2018 CORSIKA Project, corsika-project@lists.kit.edu # (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 # 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 # Licence version 3 (GPL Version 3). See file LICENSE for a full version of
# the license. # the license.
...@@ -19,6 +17,8 @@ fi ...@@ -19,6 +17,8 @@ fi
directory=$(dirname "$tracks_dat") directory=$(dirname "$tracks_dat")
hadrons_dat="$directory/hadrons.dat" hadrons_dat="$directory/hadrons.dat"
muons_dat="$directory/muons.dat" muons_dat="$directory/muons.dat"
electrons_dat="$directory/electrons.dat"
gammas_dat="$directory/gammas.dat"
if [ ! -e "$tracks_dat" ]; then if [ ! -e "$tracks_dat" ]; then
echo "$tracks_dat does not exist" >&2 echo "$tracks_dat does not exist" >&2
...@@ -28,8 +28,10 @@ fi ...@@ -28,8 +28,10 @@ fi
if [ "$muons_dat" -ot "$tracks_dat" ] || [ "$hadrons_dat" -ot "$tracks_dat" ]; then if [ "$muons_dat" -ot "$tracks_dat" ] || [ "$hadrons_dat" -ot "$tracks_dat" ]; then
echo "splitting $tracks_dat into $muons_dat and $hadrons_dat." echo "splitting $tracks_dat into $muons_dat and $hadrons_dat."
cat "$tracks_dat" | egrep '^\s+-*13\s' > "$muons_dat" cat "$tracks_dat" | egrep '^\s+-*13\s' > "$muons_dat"
cat "$tracks_dat" | egrep -v '^\s+-*13\s' > "$hadrons_dat" cat "$tracks_dat" | egrep '^\s+-*11\s' > "$electrons_dat"
fi 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" output="$2"
if [ -z "$output" ]; then if [ -z "$output" ]; then
...@@ -55,7 +57,7 @@ do for [t=0:359:1] { ...@@ -55,7 +57,7 @@ do for [t=0:359:1] {
# set output sprintf("%03d_$output", t) # set output sprintf("%03d_$output", t)
set view 80, 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 EOF
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment