Skip to content

Commit ebe0fcf

Browse files
nepeivodaRSalibuild
nepeivodaRS
and
alibuild
authored
[PWGLF] add rapidity and gap axes to histograms in UPC analysis (#8575)
Co-authored-by: ALICE Action Bot <alibuild@cern.ch>
1 parent dd6d40c commit ebe0fcf

File tree

1 file changed

+37
-37
lines changed

1 file changed

+37
-37
lines changed

‎PWGLF/Tasks/Strangeness/strange-yield-pbpb.cxx‎

Lines changed: 37 additions & 37 deletions
Original file line numberDiff line numberDiff line change
@@ -350,11 +350,11 @@ struct strangeYieldPbPb {
350350
void addKinematicQAHistograms(HistogramRegistry& histos)
351351
{
352352
const bool isCascade = (partID > 2.5) ? true : false;
353-
histos.add(Form("%s/h2dPosEtaPt", particlenames[partID].data()), "h2dPosEtaPt", kTH3F, {axisPtCoarse, axisEta, axisSelGap});
354-
histos.add(Form("%s/h2dNegEtaPt", particlenames[partID].data()), "h2dNegEtaPt", kTH3F, {axisPtCoarse, axisEta, axisSelGap});
355-
histos.add(Form("%s/h2dRapPt", particlenames[partID].data()), "h2dRapPt", kTH3F, {axisPtCoarse, axisRap, axisSelGap});
353+
histos.add(Form("%s/h3dPosEtaPt", particlenames[partID].data()), "h3dPosEtaPt", kTH3F, {axisPtCoarse, axisEta, axisSelGap});
354+
histos.add(Form("%s/h3dNegEtaPt", particlenames[partID].data()), "h3dNegEtaPt", kTH3F, {axisPtCoarse, axisEta, axisSelGap});
355+
histos.add(Form("%s/h3dRapPt", particlenames[partID].data()), "h3dRapPt", kTH3F, {axisPtCoarse, axisRap, axisSelGap});
356356
if (isCascade) {
357-
histos.add(Form("%s/h2dBachEtaPt", particlenames[partID].data()), "h2dBachEtaPt", kTH3F, {axisPtCoarse, axisEta, axisSelGap});
357+
histos.add(Form("%s/h3dBachEtaPt", particlenames[partID].data()), "h3dBachEtaPt", kTH3F, {axisPtCoarse, axisEta, axisSelGap});
358358
}
359359
}
360360

@@ -387,7 +387,7 @@ struct strangeYieldPbPb {
387387
template <int partID>
388388
void addHistograms(HistogramRegistry& histos)
389389
{
390-
histos.add(Form("%s/h5dMass", particlenames[partID].data()), "h5dMass", kTHnSparseF, {axisFT0C, axisPt, axisInvMass.at(partID), axisSelGap, axisNchInvMass});
390+
histos.add(Form("%s/h6dMass", particlenames[partID].data()), "h6dMass", kTHnSparseF, {axisFT0C, axisPt, axisInvMass.at(partID), axisSelGap, axisNchInvMass, axisRap});
391391
histos.add(Form("%s/h2dMass", particlenames[partID].data()), "h2dMass", kTH2F, {axisInvMass.at(partID), axisSelGap});
392392
if (doPlainTopoQA) {
393393
addTopoHistograms<partID>(histos);
@@ -474,11 +474,11 @@ struct strangeYieldPbPb {
474474
}
475475

476476
histos.fill(HIST(particlenames[partID]) + HIST("/h2dMass"), invMass, gap);
477-
histos.fill(HIST(particlenames[partID]) + HIST("/h5dMass"), centrality, pT, invMass, gap, coll.multAllTracksTPCOnly() + coll.multAllTracksITSTPC());
477+
histos.fill(HIST(particlenames[partID]) + HIST("/h6dMass"), centrality, pT, invMass, gap, coll.multAllTracksTPCOnly() + coll.multAllTracksITSTPC(), rapidity);
478478
if (doKienmaticQA) {
479-
histos.fill(HIST(particlenames[partID]) + HIST("/h2dPosEtaPt"), pT, cand.positiveeta(), gap);
480-
histos.fill(HIST(particlenames[partID]) + HIST("/h2dNegEtaPt"), pT, cand.negativeeta(), gap);
481-
histos.fill(HIST(particlenames[partID]) + HIST("/h2dRapPt"), pT, rapidity, gap);
479+
histos.fill(HIST(particlenames[partID]) + HIST("/h3dPosEtaPt"), pT, cand.positiveeta(), gap);
480+
histos.fill(HIST(particlenames[partID]) + HIST("/h3dNegEtaPt"), pT, cand.negativeeta(), gap);
481+
histos.fill(HIST(particlenames[partID]) + HIST("/h3dRapPt"), pT, rapidity, gap);
482482
}
483483
if (doPlainTopoQA) {
484484
histos.fill(HIST(particlenames[partID]) + HIST("/hPosDCAToPV"), cand.dcapostopv());
@@ -624,12 +624,12 @@ struct strangeYieldPbPb {
624624
}
625625
}
626626
histos.fill(HIST(particlenames[partID]) + HIST("/h2dMass"), invMass, gap);
627-
histos.fill(HIST(particlenames[partID]) + HIST("/h5dMass"), centrality, pT, invMass, gap, coll.multAllTracksTPCOnly() + coll.multAllTracksITSTPC());
627+
histos.fill(HIST(particlenames[partID]) + HIST("/h6dMass"), centrality, pT, invMass, gap, coll.multAllTracksTPCOnly() + coll.multAllTracksITSTPC(), rapidity);
628628
if (doKienmaticQA) {
629-
histos.fill(HIST(particlenames[partID]) + HIST("/h2dPosEtaPt"), pT, cand.positiveeta(), gap);
630-
histos.fill(HIST(particlenames[partID]) + HIST("/h2dNegEtaPt"), pT, cand.negativeeta(), gap);
631-
histos.fill(HIST(particlenames[partID]) + HIST("/h2dBachEtaPt"), pT, cand.bacheloreta(), gap);
632-
histos.fill(HIST(particlenames[partID]) + HIST("/h2dRapPt"), pT, rapidity, gap);
629+
histos.fill(HIST(particlenames[partID]) + HIST("/h3dPosEtaPt"), pT, cand.positiveeta(), gap);
630+
histos.fill(HIST(particlenames[partID]) + HIST("/h3dNegEtaPt"), pT, cand.negativeeta(), gap);
631+
histos.fill(HIST(particlenames[partID]) + HIST("/h3dBachEtaPt"), pT, cand.bacheloreta(), gap);
632+
histos.fill(HIST(particlenames[partID]) + HIST("/h3dRapPt"), pT, rapidity, gap);
633633
}
634634
if (doPlainTopoQA) {
635635
histos.fill(HIST(particlenames[partID]) + HIST("/hCascCosPA"), pT, cand.casccosPA(coll.posX(), coll.posY(), coll.posZ()));
@@ -842,18 +842,18 @@ struct strangeYieldPbPb {
842842
histos.get<TH1>(HIST("hEventSelection"))->GetXaxis()->SetBinLabel(16, "has UPC flag");
843843

844844
// Event QA
845-
histos.add("eventQA/hCentrality", "hCentrality", kTH1F, {axisFT0C_QA});
846-
histos.add("eventQA/hCentralityVsTracksPVeta1", "hCentralityVsTracksPVeta1", kTH2F, {axisFT0C_QA, axisNTracksPVeta1});
847-
histos.add("eventQA/hCentralityVsTracksTotalExceptITSonly", "hCentralityVsTracksTotalExceptITSonly", kTH2F, {axisFT0C_QA, axisNTracksTotalExceptITSonly});
848-
histos.add("eventQA/hOccupancy", "hOccupancy", kTH1F, {axisOccupancy});
849-
histos.add("eventQA/hCentralityVsOccupancy", "hCentralityVsOccupancy", kTH2F, {axisFT0C_QA, axisOccupancy});
850-
histos.add("eventQA/hTracksPVeta1VsTracksGlobal", "hTracksPVeta1VsTracksGlobal", kTH2F, {axisNTracksPVeta1, axisNTracksGlobal});
851-
histos.add("eventQA/hCentralityVsTracksGlobal", "hCentralityVsTracksGlobal", kTH2F, {axisFT0C_QA, axisNTracksGlobal});
845+
histos.add("eventQA/hCentrality", "hCentrality", kTH2F, {axisFT0C_QA, axisSelGap});
846+
histos.add("eventQA/hCentralityVsTracksPVeta1", "hCentralityVsTracksPVeta1", kTH3F, {axisFT0C_QA, axisNTracksPVeta1, axisSelGap});
847+
histos.add("eventQA/hCentralityVsTracksTotalExceptITSonly", "hCentralityVsTracksTotalExceptITSonly", kTH3F, {axisFT0C_QA, axisNTracksTotalExceptITSonly, axisSelGap});
848+
histos.add("eventQA/hOccupancy", "hOccupancy", kTH2F, {axisOccupancy, axisSelGap});
849+
histos.add("eventQA/hCentralityVsOccupancy", "hCentralityVsOccupancy", kTH3F, {axisFT0C_QA, axisOccupancy, axisSelGap});
850+
histos.add("eventQA/hTracksPVeta1VsTracksGlobal", "hTracksPVeta1VsTracksGlobal", kTH3F, {axisNTracksPVeta1, axisNTracksGlobal, axisSelGap});
851+
histos.add("eventQA/hCentralityVsTracksGlobal", "hCentralityVsTracksGlobal", kTH3F, {axisFT0C_QA, axisNTracksGlobal, axisSelGap});
852852
histos.add("eventQA/hGapSide", "Gap side; Entries", kTH1F, {{5, -0.5, 4.5}});
853853
histos.add("eventQA/hSelGapSide", "Selected gap side; Entries", kTH1F, {axisSelGap});
854-
histos.add("eventQA/hPosX", "Vertex position in x", kTH1F, {{100, -0.1, 0.1}});
855-
histos.add("eventQA/hPosY", "Vertex position in y", kTH1F, {{100, -0.1, 0.1}});
856-
histos.add("eventQA/hPosZ", "Vertex position in z", kTH1F, {{100, -20., 20.}});
854+
histos.add("eventQA/hPosX", "Vertex position in x", kTH2F, {{100, -0.1, 0.1}, axisSelGap});
855+
histos.add("eventQA/hPosY", "Vertex position in y", kTH2F, {{100, -0.1, 0.1}, axisSelGap});
856+
histos.add("eventQA/hPosZ", "Vertex position in z", kTH2F, {{100, -20., 20.}, axisSelGap});
857857
histos.add("eventQA/hFT0", "hFT0", kTH3F, {axisDetectors.axisFT0Aampl, axisDetectors.axisFT0Campl, axisSelGap});
858858
histos.add("eventQA/hFDD", "hFDD", kTH3F, {axisDetectors.axisFDDAampl, axisDetectors.axisFDDCampl, axisSelGap});
859859
histos.add("eventQA/hZN", "hZN", kTH3F, {axisDetectors.axisZNAampl, axisDetectors.axisZNCampl, axisSelGap});
@@ -954,20 +954,20 @@ struct strangeYieldPbPb {
954954
}
955955

956956
template <typename TCollision>
957-
void fillHistogramsQA(TCollision const& collision)
957+
void fillHistogramsQA(TCollision const& collision, int const& gap)
958958
{
959959
// QA histograms
960960
float centrality = collision.centFT0C();
961-
histos.fill(HIST("eventQA/hCentrality"), centrality);
962-
histos.fill(HIST("eventQA/hCentralityVsTracksTotalExceptITSonly"), centrality, collision.multAllTracksTPCOnly() + collision.multAllTracksITSTPC());
963-
histos.fill(HIST("eventQA/hCentralityVsTracksPVeta1"), centrality, collision.multNTracksPVeta1());
964-
histos.fill(HIST("eventQA/hOccupancy"), collision.trackOccupancyInTimeRange());
965-
histos.fill(HIST("eventQA/hCentralityVsOccupancy"), centrality, collision.trackOccupancyInTimeRange());
966-
histos.fill(HIST("eventQA/hTracksPVeta1VsTracksGlobal"), collision.multNTracksPVeta1(), collision.multNTracksGlobal());
967-
histos.fill(HIST("eventQA/hCentralityVsTracksGlobal"), centrality, collision.multNTracksGlobal());
968-
histos.fill(HIST("eventQA/hPosX"), collision.posX());
969-
histos.fill(HIST("eventQA/hPosY"), collision.posY());
970-
histos.fill(HIST("eventQA/hPosZ"), collision.posZ());
961+
histos.fill(HIST("eventQA/hCentrality"), centrality, gap);
962+
histos.fill(HIST("eventQA/hCentralityVsTracksTotalExceptITSonly"), centrality, collision.multAllTracksTPCOnly() + collision.multAllTracksITSTPC(), gap);
963+
histos.fill(HIST("eventQA/hCentralityVsTracksPVeta1"), centrality, collision.multNTracksPVeta1(), gap);
964+
histos.fill(HIST("eventQA/hOccupancy"), collision.trackOccupancyInTimeRange(), gap);
965+
histos.fill(HIST("eventQA/hCentralityVsOccupancy"), centrality, collision.trackOccupancyInTimeRange(), gap);
966+
histos.fill(HIST("eventQA/hTracksPVeta1VsTracksGlobal"), collision.multNTracksPVeta1(), collision.multNTracksGlobal(), gap);
967+
histos.fill(HIST("eventQA/hCentralityVsTracksGlobal"), centrality, collision.multNTracksGlobal(), gap);
968+
histos.fill(HIST("eventQA/hPosX"), collision.posX(), gap);
969+
histos.fill(HIST("eventQA/hPosY"), collision.posY(), gap);
970+
histos.fill(HIST("eventQA/hPosZ"), collision.posZ(), gap);
971971
}
972972

973973
template <typename TCollision>
@@ -1523,7 +1523,7 @@ struct strangeYieldPbPb {
15231523
if (studyUPConly && (selGapSide < -0.5))
15241524
return;
15251525

1526-
fillHistogramsQA(collision);
1526+
fillHistogramsQA(collision, selGapSide);
15271527

15281528
for (auto& v0 : fullV0s) {
15291529
if ((v0.v0Type() != v0cuts.v0TypeSelection) && (v0cuts.v0TypeSelection > 0))
@@ -1549,7 +1549,7 @@ struct strangeYieldPbPb {
15491549
if (studyUPConly && (selGapSide < -0.5))
15501550
return;
15511551

1552-
fillHistogramsQA(collision);
1552+
fillHistogramsQA(collision, selGapSide);
15531553

15541554
for (auto& casc : fullCascades) {
15551555
std::bitset<selNum> selMap = computeBitmapCascade(casc, collision);

0 commit comments

Comments
 (0)