From 8ed154d79ff9cf72a8ba4c9140d7c932dbeda118 Mon Sep 17 00:00:00 2001 From: jesgum Date: Mon, 12 Oct 2026 15:03:21 +0200 Subject: [PATCH] Add mult per ancestor histogram to nbd fitter --- .../Multiplicity/multGlauberNBDFitter.cxx | 18 +++++++++++++++++- .../Tools/Multiplicity/multGlauberNBDFitter.h | 15 ++++++++++++++- 2 files changed, 31 insertions(+), 2 deletions(-) diff --git a/Common/Tools/Multiplicity/multGlauberNBDFitter.cxx b/Common/Tools/Multiplicity/multGlauberNBDFitter.cxx index e3c4635533d..c5d067c7e22 100644 --- a/Common/Tools/Multiplicity/multGlauberNBDFitter.cxx +++ b/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 , 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); } diff --git a/Common/Tools/Multiplicity/multGlauberNBDFitter.h b/Common/Tools/Multiplicity/multGlauberNBDFitter.h index 7bf8c5a0346..3c5f6370722 100644 --- a/Common/Tools/Multiplicity/multGlauberNBDFitter.h +++ b/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