Add a protection for the case of not used dictionary and a setter for the recoParam
[u/mrichter/AliRoot.git] / TRD / AliTRDReconstructor.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
16 /* $Id$ */
17
18 ///////////////////////////////////////////////////////////////////////////////
19 //                                                                           //
20 // Class for TRD reconstruction                                              //
21 //                                                                           //
22 ///////////////////////////////////////////////////////////////////////////////
23
24 #include <TFile.h>
25
26 #include "AliRunLoader.h"
27 #include "AliRawReader.h"
28 #include "AliLog.h"
29 #include "AliESDTrdTrack.h"
30 #include "AliESDEvent.h"
31
32 #include "AliTRDReconstructor.h"
33 #include "AliTRDclusterizer.h"
34 #include "AliTRDtracker.h"
35 #include "AliTRDpidESD.h"
36 #include "AliTRDgtuTrack.h"
37 #include "AliTRDrawData.h"
38 #include "AliTRDdigitsManager.h"
39 #include "AliTRDtrackerV1.h"
40 #include "AliTRDrecoParam.h"
41
42 ClassImp(AliTRDReconstructor)
43
44 Bool_t AliTRDReconstructor::fgkSeedingOn  = kFALSE;
45 Int_t  AliTRDReconstructor::fgStreamLevel = 0;      // Stream (debug) level
46 AliTRDrecoParam* AliTRDReconstructor::fgRecoParam = 0x0; 
47
48
49 //_____________________________________________________________________________
50 AliTRDReconstructor::~AliTRDReconstructor()                  { 
51         if(fgRecoParam) delete fgRecoParam;
52 }
53
54
55
56 //_____________________________________________________________________________
57 void AliTRDReconstructor::ConvertDigits(AliRawReader *rawReader
58                                       , TTree *digitsTree) const
59 {
60   //
61   // Convert raw data digits into digit objects in a root tree
62   //
63
64   AliInfo("Convert raw data digits into digit objects [RawReader -> Digit TTree]");
65
66   AliTRDrawData rawData;
67   rawReader->Reset();
68   rawReader->Select("TRD");
69   AliTRDdigitsManager *manager = rawData.Raw2Digits(rawReader);
70   manager->MakeBranch(digitsTree);
71   manager->WriteDigits();
72   delete manager;
73
74 }
75
76 //_____________________________________________________________________________
77 void AliTRDReconstructor::Reconstruct(AliRawReader *rawReader
78                                     , TTree *clusterTree) const
79 {
80   //
81   // Reconstruct clusters
82   //
83
84   AliInfo("Reconstruct TRD clusters from RAW data [RawReader -> Cluster TTree]");
85
86         // TODO move it to rec.C. check TPC
87         fgRecoParam = AliTRDrecoParam::GetLowFluxParam();
88
89   rawReader->Reset();
90   rawReader->Select("TRD");
91
92   // New (fast) cluster finder
93   AliTRDclusterizer clusterer("clusterer","TRD clusterizer");
94   clusterer.OpenOutput(clusterTree);
95   clusterer.SetAddLabels(kFALSE);
96   clusterer.Raw2ClustersChamber(rawReader);
97
98 }
99
100 //_____________________________________________________________________________
101 void AliTRDReconstructor::Reconstruct(TTree *digitsTree
102                                     , TTree *clusterTree) const
103 {
104   //
105   // Reconstruct clusters
106   //
107
108   AliInfo("Reconstruct TRD clusters from Digits [Digit TTree -> Cluster TTree]");
109
110         // TODO move it to rec.C. check TPC
111         fgRecoParam = AliTRDrecoParam::GetLowFluxParam();
112
113   AliTRDclusterizer clusterer("clusterer","TRD clusterizer");
114   clusterer.OpenOutput(clusterTree);
115   clusterer.ReadDigits(digitsTree);
116   clusterer.MakeClusters();
117
118 }
119
120 //_____________________________________________________________________________
121 AliTracker *AliTRDReconstructor::CreateTracker() const
122 {
123   //
124   // Create a TRD tracker
125   //
126
127   //return new AliTRDtracker(NULL);
128
129         // TODO move it to rec.C. check TPC
130         fgRecoParam = AliTRDrecoParam::GetLowFluxParam();
131   return new AliTRDtrackerV1();
132
133 }
134
135 //_____________________________________________________________________________
136 void AliTRDReconstructor::FillESD(TTree* /*digitsTree*/
137                                 , TTree* /*clusterTree*/
138                                 , AliESDEvent* /*esd*/) const
139 {
140   //
141   // Fill ESD
142   //
143
144 }