- new option for selecting 8 dEdx slices in ESD track for TRD
[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     kDigitsConversion= BIT(0)
25     ,kWriteClusters  = BIT(1)
26     ,kSeeding        = BIT(2)
27     ,kSteerPID       = BIT(3)
28     ,kEightSlices    = BIT(4)
29     ,kWriteTracklets = BIT(5)
30     ,kDriftGas       = BIT(6)
31     ,kHLT            = BIT(7)
32     ,kCosmic         = BIT(8)
33   };
34   enum AliTRDReconstructorTask {
35     kClusterizer = 0
36     ,kTracker    = 1
37     ,kPID        = 2
38   };
39   enum AliTRDpidMethod {
40     kLQPID = 0,
41     kNNPID = 1
42   };
43   enum AliTRDdriftGas {
44     kXe = 0,
45     kAr = 1
46   };
47   enum{
48     kNNslices = 8
49    ,kLQslices = 3
50   };
51
52   AliTRDReconstructor();
53   AliTRDReconstructor(const AliTRDReconstructor &r);
54   virtual ~AliTRDReconstructor();
55   AliTRDReconstructor& operator = (const AliTRDReconstructor&)          { return *this;}
56         
57         virtual void        Init();
58
59   virtual void        ConvertDigits(AliRawReader *rawReader, TTree *digitsTree) const;
60   virtual AliTracker* CreateTracker() const;
61
62   virtual void        FillESD(AliRawReader *, TTree *clusterTree, AliESDEvent *esd) const { FillESD((TTree * )NULL, clusterTree, esd);                    }
63   virtual void        FillESD(TTree *digitsTree, TTree *clusterTree, AliESDEvent *esd) const;
64   static TClonesArray* GetClusters() {return fgClusters;}
65   Int_t               GetNdEdxSlices() const     { return GetPIDMethod() == kNNPID ? kNNslices : kLQslices;}
66   AliTRDdriftGas      GetDriftGas() const        { return fSteerParam&kDriftGas ? kAr : kXe;}
67   AliTRDpidMethod     GetPIDMethod() const       { return fSteerParam&kSteerPID ? kNNPID : kLQPID;}
68   static const AliTRDrecoParam* GetRecoParam() { return dynamic_cast<const AliTRDrecoParam*>(AliReconstructor::GetRecoParam(2)); }
69   Int_t               GetStreamLevel(AliTRDReconstructorTask task) const    { return fStreamLevel[task];} 
70   inline void         GetTCParams(Double_t *par) const;
71   virtual Bool_t      HasDigitConversion() const                   { return fSteerParam&kDigitsConversion;  };
72   Bool_t              IsWritingClusters() const  { return fSteerParam&kWriteClusters;}
73   Bool_t              IsWritingTracklets() const { return fSteerParam&kWriteTracklets;}
74   Bool_t              IsHLT() const              { return fSteerParam&kHLT;}
75   Bool_t              IsSeeding() const          { return fSteerParam&kSeeding;}
76   Bool_t              IsCosmic() const           { return fSteerParam&kCosmic;}
77   Bool_t              IsEightSlices() const      { return fSteerParam&kEightSlices;}
78
79   virtual void        Reconstruct(AliRawReader *rawReader, TTree *clusterTree) const;
80   virtual void        Reconstruct(TTree *digitsTree, TTree *clusterTree) const;
81
82   static void         SetClusters(TClonesArray *clusters) {fgClusters = clusters;}
83   void                SetOption(Option_t *opt);
84   inline void         SetTCParams(Double_t *par);
85   void                SetStreamLevel(Int_t level, AliTRDReconstructorTask task= kTracker);
86
87 private:
88   UChar_t       fStreamLevel[5];      // stream level for each reconstruction task         
89   UInt_t        fSteerParam;          // steering flags
90   Double_t      fTCParams[8];         // Tail Cancellation parameters for drift gases 
91  
92   static TClonesArray *fgClusters;    // list of clusters for local reconstructor
93
94   ClassDef(AliTRDReconstructor, 1)    //  Class for the TRD reconstruction
95
96 };
97
98 //___________________________________________________
99 inline void AliTRDReconstructor::GetTCParams(Double_t *par) const
100 {
101   if(!par) return;
102   if(GetDriftGas()==kAr) memcpy(par, &fTCParams[4], 4*sizeof(Double_t));
103   else memcpy(par, &fTCParams[0], 4*sizeof(Double_t));
104 }
105
106 //___________________________________________________
107 inline void AliTRDReconstructor::SetTCParams(Double_t *par)
108 {
109   if(!par) return;
110   memcpy(fTCParams, par, 8*sizeof(Double_t));
111 }
112
113 #endif