4 # first declare default values
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
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)
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 (set to negative to skip dumps)
27 SIMCONFIG="$ALICE_ROOT/MUON/"$MC"Config.C"
28 EMBEDWITH="" # no embedding by default
29 REALISTIC=0 # ideal simulation by default
31 # next try to see if there are options of this script that want to change the
36 while getopts "SRZX:srxzn:tg:p:d:c:e:b:" option
39 R ) RECONSTRUCTION=1;;
46 r ) RECONSTRUCTION=0;;
51 c ) SIMCONFIG=$OPTARG;;
55 p ) RECOPTIONS=$OPTARG;;
56 e ) EMBEDWITH=$OPTARG;;
61 * ) echo "Unimplemented option chosen."
67 if [ ! -n "$OUTDIR" ]; then
68 if [ "$MC" = "" ]; then
69 OUTDIR=$CURDIR"/test_out."$NEVENTS
71 OUTDIR=$CURDIR"/"$MC"_test_out."$NEVENTS
75 # look if there are some leftover options
76 shift $(($OPTIND - 1))
78 if [ $# -gt 0 ] || [ "$EXIT" -eq 1 ]; then
79 echo "ERROR : extra option not recognized"
80 echo "Usage: `basename $0` options (-SRXsrxn:tg:p:d:c:)"
81 echo " -S (-s) perform (or not) simulation (default is do it, i.e -S)"
82 echo " -R (-r) perform (or not) reconstruction (default is do it, i.e. -R)"
83 echo " -X event (-x) perform (or not) checks and dumps (default is do it for event $DUMPEVENT, i.e. -X $DUMPEVENT)"
84 echo " -Z (-z) perform reconstruction from raw data (from digits) (default is from raw data, i.e. -Z)"
85 echo " -n nevents (int) number of events to simulate (default $NEVENTS)"
86 echo " -t will use OUTDIR as a tmp directory to generate raw data "
87 echo " -g seed (uint) seed to be used in simulation (default $SEED)"
88 echo " -p recoptions (quotified string) reconstruction options to use (default \"$RECOPTIONS\")"
89 echo " -d full path to output directory (default $OUTDIR)"
90 echo " -c full path to configuration file for simulation (default $SIMCONFIG)"
91 echo " -e full path to a galice.root file relating to SDigits to be merged (embedding)"
92 echo " -b runnumber (int) make a realistic simulation using runnumber as anchor (default 0=ideal simulation)"
97 # printout the options
98 echo "sim $SIMULATION rec $RECONSTRUCTION check $CHECKS"
99 if [ "$SIMULATION" -eq 1 ]; then
100 echo "$NEVENTS events will be simulated, using the config found at $SIMCONFIG"
102 if [ -n "$EMBEDWITH" ]; then
103 echo "Will embed simulation with $EMBEDWITH"
105 if [ "$REALISTIC" -gt 0 ]; then
106 echo "Will use anchor run $REALISTIC"
108 if [ "$RECONSTRUCTION" -eq 1 ]; then
109 echo "Reconstruction options to be used : $RECOPTIONS"
110 if [ "$RAW" -eq 0 ]; then
111 echo "Will reconstruct from digits only (not from raw data)"
114 echo "Output directory will be : $OUTDIR"
116 if [ "$SIMULATION" -eq 1 ]; then
123 # Copy *ALL* the macros we need in the output directory, not to mess
124 # with our source dir in any way.
125 cp $ALICE_ROOT/MUON/.rootrc \
126 $ALICE_ROOT/MUON/rootlogon.C \
127 $ALICE_ROOT/MUON/runReconstruction.C $ALICE_ROOT/MUON/runSimulation.C \
128 $ALICE_ROOT/MUON/UpdateCDBCTPConfig.C \
129 $ALICE_ROOT/MUON/MUONefficiency.C \
130 $ALICE_ROOT/MUON/MUONTriggerEfficiency.C \
131 $ALICE_ROOT/MUON/MUONCheck.C \
136 if [ "$SLASHTMP" -eq 0 ]; then
140 mkdir ./tmp/mdc1/tags
145 chmod 777 ./tmp/mdc1/tags
147 export ALIMDC_RAWDB1=./tmp/mdc1
148 export ALIMDC_RAWDB2=./tmp/mdc2
149 export ALIMDC_TAGDB=./tmp/mdc1/tags
152 ###############################################################################
154 # Update CTP in OCDB for MUON Trigger
156 ###############################################################################
158 if [ ! -f $ALICE_ROOT/OCDB/GRP/CTP/Config/Run0_999999999_v0_s1.root ]; then
160 echo "Updating GRP CTP config ..."
162 aliroot -b > $OUTDIR/updateCDBCTPConfig.out 2>&1 << EOF
163 .L UpdateCDBCTPConfig.C++g
164 UpdateCDBCTPConfig();
171 ###############################################################################
173 # Performing SIMULATION
175 ###############################################################################
177 if [ "$SIMULATION" -eq 1 ]; then
179 echo "Running simulation ..."
181 aliroot -l -b -q runSimulation.C\($SEED,$NEVENTS,\""$SIMCONFIG"\"\,\""$EMBEDWITH"\"\,$REALISTIC\) > $OUTDIR/testSim.out 2>&1
183 mkdir $OUTDIR/$SIMDIR
185 if [ "$RAW" -eq 1 ]; then
186 if [ "$REALISTIC" -eq 0 ]; then # we can not move for realistic simulations as we need e.g. kinematics to propagate the simulated vertex to the reco.
187 echo "Moving generated files to $SIMDIR"
188 mv $OUTDIR/*QA*.root $OUTDIR/*.log $OUTDIR/$SIMDIR
189 mv $OUTDIR/MUON*.root $OUTDIR/TrackRefs*.root $OUTDIR/$SIMDIR
190 mv $OUTDIR/Kinematics*.root $OUTDIR/galice.root $OUTDIR/$SIMDIR
193 echo "Copying generated files to $SIMDIR"
194 cp $OUTDIR/*QA*.root $OUTDIR/*.log $OUTDIR/$SIMDIR
195 cp $OUTDIR/MUON*.root $OUTDIR/Kinematics*.root $OUTDIR/galice.root $OUTDIR/TrackRefs*.root $OUTDIR/$SIMDIR
198 # save geometry file in a separate directory
199 if [ "$MC" = "g3" ]; then
200 rm -fr $ALICE_ROOT/MUON/geometry
201 mkdir $ALICE_ROOT/MUON/geometry
202 cp $OUTDIR/geometry.root $ALICE_ROOT/MUON/geometry
205 # copy input geometry file in a current directory
206 if [ "$MC" = "g4" ]; then
207 cp $ALICE_ROOT/MUON/geometry/geometry.root $OUTDIR
210 cp $OUTDIR/geometry.root $OUTDIR/$SIMDIR/geometry.root
214 ###############################################################################
216 # Performing RECONSTRUCTION
218 ###############################################################################
220 if [ "$RECONSTRUCTION" -eq 1 ]; then
222 if [ "$RAW" -eq 1 ]; then
223 if [ "$REALISTIC" -eq 0 ]; then
228 if [ "$REALISTIC" -ne 0 ]; then
232 rm -f AliESD*.root *QA*.root
234 echo "Running reconstruction ..."
240 if [ -n "$EMBEDWITH" ]; then
244 if [ "$REALISTIC" -gt 0 ]; then
248 if [ "$RAW" -eq 1 ]; then
250 aliroot -l -b -q runReconstruction\.C\($SEED,\""$OUTDIR/raw.root"\",\""$RECOPTIONS"\",$RAWOCDB\) > $OUTDIR/testReco.out 2>&1
254 aliroot -l -b -q runReconstruction\.C\($SEED,\"""\",\""$RECOPTIONS"\",$RAWOCDB\) > $OUTDIR/testReco.out 2>&1
260 ###############################################################################
262 # Performing CHECKS (and dumps)
264 ###############################################################################
266 if [ "$CHECKS" -eq 1 ]; then
268 if [ -f "$OUTDIR/$SIMDIR/galice.root" ]; then
270 echo "Running efficiency ..."
272 aliroot -b > $OUTDIR/testResults.out 2>&1 << EOF
273 .L MUONefficiency.C++g
274 // no argument assumes Upsilon but MUONefficiency(443) works on Jpsi
275 MUONefficiency("$OUTDIR/$SIMDIR/galice.root");
279 if [ -f "$OUTDIR/galice.root" ]; then
281 echo "Running Trigger efficiency ..."
282 aliroot -b > $OUTDIR/testTriggerResults.out 2>&1 << EOF
283 .L MUONTriggerEfficiency.C++g
284 MUONTriggerEfficiency("$OUTDIR/$SIMDIR/galice.root", "$OUTDIR/galice.root", 1);
288 if [ -f "$OUTDIR/AliESDs.root" ]; then
290 echo "Running check ..."
291 aliroot -b > $OUTDIR/testCheck.out 2>&1 << EOF
292 gSystem->Load("libMUONevaluation");
294 MUONCheck(0, $NEVENTS-1, "$OUTDIR/$SIMDIR/galice.root", "$OUTDIR/galice.root", "$OUTDIR/AliESDs.root");
301 if [ "$DUMPEVENT" -ge 0 ]; then
303 echo "Running dumps for selected event ($DUMPEVENT) ..."
305 if [ -f "$OUTDIR/$SIMDIR/galice.root" ]; then
307 AliCDBManager* man = AliCDBManager::Instance();
308 man->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
309 AliMUONMCDataInterface mcdSim("$OUTDIR/$SIMDIR/galice.root");
310 mcdSim.DumpKine($DUMPEVENT); > $OUTDIR/dump.$DUMPEVENT.kine
311 mcdSim.DumpHits($DUMPEVENT); > $OUTDIR/dump.$DUMPEVENT.hits
312 mcdSim.DumpTrackRefs($DUMPEVENT); > $OUTDIR/dump.$DUMPEVENT.trackrefs
313 mcdSim.DumpDigits($DUMPEVENT,true); > $OUTDIR/dump.$DUMPEVENT.simdigits
314 mcdSim.DumpSDigits($DUMPEVENT,true); > $OUTDIR/dump.$DUMPEVENT.sdigits
318 echo "$OUTDIR/$SIMDIR/galice.root is not there. Skipping sim dumps"
321 if [ -f "$OUTDIR/galice.root" ]; then
323 AliCDBManager* man = AliCDBManager::Instance();
324 man->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
325 AliMUONDataInterface dRec("$OUTDIR/galice.root");
326 dRec.DumpDigits($DUMPEVENT,true); > $OUTDIR/dump.$DUMPEVENT.recdigits
327 dRec.DumpRecPoints($DUMPEVENT); > $OUTDIR/dump.$DUMPEVENT.recpoints
328 dRec.DumpTrigger($DUMPEVENT); > $OUTDIR/dump.$DUMPEVENT.trigger
332 echo "$OUTDIR/galice.root is not there. Skipping rec dumps"
339 echo "... see results in $OUTDIR"