]>
Commit | Line | Data |
---|---|---|
121a60bd | 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 | ||
3bc9d03e | 6 | /* $Id$ */ |
7 | ||
121a60bd | 8 | /////////////////////////////////////////////////////////////////////////////// |
9 | // // | |
3bc9d03e | 10 | // Class for TRD reconstruction // |
121a60bd | 11 | // // |
12 | /////////////////////////////////////////////////////////////////////////////// | |
13 | ||
121a60bd | 14 | #include "AliReconstructor.h" |
7e88424f | 15 | #include "AliDetectorRecoParam.h" |
16 | #include "AliTRDrecoParam.h" | |
0d83b3a5 | 17 | #include "AliTRDpidUtil.h" |
7e88424f | 18 | |
66f6bfd9 | 19 | class TClonesArray; |
29f95561 | 20 | class TTreeSRedirector; |
3a039a31 | 21 | class AliRawReader; |
eb38ed55 | 22 | class AliTRDReconstructor: public AliReconstructor |
23 | { | |
3a039a31 | 24 | public: |
89971837 | 25 | enum ETRDReconstructorSteer { |
00872461 | 26 | kDigitsConversion= BIT(0) |
27 | ,kWriteClusters = BIT(1) | |
28 | ,kSeeding = BIT(2) | |
29 | ,kSteerPID = BIT(3) | |
56b32fd7 | 30 | ,kEightSlices = BIT(4) |
31 | ,kWriteTracklets = BIT(5) | |
32 | ,kDriftGas = BIT(6) | |
33 | ,kHLT = BIT(7) | |
34 | ,kCosmic = BIT(8) | |
89971837 | 35 | ,kOwner = BIT(14) |
3a039a31 | 36 | }; |
89971837 | 37 | enum ETRDReconstructorTask { |
29f95561 | 38 | kRawReader = 0 |
39 | ,kClusterizer = 1 | |
40 | ,kTracker = 2 | |
41 | ,kPID = 3 | |
89971837 | 42 | ,kNtasks = 4 // number of reconsruction tasks |
3a039a31 | 43 | }; |
89971837 | 44 | enum ETRDReconstructorGas { |
45 | kXe = 0 | |
46 | ,kAr = 1 | |
29f95561 | 47 | }; |
3bc9d03e | 48 | |
3a039a31 | 49 | AliTRDReconstructor(); |
50 | AliTRDReconstructor(const AliTRDReconstructor &r); | |
48f8adf3 | 51 | virtual ~AliTRDReconstructor(); |
3a039a31 | 52 | AliTRDReconstructor& operator = (const AliTRDReconstructor&) { return *this;} |
acd241e9 | 53 | |
54 | virtual void Init(); | |
3bc9d03e | 55 | |
3a039a31 | 56 | virtual void ConvertDigits(AliRawReader *rawReader, TTree *digitsTree) const; |
57 | virtual AliTracker* CreateTracker() const; | |
89971837 | 58 | TTreeSRedirector* GetDebugStream(ETRDReconstructorTask task) const { return task < kNtasks ? fDebugStream[task] : 0x0; } |
25ca55ce | 59 | |
3a039a31 | 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; | |
66f6bfd9 | 62 | static TClonesArray* GetClusters() {return fgClusters;} |
0d83b3a5 | 63 | Int_t GetNdEdxSlices() const { return (Int_t)AliTRDpidUtil::GetNdEdxSlices(GetPIDMethod());} |
89971837 | 64 | ETRDReconstructorGas GetDriftGas() const { return fSteerParam&kDriftGas ? kAr : kXe;} |
0d83b3a5 | 65 | AliTRDpidUtil::ETRDPIDMethod GetPIDMethod() const { return fSteerParam&kSteerPID ? AliTRDpidUtil::kNN : AliTRDpidUtil::kLQ;} |
3a039a31 | 66 | static const AliTRDrecoParam* GetRecoParam() { return dynamic_cast<const AliTRDrecoParam*>(AliReconstructor::GetRecoParam(2)); } |
89971837 | 67 | Int_t GetStreamLevel(ETRDReconstructorTask task) const { return fStreamLevel[task];} |
3a039a31 | 68 | inline void GetTCParams(Double_t *par) const; |
00872461 | 69 | virtual Bool_t HasDigitConversion() const { return fSteerParam&kDigitsConversion; }; |
3a039a31 | 70 | Bool_t IsWritingClusters() const { return fSteerParam&kWriteClusters;} |
71 | Bool_t IsWritingTracklets() const { return fSteerParam&kWriteTracklets;} | |
4302c900 | 72 | Bool_t IsHLT() const { return fSteerParam&kHLT;} |
3a039a31 | 73 | Bool_t IsSeeding() const { return fSteerParam&kSeeding;} |
91834b8d | 74 | Bool_t IsCosmic() const { return fSteerParam&kCosmic;} |
56b32fd7 | 75 | Bool_t IsEightSlices() const { return fSteerParam&kEightSlices;} |
3bc9d03e | 76 | |
3a039a31 | 77 | virtual void Reconstruct(AliRawReader *rawReader, TTree *clusterTree) const; |
78 | virtual void Reconstruct(TTree *digitsTree, TTree *clusterTree) const; | |
3bc9d03e | 79 | |
66f6bfd9 | 80 | static void SetClusters(TClonesArray *clusters) {fgClusters = clusters;} |
3a039a31 | 81 | void SetOption(Option_t *opt); |
82 | inline void SetTCParams(Double_t *par); | |
89971837 | 83 | void SetStreamLevel(Int_t level, ETRDReconstructorTask task= kTracker); |
3bc9d03e | 84 | |
3a039a31 | 85 | private: |
89971837 | 86 | UChar_t fStreamLevel[kNtasks];// stream level for each reconstruction task |
29f95561 | 87 | UInt_t fSteerParam; // steering flags |
88 | Double_t fTCParams[8]; // Tail Cancellation parameters for drift gases | |
89971837 | 89 | TTreeSRedirector *fDebugStream[kNtasks];// Debug Streamer container; |
3a039a31 | 90 | |
66f6bfd9 | 91 | static TClonesArray *fgClusters; // list of clusters for local reconstructor |
3bc9d03e | 92 | |
66f6bfd9 | 93 | ClassDef(AliTRDReconstructor, 1) // Class for the TRD reconstruction |
3bc9d03e | 94 | |
121a60bd | 95 | }; |
96 | ||
3a039a31 | 97 | //___________________________________________________ |
98 | inline void AliTRDReconstructor::GetTCParams(Double_t *par) const | |
99 | { | |
100 | if(!par) return; | |
101 | if(GetDriftGas()==kAr) memcpy(par, &fTCParams[4], 4*sizeof(Double_t)); | |
102 | else memcpy(par, &fTCParams[0], 4*sizeof(Double_t)); | |
103 | } | |
104 | ||
105 | //___________________________________________________ | |
106 | inline void AliTRDReconstructor::SetTCParams(Double_t *par) | |
107 | { | |
108 | if(!par) return; | |
109 | memcpy(fTCParams, par, 8*sizeof(Double_t)); | |
110 | } | |
111 | ||
121a60bd | 112 | #endif |