JFitL1dt.sh

Go to the documentation of this file.

#!/bin/zsh
#
# \author lnauta
# 
script=${0##*/}

# ------------------------------------------------------------------------------------------
#
#                Script to fit L1dt time offsets from triggered data and MC data.
#
# ------------------------------------------------------------------------------------------

if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi
source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null

set_variable  DEBUG            0
set_variable  WORKDIR          /tmp/
set_variable  NUMBEROFFLOORS   18
set_variable  JDAQ_TIMESLICE   JDAQTimesliceSN # SN has the most statistics in data since theres no downscaling

if ( do_usage $* ); then
    usage "$script <detector file> <input file(s)> <model file(s)> [number of floors in fit]"
fi

case $# in
    4) set_variable  NUMBEROFFLOORS  $4;&
    3) set_variable  MODELFILE       $3;
       set_variable  INPUTFILE       $2;
       set_variable  DETECTOR        $1;;
    *) fatal "Wrong number of arguments.";;
esac

cp -v $DETECTOR $WORKDIR/detector.detx

# ------------------------------------------------------------------------------------------
#                                   JMonitorL1dt for data
# ------------------------------------------------------------------------------------------

timer_start

JMonitorL1dt \
    -f $INPUTFILE                 \
    -C $JDAQ_TIMESLICE            \
    -a $WORKDIR/detector.detx     \
    -o $WORKDIR/monitor_data.root \
    -d $DEBUG --!

timer_stop
timer_print

# ------------------------------------------------------------------------------------------
#                                   JMonitorL1dt for monte carlo
# ------------------------------------------------------------------------------------------

timer_start

JMonitorL1dt \
    -f $MODELFILE                  \
    -C $JDAQ_TIMESLICE             \
    -a $WORKDIR/detector.detx      \
    -o $WORKDIR/monitor_model.root \
    -d $DEBUG --!

timer_stop
timer_print

# ------------------------------------------------------------------------------------------
#                                          JFitL1dtSlices
# ------------------------------------------------------------------------------------------

timer_start
JFitL1dtSlices \
    -a $WORKDIR/detector.detx      \
    -f $WORKDIR/monitor_data.root  \
    -m $WORKDIR/monitor_model.root \
    -o $WORKDIR/slices.root        \
    -d $DEBUG --!

timer_stop
timer_print

# ------------------------------------------------------------------------------------------
#                                             JFitL1dt
# ------------------------------------------------------------------------------------------

timer_start

JFitL1dt \
    -D                        \
    -F $NUMBEROFFLOORS        \
    -a $WORKDIR/detector.detx \
    -f $WORKDIR/slices.root   \
    -o $WORKDIR/fitl1dt.root  \
    -d $DEBUG --!

timer_stop
timer_print