Jpp
JTriggerEfficiency.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 #
3 # \author mdejong
4 #
5 script=${0##*/}
6 
7 # ------------------------------------------------------------------------------------------
8 #
9 # Example script to run JTriggerEfficiency.
10 #
11 # ------------------------------------------------------------------------------------------
12 
13 if [ -z $JPP_DIR ]; then
14  echo "Variable JPP_DIR undefined."
15  exit
16 fi
17 
18 source $JPP_DIR/setenv.sh $JPP_DIR
19 
20 set_variable DEBUG ${TRIGGEREFFICIENCY_DEBUG:-1}
21 set_variable WORKDIR ./
22 set_variable INPUT_FILE $WORKDIR/sirene.root
23 set_variable OUTPUT_FILE $WORKDIR/trigger_efficiency+background.root
24 set_variable DETECTORA $JPP_DATA/km3net_reference.detx
25 set_variable DETECTORB $JPP_DATA/km3net_reference.detx
26 set_variable TRIGGER $JPP_DATA/trigger_parameters_arca.txt
27 set_variable PMT_FILE $JPP_DATA/PMT_parameters.txt
28 set_variable BACKGROUND_HZ ${TRIGGEREFFICIENCY_BACKGROUND_HZ:-${BACKGROUND_HZ:-5e3 500 50 5 0.5}}
29 set_variable SEED ${TRIGGEREFFICIENCY_SEED:-0}
30 set_variable FACTOR ${TRIGGEREFFICIENCY_FACTOR:-${FACTOR:-1.0}}
31 set_variable OPTION ${TRIGGEREFFICIENCY_TRIGGERED_EVENTS_ONLY+-O}
32 
33 if ( do_usage $* ); then
34  usage "$script [detector file [detector file [(input file)+ [output file [trigger file [PMT parameters file]]]]]]"\
35  "\nNote that if more than one input file is specified, all other arguments must be provided."
36 fi
37 
38 case $# in
39  6) set_variable PMT_FILE $6;&
40  5) set_variable TRIGGER $5;&
41  4) set_variable OUTPUT_FILE $4;&
42  3) set_variable INPUT_FILE $3;&
43  2) set_variable DETECTORB $2;&
44  1) set_variable DETECTORA $1;;
45  0) ;;
46  *) set_variable DETECTORA $argv[1]
47  set_variable DETECTORB $argv[2]
48  set_variable INPUT_FILE $argv[3,-4]
49  set_variable OUTPUT_FILE $argv[-3]
50  set_variable TRIGGER $argv[-2]
51  set_variable PMT_FILE $argv[-1];;
52 esac
53 
54 print_variable DETECTORA DETECTORB TRIGGER INPUT_FILE OUTPUT_FILE BACKGROUND_HZ PMT_FILE
55 check_input_file $DETECTORA $DETECTORB $TRIGGER $INPUT_FILE $PMT_FILE
56 check_output_file $OUTPUT_FILE
57 
58 timer_start
59 
60 JTriggerEfficiency \
61  -a ${DETECTORA} \
62  -b ${DETECTORB} \
63  -f ${INPUT_FILE} \
64  -o ${OUTPUT_FILE} \
65  -P ${PMT_FILE} \
66  -P "QE=${FACTOR}" \
67  ${OPTION} \
68  -@ ${TRIGGER} \
69  -d ${DEBUG} \
70  -S ${SEED} \
71  -B "$BACKGROUND_HZ" \
72  --!
73 
74 timer_stop
75 timer_print
76 
77 
78 JPrintMeta \
79  -f ${OUTPUT_FILE}
80 
81 JPrintTree \
82  -f ${OUTPUT_FILE}
83 
84 JTriggerMonitor \
85  -o $WORKDIR/monitor.root \
86  -f ${OUTPUT_FILE}