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