fix bug and add centQa

Author: jesgum

Common/Tasks/centralityQa.cxx

@@ -206,6 +206,7 @@ struct CentralityQa {
     } else {
       histos.add("hCentFV0A", ";FV0A centrality (%)", kTH1D, {{nBins, 0, 105.}});
       histos.add("hCentFT0M", ";FT0M centrality (%)", kTH1D, {{nBins, 0, 105.}});
+      histos.add("hCentFT0MVar3", ";FT0M centrality (%)", kTH1D, {{nBins, 0, 105.}});
       histos.add("hCentFT0A", ";FT0A centrality (%)", kTH1D, {{nBins, 0, 105.}});
       histos.add("hCentFT0C", ";FT0C centrality (%)", kTH1D, {{nBins, 0, 105.}});
       histos.add("hCentFT0CVar1", ";FT0CVar1 centrality (%)", kTH1D, {{nBins, 0, 105.}});
@@ -220,6 +221,7 @@ struct CentralityQa {
       // profiles of midrapidity multiplicity density
       histos.add("hCentProfileFV0A", ";FV0A centrality (%)", kTProfile, {{nBins, 0, 105.}});
       histos.add("hCentProfileFT0M", ";FT0M centrality (%)", kTProfile, {{nBins, 0, 105.}});
+      histos.add("hCentProfileFT0MVar3", ";FT0M centrality (%)", kTProfile, {{nBins, 0, 105.}});
       histos.add("hCentProfileFT0A", ";FT0A centrality (%)", kTProfile, {{nBins, 0, 105.}});
       histos.add("hCentProfileFT0C", ";FT0C centrality (%)", kTProfile, {{nBins, 0, 105.}});
       histos.add("hCentProfileFT0CVar1", ";FT0CVar1 centrality (%)", kTProfile, {{nBins, 0, 105.}});
@@ -233,6 +235,7 @@ struct CentralityQa {
 
       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});
+      histos.add("hMultEta05VsCentFT0MVar3", ";FT0M centrality (%); Multiplicity PV contributors (|#it{#eta}| < 0.5)", kTH2D, {{nBins, 0, 105.}, axisMultiplicityPV});
       histos.add("hMultEta05VsCentFT0A", ";FT0A centrality (%); Multiplicity PV contributors (|#it{#eta}| < 0.5)", kTH2D, {{nBins, 0, 105.}, axisMultiplicityPV});
       histos.add("hMultEta05VsCentFT0C", ";FT0C centrality (%); Multiplicity PV contributors (|#it{#eta}| < 0.5)", kTH2D, {{nBins, 0, 105.}, axisMultiplicityPV});
       histos.add("hMultEta05VsCentFT0CVar1", ";FT0CVar1 centrality (%); Multiplicity PV contributors (|#it{#eta}| < 0.5)", kTH2D, {{nBins, 0, 105.}, axisMultiplicityPV});
@@ -323,6 +326,7 @@ struct CentralityQa {
       requires { collision.centFT0C(); } ||
       requires { collision.centFT0CVariant1(); } ||
       requires { collision.centFT0CVariant2(); } ||
+      requires { collision.centFT0MVariant3(); } ||
       requires { collision.centFDDM(); } ||
       requires { collision.centNTPV(); } ||
       requires { collision.centNGlobal(); } ||
@@ -684,7 +688,19 @@ struct CentralityQa {
   }
   PROCESS_SWITCH(CentralityQa, processRun3_FT0M, "Process with Run 3 FT0M estimator", false);
 
-  void processRun3_FT0A(soa::Join<aod::Collisions, aod::EvSels, aod::MultsRun3, aod::CentFT0As>::iterator const& col, aod::BCs const&)
+  void processRun3_FT0Mvar3(soa::Join<aod::Collisions, aod::EvSels, aod::Mults, aod::CentFT0MVariant3s>::iterator const& col)
+  {
+    if (!isCollisionAccepted(col)) {
+      return;
+    }
+    LOGF(debug, "centFT0MVar3=%.0f", col.centFT0MVariant3());
+    histos.fill(HIST("hCentFT0MVar3"), col.centFT0MVariant3());
+    histos.fill(HIST("hCentProfileFT0MVar3"), col.centFT0MVariant3(), col.multNTracksPVetaHalf());
+    histos.fill(HIST("hMultEta05VsCentFT0MVar3"), col.centFT0MVariant3(), col.multNTracksPVetaHalf());
+  }
+  PROCESS_SWITCH(CentralityQa, processRun3_FT0Mvar3, "Process with Run 3 FT0M estimator", false);
+
+  void processRun3_FT0A(soa::Join<aod::Collisions, aod::EvSels, aod::Mults, aod::CentFT0As>::iterator const& col)
   {
     if (!isCollisionAccepted(col)) {
       return;

Common/Tools/Multiplicity/MultModule.h

@@ -816,9 +816,9 @@ class MultModule
         mults.multFT0AZeq = 0.0f;
       }
       if (mults.multFT0AOuter > -1.0f && std::fabs(collision.posZ()) < 15.0f && lCalibLoaded) {
-        mults.multFT0AOuter = hVtxZFT0A->Interpolate(0.0) * mults.multFT0AOuter / hVtxZFT0A->Interpolate(collision.posZ());
+        mults.multFT0AOuterZeq = hVtxZFT0A->Interpolate(0.0) * mults.multFT0AOuter / hVtxZFT0A->Interpolate(collision.posZ());
       } else {
-        mults.multFT0AOuter = 0.0f;
+        mults.multFT0AOuterZeq = 0.0f;
       }
       if (mults.multFT0C > -1.0f && std::fabs(collision.posZ()) < 15.0f && lCalibLoaded) {
         mults.multFT0CZeq = hVtxZFT0C->Interpolate(0.0) * mults.multFT0C / hVtxZFT0C->Interpolate(collision.posZ());