Adding possibility to change the tracking window in ctgTheta
authorhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 26 May 2004 14:21:50 +0000 (14:21 +0000)
committerhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 26 May 2004 14:21:50 +0000 (14:21 +0000)
TPC/AliTPCReconstructor.cxx
TPC/AliTPCReconstructor.h
TPC/AliTPCclustererMI.cxx
TPC/AliTPCtrackerMI.cxx

index 4df397b..8e284e2 100644 (file)
@@ -32,6 +32,7 @@
 
 ClassImp(AliTPCReconstructor)
 
+Double_t AliTPCReconstructor::fgCtgRange = 1.05;
 
 //_____________________________________________________________________________
 void AliTPCReconstructor::Reconstruct(AliRunLoader* runLoader) const
index 32695b3..babcf6a 100644 (file)
@@ -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
 };
 
index bd714c8..d546627 100644 (file)
@@ -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++;
index 63f3159..a8500a2 100644 (file)
@@ -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{