In case of isolation only with neutral, do not reject track-matched clusters
authorgconesab <gconesab@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 16 Mar 2012 10:52:03 +0000 (10:52 +0000)
committergconesab <gconesab@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 16 Mar 2012 10:52:03 +0000 (10:52 +0000)
PWG/CaloTrackCorrBase/AliIsolationCut.cxx

index df09da1..8b1eead 100755 (executable)
@@ -142,7 +142,7 @@ Float_t AliIsolationCut::GetCellDensity(const AliAODPWG4ParticleCorrelation * pC
     if (coneCells > 0.) 
     {
       cellDensity = (coneCells-coneCellsBad)/coneCells;
-      printf("Energy density = %f\n", cellDensity);
+      //printf("Energy density = %f\n", cellDensity);
     }
   }
 
@@ -334,8 +334,9 @@ void  AliIsolationCut::MakeIsolationCut(const TObjArray * plCTS,
       if(calo->GetID() == pCandidate->GetCaloLabel(0) || 
          calo->GetID() == pCandidate->GetCaloLabel(1)   ) continue ;      
       
-      //Skip matched clusters with tracks
-      if( pid->IsTrackMatched(calo,reader->GetCaloUtils(),reader->GetInputEvent()) ) continue ;
+      //Skip matched clusters with tracks in case of neutral+charged analysis
+      if( fPartInCone == kNeutralAndCharged && 
+          pid->IsTrackMatched(calo,reader->GetCaloUtils(),reader->GetInputEvent()) ) continue ;
     
       //Assume that come from vertex in straight line
       calo->GetMomentum(mom,reader->GetVertex(evtIndex)) ;