Jpp
Calibrate_local.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 #
3 # \author rgruiz
4 #
5 
6 version=1.0
7 script=${0##*/}
8 
9 source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
10 
11 if ( do_usage $* ); then
12  usage "$script [edit the different options within]"
13 fi
14 
15 if (( $# < 5 )); then
16  fatal "Incorrect number of arguments."
17 fi
18 
19 set_variable DAQ_FILE $1
20 set_variable DETECTOR_FILE $2
21 set_variable INSTRUMENT $3
22 set_variable STRING $4
23 set_variable OUTPUT_DIR $5
24 
25 set_variable TOP_PMTS 1
26 set_variable BOTTOM_PMTS 4
27 set_variable MAX_NEIGHBORS 4
28 set_variable ANALYSIS_LEVEL 0
29 
30 set_variable PULSE_PERIOD 6250 #in units of 16 ns
31 set_variable STAGGER_PERIOD 220 #in units of 16 ns
32 set_variable infile $DAQ_FILE
33 
34 cd ${JPP_DIR}/examples/JNanobeacon
35 
36 set_variable OUTDIR_PULSE "${OUTPUT_DIR}/out/$INSTRUMENT/JPulseFinder"
37 
38 if [ ! -d ${OUTDIR_PULSE} ]
39 then
40  mkdir -p ${OUTDIR_PULSE}
41 fi
42 
43 set_variable OUT_FILE_PULSE ${OUTDIR_PULSE}/$(basename ${DAQ_FILE} .root)_JPF.root
44 
45 if [ ! -f $OUT_FILE_PULSE ]; then
46 
47  JPulseFinder \
48  -f $DAQ_FILE \
49  -a $DETECTOR_FILE \
50  -s $STRING \
51  -p $PULSE_PERIOD \
52  -S $STAGGER_PERIOD \
53  -o $OUT_FILE_PULSE
54 else
55  echo "\nThe pulses for the file ${DAQ_FILE} already exist in file\n$OUT_FILE_PULSE\nContinue..."
56 fi
57 
58 set_variable OUTDIR_CALIBRATE "${OUTPUT_DIR}/out/$INSTRUMENT/JInterDomCal"
59 
60 if [ ! -d ${OUTDIR_CALIBRATE} ]
61 then
62  mkdir -p ${OUTDIR_CALIBRATE}
63 fi
64 
65 set_variable OUT_FILE_TXT ${OUTDIR_CALIBRATE}/$(basename ${DAQ_FILE} .root)_JNBC.txt
66 set_variable OUT_FILE_CHECKS ${OUTDIR_CALIBRATE}/$(basename ${DAQ_FILE} .root)_JNBC.root
67 set_variable OUT_FILE_DETX ${OUTDIR_CALIBRATE}/$(basename ${DAQ_FILE} .root)_JNBC.detx
68 
69 JInterDomCal \
70  -a $DETECTOR_FILE \
71  -f $OUT_FILE_PULSE \
72  -s $STRING \
73  -t $OUT_FILE_TXT \
74  -x $OUT_FILE_DETX \
75  -c $OUT_FILE_CHECKS \
76  -u $TOP_PMTS \
77  -d $BOTTOM_PMTS \
78  -n $MAX_NEIGHBORS \
79  -l $ANALYSIS_LEVEL
80 
81 report_name="$/calibration_report_${runs[i]}.pdf"
82 
83 set_variable OUT_FILE_PDF ${OUTDIR_CALIBRATE}/$(basename ${DAQ_FILE} .root).pdf
84 
85 echo "Writing calibration report in file:"
86 echo ${OUT_FILE_PDF}
87 echo "------------------------------------\n"
88 
89 python plot_report_calibration.py -ir $OUT_FILE_CHECKS -it $OUT_FILE_TXT -o $OUT_FILE_PDF