JSlewingK40.sh

Go to the documentation of this file.

#!/bin/zsh
#
#
# \author mdejong
#
version=1.0
script=${0##*/}
 
# ------------------------------------------------------------------------------------------
#
#                         Utility script to test JSlewingK40.
#
# ------------------------------------------------------------------------------------------
 
 
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     .  
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  TMAX            3.0
set_variable  DETECTOR        $WORKDIR/detector.detx
set_variable  BACKGROUND_HZ   5e3 500 50 5 0.5
set_variable  numberOfSlices  10000
set_variable  INPUT_FILE      $WORKDIR/timeslice.root
set_variable  OUTPUT_FILE     $WORKDIR/slewing\[Monte\ Carlo\].root
 
 
if (( 1 )); then
 
    detectorType=1
    numberOfLines=1
    numberOfFloors=1
    
    JDetector \
        -D 1                                   \
        -@detectorType=$detectorType           \
        -@numberOfLines=$numberOfLines         \
        -@numberOfFloors=$numberOfFloors       \
        -o $DETECTOR -V V2                     \
        -d $DEBUG    
 
fi
 
if (( 1 )); then
 
    if [ -f "$WORKDIR/timeslice.root" ]; then
        rm -i $WORKDIR/timeslice.root
    fi
 
    if [ ! -f "$WORKDIR/timeslice.root" ]; then
 
        JRandomTimesliceWriter \
            -a $DETECTOR                           \
            -o $INPUT_FILE                         \
            -B "$BACKGROUND_HZ"                    \
            -n $numberOfSlices                     \
            -P %.slewing=1                         \
            -d $DEBUG --!
 
    fi
fi
 
 
if (( 1 )); then
    
    echo "Monitoring data."
    echo "Detector: $DETECTOR"
    
    JSlewingK40 \
        -f $WORKDIR/timeslice.root                 \
        -a $DETECTOR                               \
        -o $OUTPUT_FILE                            \
        -t "20 35"                                 \
        -d $DEBUG
    
fi
 
if (( 1 )); then
 
    JScale1D \
        -f ${OUTPUT_FILE}:h0                       \
        -f ${OUTPUT_FILE}:h1                       \
        -F "1.0/GetEntries"                        \
        -O "nosw2"                                 \
        -o $WORKDIR/scale.root --!
 
fi
 
if (( 1 )); then
 
    JFit \
        -f ${WORKDIR}/scale.root:h0                \
        -o ${WORKDIR}/fit.root                     \
        -F "[0] * exp(-0.5 * (x-[1])*(x-[1]) / ([2]*[2])) + [3]"  \
        -@"p0 = GetMaximum"                        \
        -@"p1 =  0.0"                              \
        -@"p2 = 0.5*GetRMS"                        \
        -@"p3 = GetMinimum"                        \
        -d3
 
    JPlot1D \
        -f ${WORKDIR}/fit.root:h0                  \
        -> "#Deltat [ns]"                         \
        -\^ "number of hits [au]"                  \
        -o resolution.$FORMAT                      $BATCH
 
fi
 
if (( 1 )); then
 
    JPlot1D \
        -f $WORKDIR/scale.root:h1                  \
        -> "time over threshold [ns]"             \
        -\^ "number of hits [au]"                  \
        -o tot.$FORMAT                             $BATCH
 
fi
 
if (( 1 )); then
 
    JFit \
        -f ${OUTPUT_FILE}:h2                        \
        -o ${WORKDIR}/fit.root                      \
        -F "[0]*exp([1]*sqrt(x) + [2]*x) + [3]"     \
        -@"p0 = GetMaximum"                         \
        -@"p1 = -0.01"                              \
        -@"p2 = -0.05"                              \
        -@"p3 = GetMinimum"                         \
        -d3
 
    JPlot1D \
        -f ${WORKDIR}/fit.root:h2                   \
        -y "-10 +10"                                \
        -> "TOT [ns]"                              \
        -\^ "#Deltat [ns]"                          \
        -o slewing.$FORMAT                          $BATCH
 
fi