JTransitTime2D.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   DIR          $JPP_DIR/examples/JDetector  
set_variable   OUTPUT_FILE  $WORKDIR/tt2.root            


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

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

$DIR/JTransitTime2D \
    -o $OUTPUT_FILE    \
    -d $DEBUG

JFit \
    -f ${OUTPUT_FILE}:tt2                           \
    -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:tt2                        \
    -\> "time [ns]"                                 \
    -\^ "number of events [au]"                     \