JFitL1dt.sh

Go to the documentation of this file.

#!/usr/bin/env 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         ${TMPDIR:-/tmp}/            
set_variable   NUMBEROFFLOORS  18               
set_variable   DAQ_TIMESLICE   JDAQTimesliceSN  # SN has the most statistics in data since there is 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 $DAQ_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 $DAQ_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