plot-DetectorParameters2D.sh

Go to the documentation of this file.

#!/bin/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
 
 
set_variable   DEBUG           3                                  
set_variable   WORKDIR         `pwd`                              
set_variable   OUTPUT_FILE     $WORKDIR/detector_parameters.root  
set_variable   DIR             $JPP_DIR/examples/JDetector/       
set_variable   NUMBER_OF_PMTS  31                                 
set_variable:  FORMAT          GRAPHICS_FORMAT        gif
set_variable+  BATCH           GRAPHICS_BATCH         -B
 
if do_usage $*; then
    usage "$script <detector file> (detector file)+"
fi
 
if (( $# < 2 )); then
    fatal "Wrong number of parameters."
fi
 
set_variable     DETECTOR                $argv[1]
 
INPUT_FILES=($argv[1,-1])
 
$DIR/JPlotDetectorParameters2D \
    -a "$INPUT_FILES"   \
    -o  $OUTPUT_FILE    \
    -d  $DEBUG --!
 
# parameters and plot ranges
 
typeset -A Y_RANGE
 
Y_RANGE+=(t0           "-10.0 +10.0") 
 
 
eval `JPrintDetector -a $DETECTOR -O SUMMARY`
 
kill_child_processes_at_exit
 
attach getModule -a $DETECTOR
 
typeset -Z 4 STRING
typeset -Z 2 FLOOR
 
let XMIN=-0.5
let XMAX=${NUMBER_OF_PMTS}-0.5
 
for STRING in $STRINGS[*]; do
 
    echo "Creating graphics for string $STRING."
 
    set_variable  FIRST_FLOOR  1
 
    for (( FLOOR = $FIRST_FLOOR; $FLOOR <= $LAST_FLOOR; FLOOR += 1 )); do
 
        sput $STRING $FLOOR
        sget  MODULE
 
        for KEY in ${(k)Y_RANGE}; do
 
            debug "($STRING,$FLOOR) $KEY"
 
            JPlot2D \
                -f ${OUTPUT_FILE}:${MODULE}\.${KEY}\$          \
                -z $Y_RANGE[${KEY}]                            \
                -O COLZ                                        \
                -> "index"                                    \
                -< "PMT"                                      \
                -\^ "$KEY"                                     \
                -T "(${STRING},${FLOOR})"                      \
                -o ${TMPDIR:-/tmp}/${KEY}_${STRING}_${FLOOR}.$FORMAT      $BATCH
 
            JProject2D \
                -f ${OUTPUT_FILE}:${MODULE}\.${KEY}\$          \
                -P X -y "$XMIN $XMAX"                          \
                -o ${TMPDIR:-/tmp}/project.root
 
            JScale1D \
                -f ${TMPDIR:-/tmp}/project.root:\.\*                      \
                -o ${TMPDIR:-/tmp}/scale.root                             \
                -F "1.0/${NUMBER_OF_PMTS}"                     \
                -O "nosw2"
 
            JPlot1D \
                -f ${TMPDIR:-/tmp}/scale.root:\.\*                        \
                -y $Y_RANGE[${KEY}]                            \
                -> "index"                                    \
                -\^ "<$KEY>"                                   \
                -T "(${STRING},${FLOOR})"                      \
                -o ${TMPDIR:-/tmp}/P_${KEY}_${STRING}_${FLOOR}.$FORMAT    $BATCH
 
 
        done
    done
 
    for KEY in ${(k)Y_RANGE}; do
 
        montage \
            -tile 6x3                                  \
            -geometry +0+0                             \
            ${TMPDIR:-/tmp}/${KEY}_${STRING}_{01..18}.$FORMAT     \
            $WORKDIR/${KEY}_${STRING}.$FORMAT          >& /dev/null
    done
 
    for KEY in ${(k)Y_RANGE}; do
 
        montage \
            -tile 6x3                                  \
            -geometry +0+0                             \
            ${TMPDIR:-/tmp}/P_${KEY}_${STRING}_{01..18}.$FORMAT   \
            $WORKDIR/P_${KEY}_${STRING}.$FORMAT        >& /dev/null
    done
done
 
detach