seaTimeRecalibration.sh

Go to the documentation of this file.

#!/bin/zsh
#
#
# \author acreusot 
#
version=1.0
script=${0##*/}

# ------------------------------------------------------------------------------------------
#
#                         Utility script to test sea time recalibration.
#
# ------------------------------------------------------------------------------------------
# the detector file used as a starting point is supposed to correspond to the latest sea 
# time calibration.

if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi

source $JPP_DIR/setenv.sh $JPP_DIR

set_variable    DEBUG           1    

if do_usage $*; then
    usage "$script <working directory> <detector file> <data hv odd fix> <data hv even fix>"
fi

case $# in
        4)      set_variable    DATAHVEVENFIX   $4;
                set_variable    DATAHVODDFIX    $3;
                set_variable    DETX            $2;
                set_variable    WORKDIR         $1;;
        *) fatal "wrong number of arguments!";;
esac


# generic name for all files: detector files, root data files, root results files...
set_variable    BASENAME        ${DETX%.detx}

if [ ! -f "${WORKDIR}/${DETX}" ]; then
        echo "WARNING: detector file from sea calibration not found, check that the path to the file is correct or that the file has effectively been created"
fi

if (( 1 )); then
        JCalibrateK40 \
                -d $DEBUG \
                -f ${WORKDIR}/${DATAHVODDFIX} \
                -a ${WORKDIR}/${DETX} \
                -o calibrateK40_searecalibration_oddFix.root \
                -b rates \
                -C JDAQTimesliceL1 \
                -V0+1e5 \
                -M2+31
fi

if (( 1 )); then
        JCalibrateK40 \
                -d $DEBUG \
                -f ${WORKDIR}/${DATAHVEVENFIX} \
                -a ${WORKDIR}/${DETX} \
                -o calibrateK40_searecalibration_evenFix.root \
                -b rates \
                -C JDAQTimesliceL1 \
                -V0+1e5 \
                -M2+31
fi

if (( 1 )); then
        JMergeCalibrateK40 \
                -d $DEBUG \
                -f calibrateK40_searecalibration_oddFix.root \
                -o mergecalibrateK40_searecalibration_oddFix.root
fi

if (( 1 )); then
        JMergeCalibrateK40 \
                -d $DEBUG \
                -f calibrateK40_searecalibration_evenFix.root \
                -o mergecalibrateK40_searecalibration_evenFix.root
fi

# this creates the odd/even PMT maps. If the HV of odd PMTs is kept fix, then you need to
# keep fix the time offset of odd PMTs
if (( 1 )); then
        ./createPMTMap \
                -d $DEBUG \
                -a ${WORKDIR}/${DETX} \
                -o ${WORKDIR}/${BASENAME}_pmtOdd_map.txt \
                -T "odd"
fi

# create new detector files to perform sea recalibration on the 2 populations
# this 2 new files will be merged at the end of the script
if (( 1 )); then
        cp ${WORKDIR}/${DETX} ${WORKDIR}/${BASENAME}_intra_pmtOdd.detx 
        cp ${WORKDIR}/${DETX} ${WORKDIR}/${BASENAME}_intra_pmtEven.detx        
fi

# when the HV and time offset off odd PMTs are kept fix, 
# the time offset of even PMTs are estimated
if (( 1 )); then
        JFitK40 \
                -d $DEBUG \
                -f mergecalibrateK40_searecalibration_oddFix.root \
                -a ${WORKDIR}/${BASENAME}_intra_pmtEven.detx \
                -o fitK40_searecalibration_oddFix.root \
                -! "`cat ${WORKDIR}/${BASENAME}_pmtOdd_map.txt`" \
                -A \
                -w
fi

# when the HV and time offset off even PMTs are kept fix, 
# the time offset of odd PMTs are estimated
if (( 1 )); then
        JFitK40 \
                -d $DEBUG \
                -f mergecalibrateK40_searecalibration_evenFix.root \
                -a ${WORKDIR}/${BASENAME}_intra_pmtOdd.detx \
                -o fitK40_searecalibration_evenFix.root \
                -! "`cat ${WORKDIR}/${BASENAME}_pmtOdd_map.txt`" \
                -A \
                -w \
                -r
fi

# create a new detector for saving the merging of the 2 previous detector files
if (( 1 )); then
        cp ${WORKDIR}/${DETX} ${WORKDIR}/${BASENAME}_searecalibration.detx
fi

# merge the 2 detector files of odd/even PMT populations: the time offsets of odd pmts 
# are taken from the first detector file (option -a), the ones of even PMTs from the second 
# detector file (option -b) and both are saved in the third detector file (option -c).
if (( 1 )); then
        ./mergeTimeOffsets \
                -d $DEBUG \
                -a ${WORKDIR}/${BASENAME}_intra_pmtOdd.detx \
                -b ${WORKDIR}/${BASENAME}_intra_pmtEven.detx \
                -o ${WORKDIR}/${BASENAME}_searecalibration.detx \
                -T "odd"
fi