From 3a591912033af4e72c277b62fd3e585520f295ba Mon Sep 17 00:00:00 2001 From: oussama3030 Date: Mon, 31 Aug 2026 17:32:40 +0200 Subject: [PATCH] replaced two TH3D histograms with THnD --- .../Tasks/Resonances/k892hadronphotonBkg.cxx | 24 +++++++------------ PWGLF/Tasks/Strangeness/sigmaanalysis.cxx | 24 +++++++------------ 2 files changed, 16 insertions(+), 32 deletions(-) diff --git a/PWGLF/Tasks/Resonances/k892hadronphotonBkg.cxx b/PWGLF/Tasks/Resonances/k892hadronphotonBkg.cxx index 7b1eff82351..3350ce3512b 100644 --- a/PWGLF/Tasks/Resonances/k892hadronphotonBkg.cxx +++ b/PWGLF/Tasks/Resonances/k892hadronphotonBkg.cxx @@ -317,8 +317,7 @@ struct k892hadronphotonBkg { histos.add("KStarBkg/h3dRotKStarMassVsPt", "h3dRotKStarMassVsPt", kTH3D, {axisConfig.axisCentrality, axisConfig.axisPt, axisConfig.axisKStarMass}); histos.add("KStarBkg/h3dRotKStarPtVsOPAngle", "h3dRotKStarPtVsOPAngle", kTH3D, {{140, 0.f, 7.f}, axisConfig.axisPt, axisConfig.axisKStarMass}); if (doArm) { - histos.add("KStarBkg/h3dRotKStarPtVsAPAlpha", "h3dRotKStarPtVsAPAlpha", kTH3D, {axisConfig.axisAPAlpha, axisConfig.axisPt, axisConfig.axisKStarMass}); - histos.add("KStarBkg/h3dRotKStarPtVsAPQt", "h3dRotKStarPtVsAPQt", kTH3D, {axisConfig.axisAPQt, axisConfig.axisPt, axisConfig.axisKStarMass}); + histos.add("KStarBkg/h4dRotKStarPtVsAPAlphaVsAPQt", "h4dRotKStarPtVsAPAlphaVsAPQt", kTHnD, {axisConfig.axisAPAlpha, axisConfig.axisAPQt, axisConfig.axisPt, axisConfig.axisKStarMass}); } } if (kstarBkgConfig.doEvtMixing) { @@ -326,8 +325,7 @@ struct k892hadronphotonBkg { histos.add("KStarBkg/h3dMixedKStarMassVsPt", "h3dMixedKStarMassVsPt", kTH3D, {axisConfig.axisCentrality, axisConfig.axisPt, axisConfig.axisKStarMass}); histos.add("KStarBkg/h3dMixedKStarPtVsOPAngle", "h3dMixedKStarPtVsOPAngle", kTH3D, {{140, 0.f, 7.f}, axisConfig.axisPt, axisConfig.axisKStarMass}); if (doArm) { - histos.add("KStarBkg/h3dMixedKStarPtVsAPAlpha", "h3dMixedKStarPtVsAPAlpha", kTH3D, {axisConfig.axisAPAlpha, axisConfig.axisPt, axisConfig.axisKStarMass}); - histos.add("KStarBkg/h3dMixedKStarPtVsAPQt", "h3dMixedKStarPtVsAPQt", kTH3D, {axisConfig.axisAPQt, axisConfig.axisPt, axisConfig.axisKStarMass}); + histos.add("KStarBkg/h4dMixedKStarPtVsAPAlphaVsAPQt", "h4dMixedKStarPtVsAPAlphaVsAPQt", kTHnD, {axisConfig.axisAPAlpha, axisConfig.axisAPQt, axisConfig.axisPt, axisConfig.axisKStarMass}); } } @@ -341,8 +339,7 @@ struct k892hadronphotonBkg { histos.add("LambdaStarBkg/h3dRotLambdaStarMassVsPt", "h3dRotLambdaStarMassVsPt", kTH3D, {axisConfig.axisCentrality, axisConfig.axisPt, axisConfig.axisLambdaStarMass}); histos.add("LambdaStarBkg/h3dRotLambdaStarPtVsOPAngle", "h3dRotLambdaStarPtVsOPAngle", kTH3D, {{140, 0.f, 7.f}, axisConfig.axisPt, axisConfig.axisLambdaStarMass}); if (doArm) { - histos.add("LambdaStarBkg/h3dRotLambdaStarPtVsAPAlpha", "h3dRotLambdaStarPtVsAPAlpha", kTH3D, {axisConfig.axisAPAlpha, axisConfig.axisPt, axisConfig.axisLambdaStarMass}); - histos.add("LambdaStarBkg/h3dRotLambdaStarPtVsAPQt", "h3dRotLambdaStarPtVsAPQt", kTH3D, {axisConfig.axisAPQt, axisConfig.axisPt, axisConfig.axisLambdaStarMass}); + histos.add("LambdaStarBkg/h4dRotLambdaStarPtVsAPAlphaVsAPQt", "h4dRotLambdaStarPtVsAPAlphaVsAPQt", kTHnD, {axisConfig.axisAPAlpha, axisConfig.axisAPQt, axisConfig.axisPt, axisConfig.axisLambdaStarMass}); } } if (lstarBkgConfig.doEvtMixing) { @@ -350,8 +347,7 @@ struct k892hadronphotonBkg { histos.add("LambdaStarBkg/h3dMixedLambdaStarMassVsPt", "h3dMixedLambdaStarMassVsPt", kTH3D, {axisConfig.axisCentrality, axisConfig.axisPt, axisConfig.axisLambdaStarMass}); histos.add("LambdaStarBkg/h3dMixedLambdaStarPtVsOPAngle", "h3dMixedLambdaStarPtVsOPAngle", kTH3D, {{140, 0.f, 7.f}, axisConfig.axisPt, axisConfig.axisLambdaStarMass}); if (doArm) { - histos.add("LambdaStarBkg/h3dMixedLambdaStarPtVsAPAlpha", "h3dMixedLambdaStarPtVsAPAlpha", kTH3D, {axisConfig.axisAPAlpha, axisConfig.axisPt, axisConfig.axisLambdaStarMass}); - histos.add("LambdaStarBkg/h3dMixedLambdaStarPtVsAPQt", "h3dMixedLambdaStarPtVsAPQt", kTH3D, {axisConfig.axisAPQt, axisConfig.axisPt, axisConfig.axisLambdaStarMass}); + histos.add("LambdaStarBkg/h4dMixedLambdaStarPtVsAPAlphaVsAPQt", "h4dMixedLambdaStarPtVsAPAlphaVsAPQt", kTHnD, {axisConfig.axisAPAlpha, axisConfig.axisAPQt, axisConfig.axisPt, axisConfig.axisLambdaStarMass}); } } } @@ -868,16 +864,14 @@ struct k892hadronphotonBkg { histos.fill(HIST("KStarBkg/h3dRotKStarMassVsPt"), centrality, reso.Pt(), reso.M()); histos.fill(HIST("KStarBkg/h3dRotKStarPtVsOPAngle"), openAngle, reso.Pt(), reso.M()); if (doArm) { - histos.fill(HIST("KStarBkg/h3dRotKStarPtVsAPAlpha"), apAlpha, reso.Pt(), reso.M()); - histos.fill(HIST("KStarBkg/h3dRotKStarPtVsAPQt"), apQt, reso.Pt(), reso.M()); + histos.fill(HIST("KStarBkg/h4dRotKStarPtVsAPAlphaVsAPQt"), apAlpha, apQt, reso.Pt(), reso.M()); } } else { histos.fill(HIST("LambdaStarBkg/h2dRotLambdaStarMassVsPt"), reso.M(), reso.Pt()); histos.fill(HIST("LambdaStarBkg/h3dRotLambdaStarMassVsPt"), centrality, reso.Pt(), reso.M()); histos.fill(HIST("LambdaStarBkg/h3dRotLambdaStarPtVsOPAngle"), openAngle, reso.Pt(), reso.M()); if (doArm) { - histos.fill(HIST("LambdaStarBkg/h3dRotLambdaStarPtVsAPAlpha"), apAlpha, reso.Pt(), reso.M()); - histos.fill(HIST("LambdaStarBkg/h3dRotLambdaStarPtVsAPQt"), apQt, reso.Pt(), reso.M()); + histos.fill(HIST("LambdaStarBkg/h4dRotLambdaStarPtVsAPAlphaVsAPQt"), apAlpha, apQt, reso.Pt(), reso.M()); } } } @@ -947,16 +941,14 @@ struct k892hadronphotonBkg { histos.fill(HIST("KStarBkg/h3dMixedKStarMassVsPt"), centrality, pt, mass); histos.fill(HIST("KStarBkg/h3dMixedKStarPtVsOPAngle"), openAngle, pt, mass); if (doArm) { - histos.fill(HIST("KStarBkg/h3dMixedKStarPtVsAPAlpha"), apAlpha, pt, mass); - histos.fill(HIST("KStarBkg/h3dMixedKStarPtVsAPQt"), apQt, pt, mass); + histos.fill(HIST("KStarBkg/h4dMixedKStarPtVsAPAlphaVsAPQt"), apAlpha, apQt, pt, mass); } } else { histos.fill(HIST("LambdaStarBkg/h2dMixedLambdaStarMassVsPt"), mass, pt); histos.fill(HIST("LambdaStarBkg/h3dMixedLambdaStarMassVsPt"), centrality, pt, mass); histos.fill(HIST("LambdaStarBkg/h3dMixedLambdaStarPtVsOPAngle"), openAngle, pt, mass); if (doArm) { - histos.fill(HIST("LambdaStarBkg/h3dMixedLambdaStarPtVsAPAlpha"), apAlpha, pt, mass); - histos.fill(HIST("LambdaStarBkg/h3dMixedLambdaStarPtVsAPQt"), apQt, pt, mass); + histos.fill(HIST("LambdaStarBkg/h4dMixedLambdaStarPtVsAPAlphaVsAPQt"), apAlpha, apQt, pt, mass); } } } diff --git a/PWGLF/Tasks/Strangeness/sigmaanalysis.cxx b/PWGLF/Tasks/Strangeness/sigmaanalysis.cxx index e4dda52dda6..d16e113b0b5 100644 --- a/PWGLF/Tasks/Strangeness/sigmaanalysis.cxx +++ b/PWGLF/Tasks/Strangeness/sigmaanalysis.cxx @@ -445,8 +445,7 @@ struct sigmaanalysis { histos.add(histodir + "/Sigma0/h3dMass", "h3dMass", kTH3D, {axisCentrality, axisPt, axisSigmaMass}); histos.add(histodir + "/Sigma0/h3dOPAngleVsMass", "h3dOPAngleVsMass", kTH3D, {{140, 0.0f, +7.0f}, axisPt, axisSigmaMass}); if (doArm) { - histos.add(histodir + "/Sigma0/h3dAlphaVsMass", "h3dAlphaVsMass", kTH3D, {{100, -5.0f, +5.0f}, axisPt, axisSigmaMass}); - histos.add(histodir + "/Sigma0/h3dQtarmVsMass", "h3dQtarmVsMass", kTH3D, {{100, -5.0f, +5.0f}, axisPt, axisSigmaMass}); + histos.add(histodir + "/Sigma0/h4dAlphaVsQtarmVsMass", "h4dAlphaVsQtarmVsMass", kTHnD, {axisAPAlpha, axisAPQt, axisPt, axisSigmaMass}); } histos.add(histodir + "/ASigma0/hMass", "hMass", kTH1D, {axisSigmaMass}); @@ -458,8 +457,7 @@ struct sigmaanalysis { histos.add(histodir + "/ASigma0/h3dMass", "h3dMass", kTH3D, {axisCentrality, axisPt, axisSigmaMass}); histos.add(histodir + "/ASigma0/h3dOPAngleVsMass", "h3dOPAngleVsMass", kTH3D, {{140, 0.0f, +7.0f}, axisPt, axisSigmaMass}); if (doArm) { - histos.add(histodir + "/ASigma0/h3dAlphaVsMass", "h3dAlphaVsMass", kTH3D, {{100, -5.0f, +5.0f}, axisPt, axisSigmaMass}); - histos.add(histodir + "/ASigma0/h3dQtarmVsMass", "h3dQtarmVsMass", kTH3D, {{100, -5.0f, +5.0f}, axisPt, axisSigmaMass}); + histos.add(histodir + "/ASigma0/h4dAlphaVsQtarmVsMass", "h4dAlphaVsQtarmVsMass", kTHnD, {axisAPAlpha, axisAPQt, axisPt, axisSigmaMass}); } // Process MC @@ -500,8 +498,7 @@ struct sigmaanalysis { histos.add(histodir + "/MC/Sigma0/h3dMass", "h3dMass", kTH3D, {axisCentrality, axisPt, axisSigmaMass}); histos.add(histodir + "/MC/Sigma0/h3dMCProcess", "h3dMCProcess", kTH3D, {{50, -0.5f, 49.5f}, axisPt, axisSigmaMass}); if (doArm) { - histos.add(histodir + "/MC/Sigma0/h3dAlphaVsMass", "h3dAlphaVsMass", kTH3D, {{100, -5.0f, +5.0f}, axisPt, axisSigmaMass}); - histos.add(histodir + "/MC/Sigma0/h3dQtarmVsMass", "h3dQtarmVsMass", kTH3D, {{100, -5.0f, +5.0f}, axisPt, axisSigmaMass}); + histos.add(histodir + "/MC/Sigma0/h4dAlphaVsQtarmVsMass", "h4dAlphaVsQtarmVsMass", kTHnD, {axisAPAlpha, axisAPQt, axisPt, axisSigmaMass}); } histos.add(histodir + "/MC/ASigma0/hPt", "hPt", kTH1D, {axisPt}); @@ -515,8 +512,7 @@ struct sigmaanalysis { histos.add(histodir + "/MC/ASigma0/h3dMass", "h3dMass", kTH3D, {axisCentrality, axisPt, axisSigmaMass}); histos.add(histodir + "/MC/ASigma0/h3dMCProcess", "h3dMCProcess", kTH3D, {{50, -0.5f, 49.5f}, axisPt, axisSigmaMass}); if (doArm) { - histos.add(histodir + "/MC/ASigma0/h3dAlphaVsMass", "h3dAlphaVsMass", kTH3D, {{100, -5.0f, +5.0f}, axisPt, axisSigmaMass}); - histos.add(histodir + "/MC/ASigma0/h3dQtarmVsMass", "h3dQtarmVsMass", kTH3D, {{100, -5.0f, +5.0f}, axisPt, axisSigmaMass}); + histos.add(histodir + "/MC/ASigma0/h4dAlphaVsQtarmVsMass", "h4dAlphaVsQtarmVsMass", kTHnD, {axisAPAlpha, axisAPQt, axisPt, axisSigmaMass}); } histos.add(histodir + "/MC/LambdaStar/h3dMCPtvsOPAngle_Sig", "h3dMCPtvsOPAngle_Sig", kTH3D, {{140, 0.f, 7.f}, axisPt, axisSigmaMass}); @@ -1219,8 +1215,7 @@ struct sigmaanalysis { histos.fill(HIST(MainDir[mode]) + HIST("/Sigma0/h3dMass"), centrality, sigma.pt(), sigma.sigma0Mass()); histos.fill(HIST(MainDir[mode]) + HIST("/Sigma0/h3dOPAngleVsMass"), sigma.opAngle(), sigma.pt(), sigma.sigma0Mass()); if (doArm) { - histos.fill(HIST(MainDir[mode]) + HIST("/Sigma0/h3dAlphaVsMass"), sigma.lStarAlpha(), sigma.pt(), sigma.sigma0Mass()); - histos.fill(HIST(MainDir[mode]) + HIST("/Sigma0/h3dQtarmVsMass"), sigma.lStarQtarm(), sigma.pt(), sigma.sigma0Mass()); + histos.fill(HIST(MainDir[mode]) + HIST("/Sigma0/h4dAlphaVsQtarmVsMass"), sigma.lStarAlpha(), sigma.lStarQtarm(), sigma.pt(), sigma.sigma0Mass()); } } else { @@ -1243,8 +1238,7 @@ struct sigmaanalysis { histos.fill(HIST(MainDir[mode]) + HIST("/ASigma0/h3dMass"), centrality, sigma.pt(), sigma.sigma0Mass()); histos.fill(HIST(MainDir[mode]) + HIST("/ASigma0/h3dOPAngleVsMass"), sigma.opAngle(), sigma.pt(), sigma.sigma0Mass()); if (doArm) { - histos.fill(HIST(MainDir[mode]) + HIST("/ASigma0/h3dAlphaVsMass"), sigma.lStarAlpha(), sigma.pt(), sigma.sigma0Mass()); - histos.fill(HIST(MainDir[mode]) + HIST("/ASigma0/h3dQtarmVsMass"), sigma.lStarQtarm(), sigma.pt(), sigma.sigma0Mass()); + histos.fill(HIST(MainDir[mode]) + HIST("/ASigma0/h4dAlphaVsQtarmVsMass"), sigma.lStarAlpha(), sigma.lStarQtarm(), sigma.pt(), sigma.sigma0Mass()); } } @@ -1303,8 +1297,7 @@ struct sigmaanalysis { histos.fill(HIST(MainDir[mode]) + HIST("/MC/Sigma0/hMass"), sigma.sigma0Mass()); histos.fill(HIST(MainDir[mode]) + HIST("/MC/Sigma0/h3dMass"), centrality, sigma.mcpt(), sigma.sigma0Mass()); if (doArm) { - histos.fill(HIST(MainDir[mode]) + HIST("/MC/Sigma0/h3dAlphaVsMass"), sigma.lStarAlpha(), sigma.mcpt(), sigma.sigma0Mass()); - histos.fill(HIST(MainDir[mode]) + HIST("/MC/Sigma0/h3dQtarmVsMass"), sigma.lStarQtarm(), sigma.mcpt(), sigma.sigma0Mass()); + histos.fill(HIST(MainDir[mode]) + HIST("/MC/Sigma0/h4dAlphaVsQtarmVsMass"), sigma.lStarAlpha(), sigma.lStarQtarm(), sigma.mcpt(), sigma.sigma0Mass()); } histos.fill(HIST(MainDir[mode]) + HIST("/MC/Sigma0/hMCProcess"), sigma.mcprocess()); @@ -1324,8 +1317,7 @@ struct sigmaanalysis { histos.fill(HIST(MainDir[mode]) + HIST("/MC/ASigma0/hMass"), sigma.sigma0Mass()); histos.fill(HIST(MainDir[mode]) + HIST("/MC/ASigma0/h3dMass"), centrality, sigma.mcpt(), sigma.sigma0Mass()); if (doArm) { - histos.fill(HIST(MainDir[mode]) + HIST("/MC/ASigma0/h3dAlphaVsMass"), sigma.lStarAlpha(), sigma.mcpt(), sigma.sigma0Mass()); - histos.fill(HIST(MainDir[mode]) + HIST("/MC/ASigma0/h3dQtarmVsMass"), sigma.lStarQtarm(), sigma.mcpt(), sigma.sigma0Mass()); + histos.fill(HIST(MainDir[mode]) + HIST("/MC/ASigma0/h4dAlphaVsQtarmVsMass"), sigma.lStarAlpha(), sigma.lStarQtarm(), sigma.mcpt(), sigma.sigma0Mass()); } histos.fill(HIST(MainDir[mode]) + HIST("/MC/ASigma0/hMCProcess"), sigma.mcprocess());