From dd6c1857d5d72003aeb2b04ec95143e585a43449 Mon Sep 17 00:00:00 2001 From: Andrea Date: Mon, 11 Mar 2024 09:54:03 +0100 Subject: [PATCH 1/2] Adding histogram for studying single-track pt resolution. --- DPG/Tasks/AOTTrack/MonitorFilterBit.cxx | 4 ++++ 1 file changed, 4 insertions(+) diff --git a/DPG/Tasks/AOTTrack/MonitorFilterBit.cxx b/DPG/Tasks/AOTTrack/MonitorFilterBit.cxx index ec08e3cc70b..78d8e269f33 100644 --- a/DPG/Tasks/AOTTrack/MonitorFilterBit.cxx +++ b/DPG/Tasks/AOTTrack/MonitorFilterBit.cxx @@ -151,6 +151,9 @@ struct CheckFilterBit { histos.add("Tracks/RecoMCPhysPrimCollMatch/histptMCTPConlyWithClusters", "TPConlyWithClusters;#it{p}_{T}^{gen} (GeV/#it{c});#it{#eta};#it{#varphi};NclustTPC", HistType::kTHnF, {axisPt, axisEta, axisPhi, axisNclustTPC}); histos.add("Tracks/RecoMCPhysPrimCollMatch/histptTPConlyWithClusters", "TPConlyWithClusters;#it{p}_{T} (GeV/#it{c});#it{#eta};#it{#varphi};NclustTPC", HistType::kTHnF, {axisPt, axisEta, axisPhi, axisNclustTPC}); + histos.add("Tracks/RecoMCPhysPrimCollMatch/histptMCFB0resolCurv", "FB0;#it{p}_{T}^{gen} (GeV/#it{c});#it{#eta};#it{p}_{T}^{gen}/#it{p}_{T}^{reco}-1", kTH3D, {{200,0,20},{10,-1,1},{200,-0.2,0.2}}); + + histos.add("Tracks/RecoMCRad1to15cmCollMatch/histptFB0", "FB0;#it{p}_{T} (GeV/#it{c});#it{#eta};#it{#varphi}", kTH3D, {axisPt, axisEta, axisPhi}); histos.add("Tracks/RecoMCRad1to15cmCollMatch/histptFB1", "FB1;#it{p}_{T} (GeV/#it{c});#it{#eta};#it{#varphi}", kTH3D, {axisPt, axisEta, axisPhi}); histos.add("Tracks/RecoMCRad1to15cmCollMatch/histptFB2", "FB2;#it{p}_{T} (GeV/#it{c});#it{#eta};#it{#varphi}", kTH3D, {axisPt, axisEta, axisPhi}); @@ -440,6 +443,7 @@ struct CheckFilterBit { if (track.isGlobalTrack()) { histos.fill(HIST("Tracks/RecoMCPhysPrimCollMatch/histptFB0"), track.pt(), track.eta(), track.phi()); histos.fill(HIST("Tracks/RecoMCPhysPrimCollMatch/histptMCFB0"), mcparticle.pt(), track.eta(), track.phi()); + histos.fill(HIST("Tracks/RecoMCPhysPrimCollMatch/histptMCFB0resolCurv"),mcparticle.pt(), track.eta(),(mcparticle.pt()/track.pt()-1.));// did not use signed1Pt just to avoid getting the charge for MCpart; last variable is relative difference on curvature (1/pt_reco - 1/pt_gen)/(1/pt_gen) } if (track.itsChi2NCl() > 0. && track.tpcChi2NCl() < 0.) { histos.fill(HIST("Tracks/RecoMCPhysPrimCollMatch/histptITSonly"), track.pt(), track.eta(), track.phi()); From 96656b0493648d39c5c5c3b1f77369d625bc2c37 Mon Sep 17 00:00:00 2001 From: Andrea Date: Mon, 11 Mar 2024 09:58:42 +0100 Subject: [PATCH 2/2] fixing clang formatting --- DPG/Tasks/AOTTrack/MonitorFilterBit.cxx | 5 ++--- 1 file changed, 2 insertions(+), 3 deletions(-) diff --git a/DPG/Tasks/AOTTrack/MonitorFilterBit.cxx b/DPG/Tasks/AOTTrack/MonitorFilterBit.cxx index 78d8e269f33..104dfabbaf8 100644 --- a/DPG/Tasks/AOTTrack/MonitorFilterBit.cxx +++ b/DPG/Tasks/AOTTrack/MonitorFilterBit.cxx @@ -151,9 +151,8 @@ struct CheckFilterBit { histos.add("Tracks/RecoMCPhysPrimCollMatch/histptMCTPConlyWithClusters", "TPConlyWithClusters;#it{p}_{T}^{gen} (GeV/#it{c});#it{#eta};#it{#varphi};NclustTPC", HistType::kTHnF, {axisPt, axisEta, axisPhi, axisNclustTPC}); histos.add("Tracks/RecoMCPhysPrimCollMatch/histptTPConlyWithClusters", "TPConlyWithClusters;#it{p}_{T} (GeV/#it{c});#it{#eta};#it{#varphi};NclustTPC", HistType::kTHnF, {axisPt, axisEta, axisPhi, axisNclustTPC}); - histos.add("Tracks/RecoMCPhysPrimCollMatch/histptMCFB0resolCurv", "FB0;#it{p}_{T}^{gen} (GeV/#it{c});#it{#eta};#it{p}_{T}^{gen}/#it{p}_{T}^{reco}-1", kTH3D, {{200,0,20},{10,-1,1},{200,-0.2,0.2}}); + histos.add("Tracks/RecoMCPhysPrimCollMatch/histptMCFB0resolCurv", "FB0;#it{p}_{T}^{gen} (GeV/#it{c});#it{#eta};#it{p}_{T}^{gen}/#it{p}_{T}^{reco}-1", kTH3D, {{200, 0, 20}, {10, -1, 1}, {200, -0.2, 0.2}}); - histos.add("Tracks/RecoMCRad1to15cmCollMatch/histptFB0", "FB0;#it{p}_{T} (GeV/#it{c});#it{#eta};#it{#varphi}", kTH3D, {axisPt, axisEta, axisPhi}); histos.add("Tracks/RecoMCRad1to15cmCollMatch/histptFB1", "FB1;#it{p}_{T} (GeV/#it{c});#it{#eta};#it{#varphi}", kTH3D, {axisPt, axisEta, axisPhi}); histos.add("Tracks/RecoMCRad1to15cmCollMatch/histptFB2", "FB2;#it{p}_{T} (GeV/#it{c});#it{#eta};#it{#varphi}", kTH3D, {axisPt, axisEta, axisPhi}); @@ -443,7 +442,7 @@ struct CheckFilterBit { if (track.isGlobalTrack()) { histos.fill(HIST("Tracks/RecoMCPhysPrimCollMatch/histptFB0"), track.pt(), track.eta(), track.phi()); histos.fill(HIST("Tracks/RecoMCPhysPrimCollMatch/histptMCFB0"), mcparticle.pt(), track.eta(), track.phi()); - histos.fill(HIST("Tracks/RecoMCPhysPrimCollMatch/histptMCFB0resolCurv"),mcparticle.pt(), track.eta(),(mcparticle.pt()/track.pt()-1.));// did not use signed1Pt just to avoid getting the charge for MCpart; last variable is relative difference on curvature (1/pt_reco - 1/pt_gen)/(1/pt_gen) + histos.fill(HIST("Tracks/RecoMCPhysPrimCollMatch/histptMCFB0resolCurv"), mcparticle.pt(), track.eta(), (mcparticle.pt() / track.pt() - 1.)); // did not use signed1Pt just to avoid getting the charge for MCpart; last variable is relative difference on curvature (1/pt_reco - 1/pt_gen)/(1/pt_gen) } if (track.itsChi2NCl() > 0. && track.tpcChi2NCl() < 0.) { histos.fill(HIST("Tracks/RecoMCPhysPrimCollMatch/histptITSonly"), track.pt(), track.eta(), track.phi());