JMuonPath.sh

Go to the documentation of this file.

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

# ------------------------------------------------------------------------------------------
#
#                 Utility script to test JMuonPath 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   PDF               $JPP_DATA/J%p.dat
set_variable:  RATE_HZ           RECONSTRUCTION_RATE_HZ              6.0e3
set_variable   TTS_NS            2.0
set_variable   E_GEV             1e3
set_variable   ROAD_WIDTH_M      200
set_variable:  FORMAT            GRAPHICS_FORMAT        gif
set_variable+  BATCH             GRAPHICS_BATCH         -B

if do_usage $*; then
    usage "$script <detector file> <input file> <output file> <option>"\
          "\nPossible options: run, plot, clean."
fi

if (( $# != 4 )); then
    fatal "Wrong number of arguments."
fi

set_variable  DETECTOR      $1
set_variable  INPUT_FILE    $2
set_variable  OUTPUT_FILE   $3
set_variable  OPTION        $4

if [[ ! "$OUTPUT_FILE" == *"%"* ]]; then
    fatal "The output file must have the wildcard % in the name, e.g. %.root"
fi

if [[ $OPTION != "run" && $OPTION != "plot" && $OPTION != "clean" ]]; then
    fatal "Invalid option $OPTION."
fi

if [[ $OPTION == "run" ]]; then

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

    $JPP_DIR/examples/JReconstruction/JMuonPath \
        -a $DETECTOR                    \
        -P $PDF                         \
        -f ${OUTPUT_FILE/\%/mc}         \
        -o ${OUTPUT_FILE/\%/path}       \
        -N 1                            \
        -T ${TTS_NS}                    \
        -B ${RATE_HZ}                   \
        -E ${E_GEV}                     \
        -R ${ROAD_WIDTH_M}              \
        -d ${DEBUG} --!
fi

if [[ $OPTION == "plot" ]]; then

    JMuonPostfit \
        -f ${OUTPUT_FILE/\%/path}           \
        -o postfit\[path\].root             \
        -A 1003                             \
        -O N

    JPlot1D \
        -f postfit\[path\].root:hx          \
        -\> "angle [deg]"                   \
        -\^ "number of events [a.u.]"       \
        -XX                                 \
        -T ""                               \
        -o angle.$FORMAT                    $BATCH

    QUANTILES="0.33 0.5 0.66"

    JQuantiles2D \
        -f postfit\[path\].root:h2          \
        -Q "$QUANTILES"                     \
        -o quantiles\[path\].root

    JPlot1D \
        -f quantiles\[path\].root:h21       \
        -\> "number of modules"             \
        -\^ "resolution (reco muon) [deg]"  \
        -y "1e-2 1e2"                       \
        -Y                                  \
        -T ""                               \
        -o resolution.$FORMAT               $BATCH
fi

if [[ $OPTION == "clean" ]]; then
    rm -f ${OUTPUT_FILE/\%/mc}
    rm -f ${OUTPUT_FILE/\%/path}
    rm -f postfit\[path\].root
    rm -f quantiles\[path\].root
    rm -f angle.$FORMAT
    rm -f resolution.$FORMAT
fi