Add mult per ancestor histogram to nbd fitter

Author: jesgum

Common/Tools/Multiplicity/multGlauberNBDFitter.cxx

@@ -443,7 +443,20 @@ double multGlauberNBDFitter::ContinuousNBD(const double n, const double mu, cons
   return F;
 }
 
-void multGlauberNBDFitter::CalculateAvNpNc(TProfile* lNPartProf, TProfile* lNCollProf, TH2F* lNPart2DPlot, TH2F* lNColl2DPlot, TH1F* hPercentileMap, double lLoRange, double lHiRange, TH3D* lNpNcEcc, TH2F* lEcc2DPlot, TH3D* lNpNcB, TH2F* lB2DPlot, TH2F* lNancestor2DPlot, double fProbabilityCutoff)
+void multGlauberNBDFitter::CalculateAvNpNc(TProfile* lNPartProf,
+                                           TProfile* lNCollProf,
+                                           TH2F* lNPart2DPlot,
+                                           TH2F* lNColl2DPlot,
+                                           TH1F* hPercentileMap,
+                                           double lLoRange,
+                                           double lHiRange,
+                                           TH3D* lNpNcEcc,
+                                           TH2F* lEcc2DPlot,
+                                           TH3D* lNpNcB,
+                                           TH2F* lB2DPlot,
+                                           TH2F* lNancestor2DPlot,
+                                           TH2F* lMult2DPlot,
+                                           double fProbabilityCutoff)
 {
   LOG(info) << "Calculating <Npart>, <Ncoll> in centrality bins...";
   LOG(info) << "Range to calculate: " << lLoRange << " to " << lHiRange;
@@ -557,6 +570,9 @@ void multGlauberNBDFitter::CalculateAvNpNc(TProfile* lNPartProf, TProfile* lNCol
         // fill cross-check histogram with lNancestorCount at lNancestors value
         lNancestor2DPlot->Fill(lMultValueToFill, lNancestors, lProbability);
       }
+      if (lMult2DPlot) {
+        lMult2DPlot->Fill(lMultValueToFill, lMultValue, lProbability);
+      }
       if (lNPart2DPlot) {
         lNPart2DPlot->Fill(lMultValueToFill, fNpart[ibin], lProbability);
       }

Common/Tools/Multiplicity/multGlauberNBDFitter.h

@@ -131,7 +131,20 @@ class multGlauberNBDFitter : public TNamed
 
   // For estimating Npart, Ncoll in multiplicity bins
   // also viable: eccentricity, impact parameter, ancestor cross-check plot
-  void CalculateAvNpNc(TProfile* lNPartProf, TProfile* lNCollProf, TH2F* lNPart2DPlot, TH2F* lNColl2DPlot, TH1F* hPercentileMap, double lLoRange = -1, double lHiRange = -1, TH3D* lNpNcEcc = nullptr, TH2F* lEcc2DPlot = nullptr, TH3D* lNpNcB = nullptr, TH2F* lB2DPlot = nullptr, TH2F* lNancestor2DPlot = nullptr, double fProbabilityCutoff = -1);
+  void CalculateAvNpNc(TProfile* lNPartProf,
+                       TProfile* lNCollProf,
+                       TH2F* lNPart2DPlot,
+                       TH2F* lNColl2DPlot,
+                       TH1F* hPercentileMap,
+                       double lLoRange = -1,
+                       double lHiRange = -1,
+                       TH3D* lNpNcEcc = nullptr,
+                       TH2F* lEcc2DPlot = nullptr,
+                       TH3D* lNpNcB = nullptr,
+                       TH2F* lB2DPlot = nullptr,
+                       TH2F* lNancestor2DPlot = nullptr,
+                       TH2F* lMult2DPlot = nullptr,
+                       double fProbabilityCutoff = -1);
 
  private:
   // This function serves as the (analytical) NBD