From e175a6c61f665eb8ea1df6a58b59df508797ef1f Mon Sep 17 00:00:00 2001 From: spolitan Date: Thu, 11 Jun 2026 15:20:34 +0200 Subject: [PATCH 1/3] Adding th2 in tracking QA --- DPG/Tasks/AOTTrack/qaEfficiency.cxx | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/DPG/Tasks/AOTTrack/qaEfficiency.cxx b/DPG/Tasks/AOTTrack/qaEfficiency.cxx index 20abc845190..d00bce40567 100644 --- a/DPG/Tasks/AOTTrack/qaEfficiency.cxx +++ b/DPG/Tasks/AOTTrack/qaEfficiency.cxx @@ -110,6 +110,7 @@ std::array, nParticles> hPtGeneratedRecoEv; std::array, nParticles> hPtItsPrm; std::array, nParticles> hPtItsTpcPrm; std::array, nParticles> hPtTrkItsTpcPrm; +std::array, nParticles> hPtGenVsPtTrkItsTpcPrm; // histogram to store generated pt vs reco pt for primaries std::array, nParticles> hPtItsTpcTofPrm; std::array, nParticles> hPtTrkItsTpcTofPrm; std::array, nParticles> hPtGeneratedPrm; @@ -376,11 +377,11 @@ struct QaEfficiency { hPtItsPrm[histogramIndex] = histos.add(Form("MC/pdg%i/pt/prm/its", PDGs[histogramIndex]), "ITS tracks (primaries) " + tagPt, kTH1D, {axisPt}); hPtItsTpcPrm[histogramIndex] = histos.add(Form("MC/pdg%i/pt/prm/its_tpc", PDGs[histogramIndex]), "ITS-TPC tracks (primaries) " + tagPt, kTH1D, {axisPt}); hPtTrkItsTpcPrm[histogramIndex] = histos.add(Form("MC/pdg%i/pt/prm/trk/its_tpc", PDGs[histogramIndex]), "ITS-TPC tracks (reco primaries) " + tagPt, kTH1D, {axisPt}); + hPtGenVsPtTrkItsTpcPrm[histogramIndex] = histos.add(Form("MC/pdg%i/pt/prm/generated_vs_reco", PDGs[histogramIndex]), "Abs(Gen - Reco) pT vs Gen pT (primaries) " + tagPt, kTH2D, {axisPt, axisPt}); hPtItsTpcTofPrm[histogramIndex] = histos.add(Form("MC/pdg%i/pt/prm/its_tpc_tof", PDGs[histogramIndex]), "ITS-TPC-TOF tracks (primaries) " + tagPt, kTH1D, {axisPt}); hPtTrkItsTpcTofPrm[histogramIndex] = histos.add(Form("MC/pdg%i/pt/prm/trk/its_tpc_tof", PDGs[histogramIndex]), "ITS-TPC-TOF tracks (reco primaries) " + tagPt, kTH1D, {axisPt}); hPtGeneratedPrm[histogramIndex] = histos.add(Form("MC/pdg%i/pt/prm/generated", PDGs[histogramIndex]), "Generated (primaries) " + tagPt, kTH1D, {axisPt}); hPtGeneratedPrmRecoEv[histogramIndex] = histos.add(Form("MC/pdg%i/pt/prm/generated_reco_ev", PDGs[histogramIndex]), "Generated Reco Ev. " + tagPt, kTH1D, {axisPt}); - // Str hPtItsTpcStr[histogramIndex] = histos.add(Form("MC/pdg%i/pt/str/its_tpc", PDGs[histogramIndex]), "ITS-TPC tracks (from weak decays) " + tagPt, kTH1D, {axisPt}); hPtTrkItsTpcStr[histogramIndex] = histos.add(Form("MC/pdg%i/pt/str/trk/its_tpc", PDGs[histogramIndex]), "ITS-TPC tracks (reco from weak decays) " + tagPt, kTH1D, {axisPt}); @@ -1156,6 +1157,7 @@ struct QaEfficiency { if (passedITS && passedTPC) { hPtItsTpcPrm[histogramIndex]->Fill(mcParticle.pt()); hPtTrkItsTpcPrm[histogramIndex]->Fill(track.pt()); + hPtGenVsPtTrkItsTpcPrm[histogramIndex]->Fill(mcParticle.pt(), abs(track.pt() - mcParticle.pt())); hEtaItsTpcPrm[histogramIndex]->Fill(mcParticle.eta()); hEtaTrkItsTpcPrm[histogramIndex]->Fill(track.eta()); hPhiItsTpcPrm[histogramIndex]->Fill(mcParticle.phi()); From b2a15294192430c721496059f95e64545d3d2e43 Mon Sep 17 00:00:00 2001 From: spolitan Date: Fri, 12 Jun 2026 09:09:26 +0200 Subject: [PATCH 2/3] Adding th2 for pt gen vs pt reck in track QA --- DPG/Tasks/AOTTrack/qaEfficiency.cxx | 10 +++++++--- 1 file changed, 7 insertions(+), 3 deletions(-) diff --git a/DPG/Tasks/AOTTrack/qaEfficiency.cxx b/DPG/Tasks/AOTTrack/qaEfficiency.cxx index d00bce40567..7a08f60c0a9 100644 --- a/DPG/Tasks/AOTTrack/qaEfficiency.cxx +++ b/DPG/Tasks/AOTTrack/qaEfficiency.cxx @@ -110,7 +110,8 @@ std::array, nParticles> hPtGeneratedRecoEv; std::array, nParticles> hPtItsPrm; std::array, nParticles> hPtItsTpcPrm; std::array, nParticles> hPtTrkItsTpcPrm; -std::array, nParticles> hPtGenVsPtTrkItsTpcPrm; // histogram to store generated pt vs reco pt for primaries +std::array, nParticles> hDeltaPtVsPtTrkItsTpcPrm; +std::array, nParticles> hPtGenVsPtTrkItsTpcPrm; std::array, nParticles> hPtItsTpcTofPrm; std::array, nParticles> hPtTrkItsTpcTofPrm; std::array, nParticles> hPtGeneratedPrm; @@ -377,11 +378,13 @@ struct QaEfficiency { hPtItsPrm[histogramIndex] = histos.add(Form("MC/pdg%i/pt/prm/its", PDGs[histogramIndex]), "ITS tracks (primaries) " + tagPt, kTH1D, {axisPt}); hPtItsTpcPrm[histogramIndex] = histos.add(Form("MC/pdg%i/pt/prm/its_tpc", PDGs[histogramIndex]), "ITS-TPC tracks (primaries) " + tagPt, kTH1D, {axisPt}); hPtTrkItsTpcPrm[histogramIndex] = histos.add(Form("MC/pdg%i/pt/prm/trk/its_tpc", PDGs[histogramIndex]), "ITS-TPC tracks (reco primaries) " + tagPt, kTH1D, {axisPt}); - hPtGenVsPtTrkItsTpcPrm[histogramIndex] = histos.add(Form("MC/pdg%i/pt/prm/generated_vs_reco", PDGs[histogramIndex]), "Abs(Gen - Reco) pT vs Gen pT (primaries) " + tagPt, kTH2D, {axisPt, axisPt}); + hDeltaPtVsPtTrkItsTpcPrm[histogramIndex] = histos.add(Form("MC/pdg%i/pt/prm/generated_vs_reco_delta", PDGs[histogramIndex]), "Abs(Gen - Reco) pT vs Gen pT (primaries) " + tagPt, kTH2D, {axisPt, axisPt}); + hPtGenVsPtTrkItsTpcPrm[histogramIndex] = histos.add(Form("MC/pdg%i/pt/prm/generated_vs_reco", PDGs[histogramIndex]), "Reco pT vs Gen pT (primaries) " + tagPt, kTH2D, {axisPt, axisPt}); hPtItsTpcTofPrm[histogramIndex] = histos.add(Form("MC/pdg%i/pt/prm/its_tpc_tof", PDGs[histogramIndex]), "ITS-TPC-TOF tracks (primaries) " + tagPt, kTH1D, {axisPt}); hPtTrkItsTpcTofPrm[histogramIndex] = histos.add(Form("MC/pdg%i/pt/prm/trk/its_tpc_tof", PDGs[histogramIndex]), "ITS-TPC-TOF tracks (reco primaries) " + tagPt, kTH1D, {axisPt}); hPtGeneratedPrm[histogramIndex] = histos.add(Form("MC/pdg%i/pt/prm/generated", PDGs[histogramIndex]), "Generated (primaries) " + tagPt, kTH1D, {axisPt}); hPtGeneratedPrmRecoEv[histogramIndex] = histos.add(Form("MC/pdg%i/pt/prm/generated_reco_ev", PDGs[histogramIndex]), "Generated Reco Ev. " + tagPt, kTH1D, {axisPt}); + // Str hPtItsTpcStr[histogramIndex] = histos.add(Form("MC/pdg%i/pt/str/its_tpc", PDGs[histogramIndex]), "ITS-TPC tracks (from weak decays) " + tagPt, kTH1D, {axisPt}); hPtTrkItsTpcStr[histogramIndex] = histos.add(Form("MC/pdg%i/pt/str/trk/its_tpc", PDGs[histogramIndex]), "ITS-TPC tracks (reco from weak decays) " + tagPt, kTH1D, {axisPt}); @@ -1157,7 +1160,8 @@ struct QaEfficiency { if (passedITS && passedTPC) { hPtItsTpcPrm[histogramIndex]->Fill(mcParticle.pt()); hPtTrkItsTpcPrm[histogramIndex]->Fill(track.pt()); - hPtGenVsPtTrkItsTpcPrm[histogramIndex]->Fill(mcParticle.pt(), abs(track.pt() - mcParticle.pt())); + hDeltaPtVsPtTrkItsTpcPrm[histogramIndex]->Fill(mcParticle.pt(), abs(track.pt() - mcParticle.pt())); + hPtGenVsPtTrkItsTpcPrm[histogramIndex]->Fill(mcParticle.pt(), track.pt()); hEtaItsTpcPrm[histogramIndex]->Fill(mcParticle.eta()); hEtaTrkItsTpcPrm[histogramIndex]->Fill(track.eta()); hPhiItsTpcPrm[histogramIndex]->Fill(mcParticle.phi()); From eeaee9c8992361525eb3691a5e8c84942aa4e276 Mon Sep 17 00:00:00 2001 From: ALICE Action Bot Date: Fri, 12 Jun 2026 09:46:25 +0000 Subject: [PATCH 3/3] Please consider the following formatting changes --- DPG/Tasks/AOTTrack/qaEfficiency.cxx | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/DPG/Tasks/AOTTrack/qaEfficiency.cxx b/DPG/Tasks/AOTTrack/qaEfficiency.cxx index 7a08f60c0a9..17aa6734cbc 100644 --- a/DPG/Tasks/AOTTrack/qaEfficiency.cxx +++ b/DPG/Tasks/AOTTrack/qaEfficiency.cxx @@ -384,7 +384,7 @@ struct QaEfficiency { hPtTrkItsTpcTofPrm[histogramIndex] = histos.add(Form("MC/pdg%i/pt/prm/trk/its_tpc_tof", PDGs[histogramIndex]), "ITS-TPC-TOF tracks (reco primaries) " + tagPt, kTH1D, {axisPt}); hPtGeneratedPrm[histogramIndex] = histos.add(Form("MC/pdg%i/pt/prm/generated", PDGs[histogramIndex]), "Generated (primaries) " + tagPt, kTH1D, {axisPt}); hPtGeneratedPrmRecoEv[histogramIndex] = histos.add(Form("MC/pdg%i/pt/prm/generated_reco_ev", PDGs[histogramIndex]), "Generated Reco Ev. " + tagPt, kTH1D, {axisPt}); - + // Str hPtItsTpcStr[histogramIndex] = histos.add(Form("MC/pdg%i/pt/str/its_tpc", PDGs[histogramIndex]), "ITS-TPC tracks (from weak decays) " + tagPt, kTH1D, {axisPt}); hPtTrkItsTpcStr[histogramIndex] = histos.add(Form("MC/pdg%i/pt/str/trk/its_tpc", PDGs[histogramIndex]), "ITS-TPC tracks (reco from weak decays) " + tagPt, kTH1D, {axisPt});