Jpp  19.1.0
the software that should make you happy
JMuonStart.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 JMuonStart on Monte Carlo true muons.
11 #
12 # ------------------------------------------------------------------------------------------
13 
14 if [ -z $JPP_DIR ]; then
15  echo "Variable JPP_DIR undefined."
16  exit
17 fi
18 
19 source $JPP_DIR/setenv.sh $JPP_DIR
20 
21 set_variable DEBUG 2
22 set_variable WORKDIR $JPP_DATA
23 set_variable INPUT_FILE $WORKDIR/trigger_efficiency+background.root
24 set_variable INTERMEDIATE_FILE $WORKDIR/evt-mc.root
25 set_variable OUTPUT_FILE $WORKDIR/start-mc.root
26 set_variable DETECTOR $WORKDIR/km3net_reference.detx
27 set_variable PDF $WORKDIR/J%p.dat
28 set_variable: FORMAT GRAPHICS_FORMAT gif
29 set_variable+ BATCH GRAPHICS_BATCH -B
30 
31 if do_usage $*; then
32  usage "$script [detector file [input file [output file [PDF file descriptor]]]]"
33 fi
34 
35 case $# in
36  4) set_variable PDF $4;&
37  3) set_variable OUTPUT_FILE $3;&
38  2) set_variable INPUT_FILE $2;&
39  1) set_variable DETECTOR $1;;
40 esac
41 
42 if ( ! reuse_file $INTERMEDIATE_FILE ); then
43 
44  JMuonMCEvt \
45  -f $INPUT_FILE \
46  -o $INTERMEDIATE_FILE \
47  -d ${DEBUG} --!
48 fi
49 
50 if ( ! reuse_file $OUTPUT_FILE ); then
51 
52  PARAMETERS=(
53  roadWidth_m 50.0
54  R_Hz 7.0e3
55  TMin_ns -25.0
56  TMax_ns +25.0
57  Pmin1 1.0e-3
58  Pmin2 1.0e-2
59  numberOfPrefits 1)
60 
61  eval JMuonStart \
62  -a $DETECTOR \
63  -f $INTERMEDIATE_FILE \
64  -o $OUTPUT_FILE \
65  -P $PDF \
66  -@ `make_equation PARAMETERS` \
67  -d $DEBUG --!
68 fi
69 
70 JMuonPostfit \
71  -f $OUTPUT_FILE \
72  -o $WORKDIR/postfit.root \
73  -N 1 \
74  -A 5
75 
76 JPlot1D \
77  -f $WORKDIR/postfit.root:hz0 \
78  -f $WORKDIR/postfit.root:hz1 \
79  -> "#Deltaz [m]" \
80  -\^ "number of events [a.u.]" \
81  -L TR -T "" \
82  -o $WORKDIR/z.$FORMAT $BATCH
83