Skip to content

PWGHF: add information about Zvtx residual in histos for collision association probability #1046

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
vkucera merged 6 commits into from
Jul 22, 2022
Merged

PWGHF: add information about Zvtx residual in histos for collision association probability #1046

Changes from all commits
Commits
Show all changes
6 commits
Select commit Hold shift + click to select a range
177c537
PWGHF: add information about Zvtx residual in histos for collision as…
fgrosa Jul 22, 2022
65eb3df
Remove comment
fgrosa Jul 22, 2022
ea4a5ea
Add track selection
fgrosa Jul 22, 2022
438c67b
Fix track selection
fgrosa Jul 22, 2022
c01d33a
Add axis title
fgrosa Jul 22, 2022
edd885b
Add track eta
fgrosa Jul 22, 2022
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
42 changes: 24 additions & 18 deletions PWGHF/Tasks/HFMCValidation.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -273,20 +273,20 @@ struct ValidationRecLevel {
std::array<std::shared_ptr<TH1>, nCharmHadrons> histDeltaPt, histDeltaPx, histDeltaPy, histDeltaPz, histDeltaSecondaryVertexX, histDeltaSecondaryVertexY, histDeltaSecondaryVertexZ, histDeltaDecayLength;
std::array<std::array<std::array<std::shared_ptr<TH1>, 3>, 2>, nCharmHadrons> histPtDau, histEtaDau, histImpactParameterDau;
std::array<std::array<std::shared_ptr<TH1>, 2>, nCharmHadrons> histPtReco;
std::array<std::shared_ptr<TH2>, 4> histOriginTracks;
std::array<std::shared_ptr<THnSparse>, 4> histOriginTracks;

HistogramRegistry registry{"registry", {}};
void init(o2::framework::InitContext&)
{
histOriginTracks[0] = registry.add<TH2>("histOriginNonAssociatedTracks", ";origin;#it{p}_{T}^{reco} (GeV/#it{c})", HistType::kTH2F, {{4, -1.5, 2.5}, {50, 0., 10.}}); // tracks not associated to any collision
histOriginTracks[1] = registry.add<TH2>("histOriginAssociatedTracks", ";origin;#it{p}_{T}^{reco} (GeV/#it{c})", HistType::kTH2F, {{4, -1.5, 2.5}, {50, 0., 10.}}); // tracks associasted to a collision
histOriginTracks[2] = registry.add<TH2>("histOriginGoodAssociatedTracks", ";origin;#it{p}_{T}^{reco} (GeV/#it{c})", HistType::kTH2F, {{4, -1.5, 2.5}, {50, 0., 10.}}); // tracks associated to the correct collision (considering the MC collision index)
histOriginTracks[3] = registry.add<TH2>("histOriginWrongAssociatedTracks", ";origin;#it{p}_{T}^{reco} (GeV/#it{c})", HistType::kTH2F, {{4, -1.5, 2.5}, {50, 0., 10.}}); // tracks associated to the wrong collision (considering the MC collision index)
histOriginTracks[0] = registry.add<THnSparse>("histOriginNonAssociatedTracks", ";origin;#it{p}_{T}^{reco} (GeV/#it{c});#it{#eta}^{reco};#it{Z}_{vtx}^{reco}#minus#it{Z}_{vtx}^{gen} (cm)", HistType::kTHnSparseF, {{4, -1.5, 2.5}, {50, 0., 10.}, {40, -1., 1.}, {200, -1., 1.}}); // tracks not associated to any collision
histOriginTracks[1] = registry.add<THnSparse>("histOriginAssociatedTracks", ";origin;#it{p}_{T}^{reco} (GeV/#it{c});#it{#eta}^{reco};#it{Z}_{vtx}^{reco}#minus#it{Z}_{vtx}^{gen} (cm)", HistType::kTHnSparseF, {{4, -1.5, 2.5}, {50, 0., 10.}, {40, -1., 1.}, {200, -1., 1.}}); // tracks associasted to a collision
histOriginTracks[2] = registry.add<THnSparse>("histOriginGoodAssociatedTracks", ";origin;#it{p}_{T}^{reco} (GeV/#it{c});#it{#eta}^{reco};#it{Z}_{vtx}^{reco}#minus#it{Z}_{vtx}^{gen} (cm)", HistType::kTHnSparseF, {{4, -1.5, 2.5}, {50, 0., 10.}, {40, -1., 1.}, {200, -1., 1.}}); // tracks associated to the correct collision (considering the MC collision index)
histOriginTracks[3] = registry.add<THnSparse>("histOriginWrongAssociatedTracks", ";origin;#it{p}_{T}^{reco} (GeV/#it{c});#it{#eta}^{reco};#it{Z}_{vtx}^{reco}#minus#it{Z}_{vtx}^{gen} (cm)", HistType::kTHnSparseF, {{4, -1.5, 2.5}, {50, 0., 10.}, {40, -1., 1.}, {200, -1., 1.}}); // tracks associated to the wrong collision (considering the MC collision index)
for (std::size_t iHist{0}; iHist < histOriginTracks.size(); ++iHist) {
histOriginTracks[iHist]->GetXaxis()->SetBinLabel(1, "no MC particle");
histOriginTracks[iHist]->GetXaxis()->SetBinLabel(2, "no quark");
histOriginTracks[iHist]->GetXaxis()->SetBinLabel(3, "charm");
histOriginTracks[iHist]->GetXaxis()->SetBinLabel(4, "beauty");
histOriginTracks[iHist]->GetAxis(0)->SetBinLabel(1, "no MC particle");
histOriginTracks[iHist]->GetAxis(0)->SetBinLabel(2, "no quark");
histOriginTracks[iHist]->GetAxis(0)->SetBinLabel(3, "charm");
histOriginTracks[iHist]->GetAxis(0)->SetBinLabel(4, "beauty");
}
for (auto iHad = 0; iHad < nCharmHadrons; ++iHad) {
histDeltaPt[iHad] = registry.add<TH1>(Form("histDeltaPt%s", particleNames[iHad].data()), Form("Pt difference reco - MC %s; #it{p}_{T}^{reco} - #it{p}_{T}^{gen} (GeV/#it{c}); entries", labels[iHad].data()), HistType::kTH1F, {{2000, -1., 1.}});
Expand All @@ -311,27 +311,33 @@ struct ValidationRecLevel {
using HfCandProng2WithMCRec = soa::Join<aod::HfCandProng2, aod::HfCandProng2MCRec>;
using HfCandProng3WithMCRec = soa::Join<aod::HfCandProng3, aod::HfCandProng3MCRec>;
using CollisionsWithMCLabels = soa::Join<aod::Collisions, aod::McCollisionLabels>;
using TracksWithSel = soa::Join<aod::BigTracksMC, aod::TrackSelection>;

void process(HfCandProng2WithMCRec const& cand2Prongs, HfCandProng3WithMCRec const& cand3Prongs, aod::BigTracksMC const& tracks, aod::McParticles const& particlesMC, aod::McCollisions const& mcCollisions, CollisionsWithMCLabels const& collisions)
void process(HfCandProng2WithMCRec const& cand2Prongs, HfCandProng3WithMCRec const& cand3Prongs, TracksWithSel const& tracks, aod::McParticles const& particlesMC, aod::McCollisions const& mcCollisions, CollisionsWithMCLabels const& collisions)
{
// loop over tracks
for (auto& track : tracks) {
if (track.isGlobalTrackWoDCA() != (uint8_t) true) {
continue;
}
uint index = uint(track.collisionId() >= 0);
if (track.has_mcParticle()) {
auto particle = track.mcParticle(); // get corresponding MC particle to check origin
auto origin = RecoDecay::getCharmHadronOrigin(particlesMC, particle, true);
histOriginTracks[index]->Fill(origin, track.pt());
float deltaZ = -999.f;
if (index) {
auto collision = track.collision_as<CollisionsWithMCLabels>();
auto mcCollision = particle.mcCollision();
deltaZ = collision.posZ() - mcCollision.posZ();
uint index2 = 2;
if ((collision.mcCollisionId() == particle.mcCollisionId()) && std::abs(collision.posZ() - mcCollision.posZ()) < 0.02) { // 200 microns compatibility of Z vertex position also required
if (collision.mcCollisionId() == particle.mcCollisionId()) {
index2 = 1;
}
histOriginTracks[index + index2]->Fill(origin, track.pt());
histOriginTracks[index + index2]->Fill(origin, track.pt(), track.eta(), deltaZ);
}
histOriginTracks[index]->Fill(origin, track.pt(), track.eta(), deltaZ);
} else {
histOriginTracks[index]->Fill(-1., track.pt());
histOriginTracks[index]->Fill(-1.f, track.pt(), track.eta(), -999.f);
}
}

Expand Down Expand Up @@ -359,8 +365,8 @@ struct ValidationRecLevel {

if (whichHad >= 0 && whichOrigin >= 0) {
int indexParticle = 0;
if (cand2Prong.index0_as<aod::BigTracksMC>().has_mcParticle()) {
indexParticle = RecoDecay::getMother(particlesMC, cand2Prong.index0_as<aod::BigTracksMC>().mcParticle(), PDGArrayParticle[whichHad], true);
if (cand2Prong.index0_as<TracksWithSel>().has_mcParticle()) {
indexParticle = RecoDecay::getMother(particlesMC, cand2Prong.index0_as<TracksWithSel>().mcParticle(), PDGArrayParticle[whichHad], true);
}
auto mother = particlesMC.rawIteratorAt(indexParticle);
histDeltaPt[whichHad]->Fill(cand2Prong.pt() - mother.pt());
Expand Down Expand Up @@ -424,8 +430,8 @@ struct ValidationRecLevel {

if (whichHad >= 0) {
int indexParticle = 0;
if (cand3Prong.index0_as<aod::BigTracksMC>().has_mcParticle()) {
indexParticle = RecoDecay::getMother(particlesMC, cand3Prong.index0_as<aod::BigTracksMC>().mcParticle(), PDGArrayParticle[whichHad], true);
if (cand3Prong.index0_as<TracksWithSel>().has_mcParticle()) {
indexParticle = RecoDecay::getMother(particlesMC, cand3Prong.index0_as<TracksWithSel>().mcParticle(), PDGArrayParticle[whichHad], true);
}
auto mother = particlesMC.rawIteratorAt(indexParticle);
histDeltaPt[whichHad]->Fill(cand3Prong.pt() - mother.pt());
Expand Down