Jpp test-rotations-new
the software that should make you happy
Loading...
Searching...
No Matches
JMergeRunAnalyzer.sh
Go to the documentation of this file.
1#!/usr/bin/env zsh
2
3version=1.0
4script=${0##*/}
5
6source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
7
8if do_usage $*; then
9 usage "$script <input_files> <output_file(.root)> [only_directory_path]"
10fi
11
12set_variable: ONLY_DIR INPUT_DIR 0
13set_variable: WORKDIR WORKDIR `pwd`
14set_variable: FORMAT GRAPHICS_FORMAT pdf
15
16case $# in
17 3) set_variable ONLY_DIR $3;&
18 2) set_variable OUT_FILE $2
19 set_variable JRA_FILES $1;;
20 *) fatal "Wrong number of arguments"
21esac
22
23set_variable INPUT_FILES ${JRA_FILES}
24
25if [[ ${ONLY_DIR} != 0 ]]
26then
27set_variable INPUT_FILES `echo ${JRA_FILES}/*`
28fi
29
30set_variable OUT_PATH `readlink -f $(dirname ${OUT_FILE})`
31set_variable JRA_DIR ${JPP_DIR}/software/JRunAnalyzer
32set_variable FILEPREFIX "plot"
33
34if [ ! -d ${OUT_PATH} ]
35then
36 mkdir -p ${OUT_PATH}
37fi
38
39time ${JRA_DIR}/JMergeRunAnalyzer \
40 -f $INPUT_FILES \
41 -o $OUT_FILE --!
42
43#****** Create Profiling and Projection of the TH2Ds *******#
44
45JProfile2D \
46 -f ${OUT_FILE}:h_pmt_rate_distribution \
47 -f ${OUT_FILE}:h_Triggered_hits \
48 -f ${OUT_FILE}:h_Snapshot_hits \
49 -f ${OUT_FILE}:h_Number_of_overlays \
50 -o px_${OUT_FILE} \
51 -P "X" \
52 -O "g" --!
53
54JProject2D \
55 -f ${OUT_FILE}:h_Triggered_over_Snapshot_hits \
56 -o pjx_${OUT_FILE} \
57 -P "X" \
58 -y "0 0.5" --!
59
60JPlot1D \
61 -w 900x600 \
62 -f px_${OUT_FILE}:h_pmt_rate_distribution_px \
63 -o $WORKDIR/${FILEPREFIX}_pmt_rate_distribution.${FORMAT} \
64 -^ "Mean PMT rate" \
65 -B \
66 -T "%" --!
67
68convert -define pdf:use-cropbox=true -density 400 -resample 200 \
69 ${FILEPREFIX}_pmt_rate_distribution.${FORMAT} ${FILEPREFIX}_pmt_rate_distribution.png
70
71JPlot1D \
72 -w 900x600 \
73 -f px_${OUT_FILE}:h_Triggered_hits_px \
74 -f px_${OUT_FILE}:h_Triggered_hits_3dmuon_px \
75 -o $WORKDIR/${FILEPREFIX}_triggered_hits.${FORMAT} \
76 -^ "Number of triggered hits (red = JTrigger3DMuon)" \
77 -B \
78 -T "%" --!
79
80convert -define pdf:use-cropbox=true -density 400 -resample 200 \
81 ${FILEPREFIX}_triggered_hits.${FORMAT} ${FILEPREFIX}_triggered_hits.png
82
83JPlot1D \
84 -w 900x600 \
85 -f px_${OUT_FILE}:h_Snapshot_hits_px \
86 -o $WORKDIR/${FILEPREFIX}_snapshot_hits.${FORMAT} \
87 -^ "Number of snapshot hits" \
88 -B \
89 -T "%" --!
90
91convert -define pdf:use-cropbox=true -density 400 -resample 200 \
92 ${FILEPREFIX}_snapshot_hits.${FORMAT} ${FILEPREFIX}_snapshot_hits.png
93
94JPlot1D \
95 -w 900x600 \
96 -f px_${OUT_FILE}:h_Number_of_overlays_px \
97 -o $WORKDIR/${FILEPREFIX}_number_overlays.${FORMAT} \
98 -y "0 30" \
99 -^ "Number of overlays" \
100 -B \
101 -T "%" --!
102
103convert -define pdf:use-cropbox=true -density 400 -resample 200 \
104 ${FILEPREFIX}_number_overlays.${FORMAT} ${FILEPREFIX}_number_overlays.png
105
106JPlot1D \
107 -w 900x600 \
108 -f pjx_${OUT_FILE}:h_Triggered_over_Snapshot_hits_px \
109 -o $WORKDIR/${FILEPREFIX}_number_slices.${FORMAT} \
110 -^ "Number of slices" \
111 -Y \
112 -B \
113 -T "%" --!
114
115convert -define pdf:use-cropbox=true -density 400 -resample 200 \
116 ${FILEPREFIX}_number_slices.${FORMAT} ${FILEPREFIX}_number_slices.png
117
118JPlot2D \
119 -w 900x600 \
120 -f ${OUT_FILE}:h_event_duration \
121 -o $WORKDIR/${FILEPREFIX}_event_duration.${FORMAT} \
122 -Y \
123 -Z \
124 -B \
125 -O "colz" \
126 -p "57" \
127 -T "%" --!
128
129convert -define pdf:use-cropbox=true -density 400 -resample 200 \
130 ${FILEPREFIX}_event_duration.${FORMAT} ${FILEPREFIX}_event_duration.png
131
132# Check whether ghostscript is installed, and if, combine all plots
133
134if (( $+commands[gs] )); then
135 convert ${FILEPREFIX}_*.${FORMAT} ${FILEPREFIX}_JMRA.${FORMAT}
136 gs -dNOPAUSE -SDEVICE=pdfwrite -sOUTPUTFILE=${FILEPREFIX}_JMRA.${FORMAT} -dBATCH ${FILEPREFIX}_*.${FORMAT}
137 /bin/ls ${FILEPREFIX}_*.${FORMAT} | grep -v ${FILEPREFIX}_JMRA.${FORMAT} | xargs -0 -d '\n' rm
138fi
139
140rm px_${OUT_FILE}
141rm pjx_${OUT_FILE}
142
143CHECK_EXIT_CODE