]>
Commit | Line | Data |
---|---|---|
1 | #!/bin/sh | |
2 | # $Id$ | |
3 | ||
4 | # first declare default values | |
5 | ||
6 | SIMULATION=1 # will perform simulation | |
7 | RECONSTRUCTION=1 # will perform reconstruction | |
8 | RAW=1 # will reconstruct from raw data | |
9 | CHECKS=1 # will perform checks | |
10 | SLASHTMP=1 #will use /tmp to put the temporary raw data | |
11 | NEVENTS=50 # will simulate 100 events | |
12 | ||
13 | #RECOPTIONS="SAVEDIGITS NOFASTDECODERS" # reconstruction options with non-high performance decoders | |
14 | RECOPTIONS="SAVEDIGITS" # default reconstruction options | |
15 | MC="" # G3 Simulation with old Config.C | |
16 | #MC="g3" # G3 Simulation (with new config macros) | |
17 | #MC="g4" # G4 Simulation (with new config macros) | |
18 | ||
19 | OUTDIR="" | |
20 | CURDIR=`pwd` | |
21 | ||
22 | #RUN=0 # run number for OCDB access | |
23 | SEED=1234567 # random number generator seed | |
24 | SIMDIR="generated" # sub-directory where to move simulated files prior to reco | |
25 | DUMPEVENT=5 # event to be dump on files | |
26 | ||
27 | SIMCONFIG="$ALICE_ROOT/MUON/"$MC"Config.C" | |
28 | EMBEDWITH="" # no embedding by default | |
29 | ||
30 | # next try to see if there are options of this script that want to change the | |
31 | # defaults | |
32 | ||
33 | EXIT=0 | |
34 | ||
35 | while getopts "SRZX:srxzn:tg:p:d:c:e:" option | |
36 | do | |
37 | case $option in | |
38 | R ) RECONSTRUCTION=1;; | |
39 | S ) SIMULATION=1;; | |
40 | X ) | |
41 | CHECKS=1 | |
42 | DUMPEVENT=$OPTARG | |
43 | ;; | |
44 | Z ) RAW=1;; | |
45 | r ) RECONSTRUCTION=0;; | |
46 | s ) SIMULATION=0;; | |
47 | x ) CHECKS=0;; | |
48 | t ) SLASHTMP=0;; | |
49 | z ) RAW=0;; | |
50 | c ) SIMCONFIG=$OPTARG;; | |
51 | d ) OUTDIR=$OPTARG;; | |
52 | n ) NEVENTS=$OPTARG;; | |
53 | g ) SEED=$OPTARG;; | |
54 | p ) RECOPTIONS=$OPTARG;; | |
55 | e ) EMBEDWITH=$OPTARG;; | |
56 | * ) echo "Unimplemented option chosen." | |
57 | EXIT=1 | |
58 | ;; | |
59 | esac | |
60 | done | |
61 | ||
62 | if [ ! -n "$OUTDIR" ]; then | |
63 | if [ "$MC" = "" ]; then | |
64 | OUTDIR=$CURDIR"/test_out."$NEVENTS | |
65 | else | |
66 | OUTDIR=$CURDIR"/"$MC"_test_out."$NEVENTS | |
67 | fi | |
68 | fi | |
69 | ||
70 | # look if there are some leftover options | |
71 | shift $(($OPTIND - 1)) | |
72 | ||
73 | if [ $# -gt 0 ] || [ "$EXIT" -eq 1 ]; then | |
74 | echo "ERROR : extra option not recognized" | |
75 | echo "Usage: `basename $0` options (-SRXsrxn:tg:p:d:c:)" | |
76 | echo " -S (-s) perform (or not) simulation (default is do it, i.e -S)" | |
77 | echo " -R (-r) perform (or not) reconstruction (default is do it, i.e. -R)" | |
78 | echo " -X event (-x) perform (or not) checks and dumps (default is do it for event $DUMPEVENT, i.e. -X $DUMPEVENT)" | |
79 | echo " -Z (-z) perform reconstruction from raw data (from digits) (default is from raw data, i.e. -Z)" | |
80 | echo " -n nevents (int) number of events to simulate (default $NEVENTS)" | |
81 | echo " -t will use OUTDIR as a tmp directory to generate raw data " | |
82 | echo " -g seed (uint) seed to be used in simulation (default $SEED)" | |
83 | echo " -p recoptions (quotified string) reconstruction options to use (default \"$RECOPTIONS\")" | |
84 | echo " -d full path to output directory (default $OUTDIR)" | |
85 | echo " -c full path to configuration file for simulation (default $SIMCONFIG)" | |
86 | echo " -e full path to a galice.root file relating to SDigits to be merged (embedding)" | |
87 | exit 4; | |
88 | fi | |
89 | ||
90 | ||
91 | # printout the options | |
92 | echo "sim $SIMULATION rec $RECONSTRUCTION check $CHECKS" | |
93 | if [ "$SIMULATION" -eq 1 ]; then | |
94 | echo "$NEVENTS events will be simulated, using the config found at $SIMCONFIG" | |
95 | fi | |
96 | if [ -n "$EMBEDWITH" ]; then | |
97 | echo "Will embed simulation with $EMBEDWITH" | |
98 | fi | |
99 | if [ "$RECONSTRUCTION" -eq 1 ]; then | |
100 | echo "Reconstruction options to be used : $RECOPTIONS" | |
101 | if [ "$RAW" -eq 0 ]; then | |
102 | echo "Will reconstruct from digits only (not from raw data)" | |
103 | fi | |
104 | fi | |
105 | echo "Output directory will be : $OUTDIR" | |
106 | ||
107 | if [ "$SIMULATION" -eq 1 ]; then | |
108 | ||
109 | rm -fr $OUTDIR | |
110 | mkdir $OUTDIR | |
111 | ||
112 | fi | |
113 | ||
114 | # Copy *ALL* the macros we need in the output directory, not to mess | |
115 | # with our source dir in any way. | |
116 | cp $ALICE_ROOT/MUON/.rootrc \ | |
117 | $ALICE_ROOT/MUON/rootlogon.C \ | |
118 | $ALICE_ROOT/MUON/runReconstruction.C $ALICE_ROOT/MUON/runSimulation.C \ | |
119 | $ALICE_ROOT/MUON/UpdateCDBCTPConfig.C \ | |
120 | $ALICE_ROOT/MUON/MUONefficiency.C \ | |
121 | $ALICE_ROOT/MUON/MUONTriggerEfficiency.C \ | |
122 | $ALICE_ROOT/MUON/MUONCheck.C \ | |
123 | $OUTDIR | |
124 | ||
125 | cd $OUTDIR | |
126 | ||
127 | if [ "$SLASHTMP" -eq 0 ]; then | |
128 | mkdir ./tmp | |
129 | mkdir ./tmp/mdc1 | |
130 | mkdir ./tmp/mdc2 | |
131 | mkdir ./tmp/mdc1/tags | |
132 | ||
133 | chmod 777 ./tmp | |
134 | chmod 777 ./tmp/mdc1 | |
135 | chmod 777 ./tmp/mdc2 | |
136 | chmod 777 ./tmp/mdc1/tags | |
137 | ||
138 | export ALIMDC_RAWDB1=./tmp/mdc1 | |
139 | export ALIMDC_RAWDB2=./tmp/mdc2 | |
140 | export ALIMDC_TAGDB=./tmp/mdc1/tags | |
141 | fi | |
142 | ||
143 | ############################################################################### | |
144 | # | |
145 | # Update CTP in OCDB for MUON Trigger | |
146 | # | |
147 | ############################################################################### | |
148 | ||
149 | if [ ! -f $ALICE_ROOT/OCDB/GRP/CTP/Config/Run0_999999999_v0_s1.root ]; then | |
150 | ||
151 | echo "Updating GRP CTP config ..." | |
152 | ||
153 | aliroot -b > $OUTDIR/updateCDBCTPConfig.out 2>&1 << EOF | |
154 | .L UpdateCDBCTPConfig.C++g | |
155 | UpdateCDBCTPConfig(); | |
156 | .q | |
157 | EOF | |
158 | ||
159 | fi | |
160 | ||
161 | ||
162 | ############################################################################### | |
163 | # | |
164 | # Performing SIMULATION | |
165 | # | |
166 | ############################################################################### | |
167 | ||
168 | if [ "$SIMULATION" -eq 1 ]; then | |
169 | ||
170 | echo "Running simulation ..." | |
171 | ||
172 | aliroot -l -b -q runSimulation.C\($SEED,$NEVENTS,\""$SIMCONFIG"\"\,\""$EMBEDWITH"\"\) > $OUTDIR/testSim.out 2>&1 | |
173 | ||
174 | mkdir $OUTDIR/$SIMDIR | |
175 | ||
176 | if [ "$RAW" -eq 1 ]; then | |
177 | echo "Moving generated files to $SIMDIR" | |
178 | mv $OUTDIR/*QA*.root $OUTDIR/*.log $OUTDIR/$SIMDIR | |
179 | mv $OUTDIR/MUON*.root $OUTDIR/Kinematics*.root $OUTDIR/galice.root $OUTDIR/TrackRefs*.root $OUTDIR/$SIMDIR | |
180 | else | |
181 | echo "Copying generated files to $SIMDIR" | |
182 | cp $OUTDIR/*QA*.root $OUTDIR/*.log $OUTDIR/$SIMDIR | |
183 | cp $OUTDIR/MUON*.root $OUTDIR/Kinematics*.root $OUTDIR/galice.root $OUTDIR/TrackRefs*.root $OUTDIR/$SIMDIR | |
184 | fi | |
185 | ||
186 | # save geometry file in a separate directory | |
187 | if [ "$MC" = "g3" ]; then | |
188 | rm -fr $ALICE_ROOT/MUON/geometry | |
189 | mkdir $ALICE_ROOT/MUON/geometry | |
190 | cp $OUTDIR/geometry.root $ALICE_ROOT/MUON/geometry | |
191 | fi | |
192 | ||
193 | # copy input geometry file in a current directory | |
194 | if [ "$MC" = "g4" ]; then | |
195 | cp $ALICE_ROOT/MUON/geometry/geometry.root $OUTDIR | |
196 | fi | |
197 | fi | |
198 | ||
199 | ############################################################################### | |
200 | # | |
201 | # Performing RECONSTRUCTION | |
202 | # | |
203 | ############################################################################### | |
204 | ||
205 | if [ "$RECONSTRUCTION" -eq 1 ]; then | |
206 | ||
207 | if [ "$RAW" -eq 1 ]; then | |
208 | rm -f galice.root | |
209 | fi | |
210 | ||
211 | rm -f AliESD*.root *QA*.root | |
212 | ||
213 | echo "Running reconstruction ..." | |
214 | ||
215 | cd $OUTDIR | |
216 | ||
217 | BOOLEMBED=kFALSE | |
218 | ||
219 | if [ -n "$EMBEDWITH" ]; then | |
220 | BOOLEMBED=kTRUE | |
221 | fi | |
222 | ||
223 | if [ "$RAW" -eq 1 ]; then | |
224 | ||
225 | aliroot -l -b -q runReconstruction\.C\($SEED,\""$OUTDIR/raw.root"\",\""$RECOPTIONS"\",$BOOLEMBED\) > $OUTDIR/testReco.out 2>&1 | |
226 | ||
227 | else | |
228 | ||
229 | aliroot -l -b -q runReconstruction\.C\($SEED,\"""\",\""$RECOPTIONS"\",$BOOLEMBED\) > $OUTDIR/testReco.out 2>&1 | |
230 | ||
231 | fi | |
232 | ||
233 | fi | |
234 | ||
235 | ############################################################################### | |
236 | # | |
237 | # Performing CHECKS (and dumps) | |
238 | # | |
239 | ############################################################################### | |
240 | ||
241 | if [ "$CHECKS" -eq 1 ]; then | |
242 | ||
243 | if [ -f "$OUTDIR/$SIMDIR/galice.root" ]; then | |
244 | ||
245 | echo "Running efficiency ..." | |
246 | ||
247 | aliroot -b > $OUTDIR/testResults.out 2>&1 << EOF | |
248 | .L MUONefficiency.C++g | |
249 | // no argument assumes Upsilon but MUONefficiency(443) works on Jpsi | |
250 | MUONefficiency("$OUTDIR/$SIMDIR/galice.root"); | |
251 | .q | |
252 | EOF | |
253 | ||
254 | if [ -f "$OUTDIR/galice.root" ]; then | |
255 | ||
256 | echo "Running Trigger efficiency ..." | |
257 | aliroot -b > $OUTDIR/testTriggerResults.out 2>&1 << EOF | |
258 | .L MUONTriggerEfficiency.C++g | |
259 | MUONTriggerEfficiency("$OUTDIR/$SIMDIR/galice.root", "$OUTDIR/galice.root", 1); | |
260 | .q | |
261 | EOF | |
262 | ||
263 | if [ -f "$OUTDIR/AliESDs.root" ]; then | |
264 | ||
265 | echo "Running check ..." | |
266 | aliroot -b > $OUTDIR/testCheck.out 2>&1 << EOF | |
267 | gSystem->Load("libMUONevaluation"); | |
268 | .L MUONCheck.C++g | |
269 | MUONCheck(0, $NEVENTS-1, "$OUTDIR/$SIMDIR/galice.root", "$OUTDIR/galice.root", "$OUTDIR/AliESDs.root"); | |
270 | .q | |
271 | EOF | |
272 | fi | |
273 | fi | |
274 | fi | |
275 | ||
276 | echo "Running dumps for selected event ($DUMPEVENT) ..." | |
277 | ||
278 | if [ -f "$OUTDIR/$SIMDIR/galice.root" ]; then | |
279 | aliroot -l -b << EOF | |
280 | AliCDBManager* man = AliCDBManager::Instance(); | |
281 | man->SetDefaultStorage("local://$ALICE_ROOT/OCDB"); | |
282 | AliMUONMCDataInterface mcdSim("$OUTDIR/$SIMDIR/galice.root"); | |
283 | mcdSim.DumpKine($DUMPEVENT); > $OUTDIR/dump.$DUMPEVENT.kine | |
284 | mcdSim.DumpHits($DUMPEVENT); > $OUTDIR/dump.$DUMPEVENT.hits | |
285 | mcdSim.DumpTrackRefs($DUMPEVENT); > $OUTDIR/dump.$DUMPEVENT.trackrefs | |
286 | mcdSim.DumpDigits($DUMPEVENT,true); > $OUTDIR/dump.$DUMPEVENT.simdigits | |
287 | mcdSim.DumpSDigits($DUMPEVENT,true); > $OUTDIR/dump.$DUMPEVENT.sdigits | |
288 | .q | |
289 | EOF | |
290 | else | |
291 | echo "$OUTDIR/$SIMDIR/galice.root is not there. Skipping sim dumps" | |
292 | fi | |
293 | ||
294 | if [ -f "$OUTDIR/galice.root" ]; then | |
295 | aliroot -l -b << EOF | |
296 | AliCDBManager* man = AliCDBManager::Instance(); | |
297 | man->SetDefaultStorage("local://$ALICE_ROOT/OCDB"); | |
298 | AliMUONDataInterface dRec("$OUTDIR/galice.root"); | |
299 | dRec.DumpDigits($DUMPEVENT,true); > $OUTDIR/dump.$DUMPEVENT.recdigits | |
300 | dRec.DumpRecPoints($DUMPEVENT); > $OUTDIR/dump.$DUMPEVENT.recpoints | |
301 | dRec.DumpTrigger($DUMPEVENT); > $OUTDIR/dump.$DUMPEVENT.trigger | |
302 | .q | |
303 | EOF | |
304 | else | |
305 | echo "$OUTDIR/galice.root is not there. Skipping rec dumps" | |
306 | fi | |
307 | fi | |
308 | ||
309 | echo "Finished" | |
310 | echo "... see results in $OUTDIR" | |
311 | ||
312 | cd $CURDIR |