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