Version where the process for HLT and vdrift on DAQ are off(Raphaelle)

[u/mrichter/AliRoot.git] / TRD / AliTRDtrackV1.h

diff --git a/TRD/AliTRDtrackV1.h b/TRD/AliTRDtrackV1.h
index 88db8c3ddce70fc0ae35f2716a1504a30f1b730e..2065491c7ad6caa3914a589839fac92e553cc5ca 100644 (file)
--- a/TRD/AliTRDtrackV1.h
+++ b/TRD/AliTRDtrackV1.h
@@ -19,6 +19,8 @@
 #include "AliTRDseedV1.h"
 #endif
 
+template <typename Value> class TVectorT;
+typedef class TVectorT<Double_t> TVectorD;
 class AliESDtrack;
 class AliTRDcluster;
 class AliTRDReconstructor;
@@ -32,8 +34,6 @@ public:
    ,kNplane    = AliTRDgeometry::kNlayer
    ,kNcham     = AliTRDgeometry::kNstack
    ,kNsect     = AliTRDgeometry::kNsector
-   ,kNslice    =   3
-   ,kNMLPslice =   8 
    ,kMAXCLUSTERSPERTRACK = 210
   };
   
@@ -42,7 +42,8 @@ public:
     kOwner     = BIT(14)
    ,kStopped   = BIT(15) 
    ,kKink      = BIT(16) 
-   ,kPrimary   = BIT(17) 
+   ,kPrimary   = BIT(17)
+   ,kSeeder    = BIT(18)   // if set other than TPC (ITS, TRD)
   };
 
   // propagation/update error codes (up to 4 bits)
@@ -74,9 +75,11 @@ public:
   virtual ~AliTRDtrackV1();
   AliTRDtrackV1 &operator=(const AliTRDtrackV1 &ref);
   virtual void   Copy(TObject &ref) const;
- 
+  virtual TObject* Clone(const char* newname = "") const;
   Bool_t         CookPID();
-  Bool_t         CookLabel(Float_t wrong);
+  Double_t       CookTruncatedMean(const Bool_t kinvq, const Double_t mag, const Int_t charge, const Int_t kcalib, Int_t &nch, Int_t &ncls, TVectorD *Qs=NULL, TVectorD *Xs=NULL, Int_t timeBin0=-1, Int_t timeBin1=1000, Int_t tstep=1) const;
+ 
+  Int_t          CookLabel(Float_t wrong, Int_t *labs=NULL, Float_t *freq=NULL);
   AliTRDtrackV1* GetBackupTrack() const {return fBackupTrack;}
   Double_t       GetBudget(Int_t i) const { return fBudget[i];}
   AliTRDcluster* GetCluster(Int_t id);
@@ -89,6 +92,7 @@ public:
   Double_t       GetPIDsignal() const   { return 0.;}
   Double_t       GetPID(Int_t is) const { return (is >=0 && is < AliPID::kSPECIES) ? fPID[is] : -1.;}
   UChar_t        GetNumberOfTrackletsPID() const;
+  Int_t          GetNumberOfPhysicsSlices() const { return fNdEdxSlices;  };
   Double_t       GetPredictedChi2(const AliTRDseedV1 *tracklet, Double_t *cov) const;
   Double_t       GetPredictedChi2(const AliCluster* /*c*/) const                   { return 0.0; }
   Int_t          GetProlongation(Double_t xk, Double_t &y, Double_t &z) const;
@@ -101,13 +105,13 @@ public:
   AliExternalTrackParam*
                  GetTrackOut() const  { return fTrackHigh;} 
   const Int_t* GetTrackletIndexes() const { return &fTrackletIndex[0];}
-  
   Bool_t         IsEqual(const TObject *inTrack) const;
-  Bool_t         IsKink() const    { return TestBit(kKink);}
-  Bool_t         IsOwner() const   { return TestBit(kOwner);};
-  Bool_t         IsPrimary() const   { return TestBit(kPrimary);};
-  Bool_t         IsStopped() const { return TestBit(kStopped);};
+  Bool_t         IsKink() const           { return TestBit(kKink);}
+  Bool_t         IsOwner() const          { return TestBit(kOwner);};
+  Bool_t         IsPrimary() const        { return TestBit(kPrimary);};
+  Bool_t         IsStopped() const        { return TestBit(kStopped);};
   Bool_t         IsElectron() const;
+  Bool_t         IsTPCseeded() const      { return !TestBit(kSeeder);};
   inline static Bool_t IsTrackError(ETRDtrackError error, UInt_t status);
   inline static Bool_t IsLayerError(ETRDlayerError error, Int_t layer, UInt_t status);
 
@@ -122,11 +126,12 @@ public:
   void           SetESDid(Int_t id) {fESDid = id;}
   void           SetKink(Bool_t k)        { SetBit(kKink, k);}
   void           SetPrimary(Bool_t k)     { SetBit(kPrimary, k);}
+  void           SetNonTPCseeded(Bool_t k){ SetBit(kSeeder, k);}
   void           SetNumberOfClusters();
   void           SetOwner();
   void           SetPID(Short_t is, Double_t inPID){if (is >=0 && is < AliPID::kSPECIES) fPID[is]=inPID;};
   void           SetPIDquality(UChar_t /*inPIDquality*/) const {/*fPIDquality = inPIDquality*/;};
-  inline void    SetStatus(UChar_t stat, Int_t ly=-1);
+  inline void    SetErrStat(UChar_t stat, Int_t ly=-1);
   void           SetStopped(Bool_t stop) {SetBit(kStopped, stop);}
   void           SetTracklet(AliTRDseedV1 *const trklt,  Int_t index);
   void           SetTrackIn();
@@ -145,13 +150,18 @@ private:
   Double32_t   fPID[AliPID::kSPECIES]; //  PID probabilities
   Double32_t   fBudget[3];             //  Integrated material budget
   Double32_t   fDE;                    //  Integrated delta energy
+  Double32_t   fTruncatedMean;         //  Truncated mean
+  Int_t        fNchamberdEdx;          //  number of chambers used in calculating truncated mean
+  Int_t        fNclusterdEdx;          //  number of clusters used in calculating truncated mean
+  Int_t        fNdEdxSlices;           //  number of "physics slices" fill via AliTRDPIDResponse
+
   const AliTRDReconstructor *fkReconstructor;//! reconstructor link 
   AliTRDtrackV1 *fBackupTrack;         //! Backup track
   AliTRDseedV1  *fTracklet[kNplane];   //  Tracklets array defining the track
   AliExternalTrackParam *fTrackLow;    // parameters of the track which enter TRD from below (TPC) 
   AliExternalTrackParam *fTrackHigh;   // parameters of the track which enter TRD from above (HMPID, PHOS) 
 
-  ClassDef(AliTRDtrackV1, 7)          // TRD track - tracklet based
+  ClassDef(AliTRDtrackV1, 8)           // TRD track - tracklet based
 };
 
 //____________________________________________________
@@ -222,7 +232,7 @@ inline void AliTRDtrackV1::SetReconstructor(const AliTRDReconstructor *rec)
 }
 
 //____________________________________________________
-inline void AliTRDtrackV1::SetStatus(UChar_t status, Int_t ly)
+inline void AliTRDtrackV1::SetErrStat(UChar_t status, Int_t ly)
 {
   if(ly<kNplane) fStatus|=((status&0xf)<<((ly+1)*4));
   return;