TRD/AliTRDtrackletMCM.h

   1 #ifndef ALITRDTRACKLETMCM_H
   2 #define ALITRDTRACKLETMCM_H
   3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   4  * See cxx source for full Copyright notice                               */
   5
   6 /* $Id: AliTRDtrackletMCM.h 27496 2008-07-22 08:35:45Z cblume $ */
   7
   8 //-----------------------------------
   9 //
  10 // TRD tracklet word (as from FEE)
  11 // only 32-bit of information + detector ID
  12 //
  13 //----------------------------------
  14
  15 #include "AliTRDtrackletBase.h"
  16 #include "AliTRDgeometry.h"
  17 #include "AliTRDpadPlane.h"
  18
  19 class AliTRDtrackletMCM : public AliTRDtrackletBase {
  20  public:
  21   AliTRDtrackletMCM(UInt_t trackletWord = 0);
  22   AliTRDtrackletMCM(UInt_t trackletWword, Int_t hcid);
  23   AliTRDtrackletMCM(UInt_t trackletWword, Int_t hcid, Int_t rob, Int_t mcm);
  24   AliTRDtrackletMCM(const AliTRDtrackletMCM &rhs);
  25   ~AliTRDtrackletMCM();
  26
  27   // ----- Getters for contents of tracklet word -----
  28   Int_t GetYbin() const;
  29   Int_t GetdY() const;
  30   Int_t GetZbin() const { return ((fTrackletWord >> 20) & 0xf); }
  31   Int_t GetPID() const { return ((fTrackletWord >> 24) & 0xff); }
  32
  33   // ----- Getters for MCM-tracklet information -----
  34   Int_t GetMCM() const { return fMCM; }
  35   Int_t GetROB() const { return fROB; }
  36   Int_t GetLabel() const { return fLabel; }
  37
  38   // ----- Getters for offline corresponding values -----
  39   Bool_t CookPID() { return kFALSE; }
  40   Double_t GetPID(Int_t /* is */) const { return GetPID()/255.; }
  41   Int_t GetDetector() const { return fHCId / 2; }
  42   Int_t GetHCId() const { return fHCId; }
  43   Float_t GetdYdX() const { return (GetdY() * 140e-4 / 3.); }
  44   Float_t GetX() const { return fGeo->GetTime0((fHCId % 12) / 2); }
  45   Float_t GetY() const { return (GetYbin() * 160e-4); }
  46   Float_t GetZ() const { return fGeo->GetPadPlane((fHCId % 12) / 2, (fHCId / 12) % 5)->GetRowPos( 4 * (fROB / 2) + fMCM / 4); }
  47
  48   Int_t GetQ0() const { return fQ0; }
  49   Int_t GetQ1() const { return fQ1; }
  50   Int_t GetNHits() const { return fNHits; }
  51   Int_t GetNHits0() const { return fNHits0; }
  52   Int_t GetNHits1() const { return fNHits1; }
  53
  54   UInt_t GetTrackletWord() const { return fTrackletWord; }
  55   void SetTrackletWord(UInt_t trackletWord) { fTrackletWord = trackletWord; }
  56
  57   void SetDetector(Int_t id) { fHCId = 2 * id + (GetYbin() < 0 ? 0 : 1); }
  58   void SetHCId(Int_t id) { fHCId = id; }
  59   void SetMCM(Int_t mcm) { fMCM = mcm; }
  60   void SetROB(Int_t rob) { fROB = rob; }
  61   void SetLabel(Int_t label) { fLabel = label; }
  62   void SetQ0(Int_t charge) { fQ0 = charge; }
  63   void SetQ1(Int_t charge) { fQ1 = charge; }
  64   void SetNHits(Int_t nhits) { fNHits = nhits; }
  65   void SetNHits0(Int_t nhits) { fNHits0 = nhits; }
  66   void SetNHits1(Int_t nhits) { fNHits1 = nhits; }
  67
  68   void SetSlope(Float_t slope) { fSlope = slope; }
  69   void SetOffset(Float_t offset) { fOffset = offset; }
  70   void SetError(Float_t error) { fError = error; }
  71   void SetClusters(Float_t *res, Float_t *q, Int_t n);
  72
  73   Float_t GetSlope() const { return fSlope; }
  74   Float_t GetOffset() const { return fOffset; }
  75   Float_t GetError() const { return fError; }
  76   Int_t   GetNClusters() const { return fNClusters; }
  77   Float_t *GetResiduals() const { return fResiduals; }
  78   Float_t *GetClsCharges() const { return fClsCharges; }
  79
  80  protected:
  81   AliTRDgeometry *fGeo; //! TRD geometry
  82
  83   Int_t fHCId;                  // half-chamber ID (only transient)
  84   UInt_t fTrackletWord;         // tracklet word: PID | Z | deflection length | Y
  85                                 //          bits:  12   4            7          13
  86   Int_t fMCM; // MCM no. in which the tracklet was found
  87   Int_t fROB; // ROB no. on which the tracklet was found
  88
  89   Int_t fQ0; // accumulated charge in the first time window
  90   Int_t fQ1; // accumulated charge in the second time window
  91
  92   Int_t fNHits;  // no. of contributing clusters
  93   Int_t fNHits0; // no. of contributing clusters in window 0
  94   Int_t fNHits1; // no. of contributing clusters in window 1
  95
  96   Int_t fLabel; // label for MC track
  97
  98   Float_t  fSlope;            // tracklet slope
  99   Float_t  fOffset;           // tracklet offset
 100   Float_t  fError;            // tracklet error
 101   Int_t    fNClusters;        // no. of clusters
 102   Float_t *fResiduals;        //[fNClusters] cluster to tracklet residuals
 103   Float_t *fClsCharges;       //[fNClusters] cluster charge
 104
 105  private:
 106   AliTRDtrackletMCM& operator=(const AliTRDtrackletMCM &rhs);   // not implemented
 107
 108   ClassDef(AliTRDtrackletMCM, 2);
 109 };
 110
 111 #endif