Jpp
JORCAShowerPrefit.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 #
3 # \author adomi
4 #
5 version=1.0
6 script=${0##*/}
7 
8 # ------------------------------------------------------------------------------------------
9 #
10 # Utility script to test JORCAShowerPrefit.
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 set_variable DEBUG 2
24 set_variable WORKDIR $JPP_DATA
25 set_variable INPUT_FILE $WORKDIR/trigger_efficiency+background.root
26 set_variable OUTPUT_FILE $WORKDIR/shower_prefit.root
27 set_variable DETECTOR $WORKDIR/km3net_reference.detx
28 
29 if ( do_usage $* ); then
30  usage "$script [detector file [input file [output file]]]"
31 fi
32 
33 case $# in
34  3) OUTPUT_FILE=$3;&
35  2) INPUT_FILE=$2;&
36  1) DETECTOR=$1;;
37 esac
38 
39 print_variable DETECTOR INPUT_FILE OUTPUT_FILE
40 check_input_file $DETECTOR $INPUT_FILE
41 check_output_file $OUTPUT_FILE
42 
43 if (( 1 )); then
44 
45 # default input values
46 
47  set_variable TEXTRAMAX_NS 20
48  set_variable DMAX_M 50
49  set_variable SIGMA_NS 3
50  set_variable OUTLIERS 4
51  set_variable GRID_SIZE_M "-20 20"
52  set_variable GRID_STEP_M 20
53  set_variable GRID_SIZE_NS "-50 50"
54  set_variable GRID_STEP_NS 25
55 
56  timer_start
57 
58  JORCAShowerPrefit \
59  -a $DETECTOR \
60  -f $INPUT_FILE \
61  -o $OUTPUT_FILE \
62  -@ "prefit.Tmax_ns=$TEXTRAMAX_NS \
63  prefit.Dmax_m=$DMAX_M \
64  prefit.pos_grid=$GRID_SIZE_M \
65  prefit.pos_step=$GRID_STEP_M \
66  prefit.time_grid=$GRID_SIZE_NS \
67  prefit.time_step=$GRID_STEP_NS \
68  prefit.sigma_ns=$SIGMA_NS \
69  prefit.useL0=$USE_L0 \
70  prefit.numberOfOutliers=$OUTLIERS" \
71  -d ${DEBUG} --!
72 
73  timer_stop
74  timer_print
75 
76 fi