Jpp
 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 2
20 set_variable PDF $JPP_DATA/J%p.dat
21 set_variable PMT $JPP_DATA/PMT_parameters.txt
22 set_variable SUF ${JORCAREC_SUF:-jchain.orca}
23 set_variable NUMBER_OF_EVENTS ${JORCAREC_NUMBER_OF_EVENTS:-$((2**31))}
24 
25 if ( do_usage $* ); then
26  usage "$script <detector file> <input file> <output file> [PDF file descriptor [PMT parameters file]]"
27 fi
28 
29 case $# in
30  5) set_variable PMT $5;&
31  4) set_variable PDF $4;&
32  3) set_variable OUTPUT_FILE $3;
33  set_variable INPUT_FILE $2;
34  set_variable DETECTOR $1;;
35  *) fatal "Wrong number of arguments.";;
36 esac
37 
38 if [[ ! "$OUTPUT_FILE" == *"%"* ]]; then
39  fatal "The output file must have the wildcard % in the name, e.g. %.root"
40 fi
41 
42 
43 # ------------------------------------------------------------------------------------------
44 # JMuonPrefit
45 # ------------------------------------------------------------------------------------------
46 
47 PARAMETERS=(
48  sigma_ns 5
49  gridAngle_deg 5
50  useL0 1
51  numberOfOutliers 2
52  numberOfPrefits 12
53  DZMax 0.0
54  numberOfDZMax 1
55  TMaxLocal_ns 12
56  ctMin 0.0
57  roadWidth_m 50.0)
58 
59 timer_start
60 
61 eval JMuonPrefit \
62  -a $DETECTOR \
63  -f $INPUT_FILE \
64  -o ${OUTPUT_FILE/\%/jprefit} \
65  -n ${NUMBER_OF_EVENTS} \
66  -@ `make_equation PARAMETERS` \
67  -d $DEBUG --!
68 
69 CHECK_EXIT_CODE
70 
71 timer_stop
72 timer_print
73 
74 
75 # ------------------------------------------------------------------------------------------
76 # JMuonSimplex
77 # ------------------------------------------------------------------------------------------
78 
79 PARAMETERS=(
80  sigma_ns 3
81  TMaxLocal_ns 15.0
82  ctMin 0.0
83  roadWidth_m 50.0
84  numberOfPrefits 0)
85 
86 timer_start
87 
88 eval JMuonSimplex \
89  -a $DETECTOR \
90  -f ${OUTPUT_FILE/\%/jprefit} \
91  -o ${OUTPUT_FILE/\%/jsimplex} \
92  -@ `make_equation PARAMETERS` \
93  -d ${DEBUG} --!
94 
95 CHECK_EXIT_CODE
96 
97 timer_stop
98 timer_print
99 
100 
101 # ------------------------------------------------------------------------------------------
102 # JMuonStart
103 # ------------------------------------------------------------------------------------------
104 
105 PARAMETERS=(
106  roadWidth_m 50.0
107  R_Hz 7.0e3
108  TMin_ns -25.0
109  TMax_ns +25.0
110  Pmin1 1.0e-3
111  Pmin2 1.0e-2
112  numberOfPrefits 0)
113 
114 timer_start
115 
116 eval JMuonStart \
117  -a $DETECTOR \
118  -f ${OUTPUT_FILE/\%/jsimplex} \
119  -o ${OUTPUT_FILE/\%/jstart} \
120  -P $PDF \
121  -@ `make_equation PARAMETERS` \
122  -d $DEBUG --!
123 
124 CHECK_EXIT_CODE
125 
126 timer_stop
127 timer_print
128 
129 
130 # ------------------------------------------------------------------------------------------
131 # JMuonGandalf
132 # ------------------------------------------------------------------------------------------
133 
134 PARAMETERS=(
135  TTS_ns 2
136  E_GeV 10.0
137  TMin_ns -50.0
138  TMax_ns +450.0
139  ZMin_m -3.0
140  ZMax_m +9.0
141  R_Hz 7.0e3
142  roadWidth_m 50.0
143  numberOfPrefits 0)
144 
145 timer_start
146 
147 eval JMuonGandalf \
148  -a $DETECTOR \
149  -f ${OUTPUT_FILE/\%/jstart} \
150  -o ${OUTPUT_FILE/\%/jgandalf} \
151  -P $PDF \
152  -@ `make_equation PARAMETERS` \
153  -d $DEBUG --!
154 
155 CHECK_EXIT_CODE
156 
157 timer_stop
158 timer_print
159 
160 
161 # ------------------------------------------------------------------------------------------
162 # JMuonStart
163 # ------------------------------------------------------------------------------------------
164 
165 PARAMETERS=(
166  roadWidth_m 50.0
167  R_Hz 7.0e3
168  TMin_ns -25.0
169  TMax_ns +25.0
170  Pmin1 1.0e-3
171  Pmin2 1.0e-2
172  numberOfPrefits 1)
173 
174 timer_start
175 
176 eval JMuonStart \
177  -a $DETECTOR \
178  -f ${OUTPUT_FILE/\%/jgandalf} \
179  -o ${OUTPUT_FILE/\%/jstart2} \
180  -P $PDF \
181  -@ `make_equation PARAMETERS` \
182  -d $DEBUG --!
183 
184 CHECK_EXIT_CODE
185 
186 timer_stop
187 timer_print
188 
189 
190 # ------------------------------------------------------------------------------------------
191 # JMuonEnergy
192 # ------------------------------------------------------------------------------------------
193 
194 PARAMETERS=(
195  roadWidth_m 50.0
196  R_Hz 7.0e3
197  EMin_log 0.0
198  EMax_log 7.0
199  TMin_ns -25.0
200  TMax_ns +25.0
201  ZMin_m -3.0
202  resolution 0.01
203  mestimator 3
204  numberOfPrefits 1)
205 
206 timer_start
207 
208 eval JMuonEnergy \
209  -a $DETECTOR \
210  -f ${OUTPUT_FILE/\%/jstart2} \
211  -o ${OUTPUT_FILE/\%/jenergy} \
212  -P $PDF \
213  -@ `make_equation PARAMETERS` \
214  -d $DEBUG --!
215 
216 CHECK_EXIT_CODE
217 
218 timer_stop
219 timer_print
220 
221 
222 # ------------------------------------------------------------------------------------------
223 # rename outputfile and remove intermediate files
224 # ------------------------------------------------------------------------------------------
225 
226 timer_start
227 
228 mv \
229  ${OUTPUT_FILE/\%/jenergy} \
230  ${OUTPUT_FILE/\%/${SUF}} \
231  -v
232 
233 rm \
234  ${OUTPUT_FILE/\%/jprefit} \
235  ${OUTPUT_FILE/\%/jsimplex} \
236  ${OUTPUT_FILE/\%/jstart} \
237  ${OUTPUT_FILE/\%/jgandalf} \
238  ${OUTPUT_FILE/\%/jstart2} \
239  -v
240 
241 timer_stop
242 timer_print
243 
244 
245 # ------------------------------------------------------------------------------------------
246 # JConvertEvt
247 # ------------------------------------------------------------------------------------------
248 
249 timer_start
250 
251 JConvertEvt \
252  -f ${OUTPUT_FILE/\%/${SUF}} \
253  -o ${OUTPUT_FILE/\%/${SUF}.aanet} \
254  -a $DETECTOR \
255  -P $PMT \
256  -d $DEBUG --!
257 
258 timer_stop
259 timer_print
DETECTOR
set_variable DETECTOR
Definition: JCalibrateMuon.sh:31
PDF
then usage $script< detector file >< inputfile >< outputfile > fi case set_variable PDF
Definition: JCalibrateMuon.sh:28
JSYSTEM::set_variable
void set_variable(const std::string &name, const std::string &value)
Set environment variable.
Definition: JSystemToolkit.hh:379
JDATABASE::PBS::PMT
static const JPBS_t PMT(3, 4, 2, 3)
PBS of photo-multiplier tube (PMT)
exit
exit
Definition: JPizza.sh:36
f
do set_array DAQHEADER JPrintDAQHeader f
Definition: JTuneHV.sh:74
SUF
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 PMT $JPP_DATA PMT_parameters txt set_variable SUF
Definition: JARCAMuonReconstruction.sh:22
echo
then echo
Definition: JDOMDAQDriver.sh:56
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
OUTPUT_FILE
&set_variable OUTPUT_FILE
Definition: JCalibrateMuon.sh:29
o
do $DIR JTransitTime o
Definition: JTransitTime.sh:44
usage
* usage
Definition: JLigier-local.sh:33
$INPUT_FILE
then JFileTuna f $INPUT_FILE
Definition: file-Tuna.sh:66
d
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonPath.sh:45
a
then JCalibrateToT a
Definition: JTuneHV.sh:103
n
alias put_queue eval echo n
Definition: qlib.csh:19
PARAMETERS
PARAMETERS
Definition: JAcoustics.sh:52
JTOOLS::v
data_type v[N+1][M+1]
Definition: JPolint.hh:740
in
then fatal Wrong number of arguments fi set_variable DETECTOR $argv[1] set_variable INPUT_FILE $argv[2] eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY source JAcoustics sh $DETECTOR_ID typeset A TRIPODS get_tripods $WORKDIR tripod txt TRIPODS for EMITTER in
Definition: JCanberra.sh:36
script
script
Definition: JAcoustics.sh:2
INPUT_FILE
set_variable INPUT_FILE
Definition: JCalibrateMuon.sh:30
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 case set_variable NUMBER_OF_EVENTS
Definition: JRunAnalyzer.sh:17
P
then $DIR JPlotNPE PDG P
Definition: JPlotNPE-PDG.sh:60
DEBUG
#define DEBUG(A)
Message macros.
Definition: JMessage.hh:62
Generated by   doxygen 1.8.5