JRecalibrateK40.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
script=${0##*/}
 
# ------------------------------------------------------------------------------------------
#
#                      Script to run K40 time calibration following HV tuning.
#
# ------------------------------------------------------------------------------------------
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi
 
source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
 
set_variable:  DEBUG              CALIBRATION_DEBUG          1
set_variable:  WORKDIR            CALIBRATION_WORKDIR        `pwd`
set_variable:  DAQ_TIMESLICE      CALIBRATION_TIMESLICE      JDAQTimesliceL1
set_variable+  DATABASE_OPTION    CALIBRATION_USEPMTID       -U
set_variable+  MODEL              CALIBRATION_MODEL          -D
 
if do_usage $*; then
    usage "$script <detector file> <input file> <input file> <input file>"\
          "\nThe 1st and 3rd input file should correspond to a data taking run in which the HV settings are complete (e.g. \"vendor\" or \"tuned\")."\
          "\nThe 2nd input file should correspond to a run with hybrid HV settings."\
          "\nIn this, half of the PMTs should have the same HV settings as either the 1st or 3rd input file."
fi
 
if (( $# == 4 )); then
    set_variable  DETECTOR      $argv[1]
    set_array     INPUT_FILES   $argv[2,-1]
else
    fatal "Wrong number of arguments."
fi
 
JCookie.sh
 
#
# Get background option for given input file.
#
# \param   1     input file
# \return        option
#
function get_option()
{
    if (( `JPrintTree -f $1 -@ "type = JDAQSummaryslice" -k number_of_entries` > 0 )); then
        echo  rates
    else
        echo  counts
    fi
}
    
 
#
# Standard calibration.
#
# \param   1     detector
# \param   2     input file
#
function calibrate()
{
    JPrintDAQHeader -f $2  |  read  ID  RUN  FRAME_INDEX  UTC
 
    JCalibrateK40 \
        -a $1                                        \
        -f $2                                        \
        -C $DAQ_TIMESLICE                            \
        -o $WORKDIR/monitor_${RUN}.root              \
        -b `get_option $2`                           \
        -d $DEBUG --!
 
    JMergeCalibrateK40 \
        -f $WORKDIR/monitor_${RUN}.root              \
        -o $WORKDIR/merge_${RUN}.root                \
        -@ "LIVETIME_S = 10"                         \
        -d $DEBUG
 
    JConvertDetectorFormat \
        -a $1                                        \
        -o $WORKDIR/detector_${RUN}.datx
 
    JFitK40 \
        -a $WORKDIR/detector_${RUN}.datx             \
        -f $WORKDIR/merge_${RUN}.root                \
        -o $WORKDIR/fit_${RUN}.root                  \
        $MODEL                                       \
        -A                                           \
        -w                                           \
        -d $DEBUG
}
 
 
#
# Hybrid calibration.
#
# \param   1     detector
# \param   2     reference run
# \param   3     input file
#
function hybrid()
{
    JPrintDAQHeader -f $3  |  read  ID  RUN  FRAME_INDEX  UTC
 
    JTDC  \
        -D $DETECTOR_ID                              \
        -r "$2 $RUN"                                 \
        $DATABASE_OPTION                             \
        -o "$WORKDIR/TDC_${2}_${RUN}.txt"
 
    JCalibrateK40 \
        -a $1                                        \
        -f $3                                        \
        -C $DAQ_TIMESLICE                            \
        -o $WORKDIR/monitor_${RUN}.root              \
        -b `get_option $3`                           \
        -d $DEBUG --!
 
    JMergeCalibrateK40 \
        -f $WORKDIR/monitor_${RUN}.root              \
        -o $WORKDIR/merge_${RUN}.root                \
        -@ "LIVETIME_S = 10"                         \
        -d $DEBUG
 
    JConvertDetectorFormat \
        -a $1                                        \
        -o $WORKDIR/detector_${RUN}.datx
 
    JFitK40 \
        -a $WORKDIR/detector_${RUN}.datx             \
        -f $WORKDIR/merge_${RUN}.root                \
        -o $WORKDIR/fit_${RUN}.root                  \
        $MODEL                                       \
        -A                                           \
        -w                                           \
        -! $WORKDIR/TDC_${2}_${RUN}.txt              \
        -d $DEBUG
}
 
 
# Get run numbers
 
eval `JPrintDetector -a $DETECTOR -O IDENTIFIER`
 
RUNS=()
 
for INPUT_FILE in $INPUT_FILES[*]; do
 
    JPrintDAQHeader -f $INPUT_FILE  |  read  ID  RUN  FRAME_INDEX  UTC
 
    if (( $ID != $DETECTOR_ID )); then
        fatal "Detector identifier from $INPUT_FILE $ID != $DETECTOR_ID"
    fi
    
    RUNS+=($RUN)
 
done
 
 
calibrate  $DETECTOR  ${INPUT_FILES[1]}
 
hybrid     $WORKDIR/detector_${RUNS[1]}.datx  ${RUNS[1]}   ${INPUT_FILES[2]}
 
hybrid     $WORKDIR/detector_${RUNS[2]}.datx  ${RUNS[2]}   ${INPUT_FILES[3]}