expantion of recoParam (Markus) and fix bug in HLT mode
[u/mrichter/AliRoot.git] / TRD / AliTRDReconstructor.h
1 #ifndef ALITRDRECONSTRUCTOR_H
2 #define ALITRDRECONSTRUCTOR_H
3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4  * See cxx source for full Copyright notice                               */
5
6 /* $Id$ */
7
8 ///////////////////////////////////////////////////////////////////////////////
9 //                                                                           //
10 // Class for TRD reconstruction                                              //
11 //                                                                           //
12 ///////////////////////////////////////////////////////////////////////////////
13
14 #include "AliReconstructor.h"
15 #include "AliDetectorRecoParam.h"
16 #include "AliTRDrecoParam.h"
17
18 class TClonesArray;
19 class AliRawReader;
20 class AliTRDReconstructor: public AliReconstructor 
21 {
22 public:
23   enum AliTRDsteerParam {
24     kWriteClusters   = BIT(0)
25     ,kSeeding        = BIT(1)
26     ,kSteerPID       = BIT(2)
27     ,kWriteTracklets = BIT(3)
28     ,kDriftGas       = BIT(4)
29     ,kHLT            = BIT(5)
30     ,kCosmic         = BIT(6)
31   };
32   enum AliTRDReconstructorTask {
33     kClusterizer = 0
34     ,kTracker    = 1
35     ,kPID        = 2
36   };
37   enum AliTRDpidMethod {
38     kLQPID = 0,
39     kNNPID = 1
40   };
41   enum AliTRDdriftGas {
42     kXe = 0,
43     kAr = 1
44   };
45   enum{
46     kNNslices = 8
47    ,kLQslices = 3
48   };
49
50   AliTRDReconstructor();
51   AliTRDReconstructor(const AliTRDReconstructor &r);
52   virtual ~AliTRDReconstructor();
53   AliTRDReconstructor& operator = (const AliTRDReconstructor&)          { return *this;}
54         
55         virtual void        Init();
56
57   virtual void        ConvertDigits(AliRawReader *rawReader, TTree *digitsTree) const;
58   virtual AliTracker* CreateTracker() const;
59
60   virtual void        FillESD(AliRawReader *, TTree *clusterTree, AliESDEvent *esd) const { FillESD((TTree * )NULL, clusterTree, esd);                    }
61   virtual void        FillESD(TTree *digitsTree, TTree *clusterTree, AliESDEvent *esd) const;
62   static TClonesArray* GetClusters() {return fgClusters;}
63   Int_t               GetNdEdxSlices() const     { return GetPIDMethod() == kNNPID ? kNNslices : kLQslices;}
64   AliTRDdriftGas      GetDriftGas() const        { return fSteerParam&kDriftGas ? kAr : kXe;}
65   AliTRDpidMethod     GetPIDMethod() const       { return fSteerParam&kSteerPID ? kNNPID : kLQPID;}
66   static const AliTRDrecoParam* GetRecoParam() { return dynamic_cast<const AliTRDrecoParam*>(AliReconstructor::GetRecoParam(2)); }
67   Int_t               GetStreamLevel(AliTRDReconstructorTask task) const    { return fStreamLevel[task];} 
68   inline void         GetTCParams(Double_t *par) const;
69   virtual Bool_t      HasDigitConversion() const                   { return kFALSE;           };
70   Bool_t              IsWritingClusters() const  { return fSteerParam&kWriteClusters;}
71   Bool_t              IsWritingTracklets() const { return fSteerParam&kWriteTracklets;}
72   Bool_t              IsHLT() const              { return fSteerParam&kHLT;}
73   Bool_t              IsSeeding() const          { return fSteerParam&kSeeding;}
74   Bool_t              IsCosmic() const           { return fSteerParam&kCosmic;}
75
76   virtual void        Reconstruct(AliRawReader *rawReader, TTree *clusterTree) const;
77   virtual void        Reconstruct(TTree *digitsTree, TTree *clusterTree) const;
78
79   static void         SetClusters(TClonesArray *clusters) {fgClusters = clusters;}
80   void                SetOption(Option_t *opt);
81   inline void         SetTCParams(Double_t *par);
82   void                SetStreamLevel(Int_t level, AliTRDReconstructorTask task= kTracker);
83
84 private:
85   UChar_t       fStreamLevel[5];      // stream level for each reconstruction task         
86   UInt_t        fSteerParam;          // steering flags
87   Double_t      fTCParams[8];         // Tail Cancellation parameters for drift gases 
88  
89   static TClonesArray *fgClusters;    // list of clusters for local reconstructor
90
91   ClassDef(AliTRDReconstructor, 1)    //  Class for the TRD reconstruction
92
93 };
94
95 //___________________________________________________
96 inline void AliTRDReconstructor::GetTCParams(Double_t *par) const
97 {
98   if(!par) return;
99   if(GetDriftGas()==kAr) memcpy(par, &fTCParams[4], 4*sizeof(Double_t));
100   else memcpy(par, &fTCParams[0], 4*sizeof(Double_t));
101 }
102
103 //___________________________________________________
104 inline void AliTRDReconstructor::SetTCParams(Double_t *par)
105 {
106   if(!par) return;
107   memcpy(fTCParams, par, 8*sizeof(Double_t));
108 }
109
110 #endif