[PWGLF] reject the same event

PWGLF/Tasks/Strangeness/lambdalambda.cxx

@@ -1,4 +1,4 @@
 // Copyright 2019-2020 CERN and copyright holders of ALICE O2.
 // See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
 // All rights not expressly granted are reserved.
 //
@@ -9,7 +9,7 @@
 // granted to it by virtue of its status as an Intergovernmental Organization
 // or submit itself to any jurisdiction.
 
 /// \author Junlee Kim (jikim1290@gmail.com)
 
 #include <cmath>
 #include <array>
@@ -125,16 +125,16 @@
   Configurable<bool> cfgEffCor{"cfgEffCor", false, "flag to apply efficiency correction"};
   Configurable<std::string> cfgEffCorPath{"cfgEffCorPath", "", "path for pseudo efficiency correction"};
 
-  Configurable<int> cfgNoMixedEvents{"cfgNoMixedEvents", 5, "Number of mixed events per event"};
+  Configurable<int> cfgNoMixedEvents{"cfgNoMixedEvents", 10, "Number of mixed events per event"};
 
   ConfigurableAxis massAxis{"massAxis", {110, 2.22, 2.33}, "Invariant mass axis"};
   ConfigurableAxis ptAxis{"ptAxis", {VARIABLE_WIDTH, 0.2, 0.5, 1.0, 1.5, 2.0, 2.5, 3.0, 4.0, 5.0, 6.5, 8.0, 10.0, 100.0}, "Transverse momentum bins"};
-  ConfigurableAxis centAxis{"centAxis", {VARIABLE_WIDTH, 0, 5, 10, 20, 30, 40, 50, 60, 70, 80, 100}, "Centrality interval"};
+  ConfigurableAxis centAxis{"centAxis", {VARIABLE_WIDTH, 0, 10, 20, 50, 100}, "Centrality interval"};
   ConfigurableAxis vertexAxis{"vertexAxis", {10, -10, 10}, "vertex axis for mixing"};
 
   ConfigurableAxis RadiusAxis{"RadiusAxis", {100, 0, 5}, "radius of v0v0"};
   ConfigurableAxis CPAAxis{"CPAAxis", {102, -1.02, 1.02}, "CPA of v0v0"};
   ConfigurableAxis DistanceAxis{"DistanceAxis", {100, 0, 10}, "distance of v0v0"};
   ConfigurableAxis DCAAxis{"DCAAxis", {100, 0, 5}, "DCA of v0v0"};
 
   TF1* fMultPVCutLow = nullptr;
@@ -325,7 +325,7 @@
     v02mom.SetXYZ(v02.px() / v02.p(), v02.py() / v02.p(), v02.pz() / v02.p());
     ROOT::Math::XYZVector posdiff = v02pos - v01pos;
 
     float d = 1. - TMath::Power(v01mom.Dot(v02mom), 2);
     if (d < 1e-5)
       return 999;
     float t = posdiff.Dot(v01mom - v01mom.Dot(v02mom) * v02mom) / d;
@@ -338,7 +338,7 @@
   template <typename C1, typename C2, typename V01, typename V02>
   void FillHistograms(C1 const& c1, C2 const& c2, V01 const& V01s, V02 const& V02s)
   {
     for (auto& v01 : V01s) {
       auto postrack_v01 = v01.template posTrack_as<TrackCandidates>();
       auto negtrack_v01 = v01.template negTrack_as<TrackCandidates>();
 
@@ -371,7 +371,7 @@
         RecoV01 = ROOT::Math::PxPyPzMVector(v01.px(), v01.py(), v01.pz(), v01.mAntiLambda());
       }
 
       for (auto& v02 : V02s) {
         if (v01.v0Id() <= v02.v0Id() && doprocessDataSame)
           continue;
         auto postrack_v02 = v02.template posTrack_as<TrackCandidates>();
@@ -473,7 +473,7 @@
     auto tracksTuple = std::make_tuple(V0s);
     BinningTypeVertexContributor binningOnPositions{{vertexAxis, centAxis}, true};
     SameKindPair<EventCandidates, V0TrackCandidate, BinningTypeVertexContributor> pair{binningOnPositions, cfgNoMixedEvents, -1, collisions, tracksTuple, &cache};
     for (auto& [c1, tracks1, c2, tracks2] : pair) {
       if (cfgCentEst == 1) {
         centrality = c1.centFT0C();
       } else if (cfgCentEst == 2) {
@@ -483,6 +483,9 @@
         continue;
       if (!eventSelected(c2))
         continue;
+      if (c1.bcId() == c2.bcId())
+        continue;
+
       FillHistograms(c1, c2, tracks1, tracks2);
     }
   }