[PWGLF] Improve mixing buffer for kinematic mixing, best pair compatibilaty for doublephi analysis (#17203)

Author: skundu692

PWGLF/TableProducer/Resonances/doublephitable.cxx

@@ -720,6 +720,12 @@ struct doublephitable {
         const auto phi2Fit = k3Fit + k4Fit;
         const auto pairFit = phi1Fit + phi2Fit;
 
+        if (pairFit.Pt() < 6.0) {
+          continue;
+        }
+        if (pairFit.M() < 2.5 || pairFit.M() > 3.2) {
+          continue;
+        }
         // No pair-pT or pair-mass cut here.
 
         PhiPhiPairPayload pair;

PWGLF/Tasks/Resonances/doublephimeson.cxx

@@ -67,13 +67,14 @@ struct doublephimeson {
   Configurable<bool> additionalEvsel{"additionalEvsel", true, "Additional event selection"};
   Configurable<bool> isDeep{"isDeep", true, "Store deep angle"};
   Configurable<float> cutMinNsigmaTPC{"cutMinNsigmaTPC", -2.5, "nsigma cut TPC"};
-  Configurable<float> cutNsigmaTPC{"cutNsigmaTPC", 3.0, "nsigma cut TPC"};
-  Configurable<float> cutNsigmaTOF{"cutNsigmaTOF", 3.0, "nsigma cut TOF"};
+  Configurable<float> cutNsigmaTPC{"cutNsigmaTPC", 2.5, "nsigma cut TPC"};
+  Configurable<float> cutNsigmaTOF{"cutNsigmaTOF", 2.5, "nsigma cut TOF"};
   Configurable<float> momTOFCut{"momTOFCut", 1.8, "minimum pT cut for madnatory TOF"};
   Configurable<float> maxKaonPt{"maxKaonPt", 100.0, "maximum kaon pt cut"};
   Configurable<float> cfgCrossPhiLow{"cfgCrossPhiLow", 1.01, "Lower edge of phi mass window for cross-pairing (ghost) veto"};
   Configurable<float> cfgCrossPhiHigh{"cfgCrossPhiHigh", 1.03, "Upper edge of phi mass window for cross-pairing (ghost) veto"};
   Configurable<bool> useParametrized{"useParametrized", false, "Use pT dependent mass peak and width"};
+  Configurable<bool> useCrossPairRejection{"useCrossPairRejection", true, "Use cross pair phi signal compatibilaty"};
   // ------------------------------------------------------------
   // pT-dependent phi mass peak and width from single-phi BW fits
   //
@@ -337,17 +338,17 @@ struct doublephimeson {
     }
 
     if (PIDStrategy == 1003) {
-      if (ptcand < 0.5 && TOFHit != 1 && nsigmaTPC > -2.0 && nsigmaTPC < 2.0) {
+      if (ptcand < 0.5 && TOFHit != 1 && std::abs(nsigmaTPC) < cutNsigmaTPC) {
         return true;
       }
-      if (ptcand < 0.5 && TOFHit == 1 && std::sqrt(nsigmaTOF * nsigmaTOF + nsigmaTPC * nsigmaTPC) < 2.5) {
+      if (ptcand < 0.5 && TOFHit == 1 && std::sqrt(nsigmaTOF * nsigmaTOF + nsigmaTPC * nsigmaTPC) < cutNsigmaTPC) {
         return true;
       }
       if (ptcand >= 0.5) {
-        if (TOFHit != 1 && nsigmaTPC > -2.0 && nsigmaTPC < 2.0) {
+        if (TOFHit != 1 && nsigmaTPC > -2.0 && nsigmaTPC < cutNsigmaTPC) {
           return true;
         }
-        if (TOFHit == 1 && std::sqrt(nsigmaTOF * nsigmaTOF + nsigmaTPC * nsigmaTPC) < 2.5) {
+        if (TOFHit == 1 && std::sqrt(nsigmaTOF * nsigmaTOF + nsigmaTPC * nsigmaTPC) < cutNsigmaTOF) {
           return true;
         }
       }
@@ -1695,7 +1696,6 @@ struct doublephimeson {
         if (id2 <= id1) {
           continue;
         }
-
         const double kplus2pt = std::hypot(t2.phid1Px(), t2.phid1Py());
         const double kminus2pt = std::hypot(t2.phid2Px(), t2.phid2Py());
 
@@ -1736,25 +1736,20 @@ struct doublephimeson {
             t1.phid1Index() == t2.phid2Index() ||
             t1.phid2Index() == t2.phid1Index() ||
             t1.phid2Index() == t2.phid2Index()) {
-          // LOGF(info,"track share",t1.phid1Index(),t1.phid2Index(),t2.phid1Index(),t2.phid2Index());
+          LOGF(info, "track share %d %d %d %d", t1.phid1Index(), t1.phid2Index(), t2.phid1Index(), t2.phid2Index());
           continue;
         }
-        const double mCross12 = (k1p + k2m).M(); // K+ from phi1 + K- from phi2
-        const double mCross21 = (k2p + k1m).M(); // K+ from phi2 + K- from phi1
-
-        const bool cross12IsPhiLike = (mCross12 > cfgCrossPhiLow && mCross12 < cfgCrossPhiHigh);
-        const bool cross21IsPhiLike = (mCross21 > cfgCrossPhiLow && mCross21 < cfgCrossPhiHigh);
 
-        if (cross12IsPhiLike || cross21IsPhiLike) {
-          LOGF(info,
-               "Best-pairing rejected: mPhi1 = %3.4f, mPhi2 = %3.4f, mCross12 = %3.4f, mCross21 = %3.4f",
-               phi1.M(),
-               phi2.M(),
-               mCross12,
-               mCross21,
-               pair.Pt(),
-               pair.M());
-          continue;
+        auto cross12 = k1p + k2m;
+        auto cross21 = k2p + k1m;
+        bool alternativePairValid = cross12.M() > cfgCrossPhiLow && cross12.M() < cfgCrossPhiHigh && cross21.M() > cfgCrossPhiLow && cross21.M() < cfgCrossPhiHigh;
+        if (alternativePairValid) {
+          float scoreOriginal = deltaMPhiNominal(phi1.M(), phi2.M());
+          float scoreCross = deltaMPhiNominal(cross12.M(), cross21.M());
+          if (useCrossPairRejection && (scoreCross < scoreOriginal)) {
+            LOGF(info, "Best-pairing rejected: original score = %3.4f, cross scoremPhi2 = %3.4f", scoreOriginal, scoreCross);
+            continue; // another pairing of these four tracks is better
+          }
         }
         histos.fill(HIST("hPhiMass"), phi1.M(), phi2.M(), pair.Pt());
         histos.fill(HIST("hPhiMassNormalized"), getNormalizedMPhi(phi1.M(), phi1.Pt()), getNormalizedMPhi(phi2.M(), phi2.Pt()), pair.Pt());
@@ -1808,7 +1803,7 @@ struct doublephimeson {
       }
       if (pairPt > minExoticPt) {
         histos.fill(HIST("hPtCorrelation"), pairPt, ptcorr);
-        histos.fill(HIST("hMassCent"), p1.M(), p2.M(), collision.centrality());
+        // histos.fill(HIST("hMassCent"), p1.M(), p2.M(), collision.centrality());
         histos.fill(HIST("SEMassUnlike_AllVars"),
                     M,
                     pairPt,