JReconstruction/JEvD.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
script=${0##*/}
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi
 
source $JPP_DIR/setenv.sh $JPP_DIR
 
set_variable:  DEBUG             EVD_DEBUG             2
set_variable:  EDITOR            EVD_EDITOR
set_variable:  OPTION            EVD_OPTION            histogram
set_variable   PDF               $JPP_DATA/J%p.dat.gz
set_variable   DIR               $JPP_DIR/examples/JReconstruction
set_variable   EVENT_SELECTOR    event_selector.so
set_variable:  FORMAT            GRAPHICS_FORMAT       gif
set_variable+  BATCH             GRAPHICS_BATCH        -B
 
if do_usage $*; then
    usage "$script <detector file> <input file> [(calibration file)+]"
fi
 
if (( $# < 2 )); then
    fatal "Wrong number of arguments."
fi
 
set_variable DETECTOR        $argv[1]
set_variable INPUT_FILE      $argv[2]
set_array    CALIBRATION     $argv[3,-1]
 
make -C $DIR JEvD
make -f $DIR/JMakeEventSelector libs
 
if [[ -n $EDITOR ]]; then
    $EDITOR ${EVENT_SELECTOR:r}.cc &
fi
 
eval `JPrintDetector -a $DETECTOR -O HEADER`
 
echo "Set parameters for $DETECTOR_NAME."
 
if [[ "$DETECTOR_NAME" == "ORCA" ]]; then
 
    PARAMETERS=(
        TTS_ns              1.0
        E_GeV               10.0
        TMin_ns             -50.0
        TMax_ns             +450.0
        ZMin_m              -3.0
        ZMax_m              +9.0
        R_Hz                6.0e3
        roadWidth_m         50.0
        numberOfPrefits     0)
else
 
    PARAMETERS=(
        TTS_ns              2
        E_GeV               1.0e3
        TMin_ns             -50.0
        TMax_ns             +450.0
        R_Hz                6.0e3
        roadWidth_m         150.0
        numberOfPrefits     0)
fi
 
eval $DIR/JEvD \
    -f $INPUT_FILE                           \
    -a $DETECTOR                             \
    -+${^CALIBRATION}                        \
    -F $PDF                                  \
    -L `pwd`/$EVENT_SELECTOR                 \
    -@ `make_equation PARAMETERS`            \
    -A 4                                     \
    -o ${${INPUT_FILE:t}:r}.%.$FORMAT        \
    -O $OPTION                               \
    -d $DEBUG                                $BATCH