JORCAReconstruction.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
#
# \author mdejong
#
script=${0##*/}
 
# ------------------------------------------------------------------------------------------
#
#                         Script to run the JPP 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+=(
    muon.prefit.sigma_ns           5
    muon.prefit.gridAngle_deg      4
    muon.prefit.useL0              1
    muon.prefit.numberOfOutliers   2
    muon.prefit.numberOfPrefits    50
    muon.prefit.DZMax              0.0
    muon.prefit.numberOfDZMax      1
    muon.prefit.TMaxLocal_ns       12
    muon.prefit.ctMin              0.0
    muon.prefit.roadWidth_m        50.0
    muon.prefit.Qmin               1.5)
 
PARAMETERS+=(
    muon.simplex.sigma_ns         3
    muon.simplex.TMaxLocal_ns     15.0
    muon.simplex.ctMin            0.0
    muon.simplex.roadWidth_m      50.0
    muon.simplex.numberOfPrefits  0)
 
PARAMETERS+=(
    muon.start.roadWidth_m      50.0
    muon.start.R_Hz             7.0e3
    muon.start.TMin_ns          -25.0
    muon.start.TMax_ns          +25.0
    muon.start.Pmin1            1.0e-3
    muon.start.Pmin2            1.0e-2
    muon.start.numberOfPrefits  0)
 
PARAMETERS+=(
    muon.gandalf.TTS_ns          2
    muon.gandalf.E_GeV           10.0
    muon.gandalf.TMin_ns         -50.0
    muon.gandalf.TMax_ns         +450.0
    muon.gandalf.ZMin_m          -3.0
    muon.gandalf.ZMax_m          +9.0
    muon.gandalf.R_Hz            7.0e3
    muon.gandalf.roadWidth_m     50.0
    muon.gandalf.numberOfPrefits   0)
 
PARAMETERS+=(
    muon.start2.roadWidth_m      50.0
    muon.start2.R_Hz             7.0e3
    muon.start2.TMin_ns          -25.0
    muon.start2.TMax_ns          +25.0
    muon.start2.Pmin1            1.0e-3
    muon.start2.Pmin2            1.0e-2
    muon.start2.numberOfPrefits  1)
 
PARAMETERS+=(
    muon.energy.roadWidth_m      80.0
    muon.energy.R_Hz             7.0e3
    muon.energy.EMin_log         0.0
    muon.energy.EMax_log         7.0
    muon.energy.TMin_ns          -25.0
    muon.energy.TMax_ns          +25.0
    muon.energy.ZMin_m          -7.5
    muon.energy.resolution       0.01
    muon.energy.mestimator       3
    muon.energy.numberOfPrefits  1)
 
PARAMETERS+=(
    shower.prefit.numberOfPrefits     200
    shower.prefit.numberOfGrids       2
    shower.prefit.factoryLimit        20
    shower.prefit.sigma_ns             2.5
    shower.prefit.numberOfOutliers     4
    shower.prefit.TMaxLocal_ns        10
    shower.prefit.TMaxExtra_ns        80
    shower.prefit.ctMin               -0.3
    shower.prefit.DMax_m              50
    shower.prefit.pos_grid_m          10
    shower.prefit.pos_step_m          10
    shower.prefit.time_grid_ns        50
    shower.prefit.time_step_ns        50)
 
PARAMETERS+=(
    shower.simplex.numberOfPrefits    0
    shower.simplex.sigma_ns           2.5
    shower.simplex.TMaxLocal_ns       15
    shower.simplex.TMax_ns            60
    shower.simplex.TMin_ns           -60
    shower.simplex.DMax_m             50
    shower.simplex.ctMin              0.0
    shower.simplex.mestimator          1)
 
PARAMETERS+=(
    shower.position.numberOfPrefits    100
    shower.position.TMax_ns            120
    shower.position.TMin_ns           -120
    shower.position.DMax_m              80
    shower.position.Emin_GeV             1
    shower.position.Emax_GeV           200
    shower.position.En                   3
    shower.position.R_Hz             10000
    shower.position.TTS_ns               2)
 
PARAMETERS+=(
    shower.direction.numberOfPrefits             1
    shower.direction.TMax_ns                    30
    shower.direction.TMin_ns                   -30
    shower.direction.DMax_m                     60
    shower.direction.R_Hz                    10000
    shower.direction.scanAngle_deg              10
    shower.direction.Emin_GeV                   1
    shower.direction.Emax_GeV                   1000
    shower.direction.En                         100)
 
PARAMETERS+=(
    shower.fit.numberOfPrefits             4
    shower.fit.TMax_ns                    30
    shower.fit.TMin_ns                   -30
    shower.fit.DMax_m                     60
    shower.fit.DStep_m                     5
    shower.fit.R_Hz                    10000
    shower.fit.mestimator                  3)
 
eval $JPP_DIR/examples/JReconstruction/JORCAReconstruction \
    -a $DETECTOR                   \
    -+${^CALIBRATION}              \
    -T $TMAX_S                     \
    -f ${INPUT_FILE}               \
    -o ${OUTPUT_FILE}              \
    -P $PDF                        \
    -N ${NUMBER_OF_THREADS}        \
    -@ `make_equation PARAMETERS`  \
    -d $DEBUG --!