From 3ec57249069c67df92d60edcb94e03d64802da30 Mon Sep 17 00:00:00 2001 From: Mario Ciacco Date: Thu, 30 Nov 2023 15:59:01 +0100 Subject: [PATCH] add duplicate vertex studies: same ncontrib + t close to ITS roof boundary (#4030) * add duplicate vertices studies: same ncontrib + t close to ITS Rof boundary * remove unused configurable --- PWGLF/Tasks/QC/vertexQA.cxx | 12 ++++++++++++ 1 file changed, 12 insertions(+) diff --git a/PWGLF/Tasks/QC/vertexQA.cxx b/PWGLF/Tasks/QC/vertexQA.cxx index d4c19a57088..3c35b5b0ca6 100644 --- a/PWGLF/Tasks/QC/vertexQA.cxx +++ b/PWGLF/Tasks/QC/vertexQA.cxx @@ -66,6 +66,9 @@ struct vertexQA { Configurable nSigmaR{"nSigmaR", 1000., "Number of sigmas for transverse displacement of vertices"}; Configurable nSigmaT{"nSigmaT", 1000., "Number of sigmas for time of vertices"}; Configurable maxTime{"maxTime", 100000., "Maximum time difference between split vertices in ns"}; + Configurable minTimeDiff{"minTimeDiff", 14000., "Minimum time difference (ITS Rof peak)"}; + Configurable maxTimeDiff{"maxTimeDiff", 15500., "Maximum time difference (ITS Rof peak)"}; + Configurable minNcontrib{"minNcontrib", 500., "Minimum number of contributors (duplicate vertices)"}; ConfigurableAxis xVtxAxis{"xVtxBins", {100, -0.1f, 0.1f}, "Binning for the vertex x (y) in cm"}; ConfigurableAxis zVtxAxis{"zVtxBins", {100, -20.f, 20.f}, "Binning for the vertex z in cm"}; @@ -113,6 +116,9 @@ struct vertexQA { histos.add("zDistVsTDistVtxTimeSeriesHistogram", ";#Delta#it{t}_{vtx} (ns);#Delta#it{z}_{vtx} (cm)", HistType::kTH2F, {tDiffVtxAxisExtend, zDiffVtxAxis2}); histos.add("nContribTwoVtxTimeSeriesHistogram", ";#it{N}_{contrib}^{1};#it{N}_{contrib}^{2}", HistType::kTH2F, {nContribAxis, nContribAxis}); histos.add("nContribVsTDistTimeSeriesHistogram", ";#Delta#it{t}_{vtx} (ns);#Delta#it{N}_{contrib}", HistType::kTH2F, {tDiffVtxAxisExtendSigned, nContribDiffAxis}); + + histos.add("nContribITSRofTimeSeriesHistogram", ";#it{N}_{contrib}^{1};#it{N}_{contrib}^{2}", HistType::kTH2F, {nContribAxis, nContribAxis}); + histos.add("tDiffDuplicateTimeSeriesHistogram", ";#Delta#it{t}_{vtx} (ns);Entries", HistType::kTH1F, {tDiffVtxAxisExtend}); } std::deque colls; @@ -206,6 +212,12 @@ struct vertexQA { histos.fill(HIST("nContribTwoVtxTimeSeriesHistogram"), coll1.numContrib(), coll2.numContrib()); histos.fill(HIST("nContribVsTDistTimeSeriesHistogram"), -deltaT, coll2.numContrib() - coll1.numContrib()); histos.fill(HIST("nVtxTimeSeriesHistogram"), 2); + if (std::abs(deltaT) > minTimeDiff && std::abs(deltaT) < maxTimeDiff) { + histos.fill(HIST("nContribITSRofTimeSeriesHistogram"), coll1.numContrib(), coll2.numContrib()); + } + if (std::abs(coll1.numContrib() - coll2.numContrib()) < 3 * std::sqrt(coll1.numContrib() + coll2.numContrib()) && coll1.numContrib() > minNcontrib && coll2.numContrib() > minNcontrib) { + histos.fill(HIST("tDiffDuplicateTimeSeriesHistogram"), std::abs(deltaT)); + } colls.erase(id); colls.pop_front(); }