JDetector/plot-PMTParameters.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
#
#
# \author mdejong
#
version=1.0
script=${0##*/}
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi
 
source $JPP_DIR/setenv.sh $JPP_DIR
 
# parameters and plot ranges
 
typeset -A Y_RANGE
 
Y_RANGE+=(EFFICIENCY   " 0.0   1.5") 
Y_RANGE+=(QE           " 0.0   1.5") 
Y_RANGE+=(gain         " 0.0   1.5") 
Y_RANGE+=(gainSpread   " 0.0   1.0") 
Y_RANGE+=(riseTime_ns  " 0.0  15.0") 
Y_RANGE+=(TTS_ns       "-5.0  +5.0") 
Y_RANGE+=(threshold    " 0.0   1.0") 
Y_RANGE+=(slope        " 0.0  15.0") 
Y_RANGE+=(saturation   " 0.0 250.0") 
Y_RANGE+=(slewing      " 0.0   1.5") 
 
set_variable   DEBUG           2                            
set_variable   WORKDIR         `pwd`                         
set_variable   OUTPUT_FILE     $WORKDIR/pmt_parameters.root  
set_variable   DIR             $JPP_DIR/examples/JDetector/  
set_variable:  FORMAT          GRAPHICS_FORMAT        gif
set_variable+  BATCH           GRAPHICS_BATCH         -B
 
if do_usage $*; then
    usage "$script <detector file> <PMT parameters file> [parameter+]"\
          "\nThe list of optional parameters includes: ${(k)Y_RANGE}."
fi
 
if (( $# < 2 )); then
    fatal "Wrong number of parameters."
fi
 
set_variable     DETECTOR                $argv[1]
set_variable     INPUT_FILE              $argv[2]
 
if (( $# > 2 )); then
    KEYS=($argv[3,-1])
else
    KEYS=(${(k)Y_RANGE})
fi
 
$DIR/JPlotPMTParameters \
    -a $DETECTOR     \
    -P $INPUT_FILE   \
    -o $OUTPUT_FILE  \
    -d $DEBUG -A
 
eval `JPrintDetector -a $DETECTOR -O SUMMARY`
 
set_variable  FIRST_FLOOR  1
 
 
 
typeset -Z 4 STRING
typeset -Z 2 FLOOR
 
for STRING in $STRINGS[*]; do
 
    echo -n "Creating graphics for string $STRING.."
 
    MODULES=()
 
    for (( FLOOR = $FIRST_FLOOR; $FLOOR <= $LAST_FLOOR; FLOOR += 1 )); do
 
        eval MODULE=\${STRING_${(l:4::0::0:)STRING}\[$(($FLOOR))\]:--1}
 
        MODULES+=($MODULE)
 
        for KEY in $KEYS[*]; do
 
            debug "($STRING,$FLOOR) $KEY"
            
            JPlot1D \
                -f ${OUTPUT_FILE}:${MODULE}\.${KEY}\$                   \
                -y $Y_RANGE[${KEY}]                                     \
                -> "PMT"                                               \
                -\^ "$KEY"                                              \
                -T "(${STRING},${FLOOR})"                               \
                -o ${TMPDIR:-/tmp}/${KEY}_${STRING}_${FLOOR}.$FORMAT    $BATCH
        done
    done
 
    for KEY in $KEYS[*]; do
 
        JPlot1D \
            -f${OUTPUT_FILE}:${^MODULES}\.${KEY}\$                  \
            -y $Y_RANGE[${KEY}]                                     \
            -> "PMT"                                               \
            -\^ "$KEY"                                              \
            -T "${STRING}"                                          \
            -o $WORKDIR/${KEY}.1S_${STRING}.$FORMAT                 $BATCH
 
        montage \
            -tile 6x3                                           \
            -geometry +0+0                                      \
            ${TMPDIR:-/tmp}/${KEY}_${STRING}_{01..18}.$FORMAT   \
            $WORKDIR/${KEY}_${STRING}.$FORMAT                   >& /dev/null
    done
 
    echo " ->" $WORKDIR/${^KEYS}\[.1S\]_${STRING}.$FORMAT
 
    rm -f ${TMPDIR:-/tmp}/${^KEYS}_${STRING}_*.$FORMAT
done