Jpp test-rotations-new
the software that should make you happy
Loading...
Searching...
No Matches
JFitVolume.sh
Go to the documentation of this file.
1#!/usr/bin/env zsh
2#
3# \author lquinn
4#
5version=1.0
6script=${0##*/}
7
8# ------------------------------------------------------------------------------------------
9#
10# Utility script to test JVolume1D and JFitEfficiency
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 1
25set_variable WORKDIR ./
26set_variable TRIGGER_FILE $JPP_DATA/trigger_efficiency+background.root
27set_variable INPUT_FILE $JPP_DATA/start.root
28set_variable DETECTOR $JPP_DATA/km3net_reference.detx
29set_variable: FORMAT GRAPHICS_FORMAT gif
30set_variable+ BATCH GRAPHICS_BATCH -B
31
32if do_usage $*; then
33 usage "$script [detector file [trigger file [input file]]]"
34fi
35
36case $# in
37 3) set_variable INPUT_FILE $3;&
38 2) set_variable TRIGGER_FILE $2;&
39 1) set_variable DETECTOR $1;;
40esac
41
42if (( 1 )); then
43
44 JPrintTree -f $TRIGGER_FILE
45
46 echo "Check that the Monte Carlo and DAQ trees have different numbers of events. Else run JTriggerEfficiency without the -O option.\n"
47
48fi
49
50if (( 1 )); then
51
52 print_variable TRIGGER_FILE INPUT_FILE DETECTOR
53
54 $JPP_DIR/examples/JTrigger/JVolume1D \
55 -f $TRIGGER_FILE \
56 -o "volume[trigger].root" \
57 -X \
58 -d $DEBUG
59
60 $JPP_DIR/examples/JFit/JFitEfficiency \
61 -f $INPUT_FILE \
62 -o efficiency.root \
63 -a $DETECTOR \
64 -X -C \
65 -W "2 0" \
66 -Q 0.5 \
67 -d $DEBUG
68fi
69
70JOpera1D \
71 -f "volume[trigger].root:hV" \
72 -f efficiency.root:he \
73 -o "volume[reco].root" \
74 -u "Multiply" \
75 -d $DEBUG
76
77
78JPlot1D \
79 -f "volume[reco].root:Multiply" \
80 -f "volume[trigger].root:hV" \
81 -XX \
82 -> "E [GeV]" \
83 -\^ "V [km^{3}]" \
84 -LL \
85 -o volume.$FORMAT $BATCH