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