JTransitTime.sh

Go to the documentation of this file.

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

# ------------------------------------------------------------------------------------------
#
#     Utility script to plot transition time spread from JDETECTOR::JPMTAnalogueSignaProcessor.
#
# ------------------------------------------------------------------------------------------


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/tts.root
set_variable     NPE             1
set_variable     TTS_NS          2.0


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

case $# in
    1) OUTPUT_FILE=$1;;
esac


$JPP_DIR/examples/JDetector/JTransitTime \
    -o $OUTPUT_FILE    \
    -n 999999          \
    -N $NPE            \
    -T $TTS_NS         \
    -d $DEBUG

for HISTOGRAM in tts tt2 pmt; do

    JFit \
   -f ${OUTPUT_FILE}:${HISTOGRAM}                  \
        -F "[0]*exp(-0.5*(x-[1])*(x-[1])/([2]*[2]))"    \
        -@ "p0 = 0.8*GetMaximum"                        \
        -@ "p1 = 0.0"                                   \
        -@ "p2 = 2.0"                                   \
   -o $WORKDIR/fit.root
    
    JPlot1D \
   -f $WORKDIR/fit.root:${HISTOGRAM}               \
   -> "time [ns]"                                 \
   -\^ "number of events [au]"                     \
   -T "?"
done