JMuonMCEvt.sh

Go to the documentation of this file.

#!/usr/bin/env 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   PDF                $JPP_DATA/J%p.dat                            
set_variable   ENERGY_CORRECTION  energy_correction_arca.root
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;
       set_variable DETECTOR    $1;;
    *) fatal "Wrong number of arguments.";;
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
 
    $JPP_BIN/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 $JPP_BIN/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 $JPP_BIN/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 $JPP_BIN/JMuonEnergy \
        -a $DETECTOR                   \
        -f ${OUTPUT_FILE/\%/start}     \
        -o ${OUTPUT_FILE/\%/energy}    \
        -P $PDF                        \
        -E $ENERGY_CORRECTION          \
        -@ `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:E_0              \
        -f $WORKDIR/postfit\[MC\].root:E_1              \
        -f $WORKDIR/postfit\[MC\].root:E_2              \
        -\^ "number of events [a.u.]"                   \
        -> "log10(E_{fit}/GeV)"                        \
        -XX                                             \
        -o e0.$FORMAT                                   $BATCH
fi
 
if (( 1 )); then
 
    JPlot2D \
        -f $WORKDIR/postfit\[MC\].root:ExE              \
        -O COLZ                                         \
        -XX -YY                                         \
        -> "E_{true} [GeV]"                            \
        -< "E_{fit} [GeV]"                             \
        -o ee.$FORMAT                                   $BATCH
fi
 
if (( 1 )); then
 
    $JPP_BIN/JFit \
        -f $WORKDIR/postfit\[MC\].root:E_E              \
        -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:E_E                                 \
        -> "log10(E_{fit}/E_{true})"                   \
        -o er-truth.$FORMAT                             $BATCH
fi
 
if (( 1 )); then
 
    JPlot2D \
        -f $WORKDIR/postfit\[MC\].root:Va              \
        -OCOLZ                                         \
        -> "R^{2} [m^{2}]"                            \
        -< "z [m]"                                    \
        -o ha.$FORMAT                                  $BATCH
fi
    
if (( 1 )); then
 
    JPlot1D \
        -f $WORKDIR/postfit\[MC\].root:hz0             \
        -f $WORKDIR/postfit\[MC\].root:hz1             \
        -> "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