From 75afae2719b8ea82c4dcbb75f7a4e0b547946cc3 Mon Sep 17 00:00:00 2001 From: Mattia Faggin Date: Wed, 13 Jul 2022 10:26:48 +0200 Subject: [PATCH] DPG: add further correlation plots (lite). --- DPG/Tasks/qaEventTrackLite.cxx | 10 ++++++++++ 1 file changed, 10 insertions(+) diff --git a/DPG/Tasks/qaEventTrackLite.cxx b/DPG/Tasks/qaEventTrackLite.cxx index 30494831a62..cb80ce26235 100644 --- a/DPG/Tasks/qaEventTrackLite.cxx +++ b/DPG/Tasks/qaEventTrackLite.cxx @@ -50,6 +50,7 @@ struct qaEventTrackLite { // Event selections Configurable vtxZMax{"vtxZMax", 10.f, "Max VTX Z position"}; // Track selections + Configurable bHasITS{"bHasITS", false, "Select only DPG tracks with ITS (no forcing on TPC information)"}; Configurable bItsStandalone{"bItsStandalone", false, "Select only ITS standalone DPG tracks"}; Configurable bTpcOnly{"bTpcOnly", false, "Select only TPC only DPG tracks"}; Configurable bItsTpcMatched{"bItsTpcMatched", false, "Select ITS-TPC matched DPG tracks"}; @@ -117,9 +118,12 @@ struct qaEventTrackLite { // kine histograms histos.add("Tracks/VertexPositionZ", "", kTH1D, {{100, -20.f, 20.f, "Vertex Z (cm)"}}); + histos.add("Tracks/VertexPositionZvsEta", "", kTH2D, {axisEta, {100, -20.f, 20.f, "Vertex Z (cm)"}}); + histos.add("Tracks/VertexPositionZvsEtaHasITS", "", kTH2D, {axisEta, {100, -20.f, 20.f, "Vertex Z (cm)"}}); histos.add("Tracks/Kine/pt", "#it{p}_{T}", kTH1D, {axisPt}); histos.add("Tracks/Kine/eta", "#eta", kTH1D, {axisEta}); histos.add("Tracks/Kine/phi", "#phi", kTH1D, {axisPhi}); + histos.add("Tracks/Kine/phiVsEta", "#phi", kTH2D, {axisPhi, axisEta}); histos.add("Tracks/length", "track length in cm;#it{Length} (cm);", kTH1D, {{400, 0, 1000}}); const AxisSpec axisImpParRPhi{binsImpPar, "#it{d}_{r#it{#varphi}} (#cm)"}; const AxisSpec axisImpParZAxis{binsImpPar, "#it{d}_{z} (#cm)"}; @@ -192,6 +196,7 @@ struct qaEventTrackLite { Filter ptCut = o2::aod::dpgtrack::pt > ptMin; Filter etaCut = etaMin < o2::aod::dpgtrack::eta && o2::aod::dpgtrack::eta < etaMax; // Detector matching + Filter filterHasITS = (bHasITS.node() == false) || (bItsStandalone.node() == false && bTpcOnly.node() == false && bItsTpcMatched.node() == false && o2::aod::dpgtrack::hasITS == true); Filter itsStandaloneTracks = (bItsStandalone.node() == false) || (o2::aod::dpgtrack::hasITS == true && o2::aod::dpgtrack::hasTPC == false); Filter tpcOnlyTracks = (bTpcOnly.node() == false) || (o2::aod::dpgtrack::hasITS == false && o2::aod::dpgtrack::hasTPC == true); Filter itsTpcMatchedTracks = (bItsTpcMatched.node() == false) || (o2::aod::dpgtrack::hasITS == true && o2::aod::dpgtrack::hasTPC == true); @@ -247,9 +252,14 @@ struct qaEventTrackLite { histos.fill(HIST("Particles/Kine/phi"), track.phiMC()); } histos.fill(HIST("Tracks/VertexPositionZ"), track.dpgCollision().posZ()); + histos.fill(HIST("Tracks/VertexPositionZvsEta"), track.eta(), track.dpgCollision().posZ()); + if (track.hasITS()) { + histos.fill(HIST("Tracks/VertexPositionZvsEtaHasITS"), track.eta(), track.dpgCollision().posZ()); + } histos.fill(HIST("Tracks/Kine/pt"), track.pt()); histos.fill(HIST("Tracks/Kine/eta"), track.eta()); histos.fill(HIST("Tracks/Kine/phi"), track.phi()); + histos.fill(HIST("Tracks/Kine/phiVsEta"), track.phi(), track.eta()); histos.fill(HIST("Tracks/dcaXY"), track.dcaXY()); histos.fill(HIST("Tracks/dcaZ"), track.dcaZ()); histos.fill(HIST("Tracks/dcaXYvsPt"), track.pt(), track.dcaXY());