JTransitTime_HAMA-R12199.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
#
# \author mdejong
version=1.0
script=${0##*/}
 
# ------------------------------------------------------------------------------------------
#
#   Utility script to produce the transit time probability and generator for HAMA-R12199.
#
# ------------------------------------------------------------------------------------------
 
source $JPP_DIR/setenv.sh $JPP_DIR
 
set_variable   DEBUG         2
set_variable   INPUT_FILE    $JPP_DIR/examples/JDetector/R12199_delayed_pulses_hist.txt
set_variable   TYPE          1
 
set_variable:  FORMAT        GRAPHICS_FORMAT        gif
set_variable+  BATCH         GRAPHICS_BATCH         -B
 
if do_usage $*; then
    usage "$script [input file]"
fi
 
if   (( $# == 1 )); then
    set_variable INPUT_FILE  $1
elif (( $# != 0 )); then
    fatal "Wrong number of arguments."
fi
 
# convert control characters to newline
 
set_variable     TEMPORARY_FILE  ${TMPDIR:-/tmp}/__${INPUT_FILE##*/}__
 
cat $INPUT_FILE | tr '\r' '\n' > $TEMPORARY_FILE
 
$JPP_DIR/examples/JDetector/JTTS \
    -f $TEMPORARY_FILE \
    -o histogram.root  \
    -P pdf.txt         \
    -C cdf.txt         \
    -d $DEBUG
 
rm -f $TEMPORARY_FILE
 
$JPP_DIR/examples/JDetector/JTransitTime \
    -o PMT\[$TYPE\].root             \
    -O $TYPE                         \
    -d $DEBUG
 
JPlot1D \
    -f histogram.root:pmt            \
    -f histogram.root:g1             \
    -f PMT\[$TYPE\].root:tts         \
    -> "#Deltat [ns]"               \
    -y "1e-5 1e-0" -Y                \
    -T "HAMA-R12199"                 \
    -L TR                            \
    -o PMT.$FORMAT                   $BATCH