JSquid.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
script=${0##*/}
 
source $JPP_DIR/setenv.sh $JPP_DIR
 
zmodload zsh/mathfunc
 
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 <detector file> <input file>"
fi
 
if (( $# != 2 )); then
    fatal "Wrong number of arguments."
fi
 
set_variable DETECTOR    $argv[1]
set_variable INPUT_FILE  $argv[2]
 
eval `JPrintDetector -a $DETECTOR -O SUMMARY`
 
JPrintDAQHeader -f $INPUT_FILE  |  read -A DAQ_HEADER
 
set_variable  DETECTOR_ID    ${DAQ_HEADER[1]}
set_variable  RUN            ${DAQ_HEADER[2]}
 
if ( ! reuse_file $WORKDIR/trigger_processor.root ); then
 
    JTriggerProcessor \
        -a $DETECTOR                                      \
        -f $INPUT_FILE                                    \
        -C JDAQTimesliceSN                                \
        -@ "trigger3DMuon.enabled         =  1"           \
        -@ "trigger3DMuon.numberOfHits    =  3"           \
        -@ "trigger3DMuon.numberOfModules =  3"           \
        -@ "trigger3DMuon.gridAngle_deg   = 10"           \
        -@ "trigger3DMuon.TMaxExtra_ns    = 500.0e3"      \
        -@ "TMaxEvent_ns    = 500.0e3"                    \
        -@ "TMaxLocal_ns    = 20.0"                       \
        -@ "L2.numberOfHits = 4"                          \
        -@ "L2.TMaxLocal_ns = 20.0"                       \
        -@ "L2.ctMin        = 0.0"                        \
        -o $WORKDIR/trigger_processor.root                \
        -d $DEBUG
fi
 
JSquid \
    -a $DETECTOR                                          \
    -f $WORKDIR/trigger_processor.root                    \
    -o $WORKDIR/squid.root                                \
    -d $DEBUG
 
for STRING in $STRINGS[*]; do
 
    JPlot1D \
        -f $WORKDIR/squid.root:"^[ ]*${STRING}"           \
        -> "#Deltat [ns]"                                \
        -\^ "number of hits [a.u.]"                       \
        -T "$RUN"                                         \
        -N "X 505"                                        \
        -L TL                                             \
        -o $WORKDIR/T${(l:4::0::0:)STRING}.gif            $BATCH
done
 
typeset -i NX
typeset -i NY
 
let "NY = sqrt(${#STRINGS}) * 0.75 + 0.25"
let "NX = (${#STRINGS} + $NY - 1)/ $NY"
 
montage \
    -tile ${NX}x${NY}                                     \
    -geometry +0+0                                        \
    $WORKDIR/T[0-9][0-9][0-9][0-9].$FORMAT                \
    ${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}.$FORMAT >& /dev/null