Jpp
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
JORCAShowerReconstruction.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 #
3 # \author adomi
4 #
5 script=${0##*/}
6 
7 
8 # ------------------------------------------------------------------------------------------
9 #
10 # Script to run the JPP Shower Reconstruction Chain for ORCA.
11 #
12 # ------------------------------------------------------------------------------------------
13 
14 
15 if [ -z $JPP_DIR ]; then
16  echo "Variable JPP_DIR undefined."
17  exit
18 fi
19 source $JPP_DIR/setenv.sh $JPP_DIR
20 
21 set_variable DEBUG 2
22 set_variable PDF $JPP_DATA/J%p.dat
23 set_variable PMT $JPP_DATA/PMT_parameters.txt
24 set_variable SUF ${JORCAREC_SUF:-jchain.orca}
25 set_variable NUMBER_OF_EVENTS ${JORCAREC_NUMBER_OF_EVENTS:-$((2**31))}
26 
27 if ( do_usage $* ); then
28  usage "$script <detector file> <input file> <output file> [PDF file descriptor [PMT parameters file]]"
29 fi
30 
31 case $# in
32  5) set_variable PMT $5;&
33  4) set_variable PDF $4;&
34  3) set_variable OUTPUT_FILE $3;
35  set_variable INPUT_FILE $2;
36  set_variable DETECTOR $1;;
37  *) fatal "Wrong number of arguments.";;
38 esac
39 
40 if [[ ! "$OUTPUT_FILE" == *"%"* ]]; then
41  fatal "The output file must have the wildcard % in the name, e.g. %.root"
42 fi
43 
44 # ------------------------------------------------------------------------------------------
45 # JShowerPrefit
46 # ------------------------------------------------------------------------------------------
47 
48 maxOutFits=50
49 
50 # default input values
51 PARAMETERS=(
52  factoryLimit 40
53  NMaxHits 20
54  sigma_ns 4
55  numberOfOutliers 4
56  TMaxLocal_ns 120
57  ctMin -0.7
58  roadWidth_m 50
59  numberOfPrefits 50) #JShowerPointSimplex needs 50 prefits from JShowerPrefit
60 
61 timer_start
62 
63 eval JShowerPrefit \
64  -a $DETECTOR \
65  -f $INPUT_FILE \
66  -o ${OUTPUT_FILE/\%/jshprefit} \
67  -n ${NUMBER_OF_EVENTS} \
68  -@ `make_equation PARAMETERS` \
69  -d ${DEBUG} --!
70 
71 timer_stop
72 timer_print
73 
74 
75 # ------------------------------------------------------------------------------------------
76 # JShowerPointSimplex
77 # ------------------------------------------------------------------------------------------
78 
79 # default input values
80 PARAMETERS=(
81  numberOfPrefits 0 # 0 means all inputs are processed
82  numberOfOutfits 100
83  minPrefitsSize 9
84  sigma_ns 1
85  TMaxLocal_ns 15
86  TMax_ns 60
87  TMin_ns -60
88  roadWidth_m 60
89  ctMin 0.2
90  TWindow_ns "-130, 120"
91  pos_grid_m " -20, 20"
92  pos_step_m 20
93  time_grid_ns " -90, 90"
94  time_step_ns 30
95  simplex_step_m 2
96  simplex_step_ns 10
97  tukey_std_dev 20.0)
98 
99 timer_start
100 
101 eval JShowerPointSimplex \
102  -a $DETECTOR \
103  -f ${OUTPUT_FILE/\%/jshprefit} \
104  -o ${OUTPUT_FILE/\%/jshpointsimplex} \
105  -@ `make_equation PARAMETERS` \
106  -d ${DEBUG} --!
107 
108 timer_stop
109 timer_print
110 
111 rm ${OUTPUT_FILE/\%/jshprefit}
112 
113 # ------------------------------------------------------------------------------------------
114 # JShowerPositionFit
115 # ------------------------------------------------------------------------------------------
116 
117 # default input values
118 PARAMETERS=(
119  numberOfPrefits 0
120  numberOfOutfits 5
121  TMax_ns 120
122  TMin_ns -120
123  roadWidth_m 80
124  E_GeV 10
125  R_Hz 10000
126  TTS_ns 2)
127 
128 timer_start
129 
130 eval JShowerPositionFit \
131  -a $DETECTOR \
132  -f ${OUTPUT_FILE/\%/jshpointsimplex} \
133  -F $PDF \
134  -o ${OUTPUT_FILE/\%/jshpositionfit} \
135  -@ `make_equation PARAMETERS` \
136  -d ${DEBUG} --!
137 
138 timer_stop
139 timer_print
140 
141 rm ${OUTPUT_FILE/\%/jshpointsimplex}
142 
143 # ------------------------------------------------------------------------------------------
144 # JShowerEnergyPrefit
145 # ------------------------------------------------------------------------------------------
146 
147 # default input values
148 PARAMETERS=(
149  numberOfPrefits 0
150  numberOfOutfits 1
151  Emin_GeV 1
152  Emax_GeV 100
153  resolution 0.01
154  TMax_ns 500
155  TMin_ns -500
156  R_Hz 10000
157  mestimator 3
158  roadWidth_m 80)
159 
160 timer_start
161 
162 eval JShowerEnergyPrefit \
163  -a $DETECTOR \
164  -f ${OUTPUT_FILE/\%/jshpositionfit} \
165  -F $PDF \
166  -o ${OUTPUT_FILE/\%/jshenergyprefit} \
167  -@ `make_equation PARAMETERS` \
168  -d $DEBUG --!
169 
170 timer_stop
171 timer_print
172 
173 rm ${OUTPUT_FILE/\%/jshpositionfit}
174 
175 # ------------------------------------------------------------------------------------------
176 # JShowerDirectionPrefit
177 # ------------------------------------------------------------------------------------------
178 
179 # default input values
180 PARAMETERS=(
181  numberOfPrefits 0
182  numberOfOutfits 4
183  TMax_ns 30
184  TMin_ns -30
185  roadWidth_m 80
186  DMin_m 80 # can be lower to speed up the reconstruction (for online)
187  E_GeV 13
188  R_Hz 10000
189  mestimator 3 # 0=normal, 1=lorentzian, 2=linear, 3=null
190  MaxAngle_deg_lowE 40
191  scanAngle_deg_lowE 10
192  MaxAngle_deg_highE 20
193  scanAngle_deg_highE 5)
194 
195 timer_start
196 
197 eval JShowerDirectionPrefit \
198  -a $DETECTOR \
199  -f ${OUTPUT_FILE/\%/jshenergyprefit} \
200  -F $PDF \
201  -o ${OUTPUT_FILE/\%/jshdirectionprefit} \
202  -@ `make_equation PARAMETERS` \
203  -d $DEBUG --!
204 
205 timer_stop
206 timer_print
207 
208 rm ${OUTPUT_FILE/\%/jshenergyprefit}
209 
210 # ------------------------------------------------------------------------------------------
211 # JShowerFit
212 # ------------------------------------------------------------------------------------------
213 
214 # default input values
215 PARAMETERS=(
216  numberOfPrefits 4
217  TMax_ns 30
218  TMin_ns -30
219  roadWidth_m 80
220  R_Hz 10000
221  mestimator 3 # 0=normal, 1=lorentzian, 2=linear, 3=null
222  fit_step 0.05)
223 
224 timer_start
225 
226 eval JShowerFit \
227  -a $DETECTOR \
228  -f ${OUTPUT_FILE/\%/jshdirectionprefit} \
229  -F $PDF \
230  -o ${OUTPUT_FILE/\%/jshfit} \
231  -@ `make_equation PARAMETERS` \
232  -E $JPP_DIR/software/JReconstruction/energy_correction_showerfit_fullorca.txt \
233  -d $DEBUG --!
234 
235 timer_stop
236 timer_print
237 
238 rm ${OUTPUT_FILE/\%/jshdirectionprefit}
239 
240 # ------------------------------------------------------------------------------------------
241 # rename outputfile
242 # ------------------------------------------------------------------------------------------
243 
244 timer_start
245 
246 mv \
247  ${OUTPUT_FILE/\%/jshfit} \
248  ${OUTPUT_FILE/\%/${SUF}} \
249  -v
250 
251 timer_stop
252 timer_print
253 
254 
255 # ------------------------------------------------------------------------------------------
256 # JConvertEvt
257 # ------------------------------------------------------------------------------------------
258 
259 timer_start
260 
261 JConvertEvt \
262  -f ${OUTPUT_FILE/\%/${SUF}} \
263  -o ${OUTPUT_FILE/\%/${SUF}.aanet} \
264  -a $DETECTOR \
265  -P $PMT \
266  -d $DEBUG --!
267 
268 timer_stop
269 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)
lorentzian
then for mestimator in normal linear lorentzian
Definition: JMEstimator.sh:34
exit
exit
Definition: JPizza.sh:36
f
do set_array DAQHEADER JPrintDAQHeader f
Definition: JTuneHV.sh:79
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
F
then awk F
Definition: getMIIseacurrent.sh:82
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 $JPP_DIR software JCalibrate JCalibrateToT a
Definition: JTuneHV.sh:108
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
null
esac $JPP_BIN JLogger sh $LOGGER until pgrep JGetMessage</dev/null > dev null
Definition: examples/JLogger/JLogger.sh:33
E
then usage $script[input file[working directory[option]]] nWhere option can be E
Definition: JMuonPostfit.sh:37
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