JFrodo.sh

Go to the documentation of this file.

#!/bin/zsh
script=${0##*/}
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi
 
source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
 
set_variable:  DEBUG       CALIBRATION_DEBUG      1
set_variable   WORKDIR     ${TMPDIR:-/tmp}/
set_variable   FORMULA     "[0] * exp([1]*sqrt(x) + [2]*x)  +  [3]"
set_variable   XMIN        -0.5
set_variable   XMAX        255.5
set_variable:  FORMAT      GRAPHICS_FORMAT        gif
set_variable+  BATCH       GRAPHICS_BATCH         -B
 
if do_usage $*; then
    usage "$script <input file> [working directory]"
fi
 
case $# in
    2)  set_variable WORKDIR      $2;&
    1)  set_variable INPUT_FILE   $1;;
    *)  fatal "Wrong number of arguments.";;
esac
 
typeset -a X
 
X=(-0.5  3.5  6.5  9.5 12.5 15.5 18.5 
   20.5 21.5 22.5 23.5 24.5 25.5 26.5 27.5 28.5 29.5
   30.5 31.5 32.5 33.5 34.5 35.5 36.5 37.5 38.5 39.5
   40.5 42.5 44.5 46.5 48.5
   50.5 52.5 54.5 56.5 58.5
   60.5 65.5
   70.5 75.5
   80.5 85.5
   90.5 95.5
   100.5 120.5 140.5 160.5 180.5 200.5 250.5)
 
JRebin1D \
    -f ${INPUT_FILE}:hb                    \
    -x "$X[*]"                             \
    -o $WORKDIR/rebin.root                 
 
JPlot1D \
    -f $WORKDIR/rebin.root:hb              \
    -y "-0.2 +1.2"                         \
    -> "time-over-threshold [ns]"         \
    -\^ "#partialL/#partialt"              \
    -T ""                                  \
    -o $WORKDIR/gradient.$FORMAT           $BATCH
 
JFrodo \
    -f ${INPUT_FILE}                       \
    -F Gauss                               \
    -o $WORKDIR/frodo.root                 \
    -T "-20.5 +10.5"                       \
    -O LE                                  \
    -x "$XMIN $XMAX"                       \
    -w                                     \
    -d $DEBUG
 
JPlot1D \
    -f $WORKDIR/frodo.root:h0              \
    -y "-10.0 +15.0"                       \
    -> "time-over-threshold [ns]"         \
    -\^ "#Deltat [ns]"                     \
    -T "Time-slewing correction"           \
    -o $WORKDIR/slewing.$FORMAT            $BATCH