[PWGCF] Add centrality cut to the debug track source file

PWGCF/FemtoDream/Tasks/femtoDreamDebugTrack.cxx

@@ -1,4 +1,4 @@
 // Copyright 2019-2025 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.
 //
@@ -55,7 +55,7 @@
   ConfigurableAxis ConfBinmultPercentile{"ConfBinmultPercentile", {10, 0.0f, 100.0f}, "multiplicity percentile Binning"};
   ConfigurableAxis ConfBinpT{"ConfBinpT", {{240, 0, 6}}, "pT binning"};
   ConfigurableAxis ConfBineta{"ConfBineta", {{200, -1.5, 1.5}}, "eta binning"};
   ConfigurableAxis ConfBinphi{"ConfBinphi", {{200, 0, TMath::TwoPi()}}, "phi binning"};
 
   ConfigurableAxis ConfTempFitVarBins{"ConfTempFitVarBins", {300, -0.15, 0.15}, "Binning of the TempFitVar"};
   ConfigurableAxis ConfNsigmaTPCBins{"ConfNsigmaTPCBins", {1600, -8, 8}, "Binning of Nsigma TPC plot"};
@@ -66,6 +66,10 @@
   Configurable<int> ConfTempFitVarMomentum{"ConfTempFitVarMomentum", 0, "Momentum used for binning: 0 -> pt; 1 -> preco; 2 -> ptpc"};
   ConfigurableAxis ConfDummy{"ConfDummy", {1, 0, 1}, "Dummy axis for inv mass"};
 
+  Configurable<bool> ConfdoCentCut{"ConfdoCentCut", false, "Enable centrality cut"};
+  Configurable<float> ConfCentMax{"ConfCentMax", 100., "Upper limit of centrality cut"};
+  Configurable<float> ConfCentMin{"ConfCentMin", 0., "Lower limit of centrality cut"};
+
   using FemtoMCCollisions = Join<aod::FDCollisions, aod::FDMCCollLabels>;
   using FemtoMCCollision = FemtoMCCollisions::iterator;
 
@@ -110,7 +114,12 @@
   void FillDebugHistos(CollisionType& col, PartitionType& groupPartsOne)
   {
     eventHisto.fillQA<isMC>(col);
+
+    if (ConfdoCentCut.value && (col.multV0M() > ConfCentMax || col.multV0M() < ConfCentMin)) {
+      return;
+    }
+
     for (auto& part : groupPartsOne) {
       trackHisto.fillQA<isMC, true>(part, static_cast<aod::femtodreamparticle::MomentumType>(ConfTempFitVarMomentum.value), col.multNtr(), col.multV0M(), ConfOptCorrelatedPlots);
     }
   }
@@ -118,7 +127,7 @@
   /// process function when runnning over data/ Monte Carlo reconstructed only
   /// \param col subscribe to FemtoDreamCollision table
   /// \param parts subscribe to FemtoDreamParticles table
   void processData(o2::aod::FDCollision& col, FemtoFullParticles&)
   {
     auto groupPartsOne = partsOne->sliceByCached(aod::femtodreamparticle::fdCollisionId, col.globalIndex(), cache);
     FillDebugHistos<false>(col, groupPartsOne);
@@ -130,7 +139,7 @@
   /// \param col subscribe to FemtoDreamCollision table
   /// \param parts subscribe to the joined table of FemtoDreamParticles and FemtoDreamMCLabels table
   /// \param FemtoDramMCParticles subscribe to the table containing the Monte Carlo Truth information
   void processMC(FemtoMCCollision& col, o2::aod::FDMCCollisions&, FemtoFullParticlesMC& /*parts*/, aod::FDMCParticles&, aod::FDExtMCParticles&)
   {
     auto groupPartsOne = partsOneMC->sliceByCached(aod::femtodreamparticle::fdCollisionId, col.globalIndex(), cache);
     FillDebugHistos<true>(col, groupPartsOne);