JTransitTime_D0DU003MA.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
#
# \author mdejong
version=1.0
script=${0##*/}
 
# ------------------------------------------------------------------------------------------
#
#   Utility script to produce the transit time probability and generator for HAMA-R12199.
#   See https://elog.km3net.de/Calibration/329.
#
# ------------------------------------------------------------------------------------------
 
source $JPP_DIR/setenv.sh $JPP_DIR
 
set_variable   DEBUG         2
set_variable   DETID         D0DU003MA
set_variable   RUN           98
set_variable   DAQ_TIMESLICE JDAQTimeslice
PMTS=(7 15)
set_variable   TYPE          38
 
set_variable:  FORMAT        GRAPHICS_FORMAT       gif
set_variable+  BATCH         GRAPHICS_BATCH        -B
 
if do_usage $*; then
    usage "$script"
fi
 
if (( $# != 0 )); then
    fatal "Wrong number of arguments."
fi
 
expand_array RUNS
 
if [[ ! -f $DETID.datx ]]; then
    JDetectorDB  -D $DETID  -r $RUN  -o $DETID.datx
fi
 
eval `JPrintDetector -a $DETID.datx -O IDENTIFIER`
eval `JPrintDetector -a $DETID.datx -O SUMMARY`
 
set_variable  INPUT_FILE  KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}.root
 
if [[ ! -f $INPUT_FILE ]]; then
    iget /in2p3/km3net/data/raw/darkroom/KM3NeT_${(l:8::0::0:)DETECTOR_ID}/0/$INPUT_FILE
fi
 
JConvertDetectorFormat \
    -a $DETID.datx       \
    -o detector.datx
 
JPulsar \
    -a detector.datx     \
    -A                   \
    -f $INPUT_FILE       \
    -!"-1 "${^PMTS}      \
    -o pulsar.root       \
    -C $DAQ_TIMESLICE    \
    -d $DEBUG
 
JPulsar \
    -a detector.datx     \
    -f $INPUT_FILE       \
    -!"-1 "${^PMTS}      \
    -o pulsar.root       \
    -C $DAQ_TIMESLICE    \
    -d $DEBUG
 
typeset -Z 4 STRING
typeset -Z 2 FLOOR
 
for MODULE in $MODULES[*]; do
 
    getLocation -a detector.datx -M $MODULE | read STRING FLOOR
 
    if (( $FLOOR != 0 )); then
 
        JPlot1D \
            -f pulsar.root:$MODULE               \
            -x "-20.0 +50.0"                     \
            -Y                                   \
            -> "#Deltat [ns]"                   \
            -\^ "number of events [a.u.]"        \
            -T "?"                               \
            -O hist                              \
            -o M_${STRING}_${FLOOR}.$FORMAT      -B
    fi
done
 
for STRING in $STRINGS[*]; do
 
    echo "Creating graphics M_${STRING}.$FORMAT."
 
    montage \
        -tile 6x3                            \
        -geometry +0+0                       \
        M_${STRING}_{01..18}.$FORMAT         \
        M_${STRING}.$FORMAT                  >& /dev/null
done
 
JLegolas \
    -f pulsar.root                   \
    -o legolas.root                  \
    -C cdf.txt                       \
    -P pdf.txt                       \
    -d $DEBUG
 
$JPP_DIR/examples/JDetector/JTransitTime \
    -o PMT\[$TYPE\].root             \
    -O $TYPE                         \
    -d $DEBUG
 
JPlot1D \
    -f legolas.root:pmt              \
    -f PMT\[$TYPE\].root:tts         \
    -x "-30.0 +100.0"                \
    -y "2.0e-6 1.0"                  \
    -Y                               \
    -> "#Deltat [ns]"               \
    -\^ "probability"                \
    -O hist                          \
    -T "HAMA-R12199"                 \
    -L TR                            \
    -o PMT.$FORMAT                   $BATCH