Change the way how TPC pv veretx is rejected

diff --git a/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraBothEventCuts.cxx b/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraBothEventCuts.cxx
index 2b92657b2553ec4d197bd00d3877aeee9f2be424..6e5e24e005acdd97c54cb6dddd8b0aa562d87043 100644 (file)
--- a/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraBothEventCuts.cxx
+++ b/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraBothEventCuts.cxx
@@ -137,7 +137,8 @@ Bool_t AliSpectraBothEventCuts::IsSelected(AliVEvent * aod,AliSpectraBothTrackCu
 
 
 
-   const AliVVertex * vertex = fAOD->GetPrimaryVertex();//FIXME vertex is recreated             
+   const AliVVertex * vertex = fAOD->GetPrimaryVertex();//FIXME vertex is recreated    
+
   if(vertex)fHistoVtxBefSel->Fill(vertex->GetZ());
   fIsSelected =kFALSE;
   if(CheckVtx() && CheckCentralityCut() && CheckMultiplicityCut() && CheckVtxChi2perNDF())
@@ -194,18 +195,22 @@ Bool_t AliSpectraBothEventCuts::CheckVtx()
       fHistoCuts->Fill(kVtxNoEvent);
       return kFALSE;
     }
-   if(vertex->GetNContributors()<1)
+    if(vertex->GetNContributors()<1)
    {
-      fHistoCuts->Fill(kVtxNoEvent);
+
+      fHistoCuts->Fill(kZeroCont);
       return kFALSE;
 
    }
-   TString tmp(vertex->GetTitle());    
-   if(tmp.Contains("TPC"))
+       
+   TString tmp(vertex->GetTitle());
+   if(tmp.Contains("NoConstraint"))
    {
-       fHistoCuts->Fill(kVtxNoEvent);
-       return kFALSE;
-   }                           
+        fHistoCuts->Fill(kTPCasPV);
+        return kFALSE;
+   }
+
+
  // if (vertex->GetZ() > fVertexCutMin && vertex->GetZ() < fVertexCutMax)
    // {
     //  return kTRUE;

diff --git a/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraBothEventCuts.h b/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraBothEventCuts.h
index 4747f330df4c9a30fefec666bab6ddcfb03ecff4..ad932efa7319ba977512d949347139bf02ec5431 100644 (file)
--- a/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraBothEventCuts.h
+++ b/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraBothEventCuts.h
@@ -22,7 +22,7 @@ class AliSpectraBothTrackCuts;
 class AliSpectraBothEventCuts : public TNamed
 {
  public:
-  enum {  kProcessedEvents = 0,kPhysSelEvents,kAcceptedEvents, kVtxRange, kVtxCentral, kVtxNoEvent, kQVector, kNVtxCuts};
+  enum {  kProcessedEvents = 0,kPhysSelEvents,kAcceptedEvents, kVtxRange, kVtxCentral, kVtxNoEvent, kQVector,kTPCasPV,kZeroCont,kNVtxCuts};
 enum {kDoNotCheckforSDD=0,kwithSDD,kwithoutSDD};       
 
   // Constructors

diff --git a/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraBothHistoManager.cxx b/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraBothHistoManager.cxx
index c0ad82a9e6a283ef6b520ab65de4e952ede46f77..f766e228c418fdd123e14476c148ec88710077b0 100644 (file)
--- a/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraBothHistoManager.cxx
+++ b/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraBothHistoManager.cxx
@@ -317,7 +317,7 @@ Long64_t AliSpectraBothHistoManager::Merge(TCollection* list)
   return count+1;
 }
 
-
+//________________________________________________________________________________
 TH1* AliSpectraBothHistoManager::GetHistogram1D(UInt_t histoType, UInt_t particleType, UInt_t charge) {
   // GetHistogram using particle ID and histogram type
   Int_t baseId = -1;
@@ -377,7 +377,7 @@ TH1* AliSpectraBothHistoManager::GetHistogram1D(UInt_t histoType, UInt_t particl
 
   return GetHistogram(baseId);
 }
-
+//____________________________________________________________________________________________________
 TH2* AliSpectraBothHistoManager::GetHistogram2D(UInt_t histoType, UInt_t particleType, UInt_t charge){
   // returns histo based on ids, casting it to TH2*
   return (TH2*) GetHistogram1D(histoType,particleType,charge);