Skip to content

Commit 629ee67

Browse files
authored
[PWGLF] Update gen process function (#11239)
1 parent 2fc2e04 commit 629ee67

File tree

1 file changed

+30
-10
lines changed

1 file changed

+30
-10
lines changed

PWGLF/Tasks/GlobalEventProperties/heavyionMultiplicity.cxx

Lines changed: 30 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -806,22 +806,42 @@ struct HeavyionMultiplicity {
806806
}
807807
PROCESS_SWITCH(HeavyionMultiplicity, processppMonteCarlo, "process pp MC", false);
808808

809-
void processGen(soa::Join<aod::McCollisions, aod::MultsExtraMC>::iterator const& mccols, aod::McParticles const& GenParticles)
809+
void processGen(aod::McCollisions::iterator const&, aod::McParticles const& GenParticles)
810810
{
811-
histos.fill(HIST("MultBarrelEta10_vs_FT0A"), mccols.multMCNParticlesEta10(), mccols.multMCFT0A());
812-
histos.fill(HIST("MultBarrelEta10_vs_FT0C"), mccols.multMCNParticlesEta10(), mccols.multMCFT0C());
813-
histos.fill(HIST("MultBarrelEta10"), mccols.multMCNParticlesEta10());
814-
histos.fill(HIST("MultFT0A"), mccols.multMCFT0A());
815-
histos.fill(HIST("MultFT0C"), mccols.multMCFT0C());
816-
histos.fill(HIST("mult10_vs_FT0A"), mccols.multMCNParticlesEta10(), mccols.multMCFT0A());
817-
histos.fill(HIST("mult10_vs_FT0C"), mccols.multMCNParticlesEta10(), mccols.multMCFT0C());
811+
812+
int multFT0A = 0;
813+
int multFT0C = 0;
814+
int multBarrelEta10 = 0;
815+
816+
for (const auto& particle : GenParticles) {
817+
if (!isGenTrackSelected(particle)) {
818+
continue;
819+
}
820+
if (std::abs(particle.eta()) < 1.0) {
821+
multBarrelEta10++;
822+
}
823+
if (-3.3 < particle.eta() && particle.eta() < -2.1) {
824+
multFT0C++;
825+
}
826+
if (3.5 < particle.eta() && particle.eta() < 4.9) {
827+
multFT0A++;
828+
}
829+
}
830+
831+
histos.fill(HIST("MultBarrelEta10_vs_FT0A"), multBarrelEta10, multFT0A);
832+
histos.fill(HIST("MultBarrelEta10_vs_FT0C"), multBarrelEta10, multFT0C);
833+
histos.fill(HIST("MultBarrelEta10"), multBarrelEta10);
834+
histos.fill(HIST("MultFT0A"), multFT0A);
835+
histos.fill(HIST("MultFT0C"), multFT0C);
836+
histos.fill(HIST("mult10_vs_FT0A"), multBarrelEta10, multFT0A);
837+
histos.fill(HIST("mult10_vs_FT0C"), multBarrelEta10, multFT0C);
818838

819839
for (const auto& particle : GenParticles) {
820840
if (!isGenTrackSelected(particle)) {
821841
continue;
822842
}
823-
histos.fill(HIST("dndeta10_vs_FT0A"), particle.eta(), mccols.multMCFT0A());
824-
histos.fill(HIST("dndeta10_vs_FT0C"), particle.eta(), mccols.multMCFT0C());
843+
histos.fill(HIST("dndeta10_vs_FT0A"), particle.eta(), multFT0A);
844+
histos.fill(HIST("dndeta10_vs_FT0C"), particle.eta(), multFT0C);
825845
}
826846
}
827847
PROCESS_SWITCH(HeavyionMultiplicity, processGen, "process pure MC gen", false);

0 commit comments

Comments
 (0)