From: schutz <schutz@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Mon, 9 Feb 2004 17:06:36 +0000 (+0000)
Subject: updating
X-Git-Url: http://git.uio.no/git/?a=commitdiff_plain;h=80427b5a40b560ef66389384410a278376b71676;p=u%2Fmrichter%2FAliRoot.git

updating
---

diff --git a/PHOS/AnaESD.C b/PHOS/AnaESD.C
index 4eae6a19850..623e5fac4b5 100644
--- a/PHOS/AnaESD.C
+++ b/PHOS/AnaESD.C
@@ -74,17 +74,20 @@ void Ana()
     AliESDtrack * cp ; 
     for (index = 0 ; index < esd->GetNumberOfTracks() ; index++) {
       cp = esd->GetTrack(index) ;
-      ULong_t status = cp->GetStatus() ; 
-
-      // check if the tracks comes out of TRD
-      if ((status & AliESDtrack::kTRDout)==0) 
-	continue;
-      if ((status & AliESDtrack::kTRDStop)!=0) 
+      Double_t xyzAtPHOS[3] ; 
+      cp->GetOuterXYZ(xyzAtPHOS) ; 
+      // check if the track reaches PHOS
+      if ( (xyzAtPHOS[0] +  xyzAtPHOS[1] + xyzAtPHOS[2]) == 0.)
 	continue;
+      // the next check are only if we want high quality tracks 
+      //       ULong_t status = cp->GetStatus() ;  
+      //       if ((status & AliESDtrack::kTRDput)==0) 
+      // 	continue;
+      //       if ((status & AliESDtrack::kTRDStop)!=0) 
+      // 	continue;
 
       // Gets the Global coordinate of the track at the entrance of PHOS 
-      Double_t xyzAtPHOS[3] ; 
-      cp->GetOuterXYZ(xyzAtPHOS) ; 
+     
       TVector3 poscp(xyzAtPHOS[0], xyzAtPHOS[1], xyzAtPHOS[2]) ; 
       cout << "Charged particle # " << index << " pos (" 
 	   << poscp.X() << ", " << poscp.Y() << ", " <<poscp.Z() << ") : (" << poscp.Eta()