]> git.uio.no Git - u/mrichter/AliRoot.git/blob - TRD/AliTRDtrackletMCM.h
Update to store the accumulated charge of two time windows for PID
[u/mrichter/AliRoot.git] / 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 0; }
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
51   UInt_t GetTrackletWord() const { return fTrackletWord; }
52   void SetTrackletWord(UInt_t trackletWord) { fTrackletWord = trackletWord; }
53
54   void SetDetector(Int_t id) { fHCId = 2 * id + (GetYbin() < 0 ? 0 : 1); }
55   void SetHCId(Int_t id) { fHCId = id; }
56   void SetMCM(Int_t mcm) { fMCM = mcm; }
57   void SetROB(Int_t rob) { fROB = rob; }
58   void SetLabel(Int_t label) { fLabel = label; }
59   void SetQ0(Int_t charge) { fQ0 = charge; }
60   void SetQ1(Int_t charge) { fQ1 = charge; }
61
62  protected:
63   AliTRDgeometry *fGeo; //! TRD geometry
64
65   Int_t fHCId;                  // half-chamber ID (only transient)
66   UInt_t fTrackletWord;         // tracklet word: PID | Z | deflection length | Y 
67                                 //          bits:  12   4            7          13
68   Int_t fMCM; // MCM no. in which the tracklet was found
69   Int_t fROB; // ROB no. on which the tracklet was found
70
71   Int_t fQ0; // accumulated charge in the first time window
72   Int_t fQ1; // accumulated charge in the second time window
73
74   Int_t fLabel; // label for MC track
75
76  private:
77   AliTRDtrackletMCM& operator=(const AliTRDtrackletMCM &rhs);   // not implemented
78
79   ClassDef(AliTRDtrackletMCM, 2);
80 };
81
82 #endif