JCalibrateHeight.sh

Go to the documentation of this file.

#!/bin/zsh
script=${0##*/}

if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi

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

set_variable  DEBUG               ${CALIBRATE_DEBUG:-1}
set_variable  WORKDIR             `pwd`
set_variable  QUEUE_NAME          /tmp/queue
set_variable  QUEUE_SIZE          3    

APPLICATIONS=(prefit simplex gandalf)      # fit applications

if ( do_usage $* ); then
    usage "$script <detector file> (input file)+"
fi

if (( $# < 3 )); then
    fatal "Not enough arguments."
fi
    
set_variable  DETECTOR            $argv[1]

INPUT_FILES=($argv[2,-1])

if ( ! is_CCLyon ); then
    
    # local batch processor

    source qlib.sh
    
    $JPP_DIR/examples/scripts/mkqueue.sh  $QUEUE_NAME  $QUEUE_SIZE  continue

else
    fatal "The job submission should be changed for running in CC-Lyon."
fi

if [[ $DETECTOR != *\.detx ]]; then
    fatal "Invalid detector file <${DETECTOR}>."
fi

# set path

for (( i=1; i <= ${#INPUT_FILES}; i+=1 )); do

    set_variable  INPUT_FILE   ${INPUT_FILES[${i}]}

    if [[ -f `pwd`/$INPUT_FILE ]]; then
   INPUT_FILES[${i}]=`pwd`/$INPUT_FILE
    fi
done


for TIME_OFFSET_NS in +20.4; do
    for MULTIPLICATION in 0.94 0.96 0.97 0.98 0.99 1.00 1.01 1.02 1.03 1.04 1.06; do

   set_variable  OUTPUT_DIRECTORY  $WORKDIR/T=${TIME_OFFSET_NS}:M=${MULTIPLICATION}

   mkdir -p --  $OUTPUT_DIRECTORY

   # detector calibration

   if [[ ! -f $OUTPUT_DIRECTORY/detector.detx ]]; then
           
       JEditDetector \
           -a $DETECTOR                                             \
           -S "1 add $TIME_OFFSET_NS"                               \
           -S "1 mul $MULTIPLICATION"                               \
           -S "2 mul $MULTIPLICATION"                               \
           -o $OUTPUT_DIRECTORY/detector.detx --!
       
   fi
           
   for FILE_NAME in $INPUT_FILES[*]; do
               
       set_variable  APP          ${APPLICATIONS[-1]}
       set_variable  INPUT_FILE   ${FILE_NAME}
       set_variable  OUTPUT_FILE  ${OUTPUT_DIRECTORY}/${${FILE_NAME##*/}/\.root/_${APP}.root}
       set_variable  BATCH_FILE   ${OUTPUT_DIRECTORY}/${${FILE_NAME##*/}/\.root/.sh}
       
       rm -f        $BATCH_FILE
               
       echo "#!/bin/zsh"  >>  $BATCH_FILE

       if [[ ! -f ${OUPUT_FILE} ]]; then

           # sequentially processing of applications
                   
           for APP in $APPLICATIONS[*]; do
                       
               set_variable  OUTPUT_FILE  ${OUTPUT_DIRECTORY}/${${FILE_NAME##*/}/\.root/_${APP}.root}
                       
               echo "JARCA${(C)APP}.sh  $OUTPUT_DIRECTORY/detector.detx  $INPUT_FILE $OUTPUT_FILE >& /dev/null"  >>  $BATCH_FILE
                       
               set_variable  INPUT_FILE   ${OUTPUT_FILE}
                       
           done

           # cleanup
                   
           for (( i=1; i < ${#APPLICATIONS}; i+=1 )); do
                       
               set_variable  APP          ${APPLICATIONS[${i}]}
               set_variable  OUTPUT_FILE  ${OUTPUT_DIRECTORY}/${${FILE_NAME##*/}/\.root/_${APP}.root}
                       
               echo "rm -f  $OUTPUT_FILE"  >>  $BATCH_FILE
                       
           done

           echo "rm -f $BATCH_FILE"  >>  $BATCH_FILE
                   
           # submit job

           chmod +x    $BATCH_FILE
           puts_queue  $QUEUE_NAME  $BATCH_FILE

       fi
   done
    done
done