]>
Commit | Line | Data |
---|---|---|
1 | /************************************************************************** | |
2 | * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
3 | * * | |
4 | * Author: The ALICE Off-line Project. * | |
5 | * Contributors are mentioned in the code where appropriate. * | |
6 | * * | |
7 | * Permission to use, copy, modify and distribute this software and its * | |
8 | * documentation strictly for non-commercial purposes is hereby granted * | |
9 | * without fee, provided that the above copyright notice appears in all * | |
10 | * copies and that both the copyright notice and this permission notice * | |
11 | * appear in the supporting documentation. The authors make no claims * | |
12 | * about the suitability of this software for any purpose. It is * | |
13 | * provided "as is" without express or implied warranty. * | |
14 | **************************************************************************/ | |
15 | /* $Id$ */ | |
16 | ||
17 | //----------------------------------------------------------------------------- | |
18 | /// \class AliMUONReconstructor | |
19 | /// | |
20 | /// Implementation of AliReconstructor for MUON subsystem. | |
21 | /// | |
22 | /// The clustering mode and the associated parameters can be changed through the | |
23 | /// AliMUONRecoParam object set in the reconstruction macro or read from the CDB | |
24 | /// (see methods in AliMUONRecoParam.h file for details) | |
25 | /// | |
26 | /// Valid modes are : | |
27 | /// | |
28 | /// SIMPLEFIT : use the AliMUONClusterFinderSimpleFit clusterizer | |
29 | /// | |
30 | /// SIMPLEFITV3 : SIMPLEFIT with preclustering=PRECLUSTERV3 | |
31 | /// | |
32 | /// MLEM : use AliMUONClusterFinderMLEM and AliMUONPreClusterFinder for preclustering (default) | |
33 | /// MLEMV2 : MLEM with preclustering=PRECLUSTERV2 | |
34 | /// MLEMV3 : MLEM with preclustering=PRECLUSTERV3 | |
35 | /// | |
36 | /// PRECLUSTER : use only AliMUONPreClusterFinder. Only for debug as | |
37 | /// the produced clusters do not have a position, hence the tracking will not | |
38 | /// work | |
39 | /// PRECLUSTERV2 : another version of the preclustering | |
40 | /// PRECLUSTERV3 : yet another version of the preclustering | |
41 | /// | |
42 | /// COG : use AliMUONClusterFinderCOG clusterizer. Not really a production | |
43 | /// option either, as center-of-gravity is generally not a good estimate | |
44 | /// of the cluster position... | |
45 | /// | |
46 | /// PEAKCOG : COG cluster finder around local maxima | |
47 | /// PEAKFIT : fit around local maxima with up to 3 peaks, COG otherwise | |
48 | /// | |
49 | /// NOCLUSTERING : bypass completely the clustering stage | |
50 | /// | |
51 | /// ------ | |
52 | /// | |
53 | /// The behavior of the MUON reconstruction can also be changed, besides | |
54 | /// the usual methods found in AliReconstruction (e.g. to disable tracking) | |
55 | /// by using AliReconstruction::SetOption("MUON",options) | |
56 | /// where options should be a space separated string. | |
57 | /// | |
58 | /// Valid options are : | |
59 | /// | |
60 | /// SAVEDIGITS : if you want to save in the TreeD the *calibrated* digits | |
61 | /// that are used for the clustering | |
62 | /// | |
63 | /// DIGITSTOREV1 : use the V1 implementation of the digitstore | |
64 | /// DIGITSTOREV2R : use the V2R implementation of the digitstore | |
65 | /// | |
66 | /// NOLOCALRECONSTRUCTION : for debug, to disable local reconstruction (and hence | |
67 | /// "recover" old behavior) | |
68 | /// | |
69 | /// TRIGGERDISABLE : disable the treatment of MUON trigger | |
70 | /// | |
71 | /// NOFASTTRKDECODER : makes the digit maker class use the non-high performance decoder | |
72 | /// AliMUONPayloadTracker instead of AliMUONTrackerDDLDecoder. | |
73 | /// | |
74 | /// NOFASTTRGDECODER : makes the digit maker class use the non-high performance decoder | |
75 | /// AliMUONPayloadTrigger instead of AliMUONTriggerDDLDecoder. | |
76 | /// | |
77 | /// NOFASTDECODERS : makes the digit maker class use the non-high performance decoders | |
78 | /// AliMUONPayloadTracker and AliMUONPayloadTrigger. | |
79 | /// | |
80 | /// \author Laurent Aphecetche, Subatech | |
81 | //----------------------------------------------------------------------------- | |
82 | ||
83 | #include "AliMUONReconstructor.h" | |
84 | ||
85 | #include "AliMUONCalibrationData.h" | |
86 | #include "AliMUONClusterFinderCOG.h" | |
87 | #include "AliMUONClusterFinderMLEM.h" | |
88 | #include "AliMUONClusterFinderSimpleFit.h" | |
89 | #include "AliMUONClusterFinderPeakCOG.h" | |
90 | #include "AliMUONClusterFinderPeakFit.h" | |
91 | #include "AliMUONClusterStoreV1.h" | |
92 | #include "AliMUONClusterStoreV2.h" | |
93 | #include "AliMUONConstants.h" | |
94 | #include "AliMUONDigitCalibrator.h" | |
95 | #include "AliMUONDigitMaker.h" | |
96 | #include "AliMUONDigitStoreV1.h" | |
97 | #include "AliMUONDigitStoreV2R.h" | |
98 | #include "AliMUONGeometryTransformer.h" | |
99 | #include "AliMUONPreClusterFinder.h" | |
100 | #include "AliMUONPreClusterFinderV2.h" | |
101 | #include "AliMUONPreClusterFinderV3.h" | |
102 | #include "AliMUONSimpleClusterServer.h" | |
103 | #include "AliMUONTracker.h" | |
104 | #include "AliMUONTriggerCircuit.h" | |
105 | #include "AliMUONTriggerStoreV1.h" | |
106 | #include "AliMUONVClusterFinder.h" | |
107 | #include "AliMUONVClusterServer.h" | |
108 | #include "AliMUONVTrackStore.h" | |
109 | ||
110 | #include "AliMpArea.h" | |
111 | #include "AliMpCDB.h" | |
112 | #include "AliMpConstants.h" | |
113 | #include "AliMpDDLStore.h" | |
114 | #include "AliMpSegmentation.h" | |
115 | ||
116 | #include "AliRawReader.h" | |
117 | #include "AliCDBManager.h" | |
118 | #include "AliCodeTimer.h" | |
119 | #include "AliLog.h" | |
120 | ||
121 | #include <Riostream.h> | |
122 | #include <TObjArray.h> | |
123 | #include <TClonesArray.h> | |
124 | #include <TString.h> | |
125 | #include <TTree.h> | |
126 | ||
127 | /// \cond CLASSIMP | |
128 | ClassImp(AliMUONReconstructor) | |
129 | /// \endcond | |
130 | ||
131 | //_____________________________________________________________________________ | |
132 | AliMUONReconstructor::AliMUONReconstructor() : | |
133 | AliReconstructor(), | |
134 | fDigitMaker(0x0), | |
135 | fTransformer(new AliMUONGeometryTransformer()), | |
136 | fDigitStore(0x0), | |
137 | fTriggerCircuit(0x0), | |
138 | fCalibrationData(0x0), | |
139 | fDigitCalibrator(0x0), | |
140 | fClusterServer(0x0), | |
141 | fTriggerStore(0x0), | |
142 | fTrackStore(0x0), | |
143 | fClusterStore(0x0) | |
144 | { | |
145 | /// normal ctor | |
146 | ||
147 | AliDebug(1,""); | |
148 | ||
149 | // Unload and delete old mapping | |
150 | AliCDBManager::Instance()->UnloadFromCache("MUON/Calib/Mapping"); | |
151 | AliCDBManager::Instance()->UnloadFromCache("MUON/Calib/DDLStore"); | |
152 | delete AliMpDDLStore::Instance(); | |
153 | delete AliMpSegmentation::Instance(); | |
154 | ||
155 | // Load mapping | |
156 | if ( ! AliMpCDB::LoadDDLStore() ) { | |
157 | AliFatal("Could not access mapping from OCDB !"); | |
158 | } | |
159 | ||
160 | // Load geometry data | |
161 | fTransformer->LoadGeometryData(); | |
162 | ||
163 | } | |
164 | ||
165 | //_____________________________________________________________________________ | |
166 | AliMUONReconstructor::~AliMUONReconstructor() | |
167 | { | |
168 | AliDebug(1,""); | |
169 | ||
170 | /// dtor | |
171 | delete fDigitMaker; | |
172 | delete fDigitStore; | |
173 | delete fTransformer; | |
174 | delete fTriggerCircuit; | |
175 | delete fCalibrationData; | |
176 | delete fDigitCalibrator; | |
177 | delete fClusterServer; | |
178 | delete fTriggerStore; | |
179 | delete fTrackStore; | |
180 | delete fClusterStore; | |
181 | ||
182 | delete AliMpSegmentation::Instance(false); | |
183 | delete AliMpDDLStore::Instance(false); | |
184 | } | |
185 | ||
186 | //_____________________________________________________________________________ | |
187 | void | |
188 | AliMUONReconstructor::Calibrate(AliMUONVDigitStore& digitStore) const | |
189 | { | |
190 | /// Calibrate the digitStore | |
191 | if (!fDigitCalibrator) | |
192 | { | |
193 | CreateCalibrator(); | |
194 | } | |
195 | AliCodeTimerAuto(Form("%s::Calibrate(AliMUONVDigitStore*)",fDigitCalibrator->ClassName())) | |
196 | fDigitCalibrator->Calibrate(digitStore); | |
197 | } | |
198 | ||
199 | //_____________________________________________________________________________ | |
200 | void | |
201 | AliMUONReconstructor::ConvertDigits(AliRawReader* rawReader, | |
202 | AliMUONVDigitStore* digitStore, | |
203 | AliMUONVTriggerStore* triggerStore) const | |
204 | { | |
205 | /// Convert raw data into digit and trigger stores | |
206 | CreateDigitMaker(); | |
207 | ||
208 | AliCodeTimerStart(Form("%s::Raw2Digits(AliRawReader*,AliMUONVDigitStore*,AliMUONVTriggerStore*)", | |
209 | fDigitMaker->ClassName())) | |
210 | fDigitMaker->Raw2Digits(rawReader,digitStore,triggerStore); | |
211 | AliCodeTimerStop(Form("%s::Raw2Digits(AliRawReader*,AliMUONVDigitStore*,AliMUONVTriggerStore*)", | |
212 | fDigitMaker->ClassName())) | |
213 | Calibrate(*digitStore); | |
214 | } | |
215 | ||
216 | //_____________________________________________________________________________ | |
217 | void | |
218 | AliMUONReconstructor::ConvertDigits(AliRawReader* rawReader, TTree* digitsTree) const | |
219 | { | |
220 | /// convert raw data into a digit tree | |
221 | AliCodeTimerAuto("") | |
222 | ||
223 | Bool_t alone = ( TriggerStore() == 0 ); | |
224 | ||
225 | Bool_t ok = DigitStore()->Connect(*digitsTree,alone); | |
226 | if ( TriggerStore() ) | |
227 | { | |
228 | ok = ok && TriggerStore()->Connect(*digitsTree,kFALSE); | |
229 | } | |
230 | ||
231 | if (!ok) | |
232 | { | |
233 | AliError("Could not make branches on TreeD"); | |
234 | } | |
235 | else | |
236 | { | |
237 | ConvertDigits(rawReader,DigitStore(),TriggerStore()); | |
238 | AliCodeTimerStart("Fill digits") | |
239 | digitsTree->Fill(); | |
240 | AliCodeTimerStop("Fill digits") | |
241 | DigitStore()->Clear(); | |
242 | } | |
243 | } | |
244 | ||
245 | //_____________________________________________________________________________ | |
246 | void | |
247 | AliMUONReconstructor::CreateDigitMaker() const | |
248 | { | |
249 | /// Create (and create if necessary) the digit maker | |
250 | if (fDigitMaker) return; | |
251 | ||
252 | AliCodeTimerAuto("") | |
253 | ||
254 | TString option = GetOption(); | |
255 | Bool_t enableErrorLogging = kTRUE; | |
256 | Bool_t useFastTrackerDecoder = kTRUE; | |
257 | Bool_t useFastTriggerDecoder = kTRUE; | |
258 | if (option.Contains("NOFASTTRKDECODER")) | |
259 | { | |
260 | useFastTrackerDecoder = kFALSE; | |
261 | } | |
262 | if (option.Contains("NOFASTTRGDECODER")) | |
263 | { | |
264 | useFastTriggerDecoder = kFALSE; | |
265 | } | |
266 | if (option.Contains("NOFASTDECODERS")) | |
267 | { | |
268 | useFastTrackerDecoder = kFALSE; | |
269 | useFastTriggerDecoder = kFALSE; | |
270 | } | |
271 | fDigitMaker = new AliMUONDigitMaker( | |
272 | enableErrorLogging, useFastTrackerDecoder, useFastTriggerDecoder | |
273 | ); | |
274 | option.ToUpper(); | |
275 | if ( option.Contains("SAVEDIGITS" )) | |
276 | { | |
277 | fDigitMaker->SetMakeTriggerDigits(kTRUE); | |
278 | } | |
279 | } | |
280 | ||
281 | //_____________________________________________________________________________ | |
282 | void | |
283 | AliMUONReconstructor::CreateTriggerCircuit() const | |
284 | { | |
285 | /// Return (and create if necessary) the trigger circuit object | |
286 | if (fTriggerCircuit) return; | |
287 | ||
288 | AliCodeTimerAuto("") | |
289 | ||
290 | fTriggerCircuit = new AliMUONTriggerCircuit(fTransformer); | |
291 | ||
292 | } | |
293 | ||
294 | //_____________________________________________________________________________ | |
295 | AliTracker* | |
296 | AliMUONReconstructor::CreateTracker() const | |
297 | { | |
298 | /// Create the MUONTracker object | |
299 | ||
300 | CreateTriggerCircuit(); | |
301 | CreateDigitMaker(); | |
302 | CreateClusterServer(); | |
303 | ||
304 | AliMUONTracker* tracker(0x0); | |
305 | ||
306 | if ( ! GetRecoParam()->CombineClusterTrackReco() ) | |
307 | { | |
308 | tracker = new AliMUONTracker(GetRecoParam(), | |
309 | 0x0, | |
310 | *DigitStore(), | |
311 | fDigitMaker, | |
312 | fTransformer, | |
313 | fTriggerCircuit); | |
314 | } | |
315 | else | |
316 | { | |
317 | tracker = new AliMUONTracker(GetRecoParam(), | |
318 | fClusterServer, | |
319 | *DigitStore(), | |
320 | fDigitMaker, | |
321 | fTransformer, | |
322 | fTriggerCircuit); | |
323 | } | |
324 | ||
325 | ||
326 | return tracker; | |
327 | } | |
328 | ||
329 | //_____________________________________________________________________________ | |
330 | AliMUONVClusterFinder* | |
331 | AliMUONReconstructor::CreateClusterFinder(const char* clusterFinderType) | |
332 | { | |
333 | /// Create a given cluster finder instance | |
334 | ||
335 | AliCodeTimerAutoGeneral("") | |
336 | ||
337 | AliMUONVClusterFinder* clusterFinder(0x0); | |
338 | ||
339 | TString opt(clusterFinderType); | |
340 | opt.ToUpper(); | |
341 | ||
342 | if ( strstr(opt,"PRECLUSTERV2") ) | |
343 | { | |
344 | clusterFinder = new AliMUONPreClusterFinderV2; | |
345 | } | |
346 | else if ( strstr(opt,"PRECLUSTERV3") ) | |
347 | { | |
348 | clusterFinder = new AliMUONPreClusterFinderV3; | |
349 | } | |
350 | else if ( strstr(opt,"PRECLUSTER") ) | |
351 | { | |
352 | clusterFinder = new AliMUONPreClusterFinder; | |
353 | } | |
354 | else if ( strstr(opt,"PEAKCOG") ) | |
355 | { | |
356 | clusterFinder = new AliMUONClusterFinderPeakCOG(kFALSE,new AliMUONPreClusterFinder); | |
357 | } | |
358 | else if ( strstr(opt,"PEAKFIT") ) | |
359 | { | |
360 | clusterFinder = new AliMUONClusterFinderPeakFit(kFALSE,new AliMUONPreClusterFinder); | |
361 | } | |
362 | else if ( strstr(opt,"COG") ) | |
363 | { | |
364 | clusterFinder = new AliMUONClusterFinderCOG(new AliMUONPreClusterFinder); | |
365 | } | |
366 | else if ( strstr(opt,"SIMPLEFITV3") ) | |
367 | { | |
368 | clusterFinder = new AliMUONClusterFinderSimpleFit(new AliMUONClusterFinderCOG(new AliMUONPreClusterFinderV3)); | |
369 | } | |
370 | else if ( strstr(opt,"SIMPLEFIT") ) | |
371 | { | |
372 | clusterFinder = new AliMUONClusterFinderSimpleFit(new AliMUONClusterFinderCOG(new AliMUONPreClusterFinder)); | |
373 | } | |
374 | else if ( strstr(opt,"MLEM:DRAW") ) | |
375 | { | |
376 | clusterFinder = new AliMUONClusterFinderMLEM(kTRUE,new AliMUONPreClusterFinder); | |
377 | } | |
378 | else if ( strstr(opt,"MLEMV3") ) | |
379 | { | |
380 | clusterFinder = new AliMUONClusterFinderMLEM(kFALSE,new AliMUONPreClusterFinderV3); | |
381 | } | |
382 | else if ( strstr(opt,"MLEMV2") ) | |
383 | { | |
384 | clusterFinder = new AliMUONClusterFinderMLEM(kFALSE,new AliMUONPreClusterFinderV2); | |
385 | } | |
386 | else if ( strstr(opt,"MLEM") ) | |
387 | { | |
388 | clusterFinder = new AliMUONClusterFinderMLEM(kFALSE,new AliMUONPreClusterFinder); | |
389 | } | |
390 | else | |
391 | { | |
392 | AliErrorClass(Form("clustering mode \"%s\" does not exist",opt.Data())); | |
393 | return 0x0; | |
394 | } | |
395 | ||
396 | return clusterFinder; | |
397 | } | |
398 | ||
399 | //_____________________________________________________________________________ | |
400 | void | |
401 | AliMUONReconstructor::CreateClusterServer() const | |
402 | { | |
403 | /// Create cluster server | |
404 | ||
405 | if ( fClusterServer ) return; | |
406 | ||
407 | AliCodeTimerAuto(""); | |
408 | ||
409 | AliMUONVClusterFinder* clusterFinder = CreateClusterFinder(GetRecoParam()->GetClusteringMode()); | |
410 | ||
411 | if ( !clusterFinder ) return; | |
412 | ||
413 | AliInfo(Form("Will use %s for clusterizing",clusterFinder->ClassName())); | |
414 | ||
415 | fClusterServer = new AliMUONSimpleClusterServer(clusterFinder,*fTransformer); | |
416 | } | |
417 | ||
418 | //_____________________________________________________________________________ | |
419 | void | |
420 | AliMUONReconstructor::CreateCalibrator() const | |
421 | { | |
422 | /// Create the calibrator | |
423 | ||
424 | AliCodeTimerAuto("") | |
425 | ||
426 | Int_t runNumber = AliCDBManager::Instance()->GetRun(); | |
427 | ||
428 | AliInfo("Calibration will occur."); | |
429 | ||
430 | fCalibrationData = new AliMUONCalibrationData(runNumber); | |
431 | if ( !fCalibrationData->IsValid() ) | |
432 | { | |
433 | AliError("Could not retrieve calibrations !"); | |
434 | delete fCalibrationData; | |
435 | fCalibrationData = 0x0; | |
436 | return; | |
437 | } | |
438 | ||
439 | // Check that we get all the calibrations we'll need | |
440 | if ( !fCalibrationData->Pedestals() || | |
441 | !fCalibrationData->Gains() || | |
442 | !fCalibrationData->HV() ) | |
443 | { | |
444 | AliFatal("Could not access all required calibration data"); | |
445 | } | |
446 | ||
447 | TString opt(GetOption()); | |
448 | opt.ToUpper(); | |
449 | ||
450 | if ( strstr(opt,"NOSTATUSMAP") ) | |
451 | { | |
452 | AliWarning("NOSTATUSMAP is obsolete"); | |
453 | } | |
454 | ||
455 | TString calibMode = GetRecoParam()->GetCalibrationMode(); | |
456 | ||
457 | fDigitCalibrator = new AliMUONDigitCalibrator(*fCalibrationData,GetRecoParam(),calibMode.Data()); | |
458 | } | |
459 | ||
460 | //_____________________________________________________________________________ | |
461 | AliMUONVDigitStore* | |
462 | AliMUONReconstructor::DigitStore() const | |
463 | { | |
464 | /// Return (and create if necessary) the digit container | |
465 | if (!fDigitStore) | |
466 | { | |
467 | TString sopt(GetOption()); | |
468 | sopt.ToUpper(); | |
469 | ||
470 | AliInfo(Form("Options=%s",sopt.Data())); | |
471 | ||
472 | if ( sopt.Contains("DIGITSTOREV1") ) | |
473 | { | |
474 | fDigitStore = AliMUONVDigitStore::Create("AliMUONDigitStoreV1"); | |
475 | } | |
476 | else if ( sopt.Contains("DIGITSTOREV2R") ) | |
477 | { | |
478 | fDigitStore = AliMUONVDigitStore::Create("AliMUONDigitStoreV2R"); | |
479 | } | |
480 | else if ( sopt.Contains("DIGITSTOREV2S") ) | |
481 | { | |
482 | fDigitStore = AliMUONVDigitStore::Create("AliMUONDigitStoreV2S"); | |
483 | } | |
484 | ||
485 | if (!fDigitStore) fDigitStore = AliMUONVDigitStore::Create("AliMUONDigitStoreV2R"); | |
486 | ||
487 | AliInfo(Form("Will use %s to store digits during reconstruction",fDigitStore->ClassName())); | |
488 | } | |
489 | return fDigitStore; | |
490 | } | |
491 | ||
492 | //_____________________________________________________________________________ | |
493 | void | |
494 | AliMUONReconstructor::FillTreeR(AliMUONVTriggerStore* triggerStore, | |
495 | TTree& clustersTree) const | |
496 | { | |
497 | /// Write the trigger and cluster information into TreeR | |
498 | ||
499 | AliCodeTimerAuto("") | |
500 | ||
501 | AliDebug(1,""); | |
502 | ||
503 | Bool_t ok(kFALSE); | |
504 | Bool_t alone(kTRUE); // is trigger the only info in TreeR ? | |
505 | ||
506 | if ( ! GetRecoParam()->CombineClusterTrackReco() ) | |
507 | { | |
508 | alone = kFALSE; // we'll get both tracker and trigger information in TreeR | |
509 | } | |
510 | ||
511 | if ( triggerStore ) | |
512 | { | |
513 | ok = triggerStore->Connect(clustersTree,alone); | |
514 | if (!ok) | |
515 | { | |
516 | AliError("Could not create triggerStore branches in TreeR"); | |
517 | } | |
518 | } | |
519 | ||
520 | if ( !alone ) | |
521 | { | |
522 | if (!fClusterStore) | |
523 | { | |
524 | fClusterStore = new AliMUONClusterStoreV2; | |
525 | } | |
526 | ||
527 | CreateClusterServer(); | |
528 | ||
529 | TIter next(DigitStore()->CreateIterator()); | |
530 | fClusterServer->UseDigits(next); | |
531 | ||
532 | AliMpArea area; | |
533 | ||
534 | AliDebug(1,Form("Doing full clusterization in local reconstruction using %s ",fClusterServer->ClassName())); | |
535 | ||
536 | for ( Int_t i = 0; i < AliMpConstants::NofTrackingChambers(); ++i ) | |
537 | { | |
538 | if (GetRecoParam()->UseChamber(i)) | |
539 | { | |
540 | if ( ( i == 6 || i == 7 ) && GetRecoParam()->BypassSt4() ) continue; | |
541 | if ( ( i == 8 || i == 9 ) && GetRecoParam()->BypassSt5() ) continue; | |
542 | ||
543 | fClusterServer->Clusterize(i,*fClusterStore,area,GetRecoParam()); | |
544 | } | |
545 | } | |
546 | ||
547 | Bool_t cok = fClusterStore->Connect(clustersTree,alone); | |
548 | ||
549 | if (!cok) AliError("Could not connect clusterStore to clusterTree"); | |
550 | ||
551 | AliDebug(1,Form("Number of clusters found = %d",fClusterStore->GetSize())); | |
552 | ||
553 | StdoutToAliDebug(1,fClusterStore->Print()); | |
554 | } | |
555 | ||
556 | if (ok) // at least one type of branches created successfully | |
557 | { | |
558 | clustersTree.Fill(); | |
559 | } | |
560 | ||
561 | if (fClusterStore) fClusterStore->Clear(); | |
562 | } | |
563 | ||
564 | //_____________________________________________________________________________ | |
565 | Bool_t | |
566 | AliMUONReconstructor::HasDigitConversion() const | |
567 | { | |
568 | /// We *do* have digit conversion, but we might advertise it only | |
569 | /// if we want to save the digits. | |
570 | ||
571 | TString opt(GetOption()); | |
572 | opt.ToUpper(); | |
573 | if ( opt.Contains("SAVEDIGITS" ) && !opt.Contains("NOLOCALRECONSTRUCTION") ) | |
574 | { | |
575 | return kTRUE; | |
576 | } | |
577 | else | |
578 | { | |
579 | return kFALSE; | |
580 | } | |
581 | } | |
582 | ||
583 | //_____________________________________________________________________________ | |
584 | void | |
585 | AliMUONReconstructor::Reconstruct(AliRawReader* rawReader, TTree* clustersTree) const | |
586 | { | |
587 | /// This method is called by AliReconstruction if HasLocalReconstruction()==kTRUE AND | |
588 | /// HasDigitConversion()==kFALSE | |
589 | ||
590 | if ( !clustersTree ) | |
591 | { | |
592 | AliError("clustersTree is 0x0 !"); | |
593 | return; | |
594 | } | |
595 | ||
596 | ConvertDigits(rawReader,DigitStore(),TriggerStore()); | |
597 | ||
598 | FillTreeR(TriggerStore(),*clustersTree); | |
599 | } | |
600 | ||
601 | //_____________________________________________________________________________ | |
602 | void | |
603 | AliMUONReconstructor::Reconstruct(TTree* digitsTree, TTree* clustersTree) const | |
604 | { | |
605 | /// This method is called by AliReconstruction if HasLocalReconstruction()==kTRUE | |
606 | /// AND HasDigitConversion()==kTRUE | |
607 | ||
608 | AliCodeTimerAuto("") | |
609 | ||
610 | AliDebug(1,""); | |
611 | ||
612 | if (!digitsTree || !clustersTree) | |
613 | { | |
614 | AliError(Form("Tree is null : digitsTree=%p clustersTree=%p", | |
615 | digitsTree,clustersTree)); | |
616 | return; | |
617 | } | |
618 | ||
619 | if (!fDigitStore) | |
620 | { | |
621 | fDigitStore = AliMUONVDigitStore::Create(*digitsTree); | |
622 | if (!fDigitStore) | |
623 | { | |
624 | AliError(Form("Could not get DigitStore from %s",digitsTree->GetName())); | |
625 | } | |
626 | else | |
627 | { | |
628 | AliInfo(Form("Created %s from %s",fDigitStore->ClassName(),digitsTree->GetName())); | |
629 | } | |
630 | } | |
631 | if (!fTriggerStore) | |
632 | { | |
633 | fTriggerStore = AliMUONVTriggerStore::Create(*digitsTree); | |
634 | if (!fTriggerStore) | |
635 | { | |
636 | AliError(Form("Could not get TriggerStore from %s",digitsTree->GetName())); | |
637 | } | |
638 | else | |
639 | { | |
640 | AliInfo(Form("Created %s from %s",fTriggerStore->ClassName(),digitsTree->GetName())); | |
641 | } | |
642 | } | |
643 | ||
644 | if (!fTriggerStore && !fDigitStore) | |
645 | { | |
646 | AliError("No store at all. Nothing to do."); | |
647 | return; | |
648 | } | |
649 | ||
650 | // insure we start with empty stores | |
651 | if ( fDigitStore ) | |
652 | { | |
653 | fDigitStore->Clear(); | |
654 | Bool_t alone = ( fTriggerStore ? kFALSE : kTRUE ); | |
655 | Bool_t ok = fDigitStore->Connect(*digitsTree,alone); | |
656 | if (!ok) | |
657 | { | |
658 | AliError("Could not connect digitStore to digitsTree"); | |
659 | return; | |
660 | } | |
661 | } | |
662 | if ( fTriggerStore ) | |
663 | { | |
664 | fTriggerStore->Clear(); | |
665 | Bool_t alone = ( fDigitStore ? kFALSE : kTRUE ); | |
666 | Bool_t ok = fTriggerStore->Connect(*digitsTree,alone); | |
667 | if (!ok) | |
668 | { | |
669 | AliError("Could not connect triggerStore to digitsTree"); | |
670 | return; | |
671 | } | |
672 | } | |
673 | ||
674 | digitsTree->GetEvent(0); | |
675 | ||
676 | if ( fDigitStore ) | |
677 | { | |
678 | // Insure we got calibrated digits (if we reconstruct from pure simulated, | |
679 | // i.e. w/o going through raw data, this will be the case) | |
680 | TIter next(fDigitStore->CreateIterator()); | |
681 | AliMUONVDigit* digit = static_cast<AliMUONVDigit*>(next()); | |
682 | if (digit && !digit->IsCalibrated()) | |
683 | { | |
684 | Calibrate(*fDigitStore); | |
685 | } | |
686 | } | |
687 | ||
688 | FillTreeR(fTriggerStore,*clustersTree); | |
689 | } | |
690 | ||
691 | //_____________________________________________________________________________ | |
692 | AliMUONVTriggerStore* | |
693 | AliMUONReconstructor::TriggerStore() const | |
694 | { | |
695 | /// Return (and create if necessary and allowed) the trigger container | |
696 | TString sopt(GetOption()); | |
697 | sopt.ToUpper(); | |
698 | ||
699 | if (sopt.Contains("TRIGGERDISABLE")) | |
700 | { | |
701 | delete fTriggerStore; | |
702 | fTriggerStore = 0x0; | |
703 | } | |
704 | else | |
705 | { | |
706 | if (!fTriggerStore) | |
707 | { | |
708 | fTriggerStore = new AliMUONTriggerStoreV1; | |
709 | } | |
710 | } | |
711 | return fTriggerStore; | |
712 | } |