[PWGLF] timing cuts, post-processing update (#16058)

Co-authored-by: ALICE Action Bot <alibuild@cern.ch>

Author: nepeivodaRS

PWGLF/TableProducer/Strangeness/cascqaanalysis.cxx

@@ -89,8 +89,8 @@ struct Cascqaanalysis {
   ConfigurableAxis nChargedFT0MGenAxis{"nChargedFT0MGenAxis", {300, 0, 300}, "N_{FT0M, gen.}"};
   ConfigurableAxis nChargedFV0AGenAxis{"nChargedFV0AGenAxis", {300, 0, 300}, "N_{FV0A, gen.}"};
   ConfigurableAxis multNTracksAxis{"multNTracksAxis", {500, 0, 500}, "N_{tracks}"};
-  ConfigurableAxis signalFT0MAxis{"signalFT0MAxis", {10000, 0, 40000}, "FT0M amplitude"};
-  ConfigurableAxis signalFV0AAxis{"signalFV0AAxis", {10000, 0, 40000}, "FV0A amplitude"};
+  Configurable<int> signalFT0MNBins{"signalFT0MNBins", 1000, "Number of bins for FT0M amplitude QA axis"};
+  Configurable<int> signalFV0ANBins{"signalFV0ANBins", 1000, "Number of bins for FV0A amplitude QA axis"};
   ConfigurableAxis nCandidates{"nCandidates", {30, -0.5, 29.5}, "N_{cand.}"};
 
   // Event selection criteria
@@ -149,11 +149,21 @@ struct Cascqaanalysis {
     uint8_t typeFlag;
   } CollisionIndexAndType;
 
+  static constexpr unsigned int kNITSLayers = 7;
+  static constexpr float kGlobalTrackEtaMax = 0.5f;
+  static constexpr float kFT0CMinEta = -3.3f;
+  static constexpr float kFT0CMaxEta = -2.1f;
+  static constexpr float kFT0AMinEta = 3.5f;
+  static constexpr float kFT0AMaxEta = 4.9f;
+  static constexpr float kFV0AMinEta = 2.2f;
+  static constexpr float kFV0AMaxEta = 5.1f;
+  static constexpr size_t kNContributorsCorrelationSize = 2;
+
   template <typename TTrack>
   static int countITSHits(TTrack const& track)
   {
     int nHits = 0;
-    for (unsigned int i = 0; i < 7; ++i) {
+    for (unsigned int i = 0; i < kNITSLayers; ++i) {
       if (track.itsClusterMap() & (1 << i)) {
         ++nHits;
       }
@@ -225,6 +235,8 @@ struct Cascqaanalysis {
     }
 
     if (multQA) {
+      AxisSpec signalFT0MAxis = {std::max(1, static_cast<int>(signalFT0MNBins)), 0.f, 40000.f, "FT0M amplitude"};
+      AxisSpec signalFV0AAxis = {std::max(1, static_cast<int>(signalFV0ANBins)), 0.f, 40000.f, "FV0A amplitude"};
       if (isMC) {
         // Rec. lvl
         registry.add("hNchFT0Mglobal", "hNchFT0Mglobal", {HistType::kTH3D, {nChargedFT0MGenAxis, multNTracksAxis, eventTypeAxis}});
@@ -253,7 +265,7 @@ struct Cascqaanalysis {
      aod::cascdata::dcacascdaughters < dcacascdau);
 
   Partition<DauTracks> pvContribTracksIUEta1 = (nabs(aod::track::eta) < 1.0f) && ((aod::track::flags & static_cast<uint32_t>(o2::aod::track::PVContributor)) == static_cast<uint32_t>(o2::aod::track::PVContributor));
-  Partition<DauTracks> globalTracksIUEta05 = (nabs(aod::track::eta) < 0.5f) && (requireGlobalTrackInFilter());
+  Partition<DauTracks> globalTracksIUEta05 = (nabs(aod::track::eta) < kGlobalTrackEtaMax) && (requireGlobalTrackInFilter());
 
   template <class TCascTracksTo, typename TCascade>
   bool acceptCascCandidate(TCascade const& cascCand, float const& pvx, float const& pvy, float const& pvz)
@@ -289,7 +301,7 @@ struct Cascqaanalysis {
       if (pdgInfo->Charge() == 0) {
         continue;
       }
-      if (mcParticle.eta() < -3.3 || mcParticle.eta() > 4.9 || (mcParticle.eta() > -2.1 && mcParticle.eta() < 3.5)) {
+      if (mcParticle.eta() < kFT0CMinEta || mcParticle.eta() > kFT0AMaxEta || (mcParticle.eta() > kFT0CMaxEta && mcParticle.eta() < kFT0AMinEta)) {
         continue; // select on T0M Nch region
       }
       nchFT0++; // increment
@@ -313,7 +325,7 @@ struct Cascqaanalysis {
       if (pdgInfo->Charge() == 0) {
         continue;
       }
-      if (mcParticle.eta() < 2.2 || mcParticle.eta() > 5.1) {
+      if (mcParticle.eta() < kFV0AMinEta || mcParticle.eta() > kFV0AMaxEta) {
         continue; // select on V0A Nch region
       }
       nchFV0A++; // increment
@@ -695,7 +707,7 @@ struct Cascqaanalysis {
 
     registry.fill(HIST("hNchFT0MNAssocMCCollisions"), nchFT0, nAssocColl, evType);
 
-    if (numberOfContributors.size() == 2) {
+    if (numberOfContributors.size() == kNContributorsCorrelationSize) {
       std::sort(numberOfContributors.begin(), numberOfContributors.end());
       registry.fill(HIST("hNContributorsCorrelation"), numberOfContributors[0], numberOfContributors[1]);
     }