JCharge.sh

Go to the documentation of this file.

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

# ------------------------------------------------------------------------------------------
#
#     Utility script to plot charge distribution from JDETECTOR::JPMTDefaultSignalProcessor.
#
# ------------------------------------------------------------------------------------------


if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi


source $JPP_DIR/setenv.sh $JPP_DIR


set_variable   DEBUG        3                     
set_variable   WORKDIR      ${TMPDIR:-/tmp}                  
set_variable   OUTPUT_FILE  
set_variable                PMT_FILE                 "gain=0.9"  
set_variable   NPE          1                     
set_variable                X_RANGE                  "0 $((3*NPE))"  
set_variable:  FORMAT       GRAPHICS_FORMAT        gif
set_variable+  BATCH        GRAPHICS_BATCH         -B


if do_usage $*; then
    usage "$script [output file] [pmt file]"
fi

case $# in
    2) set_variable  PMT_FILE     $2;&
    1) set_variable  OUTPUT_FILE  $1;;
esac


$JPP_DIR/examples/JDetector/JCharge \
    -o $WORKDIR/charge.root  \
    -n 10000000              \
    -N $NPE                  \
    -P $PMT_FILE             \
    -d $DEBUG

for TEST in 1 2; do

    JPlot1D \
        -f $WORKDIR/charge.root:${TEST}.1     \
        -f $WORKDIR/charge.root:${TEST}.0     \
        -x $X_RANGE                           \
        -\> "Charge [npe]"                    \
        -\^ "Probability density [npe^{-1}]"  \
        -S 0.4                                \
        -o charge.$FORMAT                     $BATCH
    
    JSum1D \
        -f $WORKDIR/charge.root:${TEST}.0     \
        -f $WORKDIR/charge.root:${TEST}.1     \
        -o $WORKDIR/sum.root -B

    JPlot1D \
        -f $WORKDIR/sum.root:${TEST}.1        \
        -f $WORKDIR/sum.root:${TEST}.0        \
        -x $X_RANGE                           \
        -y "1e-3 1.1" -Y                      \
        -\> "Charge [npe]"                    \
        -\^ "Probability"                     \
        -O "]["                               \
        -S 0.4                                \
        -G XY                                 \
        -o charge.$FORMAT                     $BATCH
done

rm -f $WORKDIR/{charge,sum}.root