test-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
 
 
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