Jpp 19.3.0-rc.3
the software that should make you happy
Loading...
Searching...
No Matches
JFilter.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 test trigger algorithms.
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 DEBUG 2
25set_variable WORKDIR $JPP_DATA
26set_variable INPUT_FILE $WORKDIR/trigger_efficiency+background.root
27set_variable OUTPUT_FILE $WORKDIR/filter\%.root
28set_variable DETECTOR $WORKDIR/km3net_reference.detx
29set_variable HISTOGRAM 1
30
31
32if do_usage $*; then
33 usage "$script [detector file [input file [output file]]]"
34fi
35
36case $# in
37 3) set_variable OUTPUT_FILE $3;&
38 2) set_variable INPUT_FILE $2;&
39 1) set_variable DETECTOR $1;;
40esac
41
42
43if (( 1 )); then
44
45 check_input_file $DETECTOR $INPUT_FILE
46
47 for OPTION in A B C D E F; do
48
49 FILE="${OUTPUT_FILE/\%/[${OPTION}]}"
50
51 print_variable DETECTOR INPUT_FILE FILE
52
53 $JPP_DIR/examples/JTrigger/JFilter \
54 -a $DETECTOR \
55 -f $INPUT_FILE \
56 -o $FILE \
57 -c $OPTION \
58 -T 20 \
59 -R 200 \
60 -H $HISTOGRAM \
61 -d $DEBUG
62
63 done
64
65fi
66
67
68if (( 1 )); then
69
70 X_LABEL[1]="Number of hits"
71 X_LABEL[2]="E_{#nu} [GeV]"
72
73 X_OPTION[1]=""
74 X_OPTION[2]="-XX"
75
76 X_RANGE[1]="-x 1 50"
77 X_RANGE[2]="-x 2 7"
78
79 FILE_E=""
80 FILE_P=""
81
82 for OPTION in A B C D E F; do
83 FILE_E+=" ${OUTPUT_FILE/\%/[${OPTION}]}:he"
84 FILE_P+=" ${OUTPUT_FILE/\%/[${OPTION}]}:hp"
85 done
86
87 FILE_E=${FILE_E## }
88 FILE_P=${FILE_P## }
89
90 JPlot1D \
91 -f $FILE_E \
92 -y "0 1.1" \
93 -\^ "efficiency" \
94 -> $X_LABEL[$HISTOGRAM] \
95 $X_OPTION[$HISTOGRAM] \
96 $X_RANGE[$HISTOGRAM] \
97 -L TR
98
99
100 JPlot1D \
101 -f $FILE_P \
102 -y "0 1.1" \
103 -\^ "purity" \
104 -> $X_LABEL[$HISTOGRAM] \
105 $X_OPTION[$HISTOGRAM] \
106 $X_RANGE[$HISTOGRAM] \
107 -L TR
108
109fi