- return pad centre for z-position (before pad edge)
authorcblume <cblume@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 16 Jun 2010 15:15:45 +0000 (15:15 +0000)
committercblume <cblume@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 16 Jun 2010 15:15:45 +0000 (15:15 +0000)
- additional methods to provide ROB/MCM position for real tracklets

TRD/AliTRDtrackletMCM.h
TRD/AliTRDtrackletWord.cxx
TRD/AliTRDtrackletWord.h

index cebaf3b..118ce76 100644 (file)
@@ -43,7 +43,9 @@ class AliTRDtrackletMCM : public AliTRDtrackletBase {
   Float_t GetdYdX() const { return (GetdY() * 140e-4 / 3.); }
   Float_t GetX() const { return fGeo->GetTime0((fHCId % 12) / 2); }
   Float_t GetY() const { return (GetYbin() * 160e-4); }
-  Float_t GetZ() const { return fGeo->GetPadPlane((fHCId % 12) / 2, (fHCId / 12) % 5)->GetRowPos( 4 * (fROB / 2) + fMCM / 4); }
+  Float_t GetZ() const { return fGeo->GetPadPlane((fHCId % 12) / 2, (fHCId / 12) % 5)->GetRowPos( 4 * (fROB / 2) + fMCM / 4) - 
+      fGeo->GetPadPlane((fHCId % 12) / 2, (fHCId /12) % 5)->GetRowSize(4 * (fROB / 2) + fMCM / 4) * .5; }
+  Float_t GetLocalZ() const { return GetZ() - fGeo->GetPadPlane((fHCId % 12) / 2, (fHCId / 12) % 5)->GetRowPos(8); }
 
   Int_t GetQ0() const { return fQ0; }
   Int_t GetQ1() const { return fQ1; }
index 2083d1f..8c2102f 100644 (file)
@@ -25,6 +25,7 @@
 
 #include "AliTRDtrackletWord.h"
 #include "AliTRDgeometry.h"
+#include "AliTRDpadPlane.h"
 #include "AliLog.h"
 
 ClassImp(AliTRDtrackletWord)
@@ -80,3 +81,16 @@ Int_t AliTRDtrackletWord::GetdY() const
     return ((fTrackletWord >> 13) & 0x7f);
   }
 }
+
+Int_t AliTRDtrackletWord::GetROB() const
+{
+  return 2 * (GetZbin() / 4) + (GetY() > 0 ? 1 : 0);
+}
+
+Int_t AliTRDtrackletWord::GetMCM() const
+{
+  AliTRDpadPlane *pp = fgGeo->GetPadPlane(GetDetector());
+  return (((Int_t) ((GetY()) / pp->GetWidthIPad()) + 72) / 18) % 4
+    + 4 * (GetZbin() % 4);
+}
+
index 030e050..352c36f 100644 (file)
@@ -29,6 +29,9 @@ class AliTRDtrackletWord : public AliTRDtrackletBase {
   Int_t GetZbin() const { return ((fTrackletWord >> 20) & 0xf); }
   Int_t GetPID() const { return ((fTrackletWord >> 24) & 0xff); }
 
+  Int_t GetROB() const;
+  Int_t GetMCM() const; 
+
   // ----- Getters for offline corresponding values -----
   Bool_t CookPID() { return kFALSE; }
   Double_t GetPID(Int_t /* is */) const { return 0; }
@@ -37,7 +40,9 @@ class AliTRDtrackletWord : public AliTRDtrackletBase {
   Float_t GetdYdX() const { return (GetdY() * 140e-4 / 3.); }
   Float_t GetX() const { return fgGeo->GetTime0((fHCId%12)/2); }
   Float_t GetY() const { return (GetYbin() * 160e-4); }
-  Float_t GetZ() const { return fgGeo->GetPadPlane((fHCId % 12) / 2, (fHCId/12) % 5)->GetRowPos(GetZbin()); }
+  Float_t GetZ() const { return fgGeo->GetPadPlane((fHCId % 12) / 2, (fHCId/12) % 5)->GetRowPos(GetZbin()) -
+      fgGeo->GetPadPlane((fHCId % 12) / 2, (fHCId/12) % 5)->GetRowSize(GetZbin()); }
+  Float_t GetLocalZ() const { return GetZ() - fgGeo->GetPadPlane((fHCId % 12) / 2, (fHCId/12) % 5)->GetRowPos(8); }
 
   UInt_t GetTrackletWord() const { return fTrackletWord; }
   void SetTrackletWord(UInt_t trackletWord) { fTrackletWord = trackletWord; }
@@ -46,7 +51,7 @@ class AliTRDtrackletWord : public AliTRDtrackletBase {
   void SetHCId(Int_t id) { fHCId = id; }
 
  protected:
-  Int_t fHCId;                  //! half-chamber ID (only transient)
+  Int_t fHCId;                  // half-chamber ID
   UInt_t fTrackletWord;                // tracklet word: PID | Z | deflection length | Y 
                                //          bits:  12   4            7          13
   static AliTRDgeometry *fgGeo;  // pointer to TRD geometry for coordinate calculations
@@ -54,7 +59,7 @@ class AliTRDtrackletWord : public AliTRDtrackletBase {
  private:
   AliTRDtrackletWord& operator=(const AliTRDtrackletWord &rhs);   // not implemented
 
-  ClassDef(AliTRDtrackletWord, 1);
+  ClassDef(AliTRDtrackletWord, 2);
 };
 
 #endif