TRD/AliTRDReconstructor.h

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