Jpp  16.0.3
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
JORCAMuonReconstruction.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 #
3 # \author bofearraigh
4 #
5 script=${0##*/}
6 
7 # ------------------------------------------------------------------------------------------
8 #
9 # Script to run the JPP Track Reconstruction Chain for ORCA.
10 #
11 # ------------------------------------------------------------------------------------------
12 
13 if [ -z $JPP_DIR ]; then
14  echo "Variable JPP_DIR undefined."
15  exit
16 fi
17 source $JPP_DIR/setenv.sh $JPP_DIR
18 
19 set_variable: DEBUG RECONSTRUCTION_DEBUG 1
20 set_variable: PDF RECONSTRUCTION_PDF $JPP_DATA/J%p.dat
21 set_variable: SUFFIX RECONSTRUCTION_SUFFIX jchain.orca
22 set_variable: NUMBER_OF_EVENTS RECONSTRUCTION_NUMBER_OF_EVENTS $((2**31))
23 set_variable: TMAX_S RECONSTRUCTION_TMAX_S 100
24 
25 if do_usage $*; then
26  usage "$script <detector file> <input file> <output file> <PMT parameters file> [(calibration file)+]"
27 fi
28 
29 if (( $# < 4 )); then
30  fatal "Wrong number of arguments."
31 fi
32 
33 set_variable DETECTOR $argv[1]
34 set_variable INPUT_FILE $argv[2]
35 set_variable OUTPUT_FILE $argv[3]
36 set_variable PMT $argv[4]
37 set_array CALIBRATION $argv[5,-1]
38 
39 if [[ ! "$OUTPUT_FILE" == *"%"* ]]; then
40  fatal "The output file must have the wildcard % in the name, e.g. %.root"
41 fi
42 
43 typeset -A PARAMETERS
44 
45 # ------------------------------------------------------------------------------------------
46 # JMuonPrefit
47 # ------------------------------------------------------------------------------------------
48 
49 PARAMETERS=(
50  sigma_ns 5
51  gridAngle_deg 5
52  useL0 1
53  numberOfOutliers 2
54  numberOfPrefits 12
55  DZMax 0.0
56  numberOfDZMax 1
57  TMaxLocal_ns 12
58  ctMin 0.0
59  roadWidth_m 50.0)
60 
61 timer_start
62 
63 eval JMuonPrefit \
64  -a $DETECTOR \
65  -+${^CALIBRATION} \
66  -T $TMAX_S \
67  -f $INPUT_FILE \
68  -o ${OUTPUT_FILE/\%/jprefit} \
69  -n ${NUMBER_OF_EVENTS} \
70  -@ `make_equation PARAMETERS` \
71  -d $DEBUG --!
72 
73 CHECK_EXIT_CODE
74 
75 timer_stop
76 timer_print
77 
78 
79 # ------------------------------------------------------------------------------------------
80 # JMuonSimplex
81 # ------------------------------------------------------------------------------------------
82 
83 PARAMETERS=(
84  sigma_ns 3
85  TMaxLocal_ns 15.0
86  ctMin 0.0
87  roadWidth_m 50.0
88  numberOfPrefits 0)
89 
90 timer_start
91 
92 eval JMuonSimplex \
93  -a $DETECTOR \
94  -+${^CALIBRATION} \
95  -T $TMAX_S \
96  -f ${OUTPUT_FILE/\%/jprefit} \
97  -o ${OUTPUT_FILE/\%/jsimplex} \
98  -@ `make_equation PARAMETERS` \
99  -d ${DEBUG} --!
100 
101 CHECK_EXIT_CODE
102 
103 timer_stop
104 timer_print
105 
106 
107 # ------------------------------------------------------------------------------------------
108 # JMuonStart
109 # ------------------------------------------------------------------------------------------
110 
111 PARAMETERS=(
112  roadWidth_m 50.0
113  R_Hz 7.0e3
114  TMin_ns -25.0
115  TMax_ns +25.0
116  Pmin1 1.0e-3
117  Pmin2 1.0e-2
118  numberOfPrefits 0)
119 
120 timer_start
121 
122 eval JMuonStart \
123  -a $DETECTOR \
124  -+${^CALIBRATION} \
125  -T $TMAX_S \
126  -f ${OUTPUT_FILE/\%/jsimplex} \
127  -o ${OUTPUT_FILE/\%/jstart} \
128  -P $PDF \
129  -@ `make_equation PARAMETERS` \
130  -d $DEBUG --!
131 
132 CHECK_EXIT_CODE
133 
134 timer_stop
135 timer_print
136 
137 
138 # ------------------------------------------------------------------------------------------
139 # JMuonGandalf
140 # ------------------------------------------------------------------------------------------
141 
142 PARAMETERS=(
143  TTS_ns 2
144  E_GeV 10.0
145  TMin_ns -50.0
146  TMax_ns +450.0
147  ZMin_m -3.0
148  ZMax_m +9.0
149  R_Hz 7.0e3
150  roadWidth_m 50.0
151  numberOfPrefits 0)
152 
153 timer_start
154 
155 eval JMuonGandalf \
156  -a $DETECTOR \
157  -+${^CALIBRATION} \
158  -T $TMAX_S \
159  -f ${OUTPUT_FILE/\%/jstart} \
160  -o ${OUTPUT_FILE/\%/jgandalf} \
161  -P $PDF \
162  -@ `make_equation PARAMETERS` \
163  -d $DEBUG --!
164 
165 CHECK_EXIT_CODE
166 
167 timer_stop
168 timer_print
169 
170 
171 # ------------------------------------------------------------------------------------------
172 # JMuonStart
173 # ------------------------------------------------------------------------------------------
174 
175 PARAMETERS=(
176  roadWidth_m 50.0
177  R_Hz 7.0e3
178  TMin_ns -25.0
179  TMax_ns +25.0
180  Pmin1 1.0e-3
181  Pmin2 1.0e-2
182  numberOfPrefits 1)
183 
184 timer_start
185 
186 eval JMuonStart \
187  -a $DETECTOR \
188  -+${^CALIBRATION} \
189  -T $TMAX_S \
190  -f ${OUTPUT_FILE/\%/jgandalf} \
191  -o ${OUTPUT_FILE/\%/jstart2} \
192  -P $PDF \
193  -@ `make_equation PARAMETERS` \
194  -d $DEBUG --!
195 
196 CHECK_EXIT_CODE
197 
198 timer_stop
199 timer_print
200 
201 
202 # ------------------------------------------------------------------------------------------
203 # JMuonEnergy
204 # ------------------------------------------------------------------------------------------
205 
206 PARAMETERS=(
207  roadWidth_m 80.0
208  R_Hz 7.0e3
209  EMin_log 0.0
210  EMax_log 7.0
211  TMin_ns -25.0
212  TMax_ns +25.0
213  ZMin_m -3.0
214  resolution 0.01
215  mestimator 3
216  numberOfPrefits 1)
217 
218 timer_start
219 
220 eval JMuonEnergy \
221  -a $DETECTOR \
222  -+${^CALIBRATION} \
223  -T $TMAX_S \
224  -f ${OUTPUT_FILE/\%/jstart2} \
225  -o ${OUTPUT_FILE/\%/jenergy} \
226  -P $PDF \
227  -@ `make_equation PARAMETERS` \
228  -d $DEBUG --!
229 
230 CHECK_EXIT_CODE
231 
232 timer_stop
233 timer_print
234 
235 
236 # ------------------------------------------------------------------------------------------
237 # rename outputfile and remove intermediate files
238 # ------------------------------------------------------------------------------------------
239 
240 timer_start
241 
242 mv -v \
243  ${OUTPUT_FILE/\%/jenergy} \
244  ${OUTPUT_FILE/\%/${SUFFIX}}
245 
246 rm -vf \
247  ${OUTPUT_FILE/\%/jprefit} \
248  ${OUTPUT_FILE/\%/jsimplex} \
249  ${OUTPUT_FILE/\%/jstart} \
250  ${OUTPUT_FILE/\%/jgandalf} \
251  ${OUTPUT_FILE/\%/jstart2}
252 
253 timer_stop
254 timer_print
255 
256 
257 # ------------------------------------------------------------------------------------------
258 # JConvertEvt
259 # ------------------------------------------------------------------------------------------
260 
261 timer_start
262 
263 JConvertEvt \
264  -f ${OUTPUT_FILE/\%/${SUFFIX}} \
265  -o ${OUTPUT_FILE/\%/${SUFFIX}.aanet} \
266  -a $DETECTOR \
267  -+${^CALIBRATION} \
268  -T $TMAX_S \
269  -P $PMT \
270  -d $DEBUG --!
271 
272 timer_stop
273 timer_print
PDF
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR set_variable DEBUG set_variable PDF $JPP_DATA J p dat set_variable HISTOGRAM if do_usage *then usage $script< detector file >< inputfile >< outputfile > fi case set_variable PDF
Definition: JARCACalibrateMuon.sh:23
NUMBER_OF_EVENTS
then usage $script< input_file >< detector_file >< output_file(.root)> fi set_variable NUMBER_OF_TIMESLICES set_variable NUMBER_OF_SUMMARYSLICES set_variable NUMBER_OF_EVENTS set_variable ALL_PLOTS case set_variable NUMBER_OF_EVENTS
Definition: JRunAnalyzer.sh:18
$argv
*set_variable DETECTOR_ID $argv[1] set_array RUNS $argv[2,-2]
Definition: JDataQuality.sh:24
file
then set_variable PMT_FILE set_variable DAQ_FILE set_variable OUTPUT_FILE set_variable DETECTOR else fatal Wrong number of arguments fi set_variable RUNBYRUN file
Definition: JRandomEventRunByRun.sh:46
JDATABASE::PBS::PMT
static const JPBS_t PMT(3, 4, 2, 3)
PBS of photo-multiplier tube (PMT)
f
o $QUALITY_ROOT d $DEBUG!JPlot1D f
Definition: JDataQuality.sh:66
mestimator
mestimator
Definition: acoustics-fit.sh:494
parameters
*fatal Wrong number of arguments esac JCookie sh typeset Z DETECTOR typeset Z SOURCE_RUN typeset Z TARGET_RUN set_variable PARAMETERS_FILE $WORKDIR parameters
Definition: diff-Tuna.sh:38
exit
exit
Definition: JPizza.sh:36
echo
then echo
Definition: JDOMDAQDriver.sh:56
JTOOLS::n
const int n
Definition: JPolint.hh:676
timer_start
then timer_start(JRandomTimesliceWriter-a ${DETECTOR}-o ${PIPE}-n ${NUMBER_OF_SLICES}-P ${PMT_FILE}-d ${DEBUG}-S ${SEED}-B"$BACKGROUND_HZ"-N"$RECYCLING"--!&&for((i=0;$i!=10;++i));do;echo 'EOF' > $PIPE;usleep 100000;done &&rm-f $PIPE)&JTimesliceReprocessor-a $
Definition: JTimesliceReprocessor.sh:57
set_variable
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable
Definition: JAcoustics.sh:20
T
do set_variable OUTPUT_DIRECTORY $WORKDIR T
Definition: JCalibrateHeight.sh:61
usage
* usage
Definition: JLigier-local.sh:34
o
then JConvertDetectorFormat a $DETECTOR[1] o
Definition: JCompareDetector.sh:45
$INPUT_FILE
then JFileTuna f $INPUT_FILE
Definition: file-Tuna.sh:66
INPUT_FILE
set_variable INPUT_FILE
Definition: JARCACalibrateMuon.sh:30
OUTPUT_FILE
&set_variable OUTPUT_FILE
Definition: JARCACalibrateMuon.sh:29
d
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonPath.sh:47
$script
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR set_variable DEBUG set_variable WORKDIR if do_usage *then usage $script(input file)+" fi if (( $ fatal "Invalid number of arguments." fi JCookie.sh set_array INPUT_FILES $argv[1
a
then JCalibrateToT a
Definition: JTuneHV.sh:116
$DETECTOR
then set_variable MODULE getModule a $DETECTOR L $STRING $FLOOR JEditDetector a $DETECTOR M $MODULE add $X o $DETECTOR else echo No update of detector $DETECTOR
Definition: module-Z:fit.sh:118
DEBUG
then usage $script< input_file >< detector_file > fi set_variable OUTPUT_DIR set_variable SELECTOR JDAQTimesliceL1 set_variable DEBUG case set_variable DEBUG
Definition: JCalibratePMT.sh:21
DETECTOR
set_variable DETECTOR
Definition: JARCACalibrateMuon.sh:31
detector
do set_variable DETECTOR_TXT $WORKDIR detector
Definition: detector-XY:fit1d.sh:38
JTOOLS::v
data_type v[N+1][M+1]
Definition: JPolint.hh:756
A
source $JPP_DIR setenv csh $JPP_DIR &dev null eval JShellParser o a A
Definition: JShellParser.csh:15
script
script
Definition: JAcoustics.sh:2
P
then $DIR JPlotNPE PDG P
Definition: JPlotNPE-PDG.sh:62
if
do if[[!-f $ACOUSTICS_WORKDIR/${KEY}.txt]]
Definition: JAcoustics.sh:39
Generated by   doxygen 1.8.5