]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - TRD/AliTRDReconstructor.h
The present commit corresponds to an important change in the way the
[u/mrichter/AliRoot.git] / TRD / AliTRDReconstructor.h
index 624746f69a5c04af815566a6c5f14c3d6aaef4ee..19b0ba7c01aa927f811b7a97ccd4263cf94b84f0 100644 (file)
 ///////////////////////////////////////////////////////////////////////////////
 
 #include "AliReconstructor.h"
+#include "AliDetectorRecoParam.h"
+#include "AliTRDrecoParam.h"
+#include "AliTRDpidUtil.h"
 
+class TClonesArray;
+class TTreeSRedirector;
 class AliRawReader;
+class AliTRDReconstructor: public AliReconstructor 
+{
+public:
+  enum AliTRDsteerParam {
+    kDigitsConversion= BIT(0)
+    ,kWriteClusters  = BIT(1)
+    ,kSeeding        = BIT(2)
+    ,kSteerPID       = BIT(3)
+    ,kEightSlices    = BIT(4)
+    ,kWriteTracklets = BIT(5)
+    ,kDriftGas       = BIT(6)
+    ,kHLT            = BIT(7)
+    ,kCosmic         = BIT(8)
+  };
+  enum AliTRDReconstructorTask {
+    kRawReader    = 0
+    ,kClusterizer = 1
+    ,kTracker     = 2
+    ,kPID         = 3
+  };
+  enum AliTRDdriftGas {
+    kXe = 0,
+    kAr = 1
+  };
+  enum{
+    kOwner = BIT(14)
+  };
+
+  AliTRDReconstructor();
+  AliTRDReconstructor(const AliTRDReconstructor &r);
+  virtual ~AliTRDReconstructor();
+  AliTRDReconstructor& operator = (const AliTRDReconstructor&)          { return *this;}
+       
+       virtual void        Init();
 
-class AliTRDReconstructor: public AliReconstructor {
-
- public:
-
-  AliTRDReconstructor():AliReconstructor()                       { };
-  virtual ~AliTRDReconstructor()                                 { };
-
-  //virtual Bool_t      HasDigitConversion() const                 { return kTRUE; };
-  virtual Bool_t      HasDigitConversion() const                 { return kFALSE; };
   virtual void        ConvertDigits(AliRawReader *rawReader, TTree *digitsTree) const;
+  virtual AliTracker* CreateTracker() const;
+  TTreeSRedirector*   GetDebugStream(AliTRDReconstructorTask task) const { return task < 4 ? fDebugStream[task] : 0x0; }
 
-  virtual Bool_t      HasLocalReconstruction() const             { return kTRUE; };
-  virtual void        Reconstruct(AliRunLoader *runLoader, AliRawReader *rawReader) const;
-  virtual void        Reconstruct(AliRawReader *rawReader, TTree *clusterTree) const;
-  virtual void        Reconstruct(TTree *digitsTree, TTree *clusterTree) const;
-  virtual void        Reconstruct(AliRunLoader *runLoader) const;
-
-  virtual AliTracker *CreateTracker(AliRunLoader *runLoader) const;
-
-  virtual void        FillESD(AliRunLoader *runLoader, AliRawReader *rawReader, AliESDEvent *esd) const;
-  virtual void        FillESD(AliRawReader *rawReader, TTree *clusterTree, AliESDEvent *esd) const;
+  virtual void        FillESD(AliRawReader *, TTree *clusterTree, AliESDEvent *esd) const { FillESD((TTree * )NULL, clusterTree, esd);                    }
   virtual void        FillESD(TTree *digitsTree, TTree *clusterTree, AliESDEvent *esd) const;
-  virtual void        FillESD(AliRunLoader *runLoader, AliESDEvent *esd) const;
-
-  static  void        SetSeedingOn(Bool_t seeding)               { fgkSeedingOn  = seeding; }  
-  static  void        SetStreamLevel(Int_t level)                { fgStreamLevel = level;   }
+  static TClonesArray* GetClusters() {return fgClusters;}
+  Int_t               GetNdEdxSlices() const     { return (Int_t)AliTRDpidUtil::GetNdEdxSlices(GetPIDMethod());}
+  AliTRDdriftGas      GetDriftGas() const        { return fSteerParam&kDriftGas ? kAr : kXe;}
+  AliTRDpidUtil::ETRDPIDMethod       GetPIDMethod() const       { return fSteerParam&kSteerPID ? 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 fSteerParam&kDigitsConversion;  };
+  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              IsEightSlices() const      { return fSteerParam&kEightSlices;}
 
-  static  Bool_t      SeedingOn()                                { return fgkSeedingOn;     }
-  static  Int_t       StreamLevel()                              { return fgStreamLevel;    }
+  virtual void        Reconstruct(AliRawReader *rawReader, TTree *clusterTree) const;
+  virtual void        Reconstruct(TTree *digitsTree, TTree *clusterTree) const;
 
- private:
+  static void         SetClusters(TClonesArray *clusters) {fgClusters = clusters;}
+  void               SetOption(Option_t *opt);
+  inline void         SetTCParams(Double_t *par);
+  void                SetStreamLevel(Int_t level, AliTRDReconstructorTask task= kTracker);
 
-  static  Bool_t   fgkSeedingOn;  //  Set flag for seeding during reconstruction
-  static  Int_t    fgStreamLevel; //  Flag for streaming
+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 
+  TTreeSRedirector *fDebugStream[4];      // Debug Streamer container;
+  static TClonesArray *fgClusters;    // list of clusters for local reconstructor
 
-  ClassDef(AliTRDReconstructor,0) //  Class for the TRD reconstruction
+  ClassDef(AliTRDReconstructor, 1)    //  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