software/JReconstruction/JORCAShowerReconstruction.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
#
# \author mdejong
#
script=${0##*/}
 
 
# ------------------------------------------------------------------------------------------
#
#                         Script to run the JPP Shower Reconstruction Chain for ORCA.
#
# ------------------------------------------------------------------------------------------
 
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi
 
source $JPP_DIR/setenv.sh $JPP_DIR
 
 
source reconstruction.sh
 
set_variable:  DEBUG              RECONSTRUCTION_DEBUG                1
set_variable:  PDF                RECONSTRUCTION_PDF                  $JPP_DATA/J%p.dat.gz
set_variable:  ENERGY_CORRECTION  RECONSTRUCTION_ENERGY_CORRECTION    energy_correction_showerfit_fullorca.txt
set_variable:  NUMBER_OF_EVENTS   RECONSTRUCTION_NUMBER_OF_EVENTS     $((2**31))
set_variable:  NUMBER_OF_THREADS  RECONSTRUCTION_NUMBER_OF_THREADS    0
set_variable+  OVERWRITE          RECONSTRUCTION_OVERWRITE            1
set_variable:  TMAX_S             RECONSTRUCTION_TMAX_S               100
set_variable:  RATE_HZ            RECONSTRUCTION_RATE_HZ              7.0e3
 
if do_usage $*; then
    usage "$script <detector file> <input file> <output file> <PMT parameters file> [(calibration file)+]"
fi
 
if (( $# < 4 )); then
    fatal "Wrong number of arguments."
fi
 
set_variable DETECTOR        $argv[1]
set_variable INPUT_FILE      $argv[2]
set_variable OUTPUT_FILE     $argv[3]
set_variable PMT             $argv[4]
set_array    CALIBRATION     $argv[5,-1]
 
typeset -A PARAMETERS
 
PARAMETERS+=(
    prefit.numberOfPrefits          200
    prefit.numberOfGrids            2
    prefit.factoryLimit             20
    prefit.sigma_ns                 2.5
    prefit.numberOfOutliers         4
    prefit.TMaxLocal_ns             10.0
    prefit.TMaxExtra_ns             80.0
    prefit.ctMin                    -0.3
    prefit.DMax_m                   50.0
    prefit.pos_grid_m               10.0
    prefit.pos_step_m               10.0
    prefit.time_grid_ns             50.0
    prefit.time_step_ns             50.0)
 
PARAMETERS+=(
    simplex.numberOfPrefits         0
    simplex.numberOfPostfits        25
    simplex.Qwatershed              2
    simplex.sigma_ns                2.5
    simplex.TMaxLocal_ns            15.0
    simplex.TMax_ns                 +60.0
    simplex.TMin_ns                 -60.0
    simplex.DMax_m                  50.0
    simplex.ctMin                   0.0
    simplex.mestimator              1)
 
PARAMETERS+=(
    position.numberOfPrefits        0
    position.TMax_ns                +120.0
    position.TMin_ns                -120.0
    position.DMax_m                 80.0
    position.Emin_GeV               1.0
    position.Emax_GeV               10000.0
    position.En                     5
    position.R_Hz                   $RATE_HZ
    position.TTS_ns                 2.0)
 
PARAMETERS+=(
    direction.numberOfPrefits       1
    direction.TMax_ns               +30.0
    direction.TMin_ns               -30.0
    direction.DMax_m                60.0
    direction.R_Hz                  $RATE_HZ
    direction.scanAngle_deg         10
    direction.Emin_GeV              1
    direction.Emax_GeV              1000
    direction.En                    100)
 
PARAMETERS+=(
    fit.numberOfPrefits             4
    fit.TMax_ns                     +30.0
    fit.TMin_ns                     -30.0
    fit.DMax_m                      60.0
    fit.DStep_m                     5
    fit.R_Hz                        $RATE_HZ
    fit.mestimator                  3)
 
eval JORCAShowerReconstruction \
    -a $DETECTOR                             \
    -+${^CALIBRATION}                        \
    -T $TMAX_S                               \
    -f ${INPUT_FILE}                         \
    -o ${OUTPUT_FILE}                        \
    -F $PDF                                  \
    -n ${NUMBER_OF_EVENTS}                   \
    -N ${NUMBER_OF_THREADS}                  \
    -@ `make_equation PARAMETERS`            \
    -E $ENERGY_CORRECTION                    \
    -P $PMT                                  \
    -d $DEBUG --!
 
CHECK_EXIT_CODE