JAlignDetector.sh

Go to the documentation of this file.

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

# --------------------------------------------------------------------------------------------
#
#          Auxiliary script to align two detectors after acoustic-fit.sh procedure.
#
# --------------------------------------------------------------------------------------------

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

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

set_variable:  DEBUG       ACOUSTICS_DEBUG        2
set_variable:  WORKDIR     ACOUSTICS_WORKDIR      ${TMPDIR:-/tmp}
set_variable:  FORMAT      GRAPHICS_FORMAT        gif
set_variable+  BATCH       GRAPHICS_BATCH         -B
set_variable   SIGMA_M     0.2

if do_usage $*; then
    usage "$script <directory A>  <directory B>"\
          "\nIn each directory, there should be a file \"detector.datx\" and \"tripod.txt\","\
          "\nconform the result of script acoustic-fit.sh."
fi

if (( $# != 2 )); then
    fatal "Wrong number of arguments."
fi

set_variable   DIR_A         $1
set_variable   DIR_B         $2

JConvertDetectorFormat  -a $DIR_A/detector.datx  -o $WORKDIR/detector_a.datx
JConvertDetectorFormat  -a $DIR_B/detector.datx  -o $WORKDIR/detector_b.datx
cp                         $DIR_A/tripod.txt        $WORKDIR/tripod_a.txt
cp                         $DIR_B/tripod.txt        $WORKDIR/tripod_b.txt

eval `JPrintDetector -a $WORKDIR/detector_a.datx -O SUMMARY`


JAlignDetector \
    -a $WORKDIR/detector_a.datx          \
    -b $WORKDIR/detector_b.datx          \
    -A                                   \
    -T $WORKDIR/tripod_a.txt             \
    -s $SIGMA_M                          \
    -d $DEBUG --!

for STRING in $STRINGS[*]; do

    set_variable  MODULE  `getModule -a $WORKDIR/detector_a.datx -L "$STRING 0"`
    
    JEditDetector  -a $WORKDIR/detector_a.datx  -M "$MODULE setz 1.0"  -o $WORKDIR/detector_a.datx
    JEditDetector  -a $WORKDIR/detector_b.datx  -M "$MODULE setz 1.0"  -o $WORKDIR/detector_b.datx
    
done

echo "Output stored at $WORKDIR/detector_a.datx and $WORKDIR/tripod_a.txt."


JDrawDetector2D \
    -a $WORKDIR/detector_a.datx          \
    -a $WORKDIR/detector_b.datx          \
    -L BL                                \
    -o detector.$FORMAT                  $BATCH
    
JDrawDetector2D \
    -T $WORKDIR/tripod_a.txt             \
    -T $WORKDIR/tripod_b.txt             \
    -L BL                                \
    -o tripod.$FORMAT                    $BATCH
    
JCompareDetector \
    -a $WORKDIR/detector_a.datx          \
    -b $WORKDIR/detector_b.datx          \
    -o $WORKDIR/abc.root                 >& /dev/null

for KEY in X Y Z; do

    JPlot2D \
        -f $WORKDIR/abc.root:${KEY}      \
        -\> "string"                     \
        -\< "floor"                      \
        -T "$KEY"                        \
        -O COLZ                          \
        -o ${KEY}.$FORMAT                $BATCH

done