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         .

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

case $# in
    1) 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/GetEntries"                          \
        -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.gif

fi

if (( 1 )); then

    JPlot1D \
        -f $WORKDIR/scale.root:h1                  \
        -> "time over threshold [ns]"             \
        -\^ "number of hits [au]"                  \
        -o tot.gif

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

fi