pre-calibration.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
#
 
script=${0##*/}
 
echo JPP_DIR IS $JPP_DIR
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi
 
source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
 
set_variable   DEBUG                    2
set_variable   WORKDIR                  `pwd`
set_variable   ACOUSTICS_THREADS        32
 
if do_usage $*; then
    usage "$script [version]"\
          "\nAuxiliary script to run pre-calibrations of all detectors,"\
          "\nincluding automatic disabling of 'bad' transmissions (see JAcousticsDisable.sh)."
fi
 
if (( $# == 1 )); then
    VERSION=${1}
else
    fatal "Wrong number of arguments."
fi
 
WORKDIR=`pwd`
# -----------------------------------------------------------------------------------------------
# Run pre-calibration of given detector with optional inputs.
# The last directory corresponds to the target detector and
# the optional preceding directories to the source detectors.
#
# \param 1-N   directories
# -----------------------------------------------------------------------------------------------
fuction run()
{
    DETECTORS=($argv[1,-2])     
    DETECTOR=$argv[-1]
 
    if [[ ! -f $JPP_DIR/software/JAcoustics/pre-calibration_${DETECTOR}.sh ]]; then
        fatal "No calibration for detector $DETECTOR."
    fi
 
    DIR=$WORKDIR/${DETECTOR}.${VERSION}.A     # without disabled transmissions
 
    if [[ ! -d $DIR ]]; then
 
        mkdir -p $DIR
        pushd    $DIR
 
        JMergeDisable -f$WORKDIR/${^DETECTORS}.${VERSION}.B/disable.txt -o disable.log
 
        # replace disabled transmissions by those of input detectors
 
        sed '/cat>disable.txt<<EOF/,/EOF/ccat>disable.txt<<EOF\n# disable transmissions\nEOF'  \
            $JPP_DIR/software/JAcoustics/pre-calibration_${DETECTOR}.sh                        >  pre-calibration_${DETECTOR}.sh
        sed -i '/disable transmissions/r disable.log'                                             pre-calibration_${DETECTOR}.sh
 
        chmod +x pre-calibration_${DETECTOR}.sh
 
        ./pre-calibration_${DETECTOR}.sh  $WORKDIR/${^DETECTORS}.${VERSION}.B  >&  pre-calibration.log
 
 
        # disable transmissions
        eval `JPrintDetector -a detector.datx -O IDENTIFIER`
        JPrintDetector -a detector.datx -O IDENTIFIER
        echo Disabling the acoustics
        echo The files to Disable are KM3NeT_${(l:8::0::0:)DETECTOR_ID}_*_event.root
        JAcousticsDisable.sh  detector.datx  KM3NeT_${(l:8::0::0:)DETECTOR_ID}_*_event.root  disable.log
        echo Acoustics disabled
        popd
 
    fi  
 
    DIR=$WORKDIR/${DETECTOR}.${VERSION}.B     # with disabled transmissions
 
    if [[ ! -d $DIR ]]; then
 
        mkdir -p $DIR
        pushd    $DIR
 
        cp $WORKDIR/${DETECTOR}.${VERSION}.A/disable.log .
 
        # introduce disabled transmissions
 
        sed '/cat>disable.txt<<EOF/,/EOF/ccat>disable.txt<<EOF\n# disable transmissions\nEOF'  \
            $JPP_DIR/software/JAcoustics/pre-calibration_${DETECTOR}.sh                        >  pre-calibration_${DETECTOR}.sh
        sed -i '/disable transmissions/r disable.log'                                             pre-calibration_${DETECTOR}.sh
 
        chmod +x pre-calibration_${DETECTOR}.sh
 
        ./pre-calibration_${DETECTOR}.sh  $WORKDIR/${^DETECTORS}.${VERSION}.B  >&  pre-calibration.log
 
        popd
    fi
}
 
# -----------------------------------------------------------------------------------------------
# Run 'empty' pre-calibration of given detector with required inputs.
# The last directory corresponds to the target detector and
# the required preceding directories to the source detectors.
#
# \param 1-N   directories
# -----------------------------------------------------------------------------------------------
fuction fcc()
{
    DETECTORS=($argv[1,-2])
    DETECTOR=$argv[-1]
 
    DIR=$WORKDIR/${DETECTOR}.${VERSION}.B
 
    if [[ ! -d $DIR ]]; then
 
        mkdir -p $DIR
        pushd    $DIR
 
        $JPP_DIR/software/JAcoustics/pre-calibration_${DETECTOR}.sh  $WORKDIR/${^DETECTORS}.${VERSION}.B  >&  pre-calibration.B.log
 
        popd
    fi
}
 
# ORCA
 
run                        D_ORCA006
run  D_ORCA006             D0ORCA010
run  D0ORCA010             D1ORCA011
fcc  D0ORCA010  D1ORCA011  D0ORCA007
fcc  D0ORCA010  D1ORCA011  D0ORCA011
fcc  D0ORCA010  D1ORCA011  D1ORCA013
run  D1ORCA011             D0ORCA015
run  D0ORCA015             D1ORCA015
run  D1ORCA015             D0ORCA018
run  D1ORCA015             DBORCA018
run  DBORCA018             D0ORCA023
run  DBORCA018             D1ORCA023
run  DBORCA018             D1ORCA024
run  DBORCA018             D2ORCA024
 
## ARCA
 
run                        D0ARCA006
run  D0ARCA006             D0ARCA009
run  D0ARCA009             D0ARCA021
fcc  D0ARCA021             D0ARCA020
run  D0ARCA021             D0ARCA028
run  D0ARCA028             D0ARCA030
run  D0ARCA030             D0ARCW003
run  D0ARCA030             D1ARCA030