Add qa for new centrality tables

Author: jesgum

Common/Tasks/centralityQa.cxx

@@ -155,6 +155,8 @@ struct CentralityQa {
       histos.add("hCentNTPV", ";NTPV centrality (%)", kTH1D, {{nBins, 0, 105.}});
       histos.add("hCentNGlobal", ";NGlobal centrality (%)", kTH1D, {{nBins, 0, 105.}});
       histos.add("hCentMFT", ";MFT centrality (%)", kTH1D, {{nBins, 0, 105.}});
+      histos.add("hCentFT0MLightIonAnchorCols", ";FT0MLightIonAnchorCols centrality (%)", kTH1D, {{nBins, 0, 105.}});
+      histos.add("hCentFT0MLightIonAnchorBCs", ";FT0MLightIonAnchorBCs centrality (%)", kTH1D, {{nBins, 0, 105.}});
 
       // profiles of midrapidity multiplicity density
       histos.add("hCentProfileFV0A", ";FV0A centrality (%)", kTProfile, {{nBins, 0, 105.}});
@@ -167,6 +169,8 @@ struct CentralityQa {
       histos.add("hCentProfileNTPV", ";NTPV centrality (%)", kTProfile, {{nBins, 0, 105.}});
       histos.add("hCentProfileNGlobal", ";NGlobal centrality (%)", kTProfile, {{nBins, 0, 105.}});
       histos.add("hCentProfileMFT", ";MFT centrality (%)", kTProfile, {{nBins, 0, 105.}});
+      histos.add("hCentProfileFT0MLightIonAnchorCols", ";FT0MLightIonAnchorCols centrality (%)", kTProfile, {{nBins, 0, 105.}});
+      histos.add("hCentProfileFT0MLightIonAnchorBCs", ";FT0MLightIonAnchorBCs centrality (%)", kTProfile, {{nBins, 0, 105.}});
 
       histos.add("hMultEta05VsCentFV0A", ";FV0A centrality (%); Multiplicity PV contributors (|#it{#eta}| < 0.5)", kTH2D, {{nBins, 0, 105.}, axisMultiplicityPV});
       histos.add("hMultEta05VsCentFT0M", ";FT0M centrality (%); Multiplicity PV contributors (|#it{#eta}| < 0.5)", kTH2D, {{nBins, 0, 105.}, axisMultiplicityPV});
@@ -178,6 +182,8 @@ struct CentralityQa {
       histos.add("hMultEta05VsCentNTPV", ";NTPV centrality (%); Multiplicity PV contributors (|#it{#eta}| < 0.5)", kTH2D, {{nBins, 0, 105.}, axisMultiplicityPV});
       histos.add("hMultEta05VsCentNGlobal", ";NGlobal centrality (%); Multiplicity PV contributors (|#it{#eta}| < 0.5)", kTH2D, {{nBins, 0, 105.}, axisMultiplicityPV});
       histos.add("hMultEta05VsCentMFT", ";MFT centrality (%); Multiplicity PV contributors (|#it{#eta}| < 0.5)", kTH2D, {{nBins, 0, 105.}, axisMultiplicityPV});
+      histos.add("hMultEta05VsCentFT0MLightIonAnchorCols", ";FT0MLightIonAnchorCols centrality (%); Multiplicity PV contributors (|#it{#eta}| < 0.5)", kTH2D, {{nBins, 0, 105.}, axisMultiplicityPV});
+      histos.add("hMultEta05VsCentFT0MLightIonAnchorBCs", ";FT0MLightIonAnchorBCs centrality (%); Multiplicity PV contributors (|#it{#eta}| < 0.5)", kTH2D, {{nBins, 0, 105.}, axisMultiplicityPV});
 
       if (isMC) {
         histos.add("hMultEta05VsGenMultFV0A", ";Multiplicity FV0A; Multiplicity PV contributors (|#it{#eta}| < 0.5)", kTH2D, {axisMultiplicity, axisMultiplicityPV});
@@ -656,6 +662,28 @@ struct CentralityQa {
   }
   PROCESS_SWITCH(CentralityQa, processRun3_MFT, "Process with Run 3 MFT estimator", false);
 
+  void processRun3_FT0MLightIonAnchorCol(soa::Join<aod::Collisions, aod::EvSels, aod::Mults, aod::CentFT0MLightIonAnchorCols>::iterator const& col)
+  {
+    if (!isCollisionAccepted(col)) {
+      return;
+    }
+    histos.fill(HIST("hCentFT0MLightIonAnchorCols"), col.centFT0MLightIonAnchorCol());
+    histos.fill(HIST("hCentProfileFT0MLightIonAnchorCols"), col.centFT0MLightIonAnchorCol(), col.multNTracksPVetaHalf());
+    histos.fill(HIST("hMultEta05VsCentFT0MLightIonAnchorCols"), col.centFT0MLightIonAnchorCol(), col.multNTracksPVetaHalf());
+  }
+  PROCESS_SWITCH(CentralityQa, processRun3_FT0MLightIonAnchorCol, "Process with Run 3 FT0MLightIonAnchorCol estimator", false);
+
+  void processRun3_FT0MLightIonAnchorBC(soa::Join<aod::Collisions, aod::EvSels, aod::Mults, aod::CentFT0MLightIonAnchorBCs>::iterator const& col)
+  {
+    if (!isCollisionAccepted(col)) {
+      return;
+    }
+    histos.fill(HIST("hCentFT0MLightIonAnchorBCs"), col.centFT0MLightIonAnchorBC());
+    histos.fill(HIST("hCentProfileFT0MLightIonAnchorBCs"), col.centFT0MLightIonAnchorBC(), col.multNTracksPVetaHalf());
+    histos.fill(HIST("hMultEta05VsCentFT0MLightIonAnchorBCs"), col.centFT0MLightIonAnchorBC(), col.multNTracksPVetaHalf());
+  }
+  PROCESS_SWITCH(CentralityQa, processRun3_FT0MLightIonAnchorBC, "Process with Run 3 FT0MLightIonAnchorBC estimator", false);
+
   void processMonteCarloRun3_FV0A(soa::Join<aod::Collisions, aod::McCollisionLabels, aod::EvSels, aod::Mults, aod::CentFV0As>::iterator const& col,
                                   soa::Join<aod::McCollisions, aod::MultMCExtras> const& /*mcCollisions*/,
                                   soa::Join<aod::BCs, aod::Run3MatchedToBCSparse> const& /*bcs*/)