Jpp 19.3.0-rc.3
the software that should make you happy
Loading...
Searching...
No Matches
examples/JReconstruction/JARCAMuonReconstruction.sh
Go to the documentation of this file.
1#!/usr/bin/env zsh
2#
3# \author mdejong
4#
5script=${0##*/}
6
7# ------------------------------------------------------------------------------------------
8#
9# Script to run the JPP Track Reconstruction Chain for ARCA.
10#
11# ------------------------------------------------------------------------------------------
12
13JPP_DIR=/localstore/km3net/mjg/GIT/storage/
14
15if [ -z $JPP_DIR ]; then
16 echo "Variable JPP_DIR undefined."
17 exit
18fi
19
20source $JPP_DIR/setenv.sh $JPP_DIR
21source $JPP_DIR/externals/km3net-dataformat/definitions/reconstruction.sh
22
23set_variable: DEBUG RECONSTRUCTION_DEBUG 1
24set_variable: PDF RECONSTRUCTION_PDF $JPP_DATA/J%p.dat
25set_variable: ENERGY_CORRECTION RECONSTRUCTION_ENERGY_CORRECTION energy_correction_arca.txt
26set_variable: NUMBER_OF_EVENTS RECONSTRUCTION_NUMBER_OF_EVENTS $((2**31))
27set_variable: NUMBER_OF_THREADS RECONSTRUCTION_NUMBER_OF_THREADS 0
28set_variable+ OVERWRITE RECONSTRUCTION_OVERWRITE 1
29set_variable: TMAX_S RECONSTRUCTION_TMAX_S 100
30set_variable: RATE_HZ RECONSTRUCTION_RATE_HZ 6.0e3
31
32if do_usage $*; then
33 usage "$script <detector file> <input file> <output file> <PMT parameters file> [(calibration file)+]"
34fi
35
36if (( $# < 4 )); then
37 fatal "Wrong number of arguments."
38fi
39
40set_variable DETECTOR $argv[1]
41set_variable INPUT_FILE $argv[2]
42set_variable OUTPUT_FILE $argv[3]
43set_variable PMT $argv[4]
44set_array CALIBRATION $argv[5,-1]
45
46if (( ${#CALIBRATION} != 0 )); then
47 getMechanics.sh $DETECTOR
48fi
49
50typeset -A PARAMETERS
51
52PARAMETERS+=(
53 prefit.sigma_ns 5
54 prefit.gridAngle_deg 1
55 prefit.useL0 0
56 prefit.numberOfOutliers 3
57 prefit.numberOfPrefits 50
58 prefit.DZMax 0.0
59 prefit.numberOfDZMax 1
60 prefit.TMaxLocal_ns 18
61 prefit.ctMin 0.0
62 prefit.roadWidth_m 200.0)
63
64PARAMETERS+=(
65 simplex.sigma_ns 3
66 simplex.TMaxLocal_ns 15.0
67 simplex.ctMin 0.0
68 simplex.roadWidth_m 200.0
69 simplex.numberOfPrefits 0)
70
71PARAMETERS+=(
72 gandalf.TTS_ns 2
73 gandalf.E_GeV 1.0e3
74 gandalf.TMin_ns -50.0
75 gandalf.TMax_ns +450.0
76 gandalf.R_Hz $RATE_HZ
77 gandalf.roadWidth_m 175.0
78 gandalf.numberOfPrefits 0)
79
80PARAMETERS+=(
81 start.roadWidth_m 160.0
82 start.R_Hz $RATE_HZ
83 start.TMin_ns -25.0
84 start.TMax_ns +75.0
85 start.Pmin1 1.0e-3
86 start.Pmin2 1.0e-2
87 start.numberOfPrefits 1)
88
89PARAMETERS+=(
90 energy.roadWidth_m 300.0
91 energy.R_Hz $RATE_HZ
92 energy.EMin_log 0.0
93 energy.EMax_log 8.0
94 energy.TMin_ns -50.0
95 energy.TMax_ns +450.0
96 energy.mestimator 1
97 energy.numberOfPrefits 1)
98
99eval $JPP_DIR/examples/JReconstruction/JARCAMuonReconstruction \
100 -a $DETECTOR \
101 -+${^CALIBRATION} \
102 -T $TMAX_S \
103 -f ${INPUT_FILE} \
104 -o ${OUTPUT_FILE} \
105 -P $PDF \
106 -E $ENERGY_CORRECTION \
107 -n ${NUMBER_OF_EVENTS} \
108 -N ${NUMBER_OF_THREADS} \
109 -@ `make_equation PARAMETERS` \
110 -d $DEBUG --!
111