Jpp test-rotations-new
the software that should make you happy
Loading...
Searching...
No Matches
JMultiplicityK40.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 make coincidence rate plots.
11#
12# ------------------------------------------------------------------------------------------
13
14
15if [ -z $JPP_DIR ]; then
16 echo "Variable JPP_DIR undefined."
17 exit
18fi
19
20
21source $JPP_DIR/setenv.sh $JPP_DIR
22
23
24set_variable DIR $JPP_DIR/examples/JPhysics/
25set_variable DEBUG 2
26set_variable WORKDIR ${TMPDIR:-/tmp}/
27set_variable: FORMAT GRAPHICS_FORMAT gif
28set_variable+ BATCH GRAPHICS_BATCH -B
29
30if do_usage $*; then
31 usage "$script [working directory]"
32fi
33
34case $# in
35 1) set_variable WORKDIR $1;;
36esac
37
38if [[ ! -f $WORKDIR/k40\[KM3NeT\].root ]]; then
39
40 $DIR/JMultiplicityK40 \
41 -n 100000000 \
42 -F 4 \
43 -G 2 \
44 -D "0.216 5" \
45 -A 1 \
46 -U \
47 -o $WORKDIR/k40\[KM3NeT\].root
48fi
49
50for (( a = 0.0; $a <= 1.0; a += 0.5 )); do
51
52 set_variable LABEL `printf '%3.1f' $a`
53
54 if [[ ! -f $WORKDIR/k40\[${LABEL}\].root ]]; then
55
56 $DIR/JMultiplicityK40 \
57 -n 100000000 \
58 -F 4 \
59 -G 2 \
60 -D "0.216 5" \
61 -A 3 -a $a \
62 -U \
63 -o $WORKDIR/k40\[${LABEL}\].root
64 fi
65done
66
67INPUT_FILES=(`ls $WORKDIR/k40\[*\].root`)
68
69echo $INPUT_FILES[*]
70
71JPlot1D \
72 -f "${INPUT_FILES//\.root/.root:pmt}" \
73 -> "cos(#theta)" \
74 -\^ "acceptance" \
75 -y "0 1.6" \
76 -L TR \
77 -T "PMT" \
78 -o pmt.$FORMAT $BATCH
79
80M=2
81
82JPlot1D \
83 -f "${INPUT_FILES//\.root/.root:M.${M}}" \
84 -> "D [m]" \
85 -\^ "rate [Hz]" \
86 -x "0 5" \
87 -y "1 1e2" -Y \
88 -L TR \
89 -T "M = ${M}" \
90 -o k40.$FORMAT $BATCH