pre-calibration_D0ARCW003.sh

Go to the documentation of this file.

#!/bin/zsh
script=${0##*/}
 
# --------------------------------------------------------------------------------------------
#
# Steering script for the global-fit-of-global-fits procedure of the P0ARCW016 (233) detector.
#
# --------------------------------------------------------------------------------------------
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi
 
source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
 
set_variable  DETECTOR_ID            236
set_array     RUNS                   177-190
set_variable: MUL                    ACOUSTICS_MUL          -1.21e-2
set_variable: THREADS                ACOUSTICS_THREADS       32
set_variable  SINGULARITY            -
set_variable  WORKDIR                `pwd`
set_variable  HOMEDIR                $WORKDIR
 
if do_usage $*; then
    usage "$script <source directory> [user directory]"\
          "\nThe source directory corresponds to a pre-calibrated D0ARCA030 detector."\
          "\nThe user directory should be specified when this script is submitted in a batch queue."
fi
 
case $# in
    2)  set_variable HOMEDIR              $2;&
    1)  set_variable DIR                  ${1:a};;
    *)  fatal "Wrong number of arguments."
esac
 
source JAcousticsToolkit.sh
 
set_variable  TRIPOD                 $DIR/tripod.txt
set_variable  DETECTOR_INITIAL       detector_initial.datx
expand_array  RUNS
 
 
# Run in HOMEDIR
 
mkdir -p $HOMEDIR
cd       $HOMEDIR
 
JDetectorDB \
    -D $DETECTOR_ID       \
    -r $RUNS[1]           \
    -V ""                 \
    -WW                   \
    -o $DETECTOR_INITIAL
 
# Apply modifications.
 
JEditDetector \
    -a $DETECTOR_INITIAL     \
    -s "-1 mul $MUL"         \
    -o detector.datx --!
 
eval `JPrintDetector -a detector.datx -O SUMMARY`
 
for STRING in ${STRINGS[*]}; do
 
    set_variable  MODULE  `getModule -a detector.datx -L "$STRING 0"`
 
    JEditDetector \
        -a detector.datx         \
        -M "$MODULE setz -11.7"  \
        -o detector.datx --!
done
 
typeset -A TRIPODS
 
get_tripods  $TRIPOD          TRIPODS;  FIXED_TRIPODS=(${(k)TRIPODS[*]})
 
cp -p  $TRIPOD  tripod.txt
 
 
cat>waveform.txt<<EOF
# waveform identifier; emitter identifier
 28   9
-29   9
 25  10
-26  10
 29  13
-30  13
 34  14
-35  14
 33  16
-34  16
 23  18
-24  18
 39  19
-40  19
 45  20
-46  20
 47  21
-48  21
 41  22
-42  22
EOF
 
cat>acoustics_trigger_parameters.txt<<EOF
Q             = 0.0;
TMax_s        = 0.035;
numberOfHits  = 40;
EOF
 
cat>disable.txt<<EOF
# disable transmissions
EOF
 
cat>transmitter.txt<<EOF
# Data from xxxx_baseAcousticBeacon.csv
EOF
 
cat>hydrophone.txt<<EOF
# Data from xxxx_hydrophones.csv
EOF
 
JAcousticsEventBuilder.sh  detector.datx  $RUNS[*]
 
INPUT_FILES=(`ls KM3NeT_${(l:8::0::0:)DETECTOR_ID}_0*${^RUNS}_event.root`)
 
 
# Run in WORKDIR
 
cd  $WORKDIR
 
if [ ! $HOMEDIR -ef $WORKDIR ]; then
    cp -p $HOMEDIR/$ACOUSTICS_DETECTOR     $WORKDIR
    cp -p $HOMEDIR/${^ACOUSTICS_KEYS}.txt  $WORKDIR
    cp -p $HOMEDIR/${^ACOUSTICS_AUXS}.txt  $WORKDIR
    cp -p $HOMEDIR/${^INPUT_FILES}         $WORKDIR
fi
 
cat>script.txt<<EOF
fix  tripod  $FIXED_TRIPODS[*]
 
stage  1A  0                2  250.0e-6  10.0  20  0.5    0.5
stage  1B  0                2  250.0e-6  10.0   0  0.0005 0.2
stage  1C  ${SINGULARITY}2  0  250.0e-6  10.0  20  0.3       
stage  2A  ${SINGULARITY}2  0  100.0e-6  10.0  20  0.2    0.2
stage  2C  ${SINGULARITY}2  0  100.0e-6  10.0  20  0.2       
stage  3A  ${SINGULARITY}2  0   50.0e-6  10.0  20  0.2    0.2
stage  3C  ${SINGULARITY}2  0   50.0e-6  10.0  20  0.2       
EOF
 
JSydney.sh  detector.datx  $INPUT_FILES[*]  script.txt
 
if [ ! $HOMEDIR -ef $WORKDIR ]; then
    cp -vru * $HOMEDIR
fi