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