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    CALIBRATE_DEBUG  1
set_variable   WORKDIR  /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