Skip to content

Commit 901d207

Browse files
committed
adding QA plots
1 parent 4d7837f commit 901d207

File tree

1 file changed

+9
-0
lines changed

1 file changed

+9
-0
lines changed

PWGLF/Tasks/GlobalEventProperties/heavyionMultiplicity.cxx

Lines changed: 9 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -131,6 +131,7 @@ static constexpr TrackSelectionFlags::flagtype TrackSelectionDcaxyOnly =
131131
AxisSpec axisEvent{10, 0.5, 10.5, "#Event", "EventAxis"};
132132
AxisSpec axisVtxZ{40, -20, 20, "Vertex Z", "VzAxis"};
133133
AxisSpec axisEta{40, -2, 2, "#eta", "EtaAxis"};
134+
AxisSpec axisEtaExtended{100, -5, 5, "#eta", "EtaAxisExtended"};
134135
AxisSpec axisPhi{{0, o2::constants::math::PIQuarter, o2::constants::math::PIHalf, o2::constants::math::PIQuarter * 3., o2::constants::math::PI, o2::constants::math::PIQuarter * 5., o2::constants::math::PIHalf * 3., o2::constants::math::PIQuarter * 7., o2::constants::math::TwoPI}, "#phi", "PhiAxis"};
135136
AxisSpec axisPhi2{629, 0, o2::constants::math::TwoPI, "#phi"};
136137
AxisSpec axisCent{100, 0, 100, "#Cent"};
@@ -313,6 +314,8 @@ struct HeavyionMultiplicity {
313314
x->SetBinLabel(1, "All MC events");
314315
x->SetBinLabel(2, "MC events with reco event after event selection");
315316
x->SetBinLabel(3, "MC events with no reco events");
317+
histos.add("hgendndetaVscentGenwithNOreco", "dndeta vs impact parameter, gen events with no reco", kTH2F, {axisEtaExtended, impactParAxis});
318+
histos.add("hgendndetaVscentGenwithReco" , "dndeta vs impact parameter, gen events with at least one reco", kTH2F, {axisEtaExtended, impactParAxis});
316319
histos.add("hImpactParameterGenwithNoreco", "Impact parameter of generated MC events, with no recoevent", kTH1F, {impactParAxis});
317320
histos.add("hImpactParameterGen", "Impact parameter of generated MC events", kTH1F, {impactParAxis});
318321
histos.add("hImpactParameterRec", "Impact parameter of selected MC events", kTH1F, {impactParAxis});
@@ -949,6 +952,12 @@ struct HeavyionMultiplicity {
949952
}
950953

951954
for (const auto& particle : GenParticles) {
955+
956+
if (RecCols.size() == 0) {
957+
histos.fill(HIST("hgendndetaVscentGenwithNOreco"), particle.eta(), mcCollision.impactParameter());
958+
} else {
959+
histos.fill(HIST("hgendndetaVscentGenwithReco"), particle.eta(), mcCollision.impactParameter());
960+
}
952961

953962
if (!isGenTrackSelected(particle)) {
954963
continue;

0 commit comments

Comments
 (0)