Jpp test-rotations-new
the software that should make you happy
Loading...
Searching...
No Matches
plot-Sirene.sh
Go to the documentation of this file.
1#!/usr/bin/env zsh
2#
3# \author mdejong
4#
5version=1.0
6script=${0##*/}
7
8# ------------------------------------------------------------------------------------------
9#
10# Utility script to plot JSirene results.
11#
12# ------------------------------------------------------------------------------------------
13
14if [ -z $JPP_DIR ]; then
15 echo "Variable JPP_DIR undefined."
16 exit
17fi
18
19source $JPP_DIR/setenv.sh $JPP_DIR
20
21zmodload zsh/mathfunc
22
23set_variable: DEBUG SIRENE_DEBUG 2
24set_variable: WORKDIR SIRENE_WORKDIR ./
25set_variable: FORMAT GRAPHICS_FORMAT gif
26set_variable+ BATCH GRAPHICS_BATCH -B
27
28if do_usage $*; then
29 usage "$script (input file)+"
30fi
31
32if (( $# == 0 )); then
33 fatal "No input files."
34fi
35
36INPUT_FILES=($argv[*])
37
38typeset -A RADIATION
39
40RADIATION[EErad]="10 19"
41RADIATION[Brems]="20 29"
42RADIATION[GNrad]="30 39"
43RADIATION[DIS]="100 101"
44
45for KEY VALUE in ${(@kv)RADIATION}; do
46
47 JPlot1D \
48 -f${^INPUT_FILES}:rad \
49 -XX -PX \
50 -y "$VALUE" -Y \
51 -> "E [GeV]" \
52 -\^ "Es [GeV]" \
53 -T "$KEY" \
54 -O "][" \
55 -o $WORKDIR/rad_${KEY}.$FORMAT $BATCH
56
57done
58
59RADIATION[MCS]="0 1"
60
61for KEY VALUE in ${(@kv)RADIATION}; do
62
63 JPlot1D \
64 -f${^INPUT_FILES}:rms \
65 -XX -PX \
66 -y "$VALUE" -Y \
67 -z "1e-9 1e-2" \
68 -> "E [GeV]" \
69 -\^ "#theta^{2} [rad^{2}]" \
70 -T "$KEY" \
71 -O "][" \
72 -o $WORKDIR/rad_${KEY}.$FORMAT --! $BATCH
73
74done
75
76JPlot1D \
77 -f${^INPUT_FILES}:cpu \
78 -XX -Y \
79 -> "E [GeV]" \
80 -\^ "cpu [ms]" \
81 -T "" \
82 -o $WORKDIR/cpu.$FORMAT $BATCH