- HF can take now all kind of histograms

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

diff --git a/PWGDQ/dielectron/AliDielectronV0Cuts.cxx b/PWGDQ/dielectron/AliDielectronV0Cuts.cxx
index 0fb12a675ab2a3c7696ad5883538008fa294ce17..9b785e8cba4d332f4fc7a56b4a11e4c849fd1286 100644 (file)
--- a/PWGDQ/dielectron/AliDielectronV0Cuts.cxx
+++ b/PWGDQ/dielectron/AliDielectronV0Cuts.cxx
@@ -112,9 +112,10 @@ void AliDielectronV0Cuts::InitEvent(AliVTrack *trk)
 
   // basic quality cut, at least one of the V0 daughters has to fullfill
   AliDielectronVarCuts dauQAcuts1;
-  dauQAcuts1.AddCut(AliDielectronVarManager::kEta,          -0.9,   0.9);
-  dauQAcuts1.AddCut(AliDielectronVarManager::kNclsTPC,      50.0, 160.0);
-  dauQAcuts1.AddCut(AliDielectronVarManager::kTPCchi2Cl,     0.0,   4.0);
+  dauQAcuts1.AddCut(AliDielectronVarManager::kPt,            0.05, 100.0);
+  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.SetRequireTPCRefit(kTRUE);
@@ -129,7 +130,7 @@ void AliDielectronV0Cuts::InitEvent(AliVTrack *trk)
   if(ev->IsA() == AliESDEvent::Class()) {
     const AliESDEvent *esdev = static_cast<const AliESDEvent*>(ev);
 
-    printf("there are %d V0s in the event \n",esdev->GetNumberOfV0s());
+    //printf("there are %d V0s in the event \n",esdev->GetNumberOfV0s());
     // loop over V0s
     for (Int_t iv=0; iv<esdev->GetNumberOfV0s(); ++iv){
       AliESDv0 *v = esdev->GetV0(iv);
@@ -147,12 +148,15 @@ void AliDielectronV0Cuts::InitEvent(AliVTrack *trk)
 
       // PID default cuts
       if(fPID>=0) {
-       if( !dauPIDcuts.IsSelected(trNeg) || !dauPIDcuts.IsSelected(trPos) ) continue;
+       if( !dauPIDcuts.IsSelected(trNeg) ) continue;
+       if( !dauPIDcuts.IsSelected(trPos) ) 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;
+      // basic track cuts
+      if( !dauQAcuts2.IsSelected(trNeg) ) continue;
+      if( !dauQAcuts2.IsSelected(trPos) ) continue;
+      if( !dauQAcuts1.IsSelected(trNeg) ) continue;
+      if( !dauQAcuts1.IsSelected(trPos) ) continue;
 
       if(fMotherPdg==22) candidate.SetGammaTracks(trNeg, 11, trPos, 11);
       else candidate.SetTracks(trNeg, fNegPdg, trPos, fPosPdg);
@@ -190,12 +194,15 @@ void AliDielectronV0Cuts::InitEvent(AliVTrack *trk)
 
       // PID default cuts
       if(fPID>=0) {
-       if( !dauPIDcuts.IsSelected(trNeg) || !dauPIDcuts.IsSelected(trPos) ) continue;
+       if( !dauPIDcuts.IsSelected(trNeg) ) continue;
+       if( !dauPIDcuts.IsSelected(trPos) ) 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;
+      // basic track cuts
+      if( !dauQAcuts2.IsSelected(trNeg) ) continue;
+      if( !dauQAcuts2.IsSelected(trPos) ) continue;
+      if( !dauQAcuts1.IsSelected(trNeg) ) continue;
+      if( !dauQAcuts1.IsSelected(trPos) ) continue;
 
       AliKFVertex v0vtx = *v;
       if(fMotherPdg==22) candidate.SetGammaTracks(trNeg, 11, trPos, 11);