2 extra members for TRD dEdX (Xianguo, Marian)
authorshahoian <shahoian@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 30 Apr 2012 15:38:11 +0000 (15:38 +0000)
committershahoian <shahoian@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 30 Apr 2012 15:38:11 +0000 (15:38 +0000)
STEER/ESD/AliESDtrack.cxx
STEER/ESD/AliESDtrack.h

index 915997d..3aa9b6a 100644 (file)
@@ -229,6 +229,8 @@ AliESDtrack::AliESDtrack() :
   fTRDncls(0),
   fTRDncls0(0),
   fTRDntracklets(0),
+  fTRDNchamberdEdx(0),
+  fTRDNclusterdEdx(0),
   fTRDnSlices(0),
   fTRDslices(0x0),
   fVertexID(-2),// -2 means an orphan track 
@@ -341,6 +343,8 @@ AliESDtrack::AliESDtrack(const AliESDtrack& track):
   fTRDncls(track.fTRDncls),
   fTRDncls0(track.fTRDncls0),
   fTRDntracklets(track.fTRDntracklets),
+  fTRDNchamberdEdx(track.fTRDNchamberdEdx),
+  fTRDNclusterdEdx(track.fTRDNclusterdEdx),
   fTRDnSlices(track.fTRDnSlices),
   fTRDslices(0x0),
   fVertexID(track.fVertexID),
@@ -465,6 +469,8 @@ AliESDtrack::AliESDtrack(const AliVTrack *track) :
   fTRDncls(0),
   fTRDncls0(0),
   fTRDntracklets(0),
+  fTRDNchamberdEdx(0),
+  fTRDNclusterdEdx(0),
   fTRDnSlices(0),
   fTRDslices(0x0),
   fVertexID(-2),  // -2 means an orphan track
@@ -609,6 +615,8 @@ AliESDtrack::AliESDtrack(TParticle * part) :
   fTRDncls(0),
   fTRDncls0(0),
   fTRDntracklets(0),
+  fTRDNchamberdEdx(0),
+  fTRDNclusterdEdx(0),
   fTRDnSlices(0),
   fTRDslices(0x0),
   fVertexID(-2),  // -2 means an orphan track
@@ -913,6 +921,8 @@ AliESDtrack &AliESDtrack::operator=(const AliESDtrack &source){
     fTPCPoints[i] = source.fTPCPoints[i];  
   }
   fTRDsignal = source.fTRDsignal;
+  fTRDNchamberdEdx = source.fTRDNchamberdEdx;
+  fTRDNclusterdEdx = source.fTRDNclusterdEdx;
 
   for(int i = 0;i < kTRDnPlanes;++i){
     fTRDTimBin[i] = source.fTRDTimBin[i];   
@@ -1147,6 +1157,9 @@ void AliESDtrack::MakeMiniESDtrack(){
   fTRDncls = 0;       
   fTRDncls0 = 0;       
   fTRDsignal = 0;      
+  fTRDNchamberdEdx = 0;
+  fTRDNclusterdEdx = 0;
+
   for (Int_t i=0;i<kTRDnPlanes;i++) {
     fTRDTimBin[i]  = 0;
   }
@@ -1413,7 +1426,8 @@ Bool_t AliESDtrack::UpdateTrackParams(const AliKalmanTrack *t, ULong_t flags){
       delete [] indexTRD;
     }    
     
-    fTRDsignal=t->GetPIDsignal();
+    //commented out by Xianguo
+    //fTRDsignal=t->GetPIDsignal();
     }
     break;
   case kTRDbackup:
@@ -2415,6 +2429,8 @@ void AliESDtrack::Print(Option_t *) const {
     for(index = 0 ; index < AliPID::kSPECIES; index++) 
       printf("%f, ", p[index]) ;
       printf("\n           signal = %f\n", GetTRDsignal()) ;
+      printf("\n           NchamberdEdx = %d\n", GetTRDNchamberdEdx()) ;
+      printf("\n           NclusterdEdx = %d\n", GetTRDNclusterdEdx()) ;
   }
   if( IsOn(kTOFpid) ){
     printf("From TOF: ") ; 
index c75d647..efd80cf 100644 (file)
@@ -245,6 +245,8 @@ public:
   
   void    SetTRDpid(const Double_t *p);
   void    SetTRDsignal(Double_t sig) {fTRDsignal = sig;}
+  void    SetTRDNchamberdEdx(UChar_t nch) {fTRDNchamberdEdx = nch;}
+  void    SetTRDNclusterdEdx(UChar_t ncls){fTRDNclusterdEdx = ncls;}
          
 // A.Bercuci
   void    SetTRDntracklets(UChar_t q){fTRDntracklets = q;}
@@ -273,7 +275,8 @@ public:
   void    SetTRDTimBin(Int_t timbin, Int_t i) {fTRDTimBin[i]=timbin;}
   void    GetTRDpid(Double_t *p) const;
   Double_t GetTRDsignal() const {return fTRDsignal;}
-
+  UChar_t GetTRDNchamberdEdx() const {return fTRDNchamberdEdx;}
+  UChar_t GetTRDNclusterdEdx() const {return fTRDNclusterdEdx;}
   Char_t   GetTRDTimBin(Int_t i) const {return fTRDTimBin[i];}
   Double_t GetTRDchi2() const {return fTRDchi2;}
   UChar_t   GetTRDclusters(Int_t *idx) const;
@@ -512,6 +515,8 @@ protected:
   UChar_t fTRDncls;        // number of clusters assigned in the TRD
   UChar_t fTRDncls0;       // number of clusters assigned in the TRD before first material cross
   UChar_t fTRDntracklets;  // number of TRD tracklets used for tracking/PID
+  UChar_t fTRDNchamberdEdx;   // number of chambers used to calculated the TRD truncated mean
+  UChar_t fTRDNclusterdEdx;   // number of clusters used to calculated the TRD truncated mean
 
   Int_t fTRDnSlices;     // number of slices used for PID in the TRD
   Double32_t *fTRDslices;  //[fTRDnSlices] 
@@ -531,7 +536,7 @@ protected:
   static bool fgkOnlineMode; //! indicate the online mode to skip some of the functionality
 
   AliESDtrack & operator=(const AliESDtrack & );
-  ClassDef(AliESDtrack,64)  //ESDtrack 
+  ClassDef(AliESDtrack,65)  //ESDtrack 
 };