JMuonMCEvt.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   PDF                $JPP_DATA/J%p.dat                            
set_variable   ENERGY_CORRECTION  energy_correction_arca.txt                   
set_variable:  FORMAT             GRAPHICS_FORMAT        gif
set_variable+  BATCH              GRAPHICS_BATCH         -B

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 event gandalf start energy; do
        rm -f ${OUTPUT_FILE/\%/${APP}}
    done

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

    PARAMETERS=(
        TTS_ns            2
        E_GeV             1.0e3
        TMin_ns            -50.0
        TMax_ns           +450.0
        R_Hz              6.0e3
        roadWidth_m       175.0
        numberOfPrefits   0)

    eval JMuonGandalf \
        -a $DETECTOR                   \
        -f ${OUTPUT_FILE/\%/event}     \
        -o ${OUTPUT_FILE/\%/gandalf}   \
        -P $PDF                        \
        -@ `make_equation PARAMETERS`  \
        -d $DEBUG --!

    PARAMETERS=(
        roadWidth_m      160.0
        R_Hz             6.0e3
        TMin_ns          -25.0
        TMax_ns          +75.0
        Pmin1            1.0e-3
        Pmin2            1.0e-2
        numberOfPrefits  1)

    eval JMuonStart \
        -a $DETECTOR                   \
        -f ${OUTPUT_FILE/\%/gandalf}   \
        -o ${OUTPUT_FILE/\%/start}     \
        -P $PDF                        \
        -@ `make_equation PARAMETERS`  \
        -d $DEBUG --!

    PARAMETERS=(
        roadWidth_m      200.0
        R_Hz             6.0e3
        EMin_log         0.0
        EMax_log         8.0
        TMin_ns           -50.0
        TMax_ns          +450.0
        mestimator       1
        numberOfPrefits  1)

    eval JMuonEnergy \
        -a $DETECTOR                   \
        -f ${OUTPUT_FILE/\%/start}     \
        -o ${OUTPUT_FILE/\%/energy}    \
        -P $PDF                        \
        -E energy_correction_arca.txt  \
        -@ `make_equation PARAMETERS`  \
        -d $DEBUG --!    
fi

if (( 1 )); then

    JMuonPostfit \
        -f ${OUTPUT_FILE/\%/energy}            \
        -O LOGE                                \
        -A 4                                   \
        -o $WORKDIR/postfit\[MC\].root
fi

if (( 1 )); then

    JPlot1D \
        -f $WORKDIR/postfit\[MC\].root:hx               \
        -\> "angle [deg]"                               \
        -\^ "number of events [a.u.]"                   \
        -XX                                             \
        -L TR                                           \
        -o angle-truth.$FORMAT                          $BATCH
fi

if (( 1 )); then

    JPlot1D \
        -f $WORKDIR/postfit\[MC\].root:e0               \
        -f $WORKDIR/postfit\[MC\].root:e1               \
        -f $WORKDIR/postfit\[MC\].root:e2               \
        -\^ "number of events [a.u.]"                   \
        -\> "log10(E_{fit}/GeV)"                        \
        -XX                                             \
        -o e0.$FORMAT                                   $BATCH
fi

if (( 1 )); then

    JPlot2D \
        -f $WORKDIR/postfit\[MC\].root:ee               \
        -O COLZ                                         \
        -XX -YY                                         \
        -\> "E_{true} [GeV]"                            \
        -\< "E_{fit} [GeV]"                             \
        -o ee.$FORMAT                                   $BATCH
fi

if (( 1 )); then

    JFit \
        -f $WORKDIR/postfit\[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.$FORMAT                             $BATCH
fi

if (( 1 )); then

    JPlot2D \
        -f $WORKDIR/postfit\[MC\].root:ha              \
        -OCOLZ                                         \
        -\> "R^{2} [m^{2}]"                            \
        -\< "z [m]"                                    \
        -o ha.$FORMAT                                  $BATCH
fi
    
if (( 1 )); then

    JPlot1D \
        -f $WORKDIR/postfit\[MC\].root:hz              \
        -\> "z [m]"                                    \
        -\^ "number of events [a.u.]"                  \
        -L TR                                          \
        -o z.$FORMAT                                   $BATCH
fi

if (( 0 )); then

    JPlot1D \
        -f $WORKDIR/postfit\[MC\].root:hx              \
        -f $WORKDIR/postfit\[MC\].root:hx              \
        -\> "angle [deg]"                              \
        -\^ "number of events [a.u.]"                  \
        -XX                                            \
        -L TR                                          \
        -o angle-truth.$FORMAT                         $BATCH
fi

if (( 0 )); then

    JSum1D -f $WORKDIR/postfit\[MC\].root:hx        -o $WORKDIR/sum\[MC\].root       -N
    JSum1D -f $WORKDIR/postfit\[gandalf\].root:hx   -o $WORKDIR/sum\[gandalf\].root  -N

    JPlot1D \
        -f $WORKDIR/sum\[gandalf\].root:hx             \
        -f $WORKDIRsum\[MC\].root:hx                   \
        -\> "angle [deg]"                              \
        -\^ "integral of events"                       \
        -XX                                            \
        -L TR                                          \
        -o sum-truth.$FORMAT                           $BATCH

fi