X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=TRD%2FAliTRDcluster.h;h=e7803054e0dd816dadc7208f02d6412895456ddb;hb=b033a1bfe19fe0f2bbb3800ba85a86d1f1a4cc71;hp=c291777004deaae99cb9435cadfee31554c6e2e0;hpb=6d50f529b7ba0b8b8cea934ac14d4969bb6107a9;p=u%2Fmrichter%2FAliRoot.git diff --git a/TRD/AliTRDcluster.h b/TRD/AliTRDcluster.h index c291777004d..e7803054e0d 100644 --- a/TRD/AliTRDcluster.h +++ b/TRD/AliTRDcluster.h @@ -11,71 +11,159 @@ // // /////////////////////////////////////////////////////////////////////////////// -#include "AliCluster.h" +#include "AliCluster.h" -class AliTRDrecPoint; +class AliTRDtrackletWord; class AliTRDcluster : public AliCluster { + friend class AliHLTTRDCluster; + friend class AliTRDtrackletOflHelper; - public: +public: + enum ETRDclusterStatus { + kInChamber = BIT(16) // Out of fiducial volume of chamber (signal tails) + ,kFivePad = BIT(17) // Deconvoluted clusters + ,kLUT = BIT(18) + ,kGAUS = BIT(19) + ,kCOG = BIT(20) + ,kXcorr = BIT(21) // steer efficient vd corrections + ,kRawSignals= BIT(22) // toggle raw digits storage + }; + enum ETRDclusterMask { + kMaskedLeft = 0 + ,kMaskedCenter = 1 + ,kMaskedRight = 2 + }; AliTRDcluster(); + AliTRDcluster(Int_t det, UChar_t col, UChar_t row, UChar_t time, const Short_t *sig, UShort_t volid); + AliTRDcluster(Int_t det, Float_t q, const Float_t *pos, const Float_t *sig + , const Int_t *tracks, Char_t npads, Short_t * const signals + , UChar_t col, UChar_t row, UChar_t time + , Char_t timebin, Float_t center, UShort_t volid); + AliTRDcluster(const AliTRDtrackletWord *const tracklet, Int_t det, UShort_t volid); AliTRDcluster(const AliTRDcluster &c); - AliTRDcluster(const AliTRDrecPoint &p); + virtual ~AliTRDcluster() {}; + AliTRDcluster &operator=(const AliTRDcluster &c); - virtual void AddTrackIndex(Int_t *i); - - Int_t IsUsed() const { return (fQ < 0) ? 1 : 0; } - void Use(Int_t = 0) { fQ = -fQ; } - - Bool_t From2pad() const { return TestBit(k2pad); } - Bool_t From3pad() const { return TestBit(k3pad); } - Bool_t From4pad() const { return TestBit(k4pad); } - Bool_t From5pad() const { return TestBit(k5pad); } - Bool_t FromLarge() const { return TestBit(kLarge); } - Bool_t Isolated() const { return (TestBit(k2pad) || TestBit(k3pad)); } - - virtual void SetDetector(Int_t d) { fDetector = d; } - void SetLocalTimeBin(Int_t t) { fTimeBin = t; } - void SetQ(Float_t q) { fQ = q; } - virtual void SetX(Float_t x) { fX = x; } - void SetCenter(Float_t c) { fCenter = c; } - void SetSignals(Short_t *signals); - void Set2pad() { SetBit(k2pad); fNPads = 2; } - void Set3pad() { SetBit(k3pad); fNPads = 3; } - void Set4pad() { SetBit(k4pad); fNPads = 4; } - void Set5pad() { SetBit(k5pad); fNPads = 5; } - void SetLarge() { SetBit(kLarge); fNPads = 6; } - - virtual Int_t GetDetector() const { return fDetector; } - Int_t GetLocalTimeBin() const { return fTimeBin; } - Float_t GetQ() const { return fQ; } - virtual Float_t GetX() const { return fX; } - Int_t GetNPads() const { return fNPads; } - Float_t GetCenter() const { return fCenter; } - Float_t GetSumS() const; - - protected: - - enum { - k2pad = 0x00000001, // 2 pad cluster - k3pad = 0x00000002, // 3 pad cluster - k4pad = 0x00000004, // 4 pad cluster - k5pad = 0x00000008, // 5 pad cluster - kLarge = 0x00000016 // Large cluster - }; + virtual void AddTrackIndex(const Int_t * const i); + void Clear(Option_t *o=""); - Int_t fDetector; // TRD detector number - Float_t fX; // Local x position (first order proportional to time bin) - Char_t fTimeBin; // Time bin number within the detector - Float_t fQ; // Amplitude - Char_t fNPads; // Number of pads in cluster - Float_t fCenter; // Center of the cluster relative to the pad - Short_t fSignals[7]; // Signals in the cluster - - ClassDef(AliTRDcluster,3) // Cluster for the TRD - + Bool_t IsEqual(const TObject *otherCluster) const; + Bool_t IsInChamber() const { return TestBit(kInChamber); } + Bool_t IsMasked() const { return fClusterMasking ? kTRUE : kFALSE; } + Bool_t IsShared() const { return IsClusterShared(); } + Bool_t IsUsed() const { return IsClusterUsed(); } + Bool_t IsFivePad() const { return TestBit(kFivePad); } + inline Bool_t IsRPhiMethod(ETRDclusterStatus m) const; + inline Bool_t IsMCcluster() const; + Bool_t HasXcorr() const { return TestBit(kXcorr); } + Bool_t HasRawSignals() const { return TestBit(kRawSignals); } + + UChar_t GetPadMaskedPosition() const { return fClusterMasking & 7; } + UChar_t GetPadMaskedStatus() const { return fClusterMasking >> 3; } + Int_t GetDetector() const { return fDetector; } + Int_t GetLocalTimeBin() const { return fLocalTimeBin; } + Float_t GetQ() const { return fQ; } + Short_t GetRawQ() const { return TestBit(kRawSignals)?(Short_t)GetSumS():0;} + Int_t GetNPads() const { return fNPads; } + Float_t GetCenter() const { return fCenter; } + Int_t GetPadCol() const { return fPadCol; } + Int_t GetPadRow() const { return fPadRow; } + Int_t GetPadTime() const { return fPadTime; } + Short_t *GetSignals() { return fSignals; } + Float_t GetSumS() const; + + static Double_t GetSX(Int_t tb, Double_t z=-1); + static Double_t GetSYdrift(Int_t tb, Int_t ly=0, Double_t z=-1); + static Double_t GetSYcharge(Float_t q); + static Double_t GetSYprf(Int_t ly, Double_t center, Double_t s2); + static Double_t GetXcorr(Int_t tb, Double_t z=-1); + static Double_t GetYcorr(Int_t ly, Float_t y); + Float_t GetXloc(Double_t t0, Double_t vd + , const Double_t *const q = 0x0 + , const Double_t *const xq = 0x0 + , Double_t z = 0.2); + Float_t GetYloc(Double_t y0, Double_t s2, Double_t W, Double_t *const y1=0x0, Double_t *const y2=0x0); + + void Print(Option_t* o="") const; + + void SetLocalTimeBin(Char_t t) { fLocalTimeBin = t; } + void SetNPads(Int_t n) { fNPads = n; } + void SetPadCol(UChar_t inPadCol) { fPadCol = inPadCol; } + void SetPadRow(UChar_t inPadRow) { fPadRow = inPadRow; } + void SetPadTime(UChar_t inPadTime) { fPadTime = inPadTime; } + void SetDetector(Short_t inDetector) { fDetector = inDetector; } + void SetQ(Float_t inQ) { fQ = inQ; } + void SetClusterMasking(UChar_t inClusterMasking) { fClusterMasking = inClusterMasking; } + void SetShared(Bool_t sh = kTRUE) { SetBit(AliCluster::kShared,sh); } + void SetFivePad(Bool_t b = kTRUE) { SetBit(kFivePad,b); } + void SetInChamber(Bool_t in = kTRUE) { SetBit(kInChamber,in); } + void SetPadMaskedPosition(UChar_t position); + void SetPadMaskedStatus(UChar_t status); + void SetSigmaY2(Float_t s2, Float_t dt, Float_t exb, Float_t x0, Float_t z=-1., Float_t tgp=0.); + void SetSignals(Short_t sig[7], Bool_t raw=kTRUE){ memcpy(fSignals, sig, 7*sizeof(Short_t)); SetBit(kRawSignals, raw);} + inline void SetRPhiMethod(ETRDclusterStatus m); + void SetXcorr(Bool_t xc = kTRUE) { SetBit(kXcorr,xc); } + + void Use(Int_t u = 1) { SetBit(AliCluster::kUsed, u ? kTRUE : kFALSE); } + +protected: + + UChar_t fPadCol; // Central pad number in column direction + UChar_t fPadRow; // Central pad number in row direction + UChar_t fPadTime; // Uncalibrated time bin number + Char_t fLocalTimeBin; // T0-calibrated time bin number + UChar_t fNPads; // Number of pads in cluster + UChar_t fClusterMasking; // Bit field containing cluster status information; + Short_t fDetector; // TRD detector number + Short_t fSignals[7]; // Raw signals if HasRawSignals() returns true. Otherwise calibrated. + Float_t fQ; // Calibrated cluster charge + Float_t fCenter; // Center of the cluster relative to the pad + +private: + + Float_t GetDYcog(const Double_t *const y1=0x0, const Double_t *const y2=0x0); + Float_t GetDYlut(const Double_t *const y1=0x0, const Double_t *const y2=0x0); + Float_t GetDYgauss(Double_t sw, Double_t *const y1=0x0, Double_t *const y2=0x0); + static void FillLUT(); + + static const Int_t fgkNlut; //! Number of bins of the LUT + static Double_t *fgLUT; //! The lookup table + + ClassDef(AliTRDcluster, 7) // Cluster for the TRD + }; +//________________________________________________ +inline Bool_t AliTRDcluster::IsMCcluster() const +{ + if( GetLabel(0) == fPadRow && + GetLabel(1) == fPadCol && + GetLabel(2) == fPadTime) return kFALSE; + return kTRUE; +} + +//________________________________________________ +inline Bool_t AliTRDcluster::IsRPhiMethod(ETRDclusterStatus m) const +{ + if(m==kLUT && TestBit(kLUT)) return kTRUE; + else if(m==kGAUS && TestBit(kGAUS)) return kTRUE; + else if(m==kCOG && (!TestBit(kLUT)&&!TestBit(kGAUS))) return kTRUE; + + return kFALSE; +} + +//________________________________________________ +inline void AliTRDcluster::SetRPhiMethod(ETRDclusterStatus m) +{ + SetBit(kCOG,0);SetBit(kLUT,0);SetBit(kGAUS,0); + switch(m){ + case kCOG: SetBit(kCOG); break; + case kLUT: SetBit(kLUT); break; + case kGAUS: SetBit(kGAUS); break; + default: SetBit(kLUT); break; + } +} #endif