Update (Zaida)

[u/mrichter/AliRoot.git] / PWG3 / vertexingHF / AliAnalysisTaskSEDStarJets.cxx

diff --git a/PWG3/vertexingHF/AliAnalysisTaskSEDStarJets.cxx b/PWG3/vertexingHF/AliAnalysisTaskSEDStarJets.cxx
index d95fc4c27e17dc2fbfc586a5281a2ec130116a8a..15a388d30759a0d87f1683d33e5629ad104aef9f 100644 (file)
--- a/PWG3/vertexingHF/AliAnalysisTaskSEDStarJets.cxx
+++ b/PWG3/vertexingHF/AliAnalysisTaskSEDStarJets.cxx
@@ -234,9 +234,6 @@ void AliAnalysisTaskSEDStarJets::UserExec(Option_t *)
   }
   
   // Load the event
-  fEvents++;
-  AliDebug(2,Form("Event %d",fEvents));
-  if (fEvents%10000 ==0) AliDebug(2,Form("Event %d",fEvents));
   AliAODEvent* aodEvent = dynamic_cast<AliAODEvent*>(fInputEvent);
   
   TClonesArray *arrayVerticesHF=0;
@@ -263,6 +260,15 @@ void AliAnalysisTaskSEDStarJets::UserExec(Option_t *)
     return;
   }else AliDebug(2, Form("Found %d vertices",arrayVerticesHF->GetEntriesFast()));   
 
+  // fix for temporary bug in ESDfilter 
+  // the AODs with null vertex pointer didn't pass the PhysSel
+  if(!aodEvent->GetPrimaryVertex()) return;
+
+  fEvents++;
+  AliDebug(2,Form("Event %d",fEvents));
+  if (fEvents%10000 ==0) AliDebug(2,Form("Event %d",fEvents));
+
+
   // Simulate a jet triggered sample
   TClonesArray *arrayofJets = (TClonesArray*)aodEvent->GetJets();
   if(aodEvent->GetNJets()<=0) return;
@@ -758,7 +764,6 @@ void AliAnalysisTaskSEDStarJets::Terminate(Option_t*)
   // a query. It always runs on the client, it can be used to present
   // the results graphically or save the results to file.
   
-  Info("Terminate","");
   AliAnalysisTaskSE::Terminate();
   
   AliInfo(Form("Found %i of that MC D*->D0pi(D0->kpi) are in acceptance, in %d events",fCountDStar,fEvents));