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