Bug fix for poly fit (ChiaraB)

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

diff --git a/PWG3/vertexingHF/AliAnalysisTaskSEDStarSpectra.cxx b/PWG3/vertexingHF/AliAnalysisTaskSEDStarSpectra.cxx
index e2626fe2f20c9e7cbf13b544144b81fb1810aeb9..2fa4c59e4048d47d68eb16324c28382fded9091e 100644 (file)
--- a/PWG3/vertexingHF/AliAnalysisTaskSEDStarSpectra.cxx
+++ b/PWG3/vertexingHF/AliAnalysisTaskSEDStarSpectra.cxx
@@ -192,11 +192,6 @@ void AliAnalysisTaskSEDStarSpectra::UserExec(Option_t *)
     return;
   }
   
-  fCEvents->Fill(1);
-  // Load the event
-  fEvents++;
-  AliInfo(Form("Event %d",fEvents));
-  if (fEvents%10000 ==0) AliInfo(Form("Event %d",fEvents));
   AliAODEvent* aodEvent = dynamic_cast<AliAODEvent*>(fInputEvent);
   TClonesArray *arrayDStartoD0pi=0;
  
@@ -216,7 +211,19 @@ void AliAnalysisTaskSEDStarSpectra::UserExec(Option_t *)
   } else {
     arrayDStartoD0pi=(TClonesArray*)aodEvent->GetList()->FindObject("Dstar");
   }
+
+
+  // fix for temporary bug in ESDfilter 
+  // the AODs with null vertex pointer didn't pass the PhysSel
+  if(!aodEvent->GetPrimaryVertex() || TMath::Abs(aodEvent->GetMagneticField())<0.001) return;
  
+
+  fCEvents->Fill(1);
+  // Load the event
+  fEvents++;
+  AliInfo(Form("Event %d",fEvents));
+  if (fEvents%10000 ==0) AliInfo(Form("Event %d",fEvents));
+
   // counters for efficiencies
   Int_t icountReco = 0;
   
@@ -365,7 +372,7 @@ void AliAnalysisTaskSEDStarSpectra::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","");
+  //Info("Terminate","");
   AliAnalysisTaskSE::Terminate();
   
   fOutput = dynamic_cast<TList*> (GetOutputData(1));