JCompareDynamics.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       DYNAMICS_DEBUG                2
set_variable:  WORKDIR     DYNAMICS_WORKDIR              ./
set_variable   TSTEP_S     600.0                        
set_variable   DIR         $JPP_DIR/examples/JDynamics/  
set_variable:  TIMESTAMP   GRAPHICS_TIMESTAMP     utc
set_variable:  FORMAT      GRAPHICS_FORMAT        gif

if do_usage $*; then
    usage "$script <detector file> <calibration file 1> <calibration file 2>"
          "\nCalibration files correspond to the output of JKatoomba[.sh]."
fi

if (( $# < 3 )); then
    fatal "Wrong number of arguments."
fi

set_variable  DETECTOR            $argv[1]
set_variable  INPUT_FILE_1        $argv[2]
set_variable  INPUT_FILE_2        $argv[3]
set_variable  OUTPUT_FILE         $WORKDIR/compare_dynamic_positions.root 

eval `JPrintDetector -a $DETECTOR -O SUMMARY`

$DIR/JCompareDynamicPositions \
    -a $DETECTOR                             \
    -f "$INPUT_FILE_1"                       \
    -F "$INPUT_FILE_2"                       \
    -o $OUTPUT_FILE                          \
    -T $TSTEP_S                              \
    -d $DEBUG --!

echo "Creating graphics output in directory $WORKDIR"

JPlot1D \
    -w 1200x600                                       \
    -f $OUTPUT_FILE:H\.\*comparetx                    \
    -\^ "#DeltaTx [mrad]"                             \
    -N "X 505"                                        \
    -t "$TIMESTAMP"                                   \
    -L BR                                             \
    -T ""                                             \
    -o $WORKDIR/comparetx.$FORMAT                     -B

JPlot1D \
    -w 1200x600                                       \
    -f $OUTPUT_FILE:K\.\*comparety                    \
    -\^ "#DeltaTy [mrad]"                             \
    -N "X 505"                                        \
    -t "$TIMESTAMP"                                   \
    -L TR                                             \
    -T ""                                             \
    -o $WORKDIR/comparety.$FORMAT                     -B

JPlot1D \
    -w 1200x600                                       \
    -f $OUTPUT_FILE:M\.\*comparetx                    \
    -\> "#DeltaTx [mrad]"                             \
    -\^ "number of events [a.u.]"                     \
    -N "X 505"                                        \
    -L BR                                             \
    -T ""                                             \
    -s 100110                                         \
    -x "0 1"  \
    -o $WORKDIR/comparetx_hist.$FORMAT                 -B

JPlot1D \
    -w 1200x600                                       \
    -f $OUTPUT_FILE:N\.\*comparety                    \
    -\> "#DeltaTy [mrad]"                             \
    -\^ "number of events [a.u.]"                     \
    -N "X 505"                                        \
    -L BR                                             \
    -T ""                                             \
    -s 100110                                         \
    -x "0 1"  \
    -o $WORKDIR/comparety_hist.$FORMAT                 -B

#clean
rm compare_dynamic_positions.root