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            ${TMPDIR:-/tmp}/                        
set_variable   DIR                $JPP_DIR/examples/JDetector  
set_variable   OUTPUT_FILE        $WORKDIR/tts%.root           
set_variable   NUMBER_OF_ENTRIES  10000000                     
set_variable:  FORMAT             GRAPHICS_FORMAT        gif
set_variable+  BATCH              GRAPHICS_BATCH         -B


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

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

typeset -A TTS

for OPTION in 1 38 41 67; do

    $DIR/JTransitTime \
        -o ${OUTPUT_FILE/\%/[${OPTION}]}  \
        -n 0                              \
        -O $OPTION                        \
        -n $NUMBER_OF_ENTRIES             \
        -d $DEBUG

    JFit \
        -f ${OUTPUT_FILE/\%/[${OPTION}]}:pmt            \
        -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                            >& /dev/null

    set_array  P2  `JPrintFit -f $WORKDIR/fit.root:pmt -@p2`

    TTS[${OPTION}]=$P2[1]

    JPlot1D \
        -f ${OUTPUT_FILE/\%/[${OPTION}]}:tts            \
        -f ${OUTPUT_FILE/\%/[${OPTION}]}:pmt            \
        -y "1e-5 1.0" -Y                                \
        -\> "time [ns]"                                 \
        -\^ "probability"                               \
        -L TR                                           \
        -T ""                                           \
        -o TTS_${OPTION}.$FORMAT                        $BATCH
done

JPlot1D \
    -f ${OUTPUT_FILE/\%/[1]}:tts                    \
    -f ${OUTPUT_FILE/\%/[38]}:tts                   \
    -f ${OUTPUT_FILE/\%/[41]}:tts                   \
    -f ${OUTPUT_FILE/\%/[67]}:tts                   \
    -y "1e-5 1.0" -Y                                \
    -\> "time [ns]"                                 \
    -\^ "probability"                               \
    -L TR                                           \
    -T ""                                           \
    -o TTS.$FORMAT                                  $BATCH

for key value in ${(@kv)TTS}; do
    printf "option %3d sigma %6.3f [ns]\n" $key $value
done