Jpp
JTriggerEfficiency1D.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 test JEfficiency.
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 DEBUG 2
25 set_variable WORKDIR ./
26 set_variable INPUT_FILE $JPP_DATA/sirene.root
27 set_variable DETECTOR $JPP_DATA/km3net_reference.detx
28 set_variable TRIGGER $JPP_DATA/trigger_parameters_arca.txt
29 set_variable BACKGROUND_HZ 0
30 
31 
32 if ( do_usage $* ); then
33  usage "$script [input file [working directory]]"
34 fi
35 
36 case $# in
37  2) WORKDIR=$2;&
38  1) INPUT_FILE=$1;;
39 esac
40 
41 
42 set_variable OUTPUT_FILE $WORKDIR/trigger_efficiency.root
43 
44 if (( 1 )); then
45 
46  if [[ -f $OUTPUT_FILE ]] then
47  rm -i $OUTPUT_FILE
48  fi
49 
50  if [[ ! -f $OUTPUT_FILE ]] then
51 
52  JTriggerEfficiency \
53  -a ${DETECTOR} \
54  -f ${INPUT_FILE} \
55  -o ${OUTPUT_FILE} \
56  -@ ${TRIGGER} \
57  -d ${DEBUG} \
58  -B "$BACKGROUND_HZ" \
59  --!
60  fi
61 fi
62 
63 set_variable INPUT_FILE $OUTPUT_FILE
64 set_variable OUTPUT_FILE $WORKDIR/efficiency.root
65 
66 if (( 1 )); then
67 
68  $JPP_DIR/examples/JTrigger/JTriggerEfficiency1D \
69  -f $INPUT_FILE \
70  -o $OUTPUT_FILE \
71  -d $DEBUG
72 
73 fi
74 
75 
76 JPlot1D \
77  -f ${OUTPUT_FILE}:h1 \
78  -> "number of hits" \
79  -\^ "number of events" \
80  -o counts.gif -T ""
81 
82 
83 JPlot1D \
84  -f ${OUTPUT_FILE}:h2 \
85  -y "0 1.1" \
86  -> "number of hits" \
87  -\^ "effieciency" \
88  -o efficiency.gif -T ""