JORCAShowerReconstruction.sh

Go to the documentation of this file.

#!/bin/zsh
#
# \author adomi
#
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

set_variable:  DEBUG             RECONSTRUCTION_DEBUG                1
set_variable:  PDF               RECONSTRUCTION_PDF                  $JPP_DATA/J%p.dat
set_variable:  SUFFIX            RECONSTRUCTION_SUFFIX               jchain.orca
set_variable:  NUMBER_OF_EVENTS  RECONSTRUCTION_NUMBER_OF_EVENTS     $((2**31))
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

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 [[ ! "$OUTPUT_FILE" == *"%"* ]]; then
    fatal "The output file must have the wildcard % in the name, e.g. %.root"
fi


# ------------------------------------------------------------------------------------------
#                                       JShowerPrefit
# ------------------------------------------------------------------------------------------

maxOutFits=50

#   default input values                                                                            
PARAMETERS=(
    factoryLimit        40
    NMaxHits            20
    sigma_ns             4
    numberOfOutliers     4
    TMaxLocal_ns       120
    ctMin               -0.7
    roadWidth_m         50
    numberOfPrefits     50)  #JShowerPointSimplex needs 50 prefits from JShowerPrefit

timer_start

eval JShowerPrefit                      \
    -a $DETECTOR                        \
    -+${^CALIBRATION}                   \
    -T $TMAX_S                          \
    -f $INPUT_FILE                      \
    -o ${OUTPUT_FILE/\%/jshprefit}      \
    -n ${NUMBER_OF_EVENTS}              \
    -@ `make_equation PARAMETERS`       \
    -d ${DEBUG} --!

timer_stop
timer_print


# ------------------------------------------------------------------------------------------
#                                       JShowerPointSimplex
# ------------------------------------------------------------------------------------------

# default input values                                                                         
PARAMETERS=(
    numberOfPrefits      0  # 0 means all inputs are processed
    numberOfOutfits    100
    minPrefitsSize       9
    sigma_ns             1
    TMaxLocal_ns        15
    TMax_ns             60
    TMin_ns            -60
    roadWidth_m         60
    ctMin                0.2
    TWindow_ns       "-130, 120"
    pos_grid_m       " -20,  20"
    pos_step_m          20
    time_grid_ns     " -90,  90"
    time_step_ns        30
    simplex_step_m       2
    simplex_step_ns     10
    tukey_std_dev       20.0)

timer_start

eval JShowerPointSimplex                  \
    -a $DETECTOR                          \
    -+${^CALIBRATION}                     \
    -T $TMAX_S                            \
    -f ${OUTPUT_FILE/\%/jshprefit}        \
    -o ${OUTPUT_FILE/\%/jshpointsimplex}  \
    -@ `make_equation PARAMETERS`         \
    -d ${DEBUG} --!

timer_stop
timer_print

rm -vf  ${OUTPUT_FILE/\%/jshprefit}       

# ------------------------------------------------------------------------------------------
#                                       JShowerPositionFit
# ------------------------------------------------------------------------------------------

# default input values                                                                                  
PARAMETERS=(
    numberOfPrefits      0
    numberOfOutfits      5
    TMax_ns            120
    TMin_ns           -120
    roadWidth_m         80
    E_GeV               10
    R_Hz             10000
    TTS_ns               2)

timer_start

eval JShowerPositionFit                   \
    -a $DETECTOR                          \
    -+${^CALIBRATION}                     \
    -T $TMAX_S                            \
    -f ${OUTPUT_FILE/\%/jshpointsimplex}  \
    -F $PDF                               \
    -o ${OUTPUT_FILE/\%/jshpositionfit}   \
    -@ `make_equation PARAMETERS`         \
    -d ${DEBUG} --!

timer_stop
timer_print

rm -vf  ${OUTPUT_FILE/\%/jshpointsimplex}       

# ------------------------------------------------------------------------------------------
#                                       JShowerEnergyPrefit
# ------------------------------------------------------------------------------------------

#   default input values                                                           
PARAMETERS=(
    numberOfPrefits             0
    numberOfOutfits             1
    Emin_GeV                    1
    Emax_GeV                  100
    resolution                  0.01
    TMax_ns                   500
    TMin_ns                  -500
    R_Hz                    10000
    mestimator                  3
    roadWidth_m                80)

timer_start

eval JShowerEnergyPrefit                  \
    -a $DETECTOR                          \
    -+${^CALIBRATION}                     \
    -T $TMAX_S                            \
    -f ${OUTPUT_FILE/\%/jshpositionfit}   \
    -F $PDF                               \
    -o ${OUTPUT_FILE/\%/jshenergyprefit}  \
    -@ `make_equation PARAMETERS`         \
    -d $DEBUG --!

timer_stop
timer_print

rm -vf  ${OUTPUT_FILE/\%/jshpositionfit}

# ------------------------------------------------------------------------------------------
#                                       JShowerDirectionPrefit
# ------------------------------------------------------------------------------------------

#   default input values                                                     
PARAMETERS=(
    numberOfPrefits             0
    numberOfOutfits             4
    TMax_ns                    30
    TMin_ns                   -30
    roadWidth_m                80
    DMin_m                     80  # can be lower to speed up the reconstruction (for online)
    E_GeV                      13
    R_Hz                    10000
    mestimator                  3  # 0=normal, 1=lorentzian, 2=linear, 3=null             
    MaxAngle_deg_lowE          40
    scanAngle_deg_lowE         10
    MaxAngle_deg_highE         20
    scanAngle_deg_highE         5)

timer_start

eval JShowerDirectionPrefit                  \
    -a $DETECTOR                             \
    -+${^CALIBRATION}                        \
    -T $TMAX_S                               \
    -f ${OUTPUT_FILE/\%/jshenergyprefit}     \
    -F $PDF                                  \
    -o ${OUTPUT_FILE/\%/jshdirectionprefit}  \
    -@ `make_equation PARAMETERS`            \
    -d $DEBUG --!

timer_stop
timer_print

rm -vf  ${OUTPUT_FILE/\%/jshenergyprefit}

# ------------------------------------------------------------------------------------------
#                                       JShowerFit
# ------------------------------------------------------------------------------------------

#   default input values                                          
PARAMETERS=(
    numberOfPrefits             4
    TMax_ns                    30
    TMin_ns                   -30
    roadWidth_m                80
    R_Hz                    10000
    mestimator                  3  # 0=normal, 1=lorentzian, 2=linear, 3=null                     
    fit_step                    0.05)

timer_start

eval JShowerFit                                                                    \
    -a $DETECTOR                                                                   \
    -+${^CALIBRATION}                                                              \
    -T $TMAX_S                                                                     \
    -f ${OUTPUT_FILE/\%/jshdirectionprefit}                                        \
    -F $PDF                                                                        \
    -o ${OUTPUT_FILE/\%/jshfit}                                                    \
    -@ `make_equation PARAMETERS`                                                  \
    -E $JPP_DIR/software/JReconstruction/energy_correction_showerfit_fullorca.txt  \
    -d $DEBUG --!

CHECK_EXIT_CODE

timer_stop
timer_print

rm -vf  ${OUTPUT_FILE/\%/jshdirectionprefit}

# ------------------------------------------------------------------------------------------
#                                      rename outputfile
# ------------------------------------------------------------------------------------------

timer_start

mv \
    ${OUTPUT_FILE/\%/jshfit}                 \
    ${OUTPUT_FILE/\%/${SUFFIX}}              \
    -v

timer_stop
timer_print


# ------------------------------------------------------------------------------------------
#                                       JConvertEvt
# ------------------------------------------------------------------------------------------

timer_start

JConvertEvt \
    -f ${OUTPUT_FILE/\%/${SUFFIX}}        \
    -o ${OUTPUT_FILE/\%/${SUFFIX}.aanet}  \
    -a $DETECTOR                          \
    -P $PMT                               \
    -d $DEBUG --!

timer_stop
timer_print