examples/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 $JPP_DIR/externals/km3net-dataformat/definitions/reconstruction.sh
 
set_variable:  DEBUG             RECONSTRUCTION_DEBUG                1
set_variable:  PDF               RECONSTRUCTION_PDF                  $JPP_DATA/J%p.dat
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
 
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]
 
if (( ${#CALIBRATION} != 0 )); then
    getMechanics.sh  $DETECTOR
fi
 
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
    prefit.TMaxExtra_ns        80
    prefit.ctMin               -0.3
    prefit.DMax_m              50
    prefit.pos_grid_m          10
    prefit.pos_step_m          10
    prefit.time_grid_ns        50
    prefit.time_step_ns        50)
 
PARAMETERS+=(
    simplex.numberOfPrefits    0
    simplex.sigma_ns           2.5
    simplex.TMaxLocal_ns       15
    simplex.TMax_ns            60
    simplex.TMin_ns           -60
    simplex.DMax_m             50
    simplex.ctMin              0.0
    simplex.mestimator          1)
 
PARAMETERS+=(
    position.numberOfPrefits    100
    position.TMax_ns            120
    position.TMin_ns           -120
    position.DMax_m              80
    position.Emin_GeV             1
    position.Emax_GeV           200
    position.En                   3
    position.R_Hz             10000
    position.TTS_ns               2)
 
PARAMETERS+=(
    direction.numberOfPrefits             1
    direction.TMax_ns                    30
    direction.TMin_ns                   -30
    direction.DMax_m                     60
    direction.R_Hz                    10000
    direction.scanAngle_deg              10
    direction.Emin_GeV                   1
    direction.Emax_GeV                   1000
    direction.En                         100)
 
PARAMETERS+=(
    fit.numberOfPrefits             4
    fit.TMax_ns                    30
    fit.TMin_ns                   -30
    fit.DMax_m                     60
    fit.DStep_m                     5
    fit.R_Hz                    10000
    fit.mestimator                  3)
 
eval  $JPP_DIR/examples/JReconstruction/JORCAShowerReconstruction \
    -a $DETECTOR                             \
    -+${^CALIBRATION}                        \
    -T $TMAX_S                               \
    -f ${INPUT_FILE}                         \
    -o ${OUTPUT_FILE}                        \
    -P $PDF                                  \
    -N ${NUMBER_OF_THREADS}                  \
    -@ `make_equation PARAMETERS`            \
    -E $ENERGY_CORRECTION                    \
    -d $DEBUG --!