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