From f4912b9b2a1cd0f891cde2d4e62ffa9f08d4ab5d Mon Sep 17 00:00:00 2001 From: murioni Date: Fri, 20 Mar 2026 16:58:33 +0100 Subject: [PATCH 1/2] adding FT0M multiplicity at generated level, check on split events --- .../heavyionMultiplicity.cxx | 16 ++++++++++++++-- 1 file changed, 14 insertions(+), 2 deletions(-) diff --git a/PWGLF/Tasks/GlobalEventProperties/heavyionMultiplicity.cxx b/PWGLF/Tasks/GlobalEventProperties/heavyionMultiplicity.cxx index 5c820ca6676..dc1890d92b1 100644 --- a/PWGLF/Tasks/GlobalEventProperties/heavyionMultiplicity.cxx +++ b/PWGLF/Tasks/GlobalEventProperties/heavyionMultiplicity.cxx @@ -343,7 +343,11 @@ struct HeavyionMultiplicity { histos.add("hMultEta05vsCentrRec", "multiplicity in eta<0.5 of selected MC events vs centrality", kTH2F, {axisCent, multAxis}); histos.add("hgendndetaVsMultEta05BeforeEvtSel", "hgendndetaBeforeEvtSel vs multiplicity in eta<0.5", kTH2F, {axisEta, multAxis}); histos.add("hgendndetaVsMultEta05AfterEvtSel", "hgendndetaAfterEvtSel vs multiplicity in eta<0.5", kTH2F, {axisEta, multAxis}); - + histos.add("hImpactParameterSplit", "Impact parameter of selected and split MC events", kTH1F,{impactParAxis}); + histos.add("hMultEta05Split" , "multiplicity in eta<0.5 of selected and split MC events", kTH1F,{multAxis}); + histos.add("hMultSplit" , "multiplicity of selected and split MC events", kTH1F,{axisFt0cMult}); + histos.add("hMultvsCentrSplit" , "multiplicity of selected and split MC events vs centrality ", kTH2F,{axisCent, axisFt0cMult}); + histos.add("hMultGen", "multiplicity of generated MC events", kTH1F, {axisFt0cMult}); histos.add("hMultRec", "multiplicity of selected MC events", kTH1F, {axisFt0cMult}); histos.add("hMultvsCentrRec", "multiplicity of selected MC events vs centrality", kTH2F, {axisCent, axisFt0cMult}); @@ -463,7 +467,9 @@ struct HeavyionMultiplicity { if (isApplyCentFT0C) { cent = col.multMCFT0C(); } - if (isApplyCentFV0A) { + else if(isApplyCentFT0M) { + cent = (col.multMCFT0C() + col.multMCFT0A())/2.; + } else if (isApplyCentFV0A) { cent = col.multMCFV0A(); } return cent; @@ -990,6 +996,12 @@ struct HeavyionMultiplicity { if (std::abs(RecCol.posZ()) >= vtxRange) { continue; } + + histos.fill(HIST("hImpactParameterSplit"), mcCollision.impactParameter()); + histos.fill(HIST("hMultEta05Split"), mcCollision.multMCNParticlesEta05()); + histos.fill(HIST("hMultSplit"), selColMultMC(mcCollision)); + histos.fill(HIST("hMultvsCentrSplit"), selColCent(RecCol), selColMultMC(mcCollision)); + if (RecCol.numContrib() <= numcontributors) { continue; } else { From 084ecddabd0aaad33ddd134062723490fc9a6c2d Mon Sep 17 00:00:00 2001 From: ALICE Action Bot Date: Fri, 20 Mar 2026 16:02:48 +0000 Subject: [PATCH 2/2] Please consider the following formatting changes --- .../heavyionMultiplicity.cxx | 17 ++++++++--------- 1 file changed, 8 insertions(+), 9 deletions(-) diff --git a/PWGLF/Tasks/GlobalEventProperties/heavyionMultiplicity.cxx b/PWGLF/Tasks/GlobalEventProperties/heavyionMultiplicity.cxx index dc1890d92b1..d6f8a6df09f 100644 --- a/PWGLF/Tasks/GlobalEventProperties/heavyionMultiplicity.cxx +++ b/PWGLF/Tasks/GlobalEventProperties/heavyionMultiplicity.cxx @@ -343,11 +343,11 @@ struct HeavyionMultiplicity { histos.add("hMultEta05vsCentrRec", "multiplicity in eta<0.5 of selected MC events vs centrality", kTH2F, {axisCent, multAxis}); histos.add("hgendndetaVsMultEta05BeforeEvtSel", "hgendndetaBeforeEvtSel vs multiplicity in eta<0.5", kTH2F, {axisEta, multAxis}); histos.add("hgendndetaVsMultEta05AfterEvtSel", "hgendndetaAfterEvtSel vs multiplicity in eta<0.5", kTH2F, {axisEta, multAxis}); - histos.add("hImpactParameterSplit", "Impact parameter of selected and split MC events", kTH1F,{impactParAxis}); - histos.add("hMultEta05Split" , "multiplicity in eta<0.5 of selected and split MC events", kTH1F,{multAxis}); - histos.add("hMultSplit" , "multiplicity of selected and split MC events", kTH1F,{axisFt0cMult}); - histos.add("hMultvsCentrSplit" , "multiplicity of selected and split MC events vs centrality ", kTH2F,{axisCent, axisFt0cMult}); - + histos.add("hImpactParameterSplit", "Impact parameter of selected and split MC events", kTH1F, {impactParAxis}); + histos.add("hMultEta05Split", "multiplicity in eta<0.5 of selected and split MC events", kTH1F, {multAxis}); + histos.add("hMultSplit", "multiplicity of selected and split MC events", kTH1F, {axisFt0cMult}); + histos.add("hMultvsCentrSplit", "multiplicity of selected and split MC events vs centrality ", kTH2F, {axisCent, axisFt0cMult}); + histos.add("hMultGen", "multiplicity of generated MC events", kTH1F, {axisFt0cMult}); histos.add("hMultRec", "multiplicity of selected MC events", kTH1F, {axisFt0cMult}); histos.add("hMultvsCentrRec", "multiplicity of selected MC events vs centrality", kTH2F, {axisCent, axisFt0cMult}); @@ -466,9 +466,8 @@ struct HeavyionMultiplicity { auto cent = -1; if (isApplyCentFT0C) { cent = col.multMCFT0C(); - } - else if(isApplyCentFT0M) { - cent = (col.multMCFT0C() + col.multMCFT0A())/2.; + } else if (isApplyCentFT0M) { + cent = (col.multMCFT0C() + col.multMCFT0A()) / 2.; } else if (isApplyCentFV0A) { cent = col.multMCFV0A(); } @@ -996,7 +995,7 @@ struct HeavyionMultiplicity { if (std::abs(RecCol.posZ()) >= vtxRange) { continue; } - + histos.fill(HIST("hImpactParameterSplit"), mcCollision.impactParameter()); histos.fill(HIST("hMultEta05Split"), mcCollision.multMCNParticlesEta05()); histos.fill(HIST("hMultSplit"), selColMultMC(mcCollision));