JCalibrateStrings.sh

Go to the documentation of this file.

#!/usr/bin/env 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} -@ "type = JDAQEvent" -k number_of_entries`
 
        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}"