JCalibrateStrings.sh

Go to the documentation of this file.

#!/bin/zsh
script=${0##*/}

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

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

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

if (( $# < 2 )); then
    fatal "Not enough arguments."
fi

set_variable  DETECTOR            $argv[1]

INPUT_FILES=($argv[2,-1])

set_variable  TRIGGER_PARAMETERS  `echo "trigger3DMuon.enabled=1;TMaxEvent_ns=15000;"`
set_variable  JDAQ_TIMESLICE      JDAQTimesliceSN
set_variable  FORMULA             "[0] + ([3]*( ([1]^(2*[4])) / (([1]*[1] + (2^(1/[4])-1) * ((x-[2])*(x-[2])) ) )))^[4]"
set_variable  WORKDIR             /tmp

mkdir -p --  ${WORKDIR}

echo "Specify the list of string numbers to calibrate, separated by white spaces (example: <1 2 3>)"

read -A STRINGS

REMOVE=($STRINGS)

for STRING in "${STRINGS[@]}"; do
    
    REMOVE[$REMOVE[(i)$STRING]]=()

    RM=()
    for R in "${REMOVE[@]}"; do
        RM+=("-r ${R}-${R}")
    done

    set_variable  REDUCED_DETECTOR    reduced.detx
    
    set_variable  GRAPH_TXT           graph-time.txt
    set_variable  GRAPH_ROOT          graph-time.root
    set_variable  FIT_TXT             graph-fit.txt
    set_variable  FIT_ROOT            graph-fit.root
    
    rm -f ${WORKDIR}/${REDUCED_DETECTOR}
    rm -f ${WORKDIR}/${GRAPH_TXT}

    JEditDetector                             \
        -a ${DETECTOR}                        \
        ${RM}                                 \
        -o ${WORKDIR}/${REDUCED_DETECTOR} --!
    
    for TIME_OFFSET_NS in -200 -100 -20 -5 0 +5 +20 +100 +200; do

        JEditDetector                           \
            -a ${WORKDIR}/${REDUCED_DETECTOR}    \
            -S "$STRING add $TIME_OFFSET_NS"    \
            -o ${WORKDIR}/detector.detx --!     
    
        set_variable  OUTPUT_FILE  JTriggerProcessor.root

        rm -f ${WORKDIR}/${OUTPUT_FILE}

        JTriggerProcessor                   \
            -a ${WORKDIR}/detector.detx     \
            -f ${INPUT_FILES}               \
            -C ${JDAQ_TIMESLICE}            \
            -@ ${TRIGGER_PARAMETERS}        \
            -o ${WORKDIR}/${OUTPUT_FILE} --!

        set_variable N_TRIGGERS  `JPrintTree -f ${WORKDIR}/${OUTPUT_FILE} | awk '/KM3NETDAQ::JDAQEvent/ { print $3} '`

        printf "%5.3f %12.3f\n" ${TIME_OFFSET_NS} ${N_TRIGGERS} >> ${WORKDIR}/${GRAPH_TXT}

    done

    JGraph                                \
        -f ${WORKDIR}/${GRAPH_TXT}        \
        -o ${WORKDIR}/${GRAPH_ROOT}

    JFit                                  \
        -f ${WORKDIR}/${GRAPH_ROOT}:\.\*  \
        -o ${WORKDIR}/${FIT_ROOT}         \
        -F "$FORMULA"                     \
        -@"p0 = 0"                        \
        -@"p1 = GetRMS"                   \
        -@"p2 = GetMean"                  \
        -@"p3 = GetMaximum"               \
        -@"p4 = 1"                        \
        -d1 >> ${WORKDIR}/${FIT_TXT} 2>&1

    set_variable   OPTIMAL_OFFSET    `JTestFit -f ${WORKDIR}/${FIT_ROOT}:${GRAPH_TXT} -@ "p2 -1000 1000" -d 3 | awk '/p2/ {print $2}'`

    echo "JEditDetector -a ${DETECTOR} -S \"$STRING_ID add $OPTIMAL_OFFSET\" -o ${WORKDIR}/${OUTPUT_DETECTOR} --!"

done