From 04475328019b09a25959b498b3cc136f06c12624 Mon Sep 17 00:00:00 2001
From: jklay <jklay@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Tue, 8 Dec 2009 08:30:43 +0000
Subject: [PATCH] added method FindExtrapolationPoint

---
 EMCAL/AliEMCALTracker.cxx | 69 ++++++++++-----------------------------
 EMCAL/AliEMCALTracker.h   |  4 ++-
 2 files changed, 20 insertions(+), 53 deletions(-)

diff --git a/EMCAL/AliEMCALTracker.cxx b/EMCAL/AliEMCALTracker.cxx
index 99ec18eb523..0b1cb78679f 100644
--- a/EMCAL/AliEMCALTracker.cxx
+++ b/EMCAL/AliEMCALTracker.cxx
@@ -945,35 +945,6 @@ Int_t AliEMCALTracker::CreateMatches()
 			}
 		}
 	}
-		
-	/*
-	// loop on clusters and tracks
-	Int_t icBest;
-	Double_t dist, distBest;
-	for (it = 0; it < nTracks; it++) {
-		AliEMCALTrack *track = (AliEMCALTrack*)fTracks->At(it);
-		if (!track) continue;
-		icBest = -1;
-		distBest = fMaxDist;
-		for (ic = 0; ic < nClusters; ic++) {
-			AliEMCALMatchCluster *cluster = (AliEMCALMatchCluster*)fClusters->At(ic);
-			if (!cluster) continue;
-			dist = CheckPair(track, cluster);
-			if (dist < distBest) {
-				distBest = dist;
-				icBest = ic;
-			}
-		}
-		if (icBest >= 0) {
-			track->SetMatchedClusterIndex(icBest);
-			track->SetMatchedClusterDist(distBest);
-			count++;
-		}
-		else {
-			track->SetMatchedClusterIndex(-1);
-		}
-	}
-	*/
 	
 	return count;
 }
@@ -1022,29 +993,6 @@ Int_t AliEMCALTracker::SolveCompetitions()
 		}
 	}
 	
-	/*
-	Int_t it1, it2, nTracks = (Int_t)fTracks->GetEntries();
-	AliEMCALTrack *track1 = 0, *track2 = 0;
-	for (it1 = 0; it1 < nTracks; it1++) {
-		track1 = (AliEMCALTrack*)fTracks->At(it1);
-		if (!track1) continue;
-		if (track1->GetMatchedClusterIndex() < 0) continue;
-		for (it2 = it1+1; it2 < nTracks; it2++) {
-			track2 = (AliEMCALTrack*)fTracks->At(it2);
-			if (!track2) continue;
-			if (track2->GetMatchedClusterIndex() < 0) continue;
-			if (track1->GetMatchedClusterIndex() != track2->GetMatchedClusterIndex()) continue;
-			count++;
-			if (track1->GetMatchedClusterDist() < track2->GetMatchedClusterDist()) {
-				track2->SetMatchedClusterIndex(-1);
-			}
-			else if (track2->GetMatchedClusterDist() < track1->GetMatchedClusterDist()) {
-				track1->SetMatchedClusterIndex(-1);
-			}
-		}
-	}
-	*/
-	
 	delete [] usedC;
 	delete [] usedT;
 
@@ -1063,6 +1011,23 @@ void AliEMCALTracker::UnloadClusters()
 	
   	Clear();
 }
+//
+//------------------------------------------------------------------------------
+//
+TVector3 AliEMCALTracker::FindExtrapolationPoint(Double_t x,Double_t y,Double_t z, AliESDtrack *track)
+{
+  //Method to determine extrapolation point of track at location x,y,z
+  AliEMCALTrack *tr = new AliEMCALTrack(*track);
+  TVector3 error(-100.,-100.,-100.);
+  if (!tr->PropagateToGlobal(x,y,z, 0.0, 0.0)) {
+    return error;
+  }
+  Double_t pos[3];
+  tr->GetXYZ(pos);
+  TVector3 ExTrPos(pos[0],pos[1],pos[2]);
+  return ExTrPos;
+}
+
 //
 //------------------------------------------------------------------------------
 //
diff --git a/EMCAL/AliEMCALTracker.h b/EMCAL/AliEMCALTracker.h
index 87934cb4c4a..7d3f16b1b5c 100644
--- a/EMCAL/AliEMCALTracker.h
+++ b/EMCAL/AliEMCALTracker.h
@@ -26,6 +26,7 @@
 #include "TMath.h"
 class TList;
 class TTree;
+class TVector3;
 class TObjArray;
 class AliESDEvent;
 class AliESDCaloCluster;
@@ -66,7 +67,8 @@ public:
 	void                SetCutNTPC(Double_t value) {fCutNTPC=value;}
 	void                SetNumberOfSteps(Int_t n) {fNPropSteps=n;if(!n)SetTrackCorrectionMode("NONE");}
 	void                SetTrackCorrectionMode(Option_t *option);
-	
+        TVector3            FindExtrapolationPoint(Double_t x, Double_t y, Double_t z, AliESDtrack *track);
+
 	enum {
 		kUnmatched = -99999
 	};
-- 
2.39.3