JMCEvt.sh

Go to the documentation of this file.

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

# ------------------------------------------------------------------------------------------
#
#       Utility script to test JGandalf - JEnergy - JStart on Monte Carlo true muons.
#
# ------------------------------------------------------------------------------------------


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/trigger_efficiency+background.root 
set_variable  DETECTOR            $WORKDIR/km3net_reference.detx
set_variable  ENERGY_CORRECTION   energy_correction_arca.txt

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

case $# in
    3) set_variable WORKDIR     $3;&
    2) set_variable INPUT_FILE  $2;&
    1) set_variable DETECTOR    $1;;
esac

set_variable  OUTPUT_FILE         $WORKDIR/%-mc.root

if (( 1 )); then

    for APP in gandalf energy start; do

   if [[ -f ${OUTPUT_FILE/\%/${APP}} ]]; then
       rm -i ${OUTPUT_FILE/\%/${APP}}
   fi
    done
    
    if ( ! reuse_file ${OUTPUT_FILE/\%/event} ); then

   JMCEvt \
       -f $INPUT_FILE               \
       -o ${OUTPUT_FILE/\%/event}   \
       -d ${DEBUG} --!
    fi

    JARCAGandalf.sh $DETECTOR  ${OUTPUT_FILE/\%/event}    ${OUTPUT_FILE/\%/gandalf}
    JARCAEnergy.sh  $DETECTOR  ${OUTPUT_FILE/\%/gandalf}  ${OUTPUT_FILE/\%/energy}
    JARCAStart.sh   $DETECTOR  ${OUTPUT_FILE/\%/energy}   ${OUTPUT_FILE/\%/start}
    
fi


if (( 1 )); then

    JPostfit \
   -f ${OUTPUT_FILE/\%/start}             \
   -O LOGE                                \
   -A 5                                   \
   -o postfit\[start-mc\].root

fi


if (( 1 )); then

    JPlot1D \
        -f postfit\[start-mc\].root:hx         \
        -> "angle [deg]"                      \
        -\^ "number of events [a.u]"           \
        -XX                                    \
        -L TR                                  \
        -o angle-truth.gif

fi

if (( 1 )); then

    JPlot1D \
        -f postfit\[start-mc\].root:e0                  \
        -f postfit\[start-mc\].root:e1                  \
        -f postfit\[start-mc\].root:e2                  \
        -\^ "number of events [a.u.]"                   \
        -> "log10(E_{fit}/GeV)"                        \
        -XX

fi

if (( 1 )); then

    JPlot2D \
        -f postfit\[start-mc\].root:ee                  \
        -O COLZ                                         \
        -XX -YY                                         \
        -> "E_{true} [GeV]"                            \
        -< "E_{fit} [GeV]"                             
        
fi

if (( 1 )); then

    JFit \
        -f postfit\[start-mc\].root:er                  \
        -F "[0]*exp(-0.5*(x-[1])*(x-[1])/([2]*[2]))"    \
        -@ "p0 = 0.8*GetMaximum"                        \
        -@ "p1 = 0.0"                                   \
        -@ "p2 = 0.4"                                   \
        -x "-0.8 +0.8"

    JPlot1D \
        -f fit.root:er                                  \
        -> "log10(E_{fit}/E_{true})"                   \
        -o er-truth.gif

fi

if (( 1 )); then

    JPlot2D \
        -f postfit\[start-mc\].root:ha         \
        -OCOLZ                                 \
        -> "R^{2} [m^{2}]"                    \
        -< "z [m]"                            

fi
    
if (( 1 )); then

    JPlot1D \
        -f postfit\[start-mc\].root:hz         \
        -> "z [m]"                            \
        -\^ "number of events [a.u]"           \
        -L TR                                  \
        -o z.gif

fi

if (( 0 )); then

    JPlot1D \
        -f postfit\[start-mc\].root:hx         \
        -f postfit\[start-mc\].root:hx         \
        -> "angle [deg]"                      \
        -\^ "number of events [a.u]"           \
        -XX                                    \
        -L TR                                  \
        -o angle-truth.gif

fi

if (( 0 )); then

    JSum1D -f postfit\[start-mc\].root:hx -o sum\[start-mc\].root -N
    JSum1D -f postfit\[gandalf\].root:hx  -o sum\[gandalf\].root  -N

    JPlot1D \
        -f sum\[gandalf\].root:hx              \
        -f sum\[start-mc\].root:hx             \
        -> "angle [deg]"                      \
        -\^ "integral of events"               \
        -XX                                    \
        -L TR                                  \
        -o sum-truth.gif

fi