JEqualizer.sh

Go to the documentation of this file.

#!/bin/zsh
#
#
# \author mdejong
#
version=1.0
script=${0##*/}

if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi

source $JPP_DIR/setenv.sh $JPP_DIR

set_variable   DEBUG          2                                 
set_variable   WORKDIR        `pwd`                             
set_variable   DETECTOR       $WORKDIR/detectors/detector.detx  
set_variable   DAQ_TIMESLICE  JDAQTimesliceL0
set_variable   OPTION         ABCDEF                            

if do_usage $*; then
    usage "$script <input file> [option]"
fi

case $# in
    2)  set_variable OPTION       $argv[2];&
    1)  set_variable INPUT_FILE   $argv[1];;
    *)  fatal "Wrong number of arguments.";;
esac


JPrintDAQHeader -f $INPUT_FILE | read DETECTOR_ID RUN FRAME_INDEX TIMESLICE_START

set_variable     DIR                    $WORKDIR/$RUN

mkdir -p $DIR

set_variable     DETECTOR               $DIR/detector.detx
set_variable     PMT_FILE               $DIR/PMT.txt

if [[ ! -f $DETECTOR ]]; then
    JDetectorDB -D $DETECTOR_ID -r $RUN -o $DETECTOR
fi

# background mode

if [[ ! -f $PMT_FILE ]]; then

    # gain calibration

    set_variable     MONITOR_FILE           $DIR/monitor_tot.root
    set_variable     MERGE_FILE             $DIR/merge_tot.root
    set_variable     FIT_FILE               $DIR/fit_tot.root
    set_variable     LOG_FILE               $DIR/fit_tot.log

    JCalibrateToT \
        -a $DETECTOR                        \
        -f $INPUT_FILE                      \
        -C $DAQ_TIMESLICE                   \
        -o $MONITOR_FILE                    \
        -d $DEBUG --!

    JMergeCalibrateToT \
        -f $MONITOR_FILE                    \
        -o $MERGE_FILE                      \
        -d $DEBUG --!

    JFitToT \
        -a $DETECTOR                        \
        -f $MERGE_FILE                      \
        -o $FIT_FILE                        \
        -w                                  \
        -P $PMT_FILE                        \
        -d $DEBUG --!                       

    # QE calibration

    set_variable     MONITOR_FILE           $DIR/monitor_k40.root
    set_variable     MERGE_FILE             $DIR/merge_k40.root
    set_variable     FIT_FILE               $DIR/fit_k40.root
    set_variable     LOG_FILE               $DIR/fit_k40.log
    set_variable     DETX_FILE              $DIR/fit_k40.detx

    JCalibrateK40 \
        -a $DETECTOR                        \
        -f $INPUT_FILE                      \
        -C $DAQ_TIMESLICE                   \
        -o $MONITOR_FILE                    \
        -d $DEBUG --!

    JMergeCalibrateK40 \
        -f $MONITOR_FILE                    \
        -o $MERGE_FILE                      \
        -d $DEBUG --!
    
    cp $DETECTOR $DETX_FILE
    
    JFitK40 \
        -a $DETX_FILE                       \
        -f $MERGE_FILE                      \
        -o $FIT_FILE                        \
        -wA                                 \
        -P $PMT_FILE                        \
        -d $DEBUG --!                       

    # final corrections

    set_variable  MU            0.20
    set_variable  TOT_NS        4.0  255.0

    JEditPMTParameters \
        -P $PMT_FILE                        \
        -o $PMT_FILE                        \
        -T "$TOT_NS"                        \
        -E $MU

    exit
fi


# interactive mode

set_variable     OUTPUT_FILE            $DIR/equalizer-${OPTION}.root

if ( ! reuse_file $OUTPUT_FILE ); then

    $JPP_DIR/examples/JCalibrate/JEqualizer \
        -f $INPUT_FILE                      \
        -P $PMT_FILE                        \
        -o $OUTPUT_FILE                     \
        -r $OPTION                          \
        -d $DEBUG
fi