Jpp  master_rocky
the software that should make you happy
JMultiplicityK40.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 #
3 # \author mdejong
4 #
5 version=1.0
6 script=${0##*/}
7 
8 # ------------------------------------------------------------------------------------------
9 #
10 # Utility script to make coincidence rate plots.
11 #
12 # ------------------------------------------------------------------------------------------
13 
14 
15 if [ -z $JPP_DIR ]; then
16  echo "Variable JPP_DIR undefined."
17  exit
18 fi
19 
20 
21 source $JPP_DIR/setenv.sh $JPP_DIR
22 
23 
24 set_variable DIR $JPP_DIR/examples/JPhysics/
25 set_variable DEBUG 2
26 set_variable WORKDIR ${TMPDIR:-/tmp}/
27 set_variable: FORMAT GRAPHICS_FORMAT gif
28 set_variable+ BATCH GRAPHICS_BATCH -B
29 
30 if do_usage $*; then
31  usage "$script [working directory]"
32 fi
33 
34 case $# in
35  1) set_variable WORKDIR $1;;
36 esac
37 
38 if [[ ! -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
48 fi
49 
50 for (( 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
65 done
66 
67 INPUT_FILES=(`ls $WORKDIR/k40\[*\].root`)
68 
69 echo $INPUT_FILES[*]
70 
71 JPlot1D \
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 
80 M=2
81 
82 JPlot1D \
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