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         /tmp/

if ( do_usage $* ); then
    usage "$script [working directory]"
fi

case $# in
    1) 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              \
   -wN                               \
   -d $DEBUG

fi

if (( 1 )); then
    
    debug "Plotting results"

    JPlot2D \
   -f $WORKDIR/monitor.root:\.\*2S  \
   -OLEGO                           \
   -> "PMT pair"                   \
   -< "#DeltaT [ns]"               \
   -T ""                            \
   -o monitor.gif

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.gif

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.gif

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.gif

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.gif
   
fi