Jpp test-rotations-old
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
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: ENERGY_CORRECTION RECONSTRUCTION_ENERGY_CORRECTION energy_correction_arca.txt
24set_variable: NUMBER_OF_EVENTS RECONSTRUCTION_NUMBER_OF_EVENTS $((2**31))
25set_variable: NUMBER_OF_THREADS RECONSTRUCTION_NUMBER_OF_THREADS 0
26set_variable+ OVERWRITE RECONSTRUCTION_OVERWRITE 1
27set_variable: TMAX_S RECONSTRUCTION_TMAX_S 100
28set_variable: RATE_HZ RECONSTRUCTION_RATE_HZ 6.0e3
29
30if do_usage $*; then
31 usage "$script <detector file> <input file> <output file> <PMT parameters file> [(calibration file)+]"
32fi
33
34if (( $# < 4 )); then
35 fatal "Wrong number of arguments."
36fi
37
38set_variable DETECTOR $argv[1]
39set_variable INPUT_FILE $argv[2]
40set_variable OUTPUT_FILE $argv[3]
41set_variable PMT $argv[4]
42set_array CALIBRATION $argv[5,-1]
43
44if (( ${#CALIBRATION} != 0 )); then
45 getMechanics.sh $DETECTOR
46fi
47
48typeset -A PARAMETERS
49
50PARAMETERS+=(
51 prefit.sigma_ns 5
52 prefit.gridAngle_deg 1
53 prefit.useL0 0
54 prefit.numberOfOutliers 3
55 prefit.numberOfPrefits 50
56 prefit.DZMax 0.0
57 prefit.numberOfDZMax 1
58 prefit.TMaxLocal_ns 18
59 prefit.ctMin 0.0
60 prefit.roadWidth_m 200.0)
61
62PARAMETERS+=(
63 simplex.sigma_ns 3
64 simplex.TMaxLocal_ns 15.0
65 simplex.ctMin 0.0
66 simplex.roadWidth_m 200.0
67 simplex.numberOfPrefits 0)
68
69PARAMETERS+=(
70 gandalf.TTS_ns 2
71 gandalf.E_GeV 1.0e3
72 gandalf.TMin_ns -50.0
73 gandalf.TMax_ns +450.0
74 gandalf.R_Hz $RATE_HZ
75 gandalf.roadWidth_m 175.0
76 gandalf.numberOfPrefits 0)
77
78PARAMETERS+=(
79 start.roadWidth_m 160.0
80 start.R_Hz $RATE_HZ
81 start.TMin_ns -25.0
82 start.TMax_ns +75.0
83 start.Pmin1 1.0e-3
84 start.Pmin2 1.0e-2
85 start.numberOfPrefits 1)
86
87PARAMETERS+=(
88 energy.roadWidth_m 300.0
89 energy.R_Hz $RATE_HZ
90 energy.EMin_log 0.0
91 energy.EMax_log 8.0
92 energy.TMin_ns -50.0
93 energy.TMax_ns +450.0
94 energy.mestimator 1
95 energy.numberOfPrefits 1)
96
97eval $JPP_DIR/examples/JReconstruction/JARCAMuonReconstruction \
98 -a $DETECTOR \
99 -+${^CALIBRATION} \
100 -T $TMAX_S \
101 -f ${INPUT_FILE} \
102 -o ${OUTPUT_FILE} \
103 -P $PDF \
104 -E $ENERGY_CORRECTION \
105 -n ${NUMBER_OF_EVENTS} \
106 -N ${NUMBER_OF_THREADS} \
107 -@ `make_equation PARAMETERS` \
108 -d $DEBUG --!
109