JCalibratePMT.sh

Go to the documentation of this file.

#!/bin/zsh
#
# \author rgruiz, mdejong
#
script=${0##*/}

# ------------------------------------------------------------------------------------------
#
#                         Script to run complete PMT calibration.
#
# ------------------------------------------------------------------------------------------

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

source $JPP_DIR/setenv.sh $JPP_DIR

set_variable:  WORKDIR           CALIBRATION_WORKDIR           ./
set_variable:  DEBUG             CALIBRATION_DEBUG             2
set_variable+  OPTION            CALIBRATION_OPTION            -A
set_variable:  QE_MAX            CALIBRATION_QEMAX             2.0
set_variable:  MU                CALIBRATION_MU                0.25
set_variable:  TOT_MIN_NS        CALIBRATION_TOT_MIN_NS        0.0
set_variable:  TOT_MAX_NS        CALIBRATION_TOT_MAX_NS        35.0

if do_usage $*; then
    usage "$script <detector_file> (input file)+ <PMT parameters file>"
fi

if (( $# < 3 )); then
    fatal "Wrong number of arguments."
fi 

set_variable   DETECTOR          $argv[1]
set_array      INPUT_FILES       $argv[2,-2]
set_variable   PMT_FILE          $argv[-1]

if [[ ! -d $WORKDIR ]]; then     
    mkdir -p $WORKDIR
fi

JCookie.sh

eval `JPrintDetector -a $DETECTOR -O IDENTIFIER`
eval `JPrintDetector -a $DETECTOR -O SUMMARY`

debug "$DETECTOR $DETECTOR_ID"

if [[ "$DETECTOR_ID" != [1-9][0-9]* ]]; then
    fatal "Invalid detector identifier \"$DETECTOR_ID\"."
fi

# evaluate inputs

set_variable   RUNSETUPID
typeset -a     RANGE

RANGE[1]=99999999
RANGE[2]=0

for INPUT_FILE in $INPUT_FILES[*]; do

    JPrintDAQHeader  -f $INPUT_FILE  |  read ID RUN FRAME_INDEX UTC    

    if (( $DETECTOR_ID != $ID )); then
        fatal "Inconsistent detector identifier $DETECTOR_ID != $ID."
    fi

    JRuns  -D $DETECTOR_ID  -@ "RUN = $RUN"  -F RUNSETUPID  |  read SETUP

    debug "Run $RUN -> $SETUP."

    if (( $RUN < $RANGE[1] )); then; RANGE[1]=$RUN; fi
    if (( $RUN > $RANGE[2] )); then; RANGE[2]=$RUN; fi

    if [[ -z "$SETUP" ]]; then
        fatal "Missing setup identifier for run $RUN."
    fi
    if [[ -n "$RUNSETUPID" && $RUNSETUPID != $SETUP ]]; then
        fatal "Inconsistent run setup identifiers $RUNSETUPID != $SETUP."
    fi

    set_variable   RUNSETUPID  $SETUP
done


# data selection

for DAQ_TIMESLICE in JDAQTimesliceL1 JDAQTimeslice; do

    JPrintTree -f${^INPUT_FILES} |& eval awk \'\{ if \( \$2 \~ \"${DAQ_TIMESLICE}$\" \) \{ printf \" %d\", \$3 \}\}\'  | read -A VALUES

    if [[ -n "$VALUES" ]] && (( ${#VALUES} == ${#INPUT_FILES} )); then
        break
    fi
done

if [[ -z "$VALUES" ]] || (( ${#VALUES} != ${#INPUT_FILES} )); then
    fatal "No valid data in input file(s)."
fi


# PMT threshold

JPMTThreshold \
    -D $DETECTOR_ID        \
    -r $RANGE[1]           \
    -P $PMT_FILE           \
    --!


# PMT TTS function

JPMTTTS \
    -D $DETECTOR_ID        \
    -P $PMT_FILE           \
    --!


# PMT gain and gain spread

set_variable  OUTPUT_FILE  $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RANGE[1]}-${(l:8::0::0:)RANGE[2]}_%.root

JCalibrateToT \
    -a $DETECTOR                              \
    -f${^INPUT_FILES}                         \
    -o ${OUTPUT_FILE/\%/calibratetot}         \
    -C $DAQ_TIMESLICE                         \
    -d $DEBUG --!

CHECK_EXIT_CODE

JMergeCalibrateToT \
    -f ${OUTPUT_FILE/\%/calibratetot}         \
    -o ${OUTPUT_FILE/\%/mergetot}             \
    -d $DEBUG --!

CHECK_EXIT_CODE

JFitToT \
    -a $DETECTOR                              \
    -f ${OUTPUT_FILE/\%/mergetot}             \
    -o ${OUTPUT_FILE/\%/fittot}               \
    -P $PMT_FILE                              \
    -w                                        \
    -x "$TOT_MIN_NS $TOT_MAX_NS"              \
    -d $DEBUG --!

CHECK_EXIT_CODE


# PMT QE

JCalibrateK40 \
    -a $DETECTOR                              \
    -f${^INPUT_FILES}                         \
    -o ${OUTPUT_FILE/\%/calibratek40}         \
    -C $DAQ_TIMESLICE                         \
    -d $DEBUG --!

CHECK_EXIT_CODE

JMergeCalibrateK40 \
    -f ${OUTPUT_FILE/\%/calibratek40}         \
    -o ${OUTPUT_FILE/\%/mergek40}             \
    -d $DEBUG --!

CHECK_EXIT_CODE

JFitK40 \
    -a $DETECTOR                              \
    -f ${OUTPUT_FILE/\%/mergek40}             \
    -o ${OUTPUT_FILE/\%/fitk40}               \
    -P $PMT_FILE                              \
    -w                                        \
    $OPTION                                   \
    -d $DEBUG --!

CHECK_EXIT_CODE


# final corrections.

JEditPMTParameters \
    -P $PMT_FILE                              \
    -E $MU                                    \
    -T "$TOT_MIN_NS $TOT_MAX_NS"              \
    -Q "0.0 $QE_MAX"                          \
    -o $PMT_FILE                              \
    -d $DEBUG --!