test-PMTParameters.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        PMT_PARAMETERS_DEBUG          2
set_variable   WORKDIR      `pwd`                         
set_variable   OUTPUT_FILE  ${TMPDIR:-/tmp}/pmt_parameters.root      
set_variable   DIR          $JPP_DIR/examples/JDetector/  


if do_usage $*; then
    usage "$script <detector file> <PMT parameters file>"
fi

case $# in
    2)  set_variable DETECTOR     $1;
        set_variable INPUT_FILE   $2;;
    *)  fatal "Wrong number of parameters."
esac


# acceptance criteria

typeset -A ACCEPTANCE_WINDOW
typeset -A NUMBER_OF_OUTLIERS

ACCEPTANCE_WINDOW+=(QE           "0.8   1.2") 
ACCEPTANCE_WINDOW+=(gain         "0.8   1.2") 
ACCEPTANCE_WINDOW+=(gainSpread   "0.2   0.4") 

NUMBER_OF_OUTLIERS+=(QE           1)
NUMBER_OF_OUTLIERS+=(gain         1)
NUMBER_OF_OUTLIERS+=(gainSpread   1)


$DIR/JPlotPMTParameters \
    -a $DETECTOR     \
    -P $INPUT_FILE   \
    -o $OUTPUT_FILE  \
    -d 0 --!

for KEY in ${(k)ACCEPTANCE_WINDOW}; do

    debug "Test $KEY"

    JWindowTest1D \
        -f ${OUTPUT_FILE}:\.\*${KEY}\$    \
        -y ${ACCEPTANCE_WINDOW[${KEY}]}   \
        -N ${NUMBER_OF_OUTLIERS[${KEY}]}  \
        -d ${DEBUG} 

done