fix error in dEdx slice calculation for LQ method (Alex Wilk)

diff --git a/TRD/AliTRDReconstructor.h b/TRD/AliTRDReconstructor.h
index acfbedf6bce214905e2dd129f8dbc78cf4795ee5..6b3027c6b737774c5cde6a2b8ed344600c993ade 100644 (file)
--- a/TRD/AliTRDReconstructor.h
+++ b/TRD/AliTRDReconstructor.h
@@ -41,6 +41,10 @@ public:
     kXe = 0,
     kAr = 1
   };
+  enum{
+    kNNslices = 8
+   ,kLQslices = 3
+  };
 
   AliTRDReconstructor();
   AliTRDReconstructor(const AliTRDReconstructor &r);
@@ -76,10 +80,6 @@ public:
   void                SetStreamLevel(Int_t level, AliTRDReconstructorTask task= kTracker);
 
 private:
-  enum{
-    kNNslices = 8
-   ,kLQslices = 3
-  };
   UChar_t       fStreamLevel[5];      // stream level for each reconstruction task         
   UInt_t        fSteerParam;          // steering flags
   Double_t      fTCParams[8];         // Tail Cancellation parameters for drift gases 

diff --git a/TRD/AliTRDseedV1.cxx b/TRD/AliTRDseedV1.cxx
index 2bdf4a5e2d83bc0686526cc16db38cd6e70781a5..3495412149f935da213e18ae977dc9cd6d9fd33c 100644 (file)
--- a/TRD/AliTRDseedV1.cxx
+++ b/TRD/AliTRDseedV1.cxx
@@ -226,7 +226,8 @@ void AliTRDseedV1::CookdEdx(Int_t nslices)
     nclusters[slice]++;
   } // End of loop over clusters
 
-  if(fReconstructor->GetPIDMethod() == AliTRDReconstructor::kLQPID){
+  //if(fReconstructor->GetPIDMethod() == AliTRDReconstructor::kLQPID){
+  if(nslices == AliTRDReconstructor::kLQslices){
   // calculate mean charge per slice (only LQ PID)
     for(int is=0; is<nslices; is++){ 
       if(nclusters[is]) fdEdx[is] /= nclusters[is];