- fast fix for track matching

diff --git a/PWGGA/GammaConv/AliCaloPhotonCuts.cxx b/PWGGA/GammaConv/AliCaloPhotonCuts.cxx
index 563a784e1857c14cc9adc2ecc55444a4189f7554..9c1d725b0605060ef68732747ee7826bb349ecd7 100644 (file)
--- a/PWGGA/GammaConv/AliCaloPhotonCuts.cxx
+++ b/PWGGA/GammaConv/AliCaloPhotonCuts.cxx
@@ -520,32 +520,32 @@ Bool_t AliCaloPhotonCuts::ClusterQualityCuts(AliVCluster* cluster, AliVEvent *ev
        cutIndex++;
 
        Double_t minR = 999.0;
-       // get the minimum radius of tracks to cluster
-       if(fHistDistanceTrackToClusterBeforeQA || fHistDistanceTrackToClusterAfterQA){
-               Float_t pos[3];
-               cluster->GetPosition(pos);  // Get cluster position
-               
-               TVector3 cp(pos);
-               int NtrMatched = 0;
-               NtrMatched = cluster->GetNTracksMatched();
-               fHistNMatchedTracks->Fill(NtrMatched);
-               //loop over tracks for Jet QA
-               TList *l = event->GetList();
-               TClonesArray *tracks = dynamic_cast<TClonesArray*>(l->FindObject("Tracks"));
-               for(int itrack = 0; itrack < NtrMatched; itrack++){
-                       AliVTrack *trackcluster = static_cast<AliVTrack*>(tracks->At(itrack));
-                       if (! trackcluster) {
-                               AliError(Form("Couldn't get ESD track %d\n", itrack));
-                               continue;
-                       }
-                       Double_t dphi = -999.0;
-                       Double_t deta = -999.0;
-                       AliPicoTrack::GetEtaPhiDiff(trackcluster, cluster, dphi, deta);
-                       Double_t dr = sqrt(dphi*dphi + deta+deta);
-                       if(dr < minR)
-                               minR = dr;
-               }//loop over tracks
-       }
+//     // get the minimum radius of tracks to cluster
+//     if(fHistDistanceTrackToClusterBeforeQA || fHistDistanceTrackToClusterAfterQA){
+//             Float_t pos[3];
+//             cluster->GetPosition(pos);  // Get cluster position
+//             
+//             TVector3 cp(pos);
+//             int NtrMatched = 0;
+//             NtrMatched = cluster->GetNTracksMatched();
+//             fHistNMatchedTracks->Fill(NtrMatched);
+//             //loop over tracks for Jet QA
+//             TList *l = event->GetList();
+//             TClonesArray *tracks = dynamic_cast<TClonesArray*>(l->FindObject("Tracks"));
+//             for(int itrack = 0; itrack < NtrMatched; itrack++){
+//                     AliVTrack *trackcluster = static_cast<AliVTrack*>(tracks->At(itrack));
+//                     if (! trackcluster) {
+//                             AliError(Form("Couldn't get ESD track %d\n", itrack));
+//                             continue;
+//                     }
+//                     Double_t dphi = -999.0;
+//                     Double_t deta = -999.0;
+//                     AliPicoTrack::GetEtaPhiDiff(trackcluster, cluster, dphi, deta);
+//                     Double_t dr = sqrt(dphi*dphi + deta+deta);
+//                     if(dr < minR)
+//                             minR = dr;
+//             }//loop over tracks
+//     }
        
        // Fill Histos before Cuts
        if(fHistClusterTimevsEBeforeQA) fHistClusterTimevsEBeforeQA->Fill(cluster->GetTOF(), cluster->E());
@@ -569,38 +569,37 @@ Bool_t AliCaloPhotonCuts::ClusterQualityCuts(AliVCluster* cluster, AliVEvent *ev
 
        // Minimum distance to track
        if (fUseDistTrackToCluster){
-               Float_t pos[3];
-               cluster->GetPosition(pos);  // Get cluster position
-               TVector3 cp(pos);
-               int NtrMatched = 0;
-               NtrMatched = cluster->GetNTracksMatched();
-               
-               //loop over tracks for Jet QA
-               TList *l = event->GetList();
-               TClonesArray *tracks = dynamic_cast<TClonesArray*>(l->FindObject("Tracks"));
-
-               for(int itrack = 0; itrack < NtrMatched; itrack++){
-               AliVTrack *trackcluster = static_cast<AliVTrack*>(tracks->At(itrack));
-               if (! trackcluster) {
-                       AliError(Form("Couldn't get ESD track %d\n", itrack));
-                       continue;
-               }
-               Double_t dphi = -999.0;
-               Double_t deta = -999.0;
-               AliPicoTrack::GetEtaPhiDiff(trackcluster, cluster, dphi, deta);
-               Double_t dr = sqrt(dphi*dphi + deta+deta);
-               if(dr < fMinDistTrackToCluster){
-                       return kFALSE;
-               }
-               
-               }//loop over tracks
+//             Float_t pos[3];
+//             cluster->GetPosition(pos);  // Get cluster position
+//             TVector3 cp(pos);
+//             int NtrMatched = 0;
+//             NtrMatched = cluster->GetNTracksMatched();
+//             
+//             //loop over tracks for Jet QA
+//             TList *l = event->GetList();
+//             TClonesArray *tracks = dynamic_cast<TClonesArray*>(l->FindObject("Tracks"));
+// 
+//             for(int itrack = 0; itrack < NtrMatched; itrack++){
+//             AliVTrack *trackcluster = static_cast<AliVTrack*>(tracks->At(itrack));
+//             if (! trackcluster) {
+//                     AliError(Form("Couldn't get ESD track %d\n", itrack));
+//                     continue;
+//             }
+//             Double_t dphi = -999.0;
+//             Double_t deta = -999.0;
+//             AliPicoTrack::GetEtaPhiDiff(trackcluster, cluster, dphi, deta);
+//             Double_t dr = sqrt(dphi*dphi + deta+deta);
+//             if(dr < fMinDistTrackToCluster){
+//                     return kFALSE;
+//             }
+//             
+//             }//loop over tracks
 
-       /*
+       
                if(cluster->GetEmcCpvDistance() < fMinDistTrackToCluster){
-                               if(fHistClusterIdentificationCuts)fHistClusterIdentificationCuts->Fill(cutIndex); //2
-                               return kFALSE;
-                       }
-       */
+                       if(fHistClusterIdentificationCuts)fHistClusterIdentificationCuts->Fill(cutIndex); //2
+                       return kFALSE;
+               }
        }       
        cutIndex++;//3, next cut