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 "AliRecoParam.h"
16 #include "AliDetectorRecoParam.h"
17 #include "AliTRDpidUtil.h"
18 #include "AliTRDrecoParam.h"
19 #include "AliTRDdigitsParam.h"
20 #include "AliESDTrdTrigger.h"
23 class TTreeSRedirector;
25 class AliTRDclusterizer;
26 class AliTRDonlineTrackMatching;
27 class AliTRDReconstructor: public AliReconstructor
30 enum ETRDReconstructorSteer {
31 kDigitsConversion= BIT(0)
32 ,kWriteClusters = BIT(1)
33 ,kWriteTracklets = BIT(2)
36 ,kProcTracklets = BIT(5) // process online tracklets
38 ,kClRadialCorr = BIT(7) // toggle radial correction in clusters
40 ,kNsteer = 8 // number of tasks
42 AliTRDReconstructor();
43 virtual ~AliTRDReconstructor();
47 virtual void ConvertDigits(AliRawReader *rawReader, TTree *digitsTree) const;
48 virtual AliTracker* CreateTracker() const;
49 TTreeSRedirector* GetDebugStream(AliTRDrecoParam::ETRDReconstructionTask task) const { return task < AliTRDrecoParam::kTRDreconstructionTasks ? fDebugStream[task] : NULL; }
51 virtual void FillESD(AliRawReader *, TTree *clusterTree, AliESDEvent *esd) const { FillESD((TTree * )NULL, clusterTree, esd); }
52 virtual void FillESD(TTree *digitsTree, TTree *clusterTree, AliESDEvent *esd) const;
53 static TClonesArray* GetClusters();
54 static TClonesArray* GetTracklets(const char *trkltype = "");
55 static TClonesArray* GetTracks();
56 static Int_t GetNTimeBins() { return fgNTimeBins;}
57 Int_t GetNdEdxSlices() const { return (Int_t)AliTRDpidUtil::GetNdEdxSlices(GetPIDMethod());}
58 AliTRDpidUtil::ETRDPIDMethod GetPIDMethod() const { return GetRecoParam()->IsPIDNeuralNetwork() ? AliTRDpidUtil::kNN : AliTRDpidUtil::kLQ;}
59 static const AliTRDrecoParam* GetRecoParam() { return dynamic_cast<const AliTRDrecoParam*>(AliReconstructor::GetRecoParam(2)); }
60 static Float_t GetMinClustersInTrack() { return fgkMinClustersInTrack;}
61 static Float_t GetLabelFraction() { return fgkLabelFraction;}
62 static Double_t GetMaxChi2() { return fgkMaxChi2;}
63 static Double_t GetMaxSnp() { return fgkMaxSnp;}
64 static Double_t GetMaxStep() { return fgkMaxStep;}
65 static Double_t GetEpsilon() { return fgkEpsilon;}
67 virtual Bool_t HasDigitConversion() const { return fSteerParam&kDigitsConversion; };
68 Bool_t IsCosmic() const { return GetRecoParam()->GetEventSpecie() & AliRecoParam::kCosmic;}
69 Bool_t IsWritingClusters() const { return fSteerParam&kWriteClusters;}
70 Bool_t IsWritingTracklets() const { return fSteerParam&kWriteTracklets;}
71 Bool_t IsHLT() const { return fSteerParam&kHLT;}
72 Bool_t IsSeeding() const { return fSteerParam&kSeeding;}
73 Bool_t IsProcessingTracklets() const { return fSteerParam&kProcTracklets;}
74 Bool_t IsDebugStreaming() const { return (fSteerParam&kDebug || AliTRDReconstructor::GetStreamLevel()>0);}
75 Bool_t UseClusterRadialCorrection() const { return fSteerParam&kClRadialCorr;}
77 static void Options(UInt_t steer=0);
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 static void SetTracklets(TClonesArray *tracklets) { fgTracklets = tracklets;}
83 static void SetTracks(TClonesArray *tracks) { fgTracks = tracks;}
84 void SetOption(Option_t *opt);
85 static Int_t GetStreamLevel() { return fgStreamLevel;}
86 static void SetStreamLevel(Int_t level) { fgStreamLevel = level;}
89 AliTRDReconstructor(const AliTRDReconstructor &r); //Not implemented
90 AliTRDReconstructor& operator = (const AliTRDReconstructor&); //Not implemented
91 void ResetContainers() const;
92 static Int_t fgStreamLevel; // flag for streaming - for TRD reconstruction
94 static Char_t const *fgSteerNames[kNsteer];//! steering names
95 static Char_t const *fgSteerFlags[kNsteer];//! steering flags
96 static Char_t const *fgTaskNames[AliTRDrecoParam::kTRDreconstructionTasks]; //! tasks names
97 static Char_t const *fgTaskFlags[AliTRDrecoParam::kTRDreconstructionTasks]; //! tasks flags
98 UInt_t fSteerParam; // steering bits
99 // configuration vars for tracking
100 static const Double_t fgkMaxChi2; // Max increment in track chi2
101 static const Float_t fgkMinClustersInTrack; // Min number of clusters in track
102 static const Float_t fgkLabelFraction; // Min fraction of same label
103 static const Double_t fgkMaxSnp; // Maximal snp for tracking
104 static const Double_t fgkMaxStep; // Maximal step for tracking
105 static const Double_t fgkEpsilon; // Precision of radial coordinate
107 TTreeSRedirector *fDebugStream[AliTRDrecoParam::kTRDreconstructionTasks];// Debug Streamer container;
109 static TClonesArray *fgClusters; // list of clusters for local reconstructor
110 static TClonesArray *fgTracklets; // list of online tracklets for local reconstructor
111 static TClonesArray *fgTracks; // list of GTU tracks for local reconstructor
112 static Int_t fgNTimeBins; // number of time bins as given by the clusterizer
113 AliTRDclusterizer *fClusterizer; //! instance of TRD clusterizer
114 static AliTRDonlineTrackMatching fgOnlineTrackMatcher; // track matcher between on-line and off-line track
115 static AliESDTrdTrigger fgTriggerFlags; // L1 trigger flags
117 ClassDef(AliTRDReconstructor, 5) // Class for the TRD reconstruction