From 56b32fd7cc24d855301881554ec385801eeab3a2 Mon Sep 17 00:00:00 2001
From: abercuci <abercuci@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Tue, 30 Sep 2008 15:25:50 +0000
Subject: [PATCH] - new option for selecting 8 dEdx slices in ESD track for TRD
 - set back default 1 expo tail for Xe TC

---
 TRD/AliTRDReconstructor.cxx | 10 ++++++++--
 TRD/AliTRDReconstructor.h   | 10 ++++++----
 TRD/AliTRDrecoParam.cxx     |  2 +-
 TRD/AliTRDrecoParam.h       |  2 +-
 TRD/TRDrecLinkDef.h         |  1 +
 TRD/libTRDrec.pkg           |  1 +
 6 files changed, 18 insertions(+), 8 deletions(-)

diff --git a/TRD/AliTRDReconstructor.cxx b/TRD/AliTRDReconstructor.cxx
index 9847b091e8c..146debd82b8 100644
--- a/TRD/AliTRDReconstructor.cxx
+++ b/TRD/AliTRDReconstructor.cxx
@@ -108,6 +108,7 @@ void AliTRDReconstructor::Init(){
   printf("\tHLT         Tracking  [hlt] : %s\n", fSteerParam&kHLT?"yes":"no");
   printf("\tCosmic Reconstruction [cos] : %s\n", fSteerParam&kCosmic?"yes":"no");
   printf("\tNN PID                 [nn] : %s\n", fSteerParam&kSteerPID?"yes":"no");
+  printf("\t8 dEdx slices in ESD   [8s] : %s\n", fSteerParam&kEightSlices?"yes":"no");
   printf("\tStreaming Levels            : Clusterizer[%d] Tracker[%d] PID[%d]\n", fStreamLevel[kClusterizer], fStreamLevel[kTracker], fStreamLevel[kPID]);
 }
 
@@ -216,11 +217,12 @@ void AliTRDReconstructor::SetOption(Option_t *opt)
 // Default steer param values
 //
 // digits conversion [dc] = false
+// drift gas [ar] = false - do not update the number of exponentials in the TC !
 // write clusters [cw] = true
+// write online tracklets [tw] = false
 // track seeding (stand alone tracking) [sa] = true
 // PID method in reconstruction (NN) [nn] = true
-// write online tracklets [tw] = false
-// drift gas [ar] = false
+// 8 dEdx slices in ESD [8s] = false 
 // HLT tracking [hlt] = false
 // Cosmic Reconstruction [cos] = false
 //
@@ -247,6 +249,10 @@ void AliTRDReconstructor::SetOption(Option_t *opt)
       SETFLG(fSteerParam, kSteerPID);
       if(sopt.Contains("!")) CLRFLG(fSteerParam, kSteerPID);
       continue;
+    } else if(sopt.Contains("8s")){
+      SETFLG(fSteerParam, kEightSlices);
+      if(sopt.Contains("!")) CLRFLG(fSteerParam, kEightSlices);
+      continue;
     } else if(sopt.Contains("tw")){
       SETFLG(fSteerParam, kWriteTracklets);
       if(sopt.Contains("!")) CLRFLG(fSteerParam, kWriteTracklets);
diff --git a/TRD/AliTRDReconstructor.h b/TRD/AliTRDReconstructor.h
index 1860ec353f6..909a1faeae0 100644
--- a/TRD/AliTRDReconstructor.h
+++ b/TRD/AliTRDReconstructor.h
@@ -25,10 +25,11 @@ public:
     ,kWriteClusters  = BIT(1)
     ,kSeeding        = BIT(2)
     ,kSteerPID       = BIT(3)
-    ,kWriteTracklets = BIT(4)
-    ,kDriftGas       = BIT(5)
-    ,kHLT            = BIT(6)
-    ,kCosmic         = BIT(7)
+    ,kEightSlices    = BIT(4)
+    ,kWriteTracklets = BIT(5)
+    ,kDriftGas       = BIT(6)
+    ,kHLT            = BIT(7)
+    ,kCosmic         = BIT(8)
   };
   enum AliTRDReconstructorTask {
     kClusterizer = 0
@@ -73,6 +74,7 @@ public:
   Bool_t              IsHLT() const              { return fSteerParam&kHLT;}
   Bool_t              IsSeeding() const          { return fSteerParam&kSeeding;}
   Bool_t              IsCosmic() const           { return fSteerParam&kCosmic;}
+  Bool_t              IsEightSlices() const      { return fSteerParam&kEightSlices;}
 
   virtual void        Reconstruct(AliRawReader *rawReader, TTree *clusterTree) const;
   virtual void        Reconstruct(TTree *digitsTree, TTree *clusterTree) const;
diff --git a/TRD/AliTRDrecoParam.cxx b/TRD/AliTRDrecoParam.cxx
index 67f5bedff9c..885c4e99870 100644
--- a/TRD/AliTRDrecoParam.cxx
+++ b/TRD/AliTRDrecoParam.cxx
@@ -55,7 +55,7 @@ AliTRDrecoParam::AliTRDrecoParam()
   ,fMinLeftRightCutSigma(8.)
   ,fClusMaxThresh(4.5)
   ,fClusSigThresh(3.5)
-  ,fTCnexp(2)
+  ,fTCnexp(1)
   ,fNumberOfPresamples(0)
   ,fNumberOfPostsamples(0)
 {
diff --git a/TRD/AliTRDrecoParam.h b/TRD/AliTRDrecoParam.h
index 06a487a5541..167c3c060cd 100644
--- a/TRD/AliTRDrecoParam.h
+++ b/TRD/AliTRDrecoParam.h
@@ -36,7 +36,7 @@ public:
   Double_t GetMaxTheta() const              { return fkMaxTheta; }
   Double_t GetMaxPhi() const                { return fkMaxPhi;   }
   Double_t GetPlaneQualityThreshold() const { return fkPlaneQualityThreshold; }
-  Double_t GetPIDThreshold(Float_t /*p*/){ return 0.;}
+  Double_t GetPIDThreshold(Float_t /*p*/) const { return 0.;}
   Double_t GetRoad0y() const                { return fkRoad0y;   }
   Double_t GetRoad0z() const                { return fkRoad0z;   }
   Double_t GetRoad1y() const                { return fkRoad1y;   }
diff --git a/TRD/TRDrecLinkDef.h b/TRD/TRDrecLinkDef.h
index 825912e6a21..6c86055b926 100644
--- a/TRD/TRDrecLinkDef.h
+++ b/TRD/TRDrecLinkDef.h
@@ -19,6 +19,7 @@
 #pragma link C++ class  AliTRDseed+;
 
 #pragma link C++ class  AliTRDpidESD+;
+#pragma link C++ class  AliTRDpidUtil+;
 
 #pragma link C++ class  AliTRDReconstructor+;
 
diff --git a/TRD/libTRDrec.pkg b/TRD/libTRDrec.pkg
index c98818a3718..c8011425692 100644
--- a/TRD/libTRDrec.pkg
+++ b/TRD/libTRDrec.pkg
@@ -8,6 +8,7 @@ SRCS= AliTRDcluster.cxx \
       AliTRDtracker.cxx \
       AliTRDseed.cxx \
       AliTRDpropagationLayer.cxx \
+      AliTRDpidUtil.cxx \
       AliTRDpidESD.cxx \
       AliTRDReconstructor.cxx \
       AliTRDtrackingAnalysis.cxx \
-- 
2.43.0