JPlotData.sh

Go to the documentation of this file.

#!/bin/zsh
#
# \author lquinn
#
version=1.0
script=${0##*/}

# ------------------------------------------------------------------------------------------
#
#                 Utility script to test make basic plots of data
#
# ------------------------------------------------------------------------------------------


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           $JPP_DATA
set_variable  INPUT_FILE        $WORKDIR/gandalf.root
set_variable  OUTPUT_FILE       $WORKDIR/postfit.root
set_variable  DETECTOR          $WORKDIR/km3net_reference.detx
set_variable  APPLICATION       5                               # JStart - See JFitApplications.hh for more info
set_variable  NUMBER_OF_PREFITS 1

if ( do_usage $* ); then
    usage "$script [detector file [input file [output file]]]"
fi

case $# in
    3) OUTPUT_FILE=$3;&
    2) INPUT_FILE=$2;&
    1) DETECTOR=$1;&
esac

print_variable     INPUT_FILE OUTPUT_FILE
check_input_file   $INPUT_FILE

if (( 1 )); then

    JPlotData                        \
   -f $INPUT_FILE               \
   -o $OUTPUT_FILE              \
        -a $DETECTOR                 \
        -A $APPLICATION              \
   -N $NUMBER_OF_PREFITS        \
   -d $DEBUG  --!

fi

if (( 1 )); then

    JPlot1D \
   -f ${OUTPUT_FILE}:hz         \
   -> "cos(#theta)"            \
   -\^ "number of events [a.u]" \
   -o zenith.gif
fi

if (( 1 )); then

    JPlot2D \
   -f ${OUTPUT_FILE}:hxy   \
   -> "x [m]"             \
   -< "y [m]"             \
        -O COLZ                 \
   -o 2Dplane.gif
fi

if (( 1 )); then

    JPlot2D \
   -f ${OUTPUT_FILE}:hzR   \
   -> "R [m]"             \
   -< "z [m]"             \
        -O COLZ                 \
   -o heightVsRadius.gif
fi

if (( 1 )); then

    JPlot1D \
   -f ${OUTPUT_FILE}:hq1          \
   -> "quality"                  \
   -\^ "number of events [a.u]"   \
   -o quality.gif
fi

if (( 1 )); then

    JPlot1D \
   -f ${OUTPUT_FILE}:hq2          \
   -> "-log(L)/NDF"            \
   -\^ "number of events [a.u]" \
   -o loglikelihood.gif
fi

if (( 1 )); then

  JPlot1D \
        -f ${OUTPUT_FILE}:hmu           \
        -> "P(upgoing) - P(downgoing)" \
   -\^ "number of events [a.u]"    \
        -o atmosphericmuon.gif
fi

if (( 1 )); then

    JPlot2D \
   -f ${OUTPUT_FILE}:hmz           \
   -> "P(upgoing) - P(downgoing)" \
   -< "cos(#theta)"               \
        -O COLZ                         \
   -o atmosphericmuon2D.gif
fi

if (( 1 )); then

    JPlot2D \
   -f ${OUTPUT_FILE}:hqz           \
   -> "-log(L)/NDF"               \
   -< "cos(#theta)"               \
        -O COLZ                         \
   -o zenithVsQuality.gif
fi