JPostfit2F.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
#
# \author mdejong
#
version=1.0
script=${0##*/}
 
# ------------------------------------------------------------------------------------------
#
#                         Utility script to optimise sorting of fits.
#
# ------------------------------------------------------------------------------------------
 
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi
 
 
source $JPP_DIR/setenv.sh $JPP_DIR
 
set_variable   DEBUG         2                          
set_variable   WORKDIR       ${TMPDIR:-/tmp}/
set_variable   ANGLE_DEG     0.0                        
set_variable   JPP           $JPP_DIR/examples/JReconstruction/
set_variable:  FORMAT        GRAPHICS_FORMAT        gif
set_variable+  BATCH         GRAPHICS_BATCH         -B
 
if do_usage $*; then
    usage "$script <input file> <input file> [<angle [deg]>]"
fi
 
case $# in
    3)  set_variable   ANGLE_DEG     $3;&
    2)  set_variable  INPUT_FILE_B   $2;
        set_variable  INPUT_FILE_A   $1;;
    *)  fatal "Wrong number of arguments."
esac
 
$JPP/JPostfit2F \
    -a $INPUT_FILE_A         \
    -b $INPUT_FILE_B         \
    -o $WORKDIR/postfit.root \
    -A $ANGLE_DEG            \
    -d ${DEBUG} --!
 
JPlot1D \
    -f $WORKDIR/postfit.root:h.A           \
    -f $WORKDIR/postfit.root:h.B           \
    -> "angle [deg]"                      \
    -\^ "number of events [a.u.]"          \
    -XX                                    \
    -L TR -T ""                            \
    -o angle1D.$FORMAT                     $BATCH
 
JPlot2D \
    -f $WORKDIR/postfit.root:h2            \
    -> "quality"                          \
    -< "angle [deg]"                      \
    -O COLZ                                \
    -T ""                                  \
    -o angle2D.$FORMAT                     $BATCH
 
for PARAMETER in \
    JQUALITY                     \
    JGANDALF_BETA0_RAD           \
    JGANDALF_BETA1_RAD           \
    JGANDALF_NUMBER_OF_HITS      \
    JSTART_NPE_MIP               \
    JSTART_NPE_MIP_TOTAL         \
    JSTART_LENGTH_METRES         \
    JENERGY_ENERGY               \
    JENERGY_CHI2;                \
    do
 
    JPlot1D \
        -f "$WORKDIR/postfit.root:^\[A\].${PARAMETER}$"    \
        -f "$WORKDIR/postfit.root:^\[B\].${PARAMETER}$"    \
        -\^ "number of events [a.u.]"                      \
        -N "X 505"                                         \
        -L TR                                              \
        -T "$PARAMETER"                                    \
        -o AB.${PARAMETER:l}.$FORMAT                       $BATCH
 
    JPlot1D \
        -f "$WORKDIR/postfit.root:^\[C\].${PARAMETER}$"    \
        -f "$WORKDIR/postfit.root:^\[D\].${PARAMETER}$"    \
        -> "#Delta"                                       \
        -\^ "number of events [a.u.]"                      \
        -N "X 505"                                         \
        -L TR                                              \
        -T "$PARAMETER"                                    \
        -o CD.${PARAMETER:l}.$FORMAT                       $BATCH
done