JCalibrateK40.sh

Go to the documentation of this file.

#!/bin/zsh
#
#
# \author mdejong
#
version=1.0
script=${0##*/}
 
# ------------------------------------------------------------------------------------------
#
#                         Utility script to test JCalibrateK40/JFitK40.
#
# ------------------------------------------------------------------------------------------
 
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     ${TMPDIR:-/tmp}/  
set_variable:  FORMAT      GRAPHICS_FORMAT        gif
set_variable+  BATCH       GRAPHICS_BATCH         -B
 
if do_usage $*; then
    usage "$script [working directory]"
fi
 
case $# in
    1) set_variable  WORKDIR  $1;;
esac
 
set_variable  DETECTOR        $WORKDIR/detector.detx
set_variable  BACKGROUND_HZ   5e3 440 50 5 0.5
set_variable  QE              1.0
set_variable  numberOfSlices  6000
 
 
if (( 1 )); then
 
    set_variable  detectorType     1
    set_variable  numberOfLines    1
    set_variable  numberOfFloors   1
    
    JDetector \
        -D 1                                   \
        -@ "detectorType   = $detectorType"    \
        -@ "numberOfLines  = $numberOfLines"   \
        -@ "numberOfFloors = $numberOfFloors"  \
        -o $DETECTOR -V V2                     \
        -d $DEBUG    
 
fi
 
if (( 1 )); then
 
    debug "Generating random background."
 
    if ( ! reuse_file $WORKDIR/timeslice.root ); then
 
        JRandomTimesliceWriter \
            -a $DETECTOR                       \
            -o $WORKDIR/timeslice.root         \
            -B "$BACKGROUND_HZ"                \
            -n $numberOfSlices                 \
            -P "pmt= 101  0 QE=$QE"            \
            -P "pmt= 101  1 QE=$QE"            \
            -P "pmt= 101  2 QE=$QE"            \
            -P "pmt= 101  1 QE=$QE"            \
            -P "pmt= 101  3 QE=$QE"            \
            -P "pmt= 101  4 QE=$QE"            \
            -P "pmt= 101  5 QE=$QE"            \
            -P "pmt= 101  6 QE=$QE"            \
            -P "pmt= 101  7 QE=$QE"            \
            -P "pmt= 101  8 QE=$QE"            \
            -P "pmt= 101  9 QE=$QE"            \
            -P "pmt= 101 10 QE=$QE"            \
            -P "pmt= 101 11 QE=$QE"            \
            -P "pmt= 101 12 QE=$QE"            \
            -P "pmt= 101 11 QE=$QE"            \
            -P "pmt= 101 13 QE=$QE"            \
            -P "pmt= 101 14 QE=$QE"            \
            -P "pmt= 101 15 QE=$QE"            \
            -P "pmt= 101 16 QE=$QE"            \
            -P "pmt= 101 17 QE=$QE"            \
            -P "pmt= 101 18 QE=$QE"            \
            -P "pmt= 101 19 QE=$QE"            \
            -P "pmt= 101 20 QE=$QE"            \
            -P "pmt= 101 21 QE=$QE"            \
            -P "pmt= 101 22 QE=$QE"            \
            -P "pmt= 101 21 QE=$QE"            \
            -P "pmt= 101 23 QE=$QE"            \
            -P "pmt= 101 24 QE=$QE"            \
            -P "pmt= 101 25 QE=$QE"            \
            -P "pmt= 101 26 QE=$QE"            \
            -P "pmt= 101 27 QE=$QE"            \
            -P "pmt= 101 28 QE=$QE"            \
            -P "pmt= 101 29 QE=$QE"            \
            -P "pmt= 101 30 QE=$QE"            \
            -d $DEBUG --!
 
    fi
fi
 
if (( 1 )); then
 
    debug "Processing data."
 
    if ( ! reuse_file $WORKDIR/trigger_processor.root ); then
 
        JTriggerProcessor \
            -a $DETECTOR                       \
            -f $WORKDIR/timeslice.root         \
            -o $WORKDIR/trigger_processor.root \
            -@"ctMin=-1.0"                     \
            -@"TMaxLocal_ns=50"                \
            -@"writeL1=1"                      \
            -d $DEBUG
 
    fi
fi
 
if (( 1 )); then
    
    debug "Monitoring data."
    
    JCalibrateK40 \
        -f $WORKDIR/trigger_processor.root \
        -C JDAQTimesliceL1                 \
        -a $DETECTOR                       \
        -o $WORKDIR/monitor.root           \
        -d $DEBUG
 
    debug "Merging data."
 
    JMergeCalibrateK40 \
        -f $WORKDIR/monitor.root          \
        -o $WORKDIR/merge.root            \
        -d $DEBUG       
    
    debug "Fitting data."
    debug "Detector: $DETECTOR"
    
    JFitK40 \
        -a $DETECTOR                      \
        -f $WORKDIR/merge.root            \
        -o $WORKDIR/fit.root              \
        -w                                \
        -d $DEBUG
 
fi
 
if (( 1 )); then
    
    debug "Plotting results"
 
    JPlot2D \
        -f $WORKDIR/monitor.root:\.\*2S  \
        -OLEGO                           \
        -> "PMT pair"                   \
        -< "#DeltaT [ns]"               \
        -T ""                            \
        -o monitor.$FORMAT               $BATCH
 
fi 
 
if (( 1 )); then
    
    debug "Plotting standard deviations"
 
    JPlot2D \
        -f $WORKDIR/fit.root:\.\*2F      \
        -OLEGO                           \
        -> "PMT pair"                   \
        -< "#DeltaT [ns]"               \
        -z "-3 +3"                       \
        -T "standard deviation"          \
        -o stdev.$FORMAT                 $BATCH
 
fi 
 
if (( 1 )); then
    
    debug "Plotting results for t0"
 
    JPlot1D \
        -f $WORKDIR/fit.root:\.\*1t0     \
        -y "-10 +10"                     \
        -> "PMT"                        \
        -\^ "#DeltaT [ns]"               \
        -L TR                            \
        -T "time calibration"            \
        -o t0.$FORMAT                    $BATCH
 
fi
 
if (( 1 )); then
 
    debug "Plotting results for sigma"
 
    JPlot1D \
        -f $WORKDIR/fit.root:\.\*1sigma  \
        -y "0 5"                         \
        -> "PMT"                        \
        -\^ "#sigma [ns]"                \
        -L TR                            \
        -T "time resolution"             \
        -o sigma.$FORMAT                 $BATCH
 
fi
 
 
if (( 1 )); then
    
    debug "Plotting results for QE"
 
    JPlot1D \
        -f $WORKDIR/fit.root:\.\*1QE     \
        -y "0 2"                         \
        -> "PMT"                        \
        -\^ "QE"                         \
        -L TR                            \
        -T "QE"                          \
        -o QE.$FORMAT                    $BATCH
        
fi