7# ------------------------------------------------------------------------------------------
9# Script to run the JPP Reconstruction Chain for ORCA.
11# ------------------------------------------------------------------------------------------
13if [ -z $JPP_DIR ]; then
14 echo "Variable JPP_DIR undefined."
18source $JPP_DIR/setenv.sh $JPP_DIR
19source $JPP_DIR/externals/km3net-dataformat/definitions/reconstruction.sh
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
29 usage "$script <detector file> <input file> <output file> <PMT parameters file> [(calibration file)+]"
33 fatal "Wrong number of arguments."
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]
42if (( ${#CALIBRATION} != 0 )); then
43 getMechanics.sh $DETECTOR
47 muon.prefit.sigma_ns 5
48 muon.prefit.gridAngle_deg 4
50 muon.prefit.numberOfOutliers 2
51 muon.prefit.numberOfPrefits 50
53 muon.prefit.numberOfDZMax 1
54 muon.prefit.TMaxLocal_ns 12
56 muon.prefit.roadWidth_m 50.0
60 muon.simplex.sigma_ns 3
61 muon.simplex.TMaxLocal_ns 15.0
62 muon.simplex.ctMin 0.0
63 muon.simplex.roadWidth_m 50.0
64 muon.simplex.numberOfPrefits 0)
67 muon.start.roadWidth_m 50.0
69 muon.start.TMin_ns -25.0
70 muon.start.TMax_ns +25.0
71 muon.start.Pmin1 1.0e-3
72 muon.start.Pmin2 1.0e-2
73 muon.start.numberOfPrefits 0)
77 muon.gandalf.E_GeV 10.0
78 muon.gandalf.TMin_ns -50.0
79 muon.gandalf.TMax_ns +450.0
80 muon.gandalf.ZMin_m -3.0
81 muon.gandalf.ZMax_m +9.0
82 muon.gandalf.R_Hz 7.0e3
83 muon.gandalf.roadWidth_m 50.0
84 muon.gandalf.numberOfPrefits 0)
87 muon.start2.roadWidth_m 50.0
88 muon.start2.R_Hz 7.0e3
89 muon.start2.TMin_ns -25.0
90 muon.start2.TMax_ns +25.0
91 muon.start2.Pmin1 1.0e-3
92 muon.start2.Pmin2 1.0e-2
93 muon.start2.numberOfPrefits 1)
96 muon.energy.roadWidth_m 80.0
97 muon.energy.R_Hz 7.0e3
98 muon.energy.EMin_log 0.0
99 muon.energy.EMax_log 7.0
100 muon.energy.TMin_ns -25.0
101 muon.energy.TMax_ns +25.0
102 muon.energy.ZMin_m -7.5
103 muon.energy.resolution 0.01
104 muon.energy.mestimator 3
105 muon.energy.numberOfPrefits 1)
108 shower.prefit.numberOfPrefits 200
109 shower.prefit.numberOfGrids 2
110 shower.prefit.factoryLimit 20
111 shower.prefit.sigma_ns 2.5
112 shower.prefit.numberOfOutliers 4
113 shower.prefit.TMaxLocal_ns 10
114 shower.prefit.TMaxExtra_ns 80
115 shower.prefit.ctMin -0.3
116 shower.prefit.DMax_m 50
117 shower.prefit.pos_grid_m 10
118 shower.prefit.pos_step_m 10
119 shower.prefit.time_grid_ns 50
120 shower.prefit.time_step_ns 50)
123 shower.simplex.numberOfPrefits 0
124 shower.simplex.sigma_ns 2.5
125 shower.simplex.TMaxLocal_ns 15
126 shower.simplex.TMax_ns 60
127 shower.simplex.TMin_ns -60
128 shower.simplex.DMax_m 50
129 shower.simplex.ctMin 0.0
130 shower.simplex.mestimator 1)
133 shower.position.numberOfPrefits 100
134 shower.position.TMax_ns 120
135 shower.position.TMin_ns -120
136 shower.position.DMax_m 80
137 shower.position.Emin_GeV 1
138 shower.position.Emax_GeV 200
140 shower.position.R_Hz 10000
141 shower.position.TTS_ns 2)
144 shower.direction.numberOfPrefits 1
145 shower.direction.TMax_ns 30
146 shower.direction.TMin_ns -30
147 shower.direction.DMax_m 60
148 shower.direction.R_Hz 10000
149 shower.direction.scanAngle_deg 10
150 shower.direction.Emin_GeV 1
151 shower.direction.Emax_GeV 1000
152 shower.direction.En 100)
155 shower.fit.numberOfPrefits 4
156 shower.fit.TMax_ns 30
157 shower.fit.TMin_ns -30
160 shower.fit.R_Hz 10000
161 shower.fit.mestimator 3)
163eval $JPP_DIR/examples/JReconstruction/JORCAReconstruction \
170 -N ${NUMBER_OF_THREADS} \
171 -@ `make_equation PARAMETERS` \