Calibrate_local.sh

Go to the documentation of this file.

#!/bin/zsh
#
# \author rgruiz
#

version=1.0
script=${0##*/}

source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null

if ( do_usage $* ); then
    usage "$script [edit the different options within]"
fi

if (( $# < 5 )); then
    fatal "Incorrect number of arguments."
fi

set_variable            DAQ_FILE   $1
set_variable       DETECTOR_FILE   $2
set_variable          INSTRUMENT   $3
set_variable              STRING   $4
set_variable          OUTPUT_DIR   $5

set_variable        TOP_PMTS   1
set_variable     BOTTOM_PMTS   4
set_variable   MAX_NEIGHBORS   4
set_variable  ANALYSIS_LEVEL   0

set_variable     PULSE_PERIOD   6250 #in units of 16 ns
set_variable   STAGGER_PERIOD    220 #in units of 16 ns
set_variable           infile      $DAQ_FILE  

cd  ${JPP_DIR}/examples/JNanobeacon

set_variable   OUTDIR_PULSE   "${OUTPUT_DIR}/out/$INSTRUMENT/JPulseFinder"

if [ ! -d ${OUTDIR_PULSE} ] 
then     
    mkdir -p ${OUTDIR_PULSE} 
fi

set_variable   OUT_FILE_PULSE   ${OUTDIR_PULSE}/$(basename ${DAQ_FILE} .root)_JPF.root

if [ ! -f $OUT_FILE_PULSE ]; then

    JPulseFinder                          \
   -f $DAQ_FILE                      \
   -a $DETECTOR_FILE                 \
   -s $STRING                        \
   -p $PULSE_PERIOD                  \
   -S $STAGGER_PERIOD                \
   -o $OUT_FILE_PULSE
else      
    echo "\nThe pulses for the file ${DAQ_FILE} already exist in file\n$OUT_FILE_PULSE\nContinue..."
fi

set_variable  OUTDIR_CALIBRATE   "${OUTPUT_DIR}/out/$INSTRUMENT/JInterDomCal"

if [ ! -d ${OUTDIR_CALIBRATE} ] 
then     
    mkdir -p ${OUTDIR_CALIBRATE} 
fi

set_variable   OUT_FILE_TXT      ${OUTDIR_CALIBRATE}/$(basename ${DAQ_FILE} .root)_JNBC.txt
set_variable   OUT_FILE_CHECKS   ${OUTDIR_CALIBRATE}/$(basename ${DAQ_FILE} .root)_JNBC.root
set_variable   OUT_FILE_DETX     ${OUTDIR_CALIBRATE}/$(basename ${DAQ_FILE} .root)_JNBC.detx

JInterDomCal            \
    -a $DETECTOR_FILE   \
    -f $OUT_FILE_PULSE  \
    -s $STRING          \
    -t $OUT_FILE_TXT    \
    -x $OUT_FILE_DETX   \
    -c $OUT_FILE_CHECKS \
    -u $TOP_PMTS        \
    -d $BOTTOM_PMTS     \
    -n $MAX_NEIGHBORS   \
    -l $ANALYSIS_LEVEL

report_name="$/calibration_report_${runs[i]}.pdf"

set_variable   OUT_FILE_PDF     ${OUTDIR_CALIBRATE}/$(basename ${DAQ_FILE} .root).pdf

echo "Writing calibration report in file:" 
echo ${OUT_FILE_PDF}
echo "------------------------------------\n"

python plot_report_calibration.py -ir $OUT_FILE_CHECKS -it $OUT_FILE_TXT -o $OUT_FILE_PDF