- histos: filling with THnBase

[u/mrichter/AliRoot.git] / PWGDQ / dielectron / AliDielectronV0Cuts.cxx

diff --git a/PWGDQ/dielectron/AliDielectronV0Cuts.cxx b/PWGDQ/dielectron/AliDielectronV0Cuts.cxx
index d5dacf40214300a8f2a0d01b9543269d3914a3fc..c09b5614f144e65fce138d33c9dfad62566880c7 100644 (file)
--- a/PWGDQ/dielectron/AliDielectronV0Cuts.cxx
+++ b/PWGDQ/dielectron/AliDielectronV0Cuts.cxx
@@ -109,11 +109,12 @@ void AliDielectronV0Cuts::InitEvent(AliVTrack *trk)
 
   // basic quality cut, at least one of the V0 daughters has to fullfill
   AliDielectronVarCuts dauQAcuts1;
-  dauQAcuts1.AddCut(AliDielectronVarManager::kPt,           0.5,  1e30);
+  dauQAcuts1.AddCut(AliDielectronVarManager::kPt,           0.3,  1e30);
   dauQAcuts1.AddCut(AliDielectronVarManager::kEta,         -0.9,   0.9);
   dauQAcuts1.AddCut(AliDielectronVarManager::kNclsTPC,     50.0, 160.0);
+  dauQAcuts1.AddCut(AliDielectronVarManager::kTPCchi2Cl,    0.0,   4.0);
   AliDielectronTrackCuts dauQAcuts2;
-  dauQAcuts2.SetRequireITSRefit(kTRUE);
+  //  dauQAcuts2.SetRequireITSRefit(kTRUE);
   dauQAcuts2.SetRequireTPCRefit(kTRUE);
 
   Int_t nV0s = 0;
@@ -139,6 +140,9 @@ void AliDielectronV0Cuts::InitEvent(AliVTrack *trk)
        continue;
       }
 
+      // reject tracks with neative ID
+      if(trNeg->GetID()<0 || trPos->GetID()) continue;
+
       // at least one of the daughter has to pass basic QA cuts
       if(!(dauQAcuts1.IsSelected(trNeg) && dauQAcuts2.IsSelected(trNeg)) ||
         !(dauQAcuts1.IsSelected(trPos) && dauQAcuts2.IsSelected(trPos))  ) continue;
@@ -175,6 +179,9 @@ void AliDielectronV0Cuts::InitEvent(AliVTrack *trk)
        continue;
       }
 
+      // reject tracks with neative ID
+      if(trNeg->GetID()<0 || trPos->GetID()) continue;
+
       // at least one of the daughter has to pass basic QA cuts
       if(!(dauQAcuts1.IsSelected(trNeg) && dauQAcuts2.IsSelected(trNeg)) ||
         !(dauQAcuts1.IsSelected(trPos) && dauQAcuts2.IsSelected(trPos))  ) continue;
@@ -196,7 +203,7 @@ void AliDielectronV0Cuts::InitEvent(AliVTrack *trk)
   else
     return;
 
-  //  printf(" Number of V0s candiates found %d \n",nV0s);
+  printf(" Number of V0s candiates found %d \n",nV0s);
 
 }
 //________________________________________________________________________
@@ -210,7 +217,7 @@ Bool_t AliDielectronV0Cuts::IsSelected(TObject* track)
 
   AliVTrack *vtrack = static_cast<AliVTrack*>(track);
   InitEvent(vtrack);
-  //  printf(" track ID %d selected result %d \n",vtrack->GetID(),(fV0TrackArr.TestBitNumber(vtrack->GetID()))^fExcludeTracks);
+  //printf(" track ID %d selected result %d %d \n",vtrack->GetID(),(fV0TrackArr.TestBitNumber(vtrack->GetID())),fExcludeTracks);
   return ( (fV0TrackArr.TestBitNumber(vtrack->GetID()))^fExcludeTracks );
 }