JMergeRunAnalyzer.sh

Go to the documentation of this file.

#!/bin/zsh

version=1.0
script=${0##*/}

source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null

if do_usage $*; then
    usage "$script <input_files> <output_file(.root)> [only_directory_path]"
fi

set_variable:    ONLY_DIR    INPUT_DIR          0
set_variable:    WORKDIR     WORKDIR          `pwd`
set_variable:    FORMAT      GRAPHICS_FORMAT   pdf

case $# in
    3)  set_variable    ONLY_DIR                 $3;&
    2)  set_variable    OUT_FILE                 $2
        set_variable    JRA_FILES                $1;;
    *)  fatal "Wrong number of arguments"
esac

set_variable    INPUT_FILES  ${JRA_FILES}

if [[ ${ONLY_DIR} != 0 ]]
then
set_variable INPUT_FILES  `echo ${JRA_FILES}/*`    
fi

set_variable    OUT_PATH        `readlink -f $(dirname ${OUT_FILE})`
set_variable    JRA_DIR         ${JPP_DIR}/software/JRunAnalyzer
set_variable    FILEPREFIX      "plot"

if [ ! -d ${OUT_PATH} ] 
then     
    mkdir -p ${OUT_PATH} 
fi

time ${JRA_DIR}/JMergeRunAnalyzer  \
     -f $INPUT_FILES               \
     -o $OUT_FILE --!

#****** Create Profiling and Projection of the TH2Ds *******#

JProfile2D \
    -f ${OUT_FILE}:h_pmt_rate_distribution \
    -f ${OUT_FILE}:h_Triggered_hits        \
    -f ${OUT_FILE}:h_Snapshot_hits         \
    -f ${OUT_FILE}:h_Number_of_overlays    \
    -o px_${OUT_FILE}                      \
    -P "X"                                 \
    -O "g" --!

JProject2D \
    -f ${OUT_FILE}:h_Triggered_over_Snapshot_hits \
    -o pjx_${OUT_FILE}                            \
    -P "X"                                        \
    -y "0 0.5" --!

JPlot1D \
    -w 900x600                                                  \
    -f px_${OUT_FILE}:h_pmt_rate_distribution_px                \
    -o $WORKDIR/${FILEPREFIX}_pmt_rate_distribution.${FORMAT}   \
    -^ "Mean PMT rate"                                          \
    -B                                                          \
    -T "%" --!

convert -define pdf:use-cropbox=true -density 400 -resample 200 \
        ${FILEPREFIX}_pmt_rate_distribution.${FORMAT} ${FILEPREFIX}_pmt_rate_distribution.png

JPlot1D \
    -w 900x600                                             \
    -f px_${OUT_FILE}:h_Triggered_hits_px                  \
    -f px_${OUT_FILE}:h_Triggered_hits_3dmuon_px           \
    -o $WORKDIR/${FILEPREFIX}_triggered_hits.${FORMAT}     \
    -^ "Number of triggered hits (red = JTrigger3DMuon)"   \
    -B                                                     \
    -T "%" --!

convert -define pdf:use-cropbox=true -density 400 -resample 200 \
        ${FILEPREFIX}_triggered_hits.${FORMAT} ${FILEPREFIX}_triggered_hits.png

JPlot1D \
    -w 900x600                                         \
    -f px_${OUT_FILE}:h_Snapshot_hits_px               \
    -o $WORKDIR/${FILEPREFIX}_snapshot_hits.${FORMAT}  \
    -^ "Number of snapshot hits"                       \
    -B                                                 \
    -T "%" --!

convert -define pdf:use-cropbox=true -density 400 -resample 200 \
        ${FILEPREFIX}_snapshot_hits.${FORMAT} ${FILEPREFIX}_snapshot_hits.png

JPlot1D \
    -w 900x600                                          \
    -f px_${OUT_FILE}:h_Number_of_overlays_px           \
    -o $WORKDIR/${FILEPREFIX}_number_overlays.${FORMAT} \
    -y "0 30"                                           \
    -^ "Number of overlays"                             \
    -B                                                  \
    -T "%" --!

convert -define pdf:use-cropbox=true -density 400 -resample 200 \
        ${FILEPREFIX}_number_overlays.${FORMAT} ${FILEPREFIX}_number_overlays.png

JPlot1D \
    -w 900x600                                           \
    -f pjx_${OUT_FILE}:h_Triggered_over_Snapshot_hits_px \
    -o $WORKDIR/${FILEPREFIX}_number_slices.${FORMAT}    \
    -^ "Number of slices"                                \
    -Y                                                   \
    -B                                                   \
    -T "%" --!

convert -define pdf:use-cropbox=true -density 400 -resample 200 \
        ${FILEPREFIX}_number_slices.${FORMAT} ${FILEPREFIX}_number_slices.png

JPlot2D \
    -w 900x600                                          \
    -f ${OUT_FILE}:h_event_duration                     \
    -o $WORKDIR/${FILEPREFIX}_event_duration.${FORMAT}  \
    -Y                                                  \
    -Z                                                  \
    -B                                                  \
    -O "colz"                                           \
    -p "57"                                             \
    -T "%" --!

convert -define pdf:use-cropbox=true -density 400 -resample 200 \
        ${FILEPREFIX}_event_duration.${FORMAT} ${FILEPREFIX}_event_duration.png

# Check whether ghostscript is installed, and if, combine all plots

if (( $+commands[gs] )); then
     convert ${FILEPREFIX}_*.${FORMAT} ${FILEPREFIX}_JMRA.${FORMAT}
     gs -dNOPAUSE -SDEVICE=pdfwrite -sOUTPUTFILE=${FILEPREFIX}_JMRA.${FORMAT} -dBATCH ${FILEPREFIX}_*.${FORMAT}
     /bin/ls ${FILEPREFIX}_*.${FORMAT} | grep -v ${FILEPREFIX}_JMRA.${FORMAT} | xargs -0 -d '\n' rm
fi

rm px_${OUT_FILE}
rm pjx_${OUT_FILE}

CHECK_EXIT_CODE