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         /tmp
set_variable     OUTPUT_FILE     $WORKDIR/charge.root
set_variable     PMT_FILE        "gain=0.9"
set_variable     NPE             1 
set_variable     xrange          "0 $((3*NPE))"


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

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


if (( 1 )); then

    $JPP_DIR/examples/JDetector/JCharge \
        -o $OUTPUT_FILE    \
        -n 999999          \
        -N $NPE            \
        -P $PMT_FILE       \
        -d $DEBUG

fi

if (( 0 )); then
    
    JPlot1D \
        -f ${OUTPUT_FILE}:h0                 \
        -f ${OUTPUT_FILE}:h1                 \
        -x $xrange                           \
        -\> "Charge [npe]"                   \
        -\^ "Probability density [npe^{-1}]" \
        -S 0.4                               \
        -y "1e-3 1.1" --                 

fi

if (( 0 )); then
    
    JSum1D \
        -f ${OUTPUT_FILE}:h0            \
        -f ${OUTPUT_FILE}:h1            \
        -o $WORKDIR/sum.root

    JPlot1D \
        -f $WORKDIR/sum.root:h0         \
        -f $WORKDIR/sum.root:h1         \
        -x $xrange                      \
        -\> "Charge [npe]"              \
        -\^ "Probability"    \
        -S 0.4                          \
        -G XY                           \
        -y "1e-3 1.1" -Y --               

fi