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          /tmp/
set_variable  JDAQ_TIMESLICE   JDAQTimeslice

if ( do_usage $* ); then
    usage "$script <detector file> <input file> [time slice [working directory]]"
fi

case $# in
    4) set_variable  WORKDIR         $4;&
    3) set_variable  JDAQ_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 $JDAQ_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.gif