]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - TPC/AliTPCseed.h
newline added at the end of file
[u/mrichter/AliRoot.git] / TPC / AliTPCseed.h
index 283ad16f5ad5f5512a95cd3f5926ea5d62f67555..2184e8aa8f34737198c9a77e8bdcd7d303ce5320 100644 (file)
@@ -14,6 +14,7 @@
 //-------------------------------------------------------
 
 #include <TError.h>
+#include <TBits.h>
 
 #include "AliTPCtrack.h"
 #include "AliComplexCluster.h"
@@ -28,7 +29,6 @@ class AliESD;
 class TClonesArray;
 
 class AliTPCseed : public AliTPCtrack {
-  friend class AliTPCtrackerMI;
   public:  
      AliTPCseed();
      virtual ~AliTPCseed();
@@ -36,7 +36,8 @@ class AliTPCseed : public AliTPCtrack {
      AliTPCseed(const AliTPCseed &s, Bool_t clusterOwner = kFALSE);
      //AliTPCseed(const AliTPCseed &t, Double_t a);
      AliTPCseed(Double_t xr, Double_t alpha, const Double_t xx[5], 
-                const Double_t cc[15], Int_t i);     
+                const Double_t cc[15], Int_t i);   
+     AliTPCseed &operator = (const AliTPCseed & param);  
      Int_t Compare(const TObject *o) const;
      void Reset(Bool_t all = kTRUE);
      Int_t GetProlongation(Double_t xr, Double_t &y, Double_t & z) const;
@@ -44,6 +45,7 @@ class AliTPCseed : public AliTPCtrack {
      virtual Bool_t Update(const AliCluster* c2, Double_t chi2, Int_t i);
      AliTPCTrackerPoint * GetTrackPoint(Int_t i);
      AliTPCclusterMI * GetClusterFast(Int_t irow){ return fClusterPointer[irow];}
+     void SetClusterPointer(Int_t irow, AliTPCclusterMI* cl) {fClusterPointer[irow]=cl;}
      void RebuildSeed(); // rebuild seed to be ready for storing
      Double_t GetDensityFirst(Int_t n);
      Double_t GetSigma2C() const {
@@ -70,13 +72,71 @@ class AliTPCseed : public AliTPCtrack {
      void SetErrorY2(Float_t sy2){fErrorY2=sy2;}
      void SetErrorZ2(Float_t sz2){fErrorZ2=sz2;}
      Float_t  CookdEdx(Double_t low=0.05, Double_t up=0.70, Int_t i1=0, Int_t i2=159, Bool_t onlyused = kFALSE);
+     Float_t  CookShape(Int_t type);
      void CookPID();
      Double_t Bethe(Double_t bg);     // return bethe-bloch
+     Double_t BetheMass(Double_t mass=0.1057);     // return bethe-bloch
      //     void CookdEdx2(Double_t low=0.05, Double_t up=0.70);
      Bool_t IsActive() const { return !(fRemoval);}
      void Desactivate(Int_t reason){ fRemoval = reason;} 
+     AliTPCclusterMI* GetClusterPointer(Int_t i) const {return fClusterPointer[i];}
+     Int_t GetSector() const {return fSector;}
+     Float_t GetCurrentSigmaY2() const {return fCurrentSigmaY2;}
+     Float_t GetCurrentSigmaZ2() const {return fCurrentSigmaZ2;}
+     Int_t GetRelativeSector() const {return fRelativeSector;}
+     Char_t GetCircular() const {return fCircular;}
+
+     void SetCurrentSigmaY2(Float_t s) {fCurrentSigmaY2=s;}
+     void SetCurrentSigmaZ2(Float_t s) {fCurrentSigmaZ2=s;}
+     void SetRelativeSector(Int_t r) {fRelativeSector=r;}
+     void SetCircular(Char_t c) {fCircular=c;}
+     void SetIsSeeding(Bool_t s) {fIsSeeding=s;}
+     void SetSeedType(Int_t s) {fSeedType=s;}
+     void SetSeed1(Int_t s) {fSeed1=s;}
+     void SetSeed2(Int_t s) {fSeed2=s;}
+     void SetESD(AliESDtrack* esd) {fEsd=esd;}
+     void SetBSigned(Bool_t s) {fBSigned=s;}
+     void SetSort(Int_t s) {fSort=s;}
+     void SetOverlapLabel(Int_t i, Int_t l) {fOverlapLabels[i]=l;}
+     void SetCurrentCluster(AliTPCclusterMI* cl) {fCurrentCluster=cl;}
+     void SetNoCluster(Int_t n) {fNoCluster=n;}
+     void SetRow(Int_t n) {fRow=n;}
+     void SetSector(Int_t n) {fSector=n;}
+     void SetCurrentClusterIndex1(Int_t n) {fCurrentClusterIndex1=n;}
+     void SetInDead(Bool_t s) {fInDead=s;}
+     void SetPoints(TClonesArray* p) {fPoints=p;}
+     void SetEPoints(TClonesArray* p) {fEPoints=p;}
+
+     Double_t TPCrPID(Int_t i) const {return fTPCr[i];}
+     Double_t* TPCrPIDs() {return fTPCr;}
+     Bool_t GetIsSeeding() const {return fIsSeeding;}
+     Int_t GetSeedType() const {return fSeedType;}
+     Int_t GetSeed1() const {return fSeed1;}
+     Int_t GetSeed2() const {return fSeed2;}
+     AliESDtrack* GetESD() {return fEsd;}
+     Float_t GetSDEDX(Int_t i) const {return fSDEDX[i];}
+     Int_t GetNCDEDX(Int_t i) const {return fNCDEDX[i];}
+     Bool_t GetBSigned() const {return fBSigned;}
+     Int_t GetSort() const {return fSort;}
+     Int_t GetOverlapLabel(Int_t i) const {return fOverlapLabels[i];}
+     AliTPCclusterMI* GetCurrentCluster() const {return fCurrentCluster;}
+     Int_t GetNoCluster() const {return fNoCluster;}
+     Int_t GetRow() const {return fRow;}
+     Int_t GetCurrentClusterIndex1() const {return fCurrentClusterIndex1;}
+     Bool_t GetInDead() const {return fInDead;}
+     Float_t GetErrorY2() const {return fErrorY2;}
+     Float_t GetErrorZ2() const {return fErrorZ2;}
      //
      //
+     void SetClusterMapBit(int ibit, Bool_t state);
+     Bool_t GetClusterMapBit(int ibit);
+     void SetSharedMapBit(int ibit, Bool_t state);
+     Bool_t GetSharedMapBit(int ibit);
+     const TBits& GetClusterMap() const { return fClusterMap; };
+     const TBits& GetSharedMap() const { return fSharedMap; };
+
+     Float_t  CookdEdxNorm(Double_t low=0.05, Double_t up=0.70, Int_t type=0, Int_t i1=0, Int_t i2=159);
+
  private:
      //     AliTPCseed & operator = (const AliTPCseed &)
      //  {::Fatal("= operator","Not Implemented\n");return *this;}
@@ -86,9 +146,9 @@ class AliTPCseed : public AliTPCtrack {
      TClonesArray * fPoints;              //!array with points along the track
      TClonesArray * fEPoints;             //! array with exact points - calculated in special macro not used in tracking
      //---CURRENT VALUES
-     Int_t fRow;                 //!current row number  
-     Int_t fSector;              //!current sector number
-     Int_t fRelativeSector;      //! index of current relative sector
+     Int_t fRow;                 // current row number  
+     Int_t fSector;              // current sector number
+     Int_t fRelativeSector;      // index of current relative sector
      Float_t fCurrentSigmaY2;    //!expected current cluster sigma Y
      Float_t fCurrentSigmaZ2;    //!expected current cluster sigma Z
      Float_t fErrorY2;           //!sigma of current cluster 
@@ -114,7 +174,9 @@ class AliTPCseed : public AliTPCtrack {
      Float_t fMAngular;           // mean angular factor
      Char_t   fCircular;           // indicates curlin track
      AliTPCTrackerPoint  fTrackPoints[160];  //track points - array track points
-     ClassDef(AliTPCseed,1)  
+     TBits   fClusterMap;       // bit is 1 if track has a hit on padrow
+     TBits   fSharedMap;        // bit is 1 if track shares a hit on padrow
+     ClassDef(AliTPCseed,3)  
 };