examples/JReconstruction/JORCAMuonReconstruction.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
#
# \author mdejong
#
script=${0##*/}
 
# ------------------------------------------------------------------------------------------
#
#                         Script to run the JPP Track 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:  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
 
PARAMETERS+=(
    prefit.sigma_ns           5
    prefit.gridAngle_deg      4
    prefit.useL0              1
    prefit.numberOfOutliers   2
    prefit.numberOfPrefits    50
    prefit.DZMax              0.0
    prefit.numberOfDZMax      1
    prefit.TMaxLocal_ns       12
    prefit.ctMin              0.0
    prefit.roadWidth_m        50.0
    prefit.Qmin               1.5)
 
PARAMETERS+=(
    simplex.sigma_ns         3
    simplex.TMaxLocal_ns     15.0
    simplex.ctMin            0.0
    simplex.roadWidth_m      50.0
    simplex.numberOfPrefits  0)
 
PARAMETERS+=(
    start.roadWidth_m      50.0
    start.R_Hz             7.0e3
    start.TMin_ns          -25.0
    start.TMax_ns          +25.0
    start.Pmin1            1.0e-3
    start.Pmin2            1.0e-2
    start.numberOfPrefits  0)
 
PARAMETERS+=(
    gandalf.TTS_ns          2
    gandalf.E_GeV           10.0
    gandalf.TMin_ns         -50.0
    gandalf.TMax_ns         +450.0
    gandalf.ZMin_m          -3.0
    gandalf.ZMax_m          +9.0
    gandalf.R_Hz            7.0e3
    gandalf.roadWidth_m     50.0
    gandalf.numberOfPrefits   0)
 
PARAMETERS+=(
    start2.roadWidth_m      50.0
    start2.R_Hz             7.0e3
    start2.TMin_ns          -25.0
    start2.TMax_ns          +25.0
    start2.Pmin1            1.0e-3
    start2.Pmin2            1.0e-2
    start2.numberOfPrefits  1)
 
PARAMETERS+=(
    energy.roadWidth_m      80.0
    energy.R_Hz             7.0e3
    energy.EMin_log         0.0
    energy.EMax_log         7.0
    energy.TMin_ns          -25.0
    energy.TMax_ns          +25.0
    energy.ZMin_m          -7.5
    energy.resolution       0.01
    energy.mestimator       3
    energy.numberOfPrefits  1)
 
eval $JPP_DIR/examples/JReconstruction/JORCAMuonReconstruction \
    -a $DETECTOR                   \
    -+${^CALIBRATION}              \
    -T $TMAX_S                     \
    -f ${INPUT_FILE}               \
    -o ${OUTPUT_FILE}              \
    -P $PDF                        \
    -N ${NUMBER_OF_THREADS}        \
    -@ `make_equation PARAMETERS`  \
    -d $DEBUG --!