@@ -37,6 +37,11 @@ struct otfParticlePrediction {
3737 const AxisSpec axisPt{binsPt, " #it{p}_{T} (GeV/#it{c})" };
3838
3939 histos.add <TH1>(" collisions/generated" , " collisions" , kTH1D , {{1 , -0.5 , 0.5 }});
40+ histos.add <TH2>(" particles/generated/pi0" , " pi0" , kTH2D , {axisPt, axisEta});
41+ histos.add <TH2>(" particles/generated/eta" , " eta" , kTH2D , {axisPt, axisEta});
42+ histos.add <TH2>(" particles/generated/etaP" , " etaP" , kTH2D , {axisPt, axisEta});
43+ histos.add <TH2>(" particles/generated/rho" , " rho" , kTH2D , {axisPt, axisEta});
44+ histos.add <TH2>(" particles/generated/omega" , " omega" , kTH2D , {axisPt, axisEta});
4045 histos.add <TH2>(" particles/generated/phi" , " phi" , kTH2D , {axisPt, axisEta});
4146 }
4247
@@ -48,9 +53,34 @@ struct otfParticlePrediction {
4853 }
4954
5055 for (const auto & particle : mcParticles) {
51- if (std::abs (particle.pdgCode ()) == 333 ) // phi
52- {
53- histos.fill (HIST (" particles/generated/phi" ), particle.pt (), particle.eta ());
56+ auto pdg = std::abs (particle.pdgCode ());
57+ if (pdg < 100 ) {
58+ continue ;
59+ }
60+ if (pdg > 1000 ) {
61+ continue ;
62+ }
63+ switch (pdg) {
64+ case 111 :
65+ histos.fill (HIST (" particles/generated/pi0" ), particle.pt (), particle.eta ());
66+ break ;
67+ case 221 :
68+ histos.fill (HIST (" particles/generated/eta" ), particle.pt (), particle.eta ());
69+ break ;
70+ case 331 :
71+ histos.fill (HIST (" particles/generated/etaP" ), particle.pt (), particle.eta ());
72+ break ;
73+ case 223 :
74+ histos.fill (HIST (" particles/generated/omega" ), particle.pt (), particle.eta ());
75+ break ;
76+ case 113 :
77+ histos.fill (HIST (" particles/generated/rho" ), particle.pt (), particle.eta ());
78+ break ;
79+ case 333 :
80+ histos.fill (HIST (" particles/generated/phi" ), particle.pt (), particle.eta ());
81+ break ;
82+ default :
83+ break ;
5484 }
5585 }
5686 }
0 commit comments