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>+ <string number>+"
fi

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

set_variable INDEX_STRING_START -1

while [[ ${${argv[${START_INDEX_STRING}]}##*.} != "root" ]]; do
    set_variable START_INDEX_STRING     $(( START_INDEX_STRING - 1 ))
done

set_variable DETECTOR    ${argv[1]}
set_array INPUT_FILES    ${argv[2, $(( START_INDEX_STRING - 1 ))     ]}
set_array STRINGS        ${argv[   $(( START_INDEX_STRING     )), -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             ${TMPDIR:-/tmp}

mkdir -p --  ${WORKDIR}

REMOVE=($STRINGS)

OUTPUT_OFFSETS=""

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}'`

    JEditDetector                           \
        -a ${DETECTOR}                      \
        -s "$STRING add $OPTIMAL_OFFSET"    \
        -o ${DETECTOR} --!      

    OUTPUT_OFFSETS="${OUTPUT_OFFSETS} JEditDetector -a ${DETECTOR} -S \"$STRING add $OPTIMAL_OFFSET\" -o ${WORKDIR}/${OUTPUT_DETECTOR} --! \n" 
    
done

echo -e "${OUTPUT_OFFSETS}"