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 */
8 ///////////////////////////////////////////////////////////////////////////////
10 // Class for TRD reconstruction //
12 ///////////////////////////////////////////////////////////////////////////////
14 #include "AliReconstructor.h"
15 #include "AliDetectorRecoParam.h"
16 #include "AliTRDrecoParam.h"
17 #include "AliTRDpidUtil.h"
20 class TTreeSRedirector;
22 class AliTRDReconstructor: public AliReconstructor
25 enum AliTRDsteerParam {
26 kDigitsConversion= BIT(0)
27 ,kWriteClusters = BIT(1)
30 ,kEightSlices = BIT(4)
31 ,kWriteTracklets = BIT(5)
36 enum AliTRDReconstructorTask {
50 AliTRDReconstructor();
51 AliTRDReconstructor(const AliTRDReconstructor &r);
52 virtual ~AliTRDReconstructor();
53 AliTRDReconstructor& operator = (const AliTRDReconstructor&) { return *this;}
57 virtual void ConvertDigits(AliRawReader *rawReader, TTree *digitsTree) const;
58 virtual AliTracker* CreateTracker() const;
59 TTreeSRedirector* GetDebugStream(AliTRDReconstructorTask task) const { return task < 4 ? fDebugStream[task] : 0x0; }
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 (Int_t)AliTRDpidUtil::GetNdEdxSlices(GetPIDMethod());}
65 AliTRDdriftGas GetDriftGas() const { return fSteerParam&kDriftGas ? kAr : kXe;}
66 AliTRDpidUtil::ETRDPIDMethod GetPIDMethod() const { return fSteerParam&kSteerPID ? AliTRDpidUtil::kNN : AliTRDpidUtil::kLQ;}
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 Bool_t IsEightSlices() const { return fSteerParam&kEightSlices;}
78 virtual void Reconstruct(AliRawReader *rawReader, TTree *clusterTree) const;
79 virtual void Reconstruct(TTree *digitsTree, TTree *clusterTree) const;
81 static void SetClusters(TClonesArray *clusters) {fgClusters = clusters;}
82 void SetOption(Option_t *opt);
83 inline void SetTCParams(Double_t *par);
84 void SetStreamLevel(Int_t level, AliTRDReconstructorTask task= kTracker);
87 UChar_t fStreamLevel[5]; // stream level for each reconstruction task
88 UInt_t fSteerParam; // steering flags
89 Double_t fTCParams[8]; // Tail Cancellation parameters for drift gases
90 TTreeSRedirector *fDebugStream[4]; // Debug Streamer container;
92 static TClonesArray *fgClusters; // list of clusters for local reconstructor
94 ClassDef(AliTRDReconstructor, 1) // Class for the TRD reconstruction
98 //___________________________________________________
99 inline void AliTRDReconstructor::GetTCParams(Double_t *par) const
102 if(GetDriftGas()==kAr) memcpy(par, &fTCParams[4], 4*sizeof(Double_t));
103 else memcpy(par, &fTCParams[0], 4*sizeof(Double_t));
106 //___________________________________________________
107 inline void AliTRDReconstructor::SetTCParams(Double_t *par)
110 memcpy(fTCParams, par, 8*sizeof(Double_t));