///////////////////////////////////////////////////////////////////////////////
#include "AliReconstructor.h"
+#include "AliRecoParam.h"
#include "AliDetectorRecoParam.h"
+#include "AliTRDpidUtil.h"
#include "AliTRDrecoParam.h"
+#include "AliTRDdigitsParam.h"
class TClonesArray;
+class TTreeSRedirector;
class AliRawReader;
+class AliTRDclusterizer;
class AliTRDReconstructor: public AliReconstructor
{
public:
- enum AliTRDsteerParam {
- kWriteClusters = BIT(0)
- ,kSeeding = BIT(1)
- ,kSteerPID = BIT(2)
- ,kWriteTracklets = BIT(3)
- ,kDriftGas = BIT(4)
- ,kHLT = BIT(5)
- ,kCosmic = BIT(6)
+ enum ETRDReconstructorSteer {
+ kDigitsConversion= BIT(0)
+ ,kWriteClusters = BIT(1)
+ ,kWriteTracklets = BIT(2)
+ ,kSeeding = BIT(3)
+ ,kHLT = BIT(4)
+ ,kProcTracklets = BIT(5) // process online tracklets
+ ,kDebug = BIT(6)
+ ,kClRadialCorr = BIT(7) // toggle radial correction in clusters
+ ,kOwner = BIT(8)
+ ,kNsteer = 8 // number of tasks
};
- enum AliTRDReconstructorTask {
- kClusterizer = 0
- ,kTracker = 1
- ,kPID = 2
- };
- enum AliTRDpidMethod {
- kLQPID = 0,
- kNNPID = 1
- };
- enum AliTRDdriftGas {
- kXe = 0,
- kAr = 1
- };
- enum{
- kNNslices = 8
- ,kLQslices = 3
- };
-
AliTRDReconstructor();
- AliTRDReconstructor(const AliTRDReconstructor &r);
virtual ~AliTRDReconstructor();
- AliTRDReconstructor& operator = (const AliTRDReconstructor&) { return *this;}
virtual void Init();
virtual void ConvertDigits(AliRawReader *rawReader, TTree *digitsTree) const;
virtual AliTracker* CreateTracker() const;
+ TTreeSRedirector* GetDebugStream(AliTRDrecoParam::ETRDReconstructionTask task) const { return task < AliTRDrecoParam::kTRDreconstructionTasks ? fDebugStream[task] : NULL; }
virtual void FillESD(AliRawReader *, TTree *clusterTree, AliESDEvent *esd) const { FillESD((TTree * )NULL, clusterTree, esd); }
virtual void FillESD(TTree *digitsTree, TTree *clusterTree, AliESDEvent *esd) const;
static TClonesArray* GetClusters() {return fgClusters;}
- Int_t GetNdEdxSlices() const { return GetPIDMethod() == kNNPID ? kNNslices : kLQslices;}
- AliTRDdriftGas GetDriftGas() const { return fSteerParam&kDriftGas ? kAr : kXe;}
- AliTRDpidMethod GetPIDMethod() const { return fSteerParam&kSteerPID ? kNNPID : kLQPID;}
+ static TClonesArray* GetTracklets() {return fgTracklets;}
+ static Int_t GetNTimeBins() {return fgNTimeBins;}
+ Int_t GetNdEdxSlices() const { return (Int_t)AliTRDpidUtil::GetNdEdxSlices(GetPIDMethod());}
+ AliTRDpidUtil::ETRDPIDMethod GetPIDMethod() const { return GetRecoParam()->IsPIDNeuralNetwork() ? AliTRDpidUtil::kNN : AliTRDpidUtil::kLQ;}
static const AliTRDrecoParam* GetRecoParam() { return dynamic_cast<const AliTRDrecoParam*>(AliReconstructor::GetRecoParam(2)); }
- Int_t GetStreamLevel(AliTRDReconstructorTask task) const { return fStreamLevel[task];}
- inline void GetTCParams(Double_t *par) const;
- virtual Bool_t HasDigitConversion() const { return kFALSE; };
+
+ virtual Bool_t HasDigitConversion() const { return fSteerParam&kDigitsConversion; };
+ Bool_t IsCosmic() const { return GetRecoParam()->GetEventSpecie() & AliRecoParam::kCosmic;}
Bool_t IsWritingClusters() const { return fSteerParam&kWriteClusters;}
Bool_t IsWritingTracklets() const { return fSteerParam&kWriteTracklets;}
Bool_t IsHLT() const { return fSteerParam&kHLT;}
Bool_t IsSeeding() const { return fSteerParam&kSeeding;}
- Bool_t IsCosmic() const { return fSteerParam&kCosmic;}
+ Bool_t IsProcessingTracklets() const { return fSteerParam&kProcTracklets;}
+ Bool_t IsDebugStreaming() const { return fSteerParam&kDebug;}
+ Bool_t UseClusterRadialCorrection() const { return fSteerParam&kClRadialCorr;}
+ static void Options(UInt_t steer=0);
virtual void Reconstruct(AliRawReader *rawReader, TTree *clusterTree) const;
virtual void Reconstruct(TTree *digitsTree, TTree *clusterTree) const;
- static void SetClusters(TClonesArray *clusters) {fgClusters = clusters;}
+ static void SetClusters(TClonesArray *clusters) { fgClusters = clusters;}
+ static void SetTracklets(TClonesArray *tracklets) { fgTracklets = tracklets;}
void SetOption(Option_t *opt);
- inline void SetTCParams(Double_t *par);
- void SetStreamLevel(Int_t level, AliTRDReconstructorTask task= kTracker);
private:
- UChar_t fStreamLevel[5]; // stream level for each reconstruction task
- UInt_t fSteerParam; // steering flags
- Double_t fTCParams[8]; // Tail Cancellation parameters for drift gases
+ AliTRDReconstructor(const AliTRDReconstructor &r); //Not implemented
+ AliTRDReconstructor& operator = (const AliTRDReconstructor&); //Not implemented
+
+ static Char_t const *fgSteerNames[kNsteer];//! steering names
+ static Char_t const *fgSteerFlags[kNsteer];//! steering flags
+ static Char_t const *fgTaskNames[AliTRDrecoParam::kTRDreconstructionTasks]; //! tasks names
+ static Char_t const *fgTaskFlags[AliTRDrecoParam::kTRDreconstructionTasks]; //! tasks flags
+ UInt_t fSteerParam; // steering bits
+ TTreeSRedirector *fDebugStream[AliTRDrecoParam::kTRDreconstructionTasks];// Debug Streamer container;
- static TClonesArray *fgClusters; // list of clusters for local reconstructor
+ static TClonesArray *fgClusters; // list of clusters for local reconstructor
+ static TClonesArray *fgTracklets; // list of online tracklets for local reconstructor
+ static Int_t fgNTimeBins; // number of time bins as given by the clusterizer
+ AliTRDclusterizer *fClusterizer; //! instance of TRD clusterizer
- ClassDef(AliTRDReconstructor, 1) // Class for the TRD reconstruction
+ ClassDef(AliTRDReconstructor, 4) // Class for the TRD reconstruction
};
-//___________________________________________________
-inline void AliTRDReconstructor::GetTCParams(Double_t *par) const
-{
- if(!par) return;
- if(GetDriftGas()==kAr) memcpy(par, &fTCParams[4], 4*sizeof(Double_t));
- else memcpy(par, &fTCParams[0], 4*sizeof(Double_t));
-}
-//___________________________________________________
-inline void AliTRDReconstructor::SetTCParams(Double_t *par)
-{
- if(!par) return;
- memcpy(fTCParams, par, 8*sizeof(Double_t));
-}
#endif
+