diff --git a/PWGLF/Tasks/GlobalEventProperties/heavyionMultiplicity.cxx b/PWGLF/Tasks/GlobalEventProperties/heavyionMultiplicity.cxx index 7c91e73f9a9..9a496f78016 100644 --- a/PWGLF/Tasks/GlobalEventProperties/heavyionMultiplicity.cxx +++ b/PWGLF/Tasks/GlobalEventProperties/heavyionMultiplicity.cxx @@ -311,7 +311,9 @@ struct HeavyionMultiplicity { auto hstat = histos.get(HIST("MCEventHist")); auto* x = hstat->GetXaxis(); x->SetBinLabel(1, "All MC events"); - x->SetBinLabel(2, "MC events with atleast one reco event"); + x->SetBinLabel(2, "MC events with reco event after event selection"); + x->SetBinLabel(3, "MC events with no reco events"); + histos.add("hImpactParameterGenwithNoreco", "Impact parameter of generated MC events, with no recoevent", kTH1F, {impactParAxis}); histos.add("hImpactParameterGen", "Impact parameter of generated MC events", kTH1F, {impactParAxis}); histos.add("hImpactParameterRec", "Impact parameter of selected MC events", kTH1F, {impactParAxis}); histos.add("hImpactParvsCentrRec", "Impact parameter of selected MC events vs centrality", kTH2F, {axisCent, impactParAxis}); @@ -915,6 +917,11 @@ struct HeavyionMultiplicity { histos.fill(HIST("MCEventHist"), 1); histos.fill(HIST("hImpactParameterGen"), mcCollision.impactParameter()); + if (RecCols.size() == 0) { + histos.fill(HIST("MCEventHist"), 3); + histos.fill(HIST("hImpactParameterGenwithNoreco"), mcCollision.impactParameter()); + } + bool atLeastOne = false; auto centrality = -999.; auto numcontributors = -999;