[PWGLF] Add centrality info in derfide data table for double phi (AliceO2Group#16782)

Author: skundu692

PWGLF/DataModel/ReducedDoublePhiTables.h

@@ -27,6 +27,7 @@ namespace redphievent
 {
 DECLARE_SOA_COLUMN(NumPos, numPos, int); //! Number of positive Kaon
 DECLARE_SOA_COLUMN(NumNeg, numNeg, int); //! Number of negative Kaon
+DECLARE_SOA_COLUMN(Centrality, centrality, float); //! Number of negative Kaon
 } // namespace redphievent
 DECLARE_SOA_TABLE(RedPhiEvents, "AOD", "REDPHIEVENT",
                   o2::soa::Index<>,
@@ -36,7 +37,8 @@ DECLARE_SOA_TABLE(RedPhiEvents, "AOD", "REDPHIEVENT",
                   collision::PosZ,
                   collision::NumContrib,
                   redphievent::NumPos,
-                  redphievent::NumNeg);
+                  redphievent::NumNeg,
+                  redphievent::Centrality);
 using RedPhiEvent = RedPhiEvents::iterator;
 
 namespace phitrack

PWGLF/TableProducer/Resonances/doublephitable.cxx

@@ -17,6 +17,7 @@
 #include "PWGLF/DataModel/ReducedDoublePhiTables.h"
 
 #include "Common/CCDB/EventSelectionParams.h"
+#include "Common/DataModel/Centrality.h"
 #include "Common/DataModel/EventSelection.h"
 #include "Common/DataModel/Multiplicity.h"
 #include "Common/DataModel/PIDResponseITS.h"
@@ -83,7 +84,7 @@ struct doublephitable {
   Filter DCAcutFilter = (nabs(aod::track::dcaXY) < cfgCutDCAxy) && (nabs(aod::track::dcaZ) < cfgCutDCAz);
   Filter PIDcutFilter = nabs(aod::pidtpc::tpcNSigmaKa) < nsigmaCutTPC;
 
-  using EventCandidates = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::Mults>>;
+  using EventCandidates = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::Mults, aod::CentFT0Ms>>;
   using TrackCandidates = soa::Filtered<soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, aod::TrackSelection, aod::pidTOFbeta, aod::pidTPCFullKa, aod::pidTOFFullKa>>;
 
   SliceCache cache;
@@ -165,7 +166,7 @@ struct doublephitable {
 
     int Npostrack = 0;
     int Nnegtrack = 0;
-
+    float centrality = collision.centFT0M();
     if (collision.sel8() && collision.selection_bit(aod::evsel::kNoTimeFrameBorder) && collision.selection_bit(aod::evsel::kNoITSROFrameBorder) && collision.selection_bit(aod::evsel::kNoSameBunchPileup) && collision.selection_bit(aod::evsel::kIsGoodZvtxFT0vsPV)) {
       auto posThisColl = posTracks->sliceByCached(aod::track::collisionId, collision.globalIndex(), cache);
       auto negThisColl = negTracks->sliceByCached(aod::track::collisionId, collision.globalIndex(), cache);
@@ -258,7 +259,7 @@ struct doublephitable {
     if (keepEventDoublePhi && numberPhi > 1 && Npostrack > 1 && Nnegtrack > 1 && (phiresonance.size() == phiresonanced1.size()) && (phiresonance.size() == phiresonanced2.size())) {
       qaRegistry.fill(HIST("hEventstat"), 1.5);
       /////////// Fill collision table///////////////
-      redPhiEvents(bc.globalBC(), currentRunNumber, bc.timestamp(), collision.posZ(), collision.numContrib(), Npostrack, Nnegtrack);
+      redPhiEvents(bc.globalBC(), currentRunNumber, bc.timestamp(), collision.posZ(), collision.numContrib(), Npostrack, Nnegtrack, centrality);
       auto indexEvent = redPhiEvents.lastIndex();
       //// Fill track table for Phi//////////////////
       for (auto if1 = phiresonance.begin(); if1 != phiresonance.end(); ++if1) {

PWGLF/Tasks/Resonances/doublephimeson.cxx

@@ -181,6 +181,7 @@ struct doublephimeson {
     histos.add("hDeltaRkaonplus", "hDeltaRkaonplus", kTH1F, {{800, 0.0, 8.0}});
     histos.add("hDeltaRkaonminus", "hDeltaRkaonminus", kTH1F, {{800, 0.0, 8.0}});
     histos.add("hPtCorrelation", "hPtCorrelation", kTH2F, {{400, 0.0, 40.0}, {5000, 0.0, 100.0}});
+    histos.add("hPtCent", "hPtCent", kTH2F, {{100, 0.0, 100.0}, {100, 0.0, 100.0}});
     const AxisSpec thnAxisdeltapt{configThnAxisDeltaPt, "Delta pt"};
     const AxisSpec thnAxisdaughterpt{configThnAxisDaughterPt, "Daughter pt"};
     const AxisSpec thnAxisInvMass{configThnAxisInvMass, "#it{M} (GeV/#it{c}^{2})"};
@@ -1780,6 +1781,7 @@ struct doublephimeson {
       }
       if (pairPt > minExoticPt) {
         histos.fill(HIST("hPtCorrelation"), pairPt, ptcorr);
+        histos.fill(HIST("hPtCent"), pairPt, collision.centrality());
         histos.fill(HIST("SEMassUnlike_AllVars"),
                     M,
                     pairPt,