From d7a1155515f6473e31d47feb1d2c09bd2b63aa11 Mon Sep 17 00:00:00 2001
From: hristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Wed, 26 May 2004 14:21:50 +0000
Subject: [PATCH] Adding possibility to change the tracking window in ctgTheta

---
 TPC/AliTPCReconstructor.cxx |  1 +
 TPC/AliTPCReconstructor.h   |  5 +++++
 TPC/AliTPCclustererMI.cxx   |  3 ++-
 TPC/AliTPCtrackerMI.cxx     | 10 ++++++----
 4 files changed, 14 insertions(+), 5 deletions(-)

diff --git a/TPC/AliTPCReconstructor.cxx b/TPC/AliTPCReconstructor.cxx
index 4df397b01b7..8e284e2d290 100644
--- a/TPC/AliTPCReconstructor.cxx
+++ b/TPC/AliTPCReconstructor.cxx
@@ -32,6 +32,7 @@
 
 ClassImp(AliTPCReconstructor)
 
+Double_t AliTPCReconstructor::fgCtgRange = 1.05;
 
 //_____________________________________________________________________________
 void AliTPCReconstructor::Reconstruct(AliRunLoader* runLoader) const
diff --git a/TPC/AliTPCReconstructor.h b/TPC/AliTPCReconstructor.h
index 32695b33e04..babcf6aa9f2 100644
--- a/TPC/AliTPCReconstructor.h
+++ b/TPC/AliTPCReconstructor.h
@@ -19,9 +19,14 @@ public:
   virtual AliTracker*  CreateTracker(AliRunLoader* runLoader) const;
   virtual void         FillESD(AliRunLoader* runLoader, AliESD* esd) const;
 
+  static void SetCtgRange(Double_t ctgRange = 1.05) {fgCtgRange = ctgRange;}
+  static Double_t GetCtgRange(){ return fgCtgRange;}
+
 private:
   AliTPCParam*         GetTPCParam(AliRunLoader* runLoader) const;
 
+  static Double_t fgCtgRange; //! +-fCtgRange is the ctg(Theta) window used for clusterization and tracking (MI) 
+
   ClassDef(AliTPCReconstructor, 0)   // class for the TPC reconstruction
 };
 
diff --git a/TPC/AliTPCclustererMI.cxx b/TPC/AliTPCclustererMI.cxx
index bd714c82e5e..d546627eff1 100644
--- a/TPC/AliTPCclustererMI.cxx
+++ b/TPC/AliTPCclustererMI.cxx
@@ -21,6 +21,7 @@
 //   Origin: Marian Ivanov 
 //-------------------------------------------------------
 
+#include "AliTPCReconstructor.h"
 #include "AliTPCclustererMI.h"
 #include "AliTPCclusterMI.h"
 #include <TObjArray.h>
@@ -707,7 +708,7 @@ void AliTPCclustererMI::FindClusters()
   //first loop - for "gold cluster" 
   fLoop=1;
   Int_t *b=&fBins[-1]+2*fMaxTime;
-  Int_t crtime = Int_t((fParam->GetZLength()-1.05*fRx)/fZWidth-5);
+  Int_t crtime = Int_t((fParam->GetZLength()-AliTPCReconstructor::GetCtgRange()*fRx)/fZWidth-5);
 
   for (Int_t i=2*fMaxTime; i<fMaxBin-2*fMaxTime; i++) {
     b++;
diff --git a/TPC/AliTPCtrackerMI.cxx b/TPC/AliTPCtrackerMI.cxx
index 63f31594422..a8500a27dec 100644
--- a/TPC/AliTPCtrackerMI.cxx
+++ b/TPC/AliTPCtrackerMI.cxx
@@ -49,6 +49,8 @@
 #include "AliTPCreco.h" 
 #include "AliTPCtrackerMI.h"
 #include "TStopwatch.h"
+
+#include "AliTPCReconstructor.h"
 //
 
 ClassImp(AliTPCseed)
@@ -1539,7 +1541,7 @@ Int_t AliTPCtrackerMI::FollowToNext(AliTPCseed& t, Int_t nr) {
   } 
   else
     {
-      if (TMath::Abs(z)<(1.05*x+10)) t.fNFoundable++;
+      if (TMath::Abs(z)<(AliTPCReconstructor::GetCtgRange()*x+10)) t.fNFoundable++;
       else
 	return 0;
     }   
@@ -1636,7 +1638,7 @@ Int_t AliTPCtrackerMI::FollowToNextFast(AliTPCseed& t, Int_t nr) {
   } 
   else
     {
-      if (TMath::Abs(z)>(1.05*x+10)) t.SetClusterIndex2(row,-1);
+      if (TMath::Abs(z)>(AliTPCReconstructor::GetCtgRange()*x+10)) t.SetClusterIndex2(row,-1);
     }   
   //calculate 
   
@@ -1719,7 +1721,7 @@ Int_t AliTPCtrackerMI::UpdateClusters(AliTPCseed& t,  Int_t nr) {
   } 
   else
     {
-      if (TMath::Abs(t.GetZ())<(1.05*t.GetX()+10)) t.fNFoundable++;
+      if (TMath::Abs(t.GetZ())<(AliTPCReconstructor::GetCtgRange()*t.GetX()+10)) t.fNFoundable++;
       else
 	return 0;      
     }
@@ -5264,7 +5266,7 @@ void AliTPCseed::CookdEdx(Double_t low, Double_t up,Int_t i1, Int_t i2, Bool_t o
 	  //  TMath::Abs( Int_t(iz) - iz + 0.5);
 	  //ampc *= 1.15*(1-0.3*dy);
 	  //ampc *= 1.15*(1-0.3*dz);
-	  //	  Float_t zfactor = (1.05-0.0004*TMath::Abs(point->GetCPoint().GetZ()));
+	  //	  Float_t zfactor = (AliTPCReconstructor::GetCtgRange()-0.0004*TMath::Abs(point->GetCPoint().GetZ()));
 	  //ampc               *=zfactor; 
 	}
 	else{ 
-- 
2.31.1