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