Jpp 19.3.0-rc.3
the software that should make you happy
Loading...
Searching...
No Matches
examples/JReconstruction/JORCAMuonReconstruction.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 ORCA.
10#
11# ------------------------------------------------------------------------------------------
12
13if [ -z $JPP_DIR ]; then
14 echo "Variable JPP_DIR undefined."
15 exit
16fi
17
18source $JPP_DIR/setenv.sh $JPP_DIR
19source $JPP_DIR/externals/km3net-dataformat/definitions/reconstruction.sh
20
21set_variable: DEBUG RECONSTRUCTION_DEBUG 1
22set_variable: PDF RECONSTRUCTION_PDF $JPP_DATA/J%p.dat
23set_variable: NUMBER_OF_EVENTS RECONSTRUCTION_NUMBER_OF_EVENTS $((2**31))
24set_variable: NUMBER_OF_THREADS RECONSTRUCTION_NUMBER_OF_THREADS 0
25set_variable+ OVERWRITE RECONSTRUCTION_OVERWRITE 1
26set_variable: TMAX_S RECONSTRUCTION_TMAX_S 100
27
28if do_usage $*; then
29 usage "$script <detector file> <input file> <output file> <PMT parameters file> [(calibration file)+]"
30fi
31
32if (( $# < 4 )); then
33 fatal "Wrong number of arguments."
34fi
35
36set_variable DETECTOR $argv[1]
37set_variable INPUT_FILE $argv[2]
38set_variable OUTPUT_FILE $argv[3]
39set_variable PMT $argv[4]
40set_array CALIBRATION $argv[5,-1]
41
42if (( ${#CALIBRATION} != 0 )); then
43 getMechanics.sh $DETECTOR
44fi
45
46PARAMETERS+=(
47 prefit.sigma_ns 5
48 prefit.gridAngle_deg 4
49 prefit.useL0 1
50 prefit.numberOfOutliers 2
51 prefit.numberOfPrefits 50
52 prefit.DZMax 0.0
53 prefit.numberOfDZMax 1
54 prefit.TMaxLocal_ns 12
55 prefit.ctMin 0.0
56 prefit.roadWidth_m 50.0
57 prefit.Qmin 1.5)
58
59PARAMETERS+=(
60 simplex.sigma_ns 3
61 simplex.TMaxLocal_ns 15.0
62 simplex.ctMin 0.0
63 simplex.roadWidth_m 50.0
64 simplex.numberOfPrefits 0)
65
66PARAMETERS+=(
67 start.roadWidth_m 50.0
68 start.R_Hz 7.0e3
69 start.TMin_ns -25.0
70 start.TMax_ns +25.0
71 start.Pmin1 1.0e-3
72 start.Pmin2 1.0e-2
73 start.numberOfPrefits 0)
74
75PARAMETERS+=(
76 gandalf.TTS_ns 2
77 gandalf.E_GeV 10.0
78 gandalf.TMin_ns -50.0
79 gandalf.TMax_ns +450.0
80 gandalf.ZMin_m -3.0
81 gandalf.ZMax_m +9.0
82 gandalf.R_Hz 7.0e3
83 gandalf.roadWidth_m 50.0
84 gandalf.numberOfPrefits 0)
85
86PARAMETERS+=(
87 start2.roadWidth_m 50.0
88 start2.R_Hz 7.0e3
89 start2.TMin_ns -25.0
90 start2.TMax_ns +25.0
91 start2.Pmin1 1.0e-3
92 start2.Pmin2 1.0e-2
93 start2.numberOfPrefits 1)
94
95PARAMETERS+=(
96 energy.roadWidth_m 80.0
97 energy.R_Hz 7.0e3
98 energy.EMin_log 0.0
99 energy.EMax_log 7.0
100 energy.TMin_ns -25.0
101 energy.TMax_ns +25.0
102 energy.ZMin_m -7.5
103 energy.resolution 0.01
104 energy.mestimator 3
105 energy.numberOfPrefits 1)
106
107eval $JPP_DIR/examples/JReconstruction/JORCAMuonReconstruction \
108 -a $DETECTOR \
109 -+${^CALIBRATION} \
110 -T $TMAX_S \
111 -f ${INPUT_FILE} \
112 -o ${OUTPUT_FILE} \
113 -P $PDF \
114 -N ${NUMBER_OF_THREADS} \
115 -@ `make_equation PARAMETERS` \
116 -d $DEBUG --!