software/JReconstruction/JARCAMuonReconstruction.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 ARCA.
#
# ------------------------------------------------------------------------------------------
 
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:  OPTION             RECONSTRUCTION_OPTION               --!
set_variable:  PDF                RECONSTRUCTION_PDF                  $JPP_DATA/J%p.dat
set_variable:  ENERGY_CORRECTION  RECONSTRUCTION_ENERGY_CORRECTION    energy_correction_arca.root
set_variable:  NUMBER_OF_EVENTS   RECONSTRUCTION_NUMBER_OF_EVENTS     $((2**31))
set_variable:  NUMBER_OF_THREADS  RECONSTRUCTION_NUMBER_OF_THREADS    0
set_variable:  TMAX_S             RECONSTRUCTION_TMAX_S               100
set_variable:  RATE_HZ            RECONSTRUCTION_RATE_HZ              6.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]
 
if (( ${#CALIBRATION} != 0 )); then
    getMechanics.sh  $DETECTOR
fi
 
typeset -A PARAMETERS
 
PARAMETERS+=(
    prefit.sigma_ns             5.0
    prefit.gridAngle_deg        1.0
    prefit.useL0                0
    prefit.numberOfOutliers     3
    prefit.numberOfPrefits      60
    prefit.numberOfPostfits     20
    prefit.TMaxLocal_ns         18.0
    prefit.ctMin                0.0
    prefit.roadWidth_m          200.0)
 
PARAMETERS+=(
    simplex.sigma_ns            3.0
    simplex.TMaxLocal_ns        15.0
    simplex.ctMin               0.0
    simplex.roadWidth_m         200.0
    simplex.numberOfPrefits     0)
 
PARAMETERS+=(
    gandalf.TTS_ns              1.0
    gandalf.E_GeV               1.0e3
    gandalf.TMin_ns             -50.0
    gandalf.TMax_ns             +450.0
    gandalf.R_Hz                $RATE_HZ
    gandalf.roadWidth_m         175.0
    gandalf.cosLR               0.7
    gandalf.numberOfPrefits     0)
 
PARAMETERS+=(
    start.roadWidth_m           160.0
    start.R_Hz                  $RATE_HZ
    start.TMin_ns               -25.0
    start.TMax_ns               +75.0
    start.Pmin                  1.0e-8
    start.Pmin1                 2.0e-4
    start.Pmin2                 1.0e-2
    start.numberOfPrefits       1
    start.numberOfPostfits      0)
 
PARAMETERS+=(
    energy.roadWidth_m          300.0
    energy.R_Hz                 $RATE_HZ
    energy.EMin_log             0.0
    energy.EMax_log             8.0
    energy.TMin_ns              -50.0
    energy.TMax_ns              +450.0
    energy.mestimator           1
    energy.numberOfPrefits      0)
 
PARAMETERS+=(
    features.TTS_ns             1.0
    features.E_GeV              1.0e3
    features.TMin_ns            -50.0
    features.TMax_ns            +450.0
    features.ZMin_m             -3.0
    features.ZMax_m             +9.0
    features.R_Hz               $RATE_HZ
    features.roadWidth_m        175.0
    features.numberOfPrefits    0)
 
eval JARCAMuonReconstruction \
    -a $DETECTOR                   \
    -+${^CALIBRATION}              \
    -T $TMAX_S                     \
    -f ${INPUT_FILE}               \
    -o ${OUTPUT_FILE}              \
    -F $PDF                        \
    -E $ENERGY_CORRECTION          \
    -n ${NUMBER_OF_EVENTS}         \
    -N ${NUMBER_OF_THREADS}        \
    -@ `make_equation PARAMETERS`  \
    -P $PMT                        \
    -d $DEBUG                      \
    $OPTION
 
CHECK_EXIT_CODE