JLegolas.sh

Go to the documentation of this file.

#!/bin/zsh
 
version=1.0
script=${0##*/}
 
# ------------------------------------------------------------------------------------------
#
#             Utility script to create PDF and CDF of transition-time distribution.
#
# ------------------------------------------------------------------------------------------
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi
 
source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
 
set_variable   DEBUG           2              
set_variable   WORKDIR         ${TMPDIR:-/tmp}/          
set_variable   DAQ_TIMESLICE   JDAQTimeslice
set_variable:  FORMAT          GRAPHICS_FORMAT       gif
set_variable+  BATCH           GRAPHICS_BATCH        -B
 
if do_usage $*; then
    usage "$script <detector file> <input file> [time slice [working directory]]"\
          "\nUtility script to create PDF and CDF of transition-time distribution."
fi
 
case $# in
    4) set_variable  WORKDIR         $4;&
    3) set_variable  DAQ_TIMESLICE   $3;&
    2) set_variable  INPUT_FILE      $2;
       set_variable  DETECTOR        $1;;
    *) fatal "Wrong number of arguments."
esac
 
cp -v $DETECTOR $WORKDIR/detector.detx
 
for (( i = 0; $i != 2; ++i )); do
 
    timer_start
 
    JPulsar \
        -f $INPUT_FILE                \
        -C $DAQ_TIMESLICE             \
        -a $WORKDIR/detector.detx     \
        -A                            \
        -! "-1  7"                    \
        -! "-1 15"                    \
        -o $WORKDIR/pulsar.root       \
        -d $DEBUG --!
 
    timer_stop
    timer_print
 
done
 
 
JLegolas \
    -f $WORKDIR/pulsar.root       \
    -P $WORKDIR/pdf.txt           \
    -C $WORKDIR/cdf.txt           \
    -o $WORKDIR/legolas.root      \
    -d $DEBUG --!
 
JPlot1D \
    -f $WORKDIR/legolas.root:pmt  \
    -f $WORKDIR/legolas.root:g1   \
    -> "#Deltat [ns]"            \
    -CC                           \
    -y "1e-6 1e-0" -Y             \
    -T "PMT"                      \
    -o $WORKDIR/tts.$FORMAT       $BATCH