JRA_projections.sh

Go to the documentation of this file.

#!/bin/bash

#     _________             ___              _                    
#    |_  | ___ \           / _ \            | |                   
#      | | |_/ /   _ _ __ / /_\ \_ __   __ _| |_   _ _______ _ __ 
#      | |    / | | | '_ \|  _  | '_ \ / _` | | | | |_  / _ \ '__|
#  /\__/ / |\ \ |_| | | | | | | | | | | (_| | | |_| |/ /  __/ |   
#  \____/\_| \_\__,_|_| |_\_| |_/_| |_|\__,_|_|\__, /___\___|_|   
#                                               __/ |             
#                                              |___/              

script=${0##*/}

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

if ( do_usage $* ); then
    usage "$script [jra file string floor out dir]"
fi

set_variable WORKDIR ${PWD}

set_variable JRA_FILE ${WORKDIR}/out/JRA_4237.root

set_variable OUT_DIR ${WORKDIR}/out/JRA_Project

set_variable STRING 1

set_variable FLOOR 4

set_variable PMT 16

case $# in 
    
    4)
   OUT_DIR=$4
   ;&
    3)
   FLOOR=$3
   ;&
    2)
   STRING=$2
   ;;
    1)
   JRA_FILE=$1
   ;;
    0)
   :
   ;;
    *)
   echo "Too many arguments"
   exit
   ;;
esac

#######################################
#Variables that we will need to define#
#######################################

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

# Parse the detector file name from the complete path, and format it so that it compiles in latex.

####################################################################
# Use JFitToT to fit the ToT distributions for each PMT in the DOM #
####################################################################

if [[ ${JRA_FILE} =~ 'JRA_'(.*)(.root) ]]; then
RUN_NUMBER="${BASH_REMATCH[1]}"
fi

N_PMTS=31

SLICES=""

for((j=0;j<N_PMTS;j++)); do
    SLICES=$SLICES" $(echo ${j} - 0.5 | bc -l | sed -e 's/^\./0./' -e 's/^-\./-0./') $(echo ${j} + 0.5 | bc)"
done

JProject2D -f "${JRA_FILE}:Detector/DU${STRING}/F${FLOOR}/h_pmt_rate_distributions_TimesliceL1" -P x -y "${SLICES}" -o "${OUT_DIR}/JRA_proj_${RUN_NUMBER}_S${STRING}F${FLOOR}_TSL1.root"

JProject2D -f "${JRA_FILE}:Detector/DU${STRING}/F${FLOOR}/h_pmt_rate_distributions_Summaryslice" -P x -y "${SLICES}" -o "${OUT_DIR}/JRA_proj_${RUN_NUMBER}_S${STRING}F${FLOOR}_Summary.root"

JProject2D -f "${JRA_FILE}:Detector/DU${STRING}/F${FLOOR}/h_fifo_almost_full" -P x -y "${SLICES}" -o "${OUT_DIR}/JRA_proj_${RUN_NUMBER}_S${STRING}F${FLOOR}_FIFO.root"

JProject2D -f "${JRA_FILE}:Detector/DU${STRING}/F${FLOOR}/h_hrv" -P x -y "${SLICES}" -o "${OUT_DIR}/JRA_proj_${RUN_NUMBER}_S${STRING}F${FLOOR}_HRV.root"

JProject2D -f "${JRA_FILE}:Detector/DU${STRING}/F${FLOOR}/h_rate_vs_time_Summaryslice" -P x -y "${SLICES}" -o "${OUT_DIR}/JRA_proj_${RUN_NUMBER}_S${STRING}F${FLOOR}_summary_rate_vs_t.root"

JProject2D -f "${JRA_FILE}:Detector/DU${STRING}/F${FLOOR}/h_rate_vs_time_TimesliceL1" -P x -y "${SLICES}" -o "${OUT_DIR}/JRA_proj_${RUN_NUMBER}_S${STRING}F${FLOOR}_L1_rate_vs_t.root"


JPlot1D -f "${OUT_DIR}/JRA_proj_${RUN_NUMBER}_S${STRING}F${FLOOR}_TSL1.root:h_pmt_rate_distributions_TimesliceL1_px\\[$(echo ${PMT} - 0.5 | bc -l | sed -e 's/^\./0./' -e 's/^-\./-0./')" -x "1e-3 1e1" -X -y "1e-1 1e6"  -O "hist" -Y -G "[X][Y]" -T "TimesliceL1 S${STRING}F${FLOOR}PMT${PMT}" -o "${OUT_DIR}/rateL1_PMT_${PMT}.png" -B

JPlot1D -f "${OUT_DIR}/JRA_proj_${RUN_NUMBER}_S${STRING}F${FLOOR}_Summary.root:h_pmt_rate_distributions_Summaryslice_px\\[$(echo ${PMT} - 0.5 | bc -l | sed -e 's/^\./0./' -e 's/^-\./-0./')" -x "1e-3 1e3" -X -y "1e-1 1e6" -Y  -O "hist" -G "[X][Y]" -T "Summaryslice S${STRING}F${FLOOR}PMT${PMT}" -o "${OUT_DIR}/rateSummary_PMT_${PMT}.png" -B

JPlot1D -f "${OUT_DIR}/JRA_proj_${RUN_NUMBER}_S${STRING}F${FLOOR}_summary_rate_vs_t.root:h_rate_vs_time_Summaryslice_px\\[$(echo ${PMT} - 0.5 | bc -l | sed -e 's/^\./0./' -e 's/^-\./-0./')" -y "1e-1 1e6" -Y -G "[X][Y]" -T "Summaryslice S${STRING}F${FLOOR}PMT${PMT}"  -O "hist" -o "${OUT_DIR}/rateSummary_t_PMT_${PMT}.png" -B

JPlot1D -f "${OUT_DIR}/JRA_proj_${RUN_NUMBER}_S${STRING}F${FLOOR}_L1_rate_vs_t.root:h_rate_vs_time_TimesliceL1_px\\[$(echo ${PMT} - 0.5 | bc -l | sed -e 's/^\./0./' -e 's/^-\./-0./')" -y "1e-1 1e6" -Y -G "[X][Y]" -T "TimesliceL1 S${STRING}F${FLOOR}PMT${PMT}"  -O "hist" -o "${OUT_DIR}/rateL1_t_PMT_${PMT}.png" -B

JPlot1D -f "${OUT_DIR}/JRA_proj_${RUN_NUMBER}_S${STRING}F${FLOOR}_FIFO.root:h_fifo_almost_full_px\\[$(echo ${PMT} - 0.5 | bc -l | sed -e 's/^\./0./' -e 's/^-\./-0./')" -y "0 2" -T "FIFO S${STRING}F${FLOOR}PMT${PMT}" -O "hist" -o "${OUT_DIR}/FIFO_${PMT}.png" -B

JPlot1D -f "${OUT_DIR}/JRA_proj_${RUN_NUMBER}_S${STRING}F${FLOOR}_HRV.root:h_hrv_px\\[$(echo ${PMT} - 0.5 | bc -l | sed -e 's/^\./0./' -e 's/^-\./-0./')" -y "0 2" -T "HRV S${STRING}F${FLOOR}PMT${PMT}"  -O "hist" -o "${OUT_DIR}/HRV_${PMT}.png" -B