JEqualizer.sh

Go to the documentation of this file.

#!/usr/bin/env 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