examples/JCalibrate/JHobbit.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

set_variable:  DEBUG       CALIBRATION_DEBUG  1
set_variable   WORKDIR     ${TMPDIR:-/tmp}/
set_variable   QUALITY     "0.9   1.1"
set_variable   CHI2        "0.0   2.0"

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])

typeset -A    X_RANGE

#X_RANGE+=(Gauss         " -5.0  +5.0")
#X_RANGE+=(Landau        " -5.0 +10.0")
X_RANGE+=(EMG           "-10.0 +50.0")
#X_RANGE+=(BreitWigner   "-15.0 +35.0")

if [[ $DETECTOR != *\.detx ]]; then
    fatal "Invalid detector file <${DETECTOR}>."
fi

set_variable   DETECTOR_NEW        ${DETECTOR/.detx/_hobbit.detx}

cp $DETECTOR  $DETECTOR_NEW

for FUNCTION T_NS in ${(kv)X_RANGE}; do
    
    JHobbit \
        -f "$INPUT_FILES[*]"                        \
        -a $DETECTOR_NEW                            \
        -A                                          \
        -o $WORKDIR/hobbit\[${FUNCTION}\].root      \
        -F $FUNCTION                                \
        -T "$T_NS"                                  \
        -d $DEBUG

    JTestRange1D \
        -f $WORKDIR/hobbit\[${FUNCTION}\].root:hq  \
        -y "$QUALITY"                              \
        -N 0                                       \
        -d $DEBUG

    CHECK_EXIT_CODE
            
    JTestRange1D \
        -f $WORKDIR/hobbit\[${FUNCTION}\].root:hc  \
        -y "$CHI2"                                 \
        -N 0                                       \
        -d $DEBUG

    CHECK_EXIT_CODE
            
done