From a4be99ba2f64bf11a9ae04621091a66cea99c421 Mon Sep 17 00:00:00 2001 From: altsybee Date: Wed, 29 Oct 2025 12:38:51 +0100 Subject: [PATCH] add more QA histos for MC --- DPG/Tasks/AOTEvent/lightIonsEvSelQa.cxx | 26 ++++++++++++++++++++++--- 1 file changed, 23 insertions(+), 3 deletions(-) diff --git a/DPG/Tasks/AOTEvent/lightIonsEvSelQa.cxx b/DPG/Tasks/AOTEvent/lightIonsEvSelQa.cxx index 03ef02d3722..e87ca9c73f3 100644 --- a/DPG/Tasks/AOTEvent/lightIonsEvSelQa.cxx +++ b/DPG/Tasks/AOTEvent/lightIonsEvSelQa.cxx @@ -655,6 +655,12 @@ struct LightIonsEvSelQa { histos.add("MCnonTVX/hMCdataVzDiff", "", kTH2F, {axisNcontrib, axisVtxZdiff}); histos.add("MCnonTVX/hMCdataBcDiffVsMult", "", kTH2F, {axisNcontrib, axisBcDiff}); histos.add("MCnonTVX/hMCdataFoundBcDiffVsMult", "", kTH2F, {axisNcontrib, axisBcDiff}); + + // + histos.add("MC_not_TF_ROF_borders/hNcontribColFromData", "", kTH1F, {axisNcontrib}); + histos.add("MC_not_TF_ROF_borders/hNcontribAccFromData", "", kTH1F, {axisNcontrib}); + histos.add("MC_not_TF_ROF_borders/hNcontribColFromData_foundBcDiff0", "", kTH1F, {axisNcontrib}); + histos.add("MC_not_TF_ROF_borders/hNcontribAccFromData_foundBcDiff0", "", kTH1F, {axisNcontrib}); } Preslice perCollision = aod::track::collisionId; @@ -946,9 +952,11 @@ struct LightIonsEvSelQa { nGlobalTracksPV++; } // end of track loop - histos.fill(HIST("hNcontribColFromData"), nPVtracks); - if (col.selection_bit(kIsTriggerTVX)) - histos.fill(HIST("hNcontribAccFromData"), nPVtracks); + if (col.selection_bit(kNoTimeFrameBorder) && col.selection_bit(kNoITSROFrameBorder)) { + histos.fill(HIST("hNcontribColFromData"), nPVtracks); + if (col.selection_bit(kIsTriggerTVX)) + histos.fill(HIST("hNcontribAccFromData"), nPVtracks); + } bool hasFT0 = foundBC.has_ft0(); bool hasFV0A = foundBC.has_fv0a(); @@ -1867,6 +1875,18 @@ struct LightIonsEvSelQa { histos.fill(HIST("MCnonTVX/hMCdataBcDiffVsMult"), col.numContrib(), bcDiff); histos.fill(HIST("MCnonTVX/hMCdataFoundBcDiffVsMult"), col.numContrib(), foundBcDiff); } + + if (col.selection_bit(kNoTimeFrameBorder) && col.selection_bit(kNoITSROFrameBorder)) { + histos.fill(HIST("MC_not_TF_ROF_borders/hNcontribColFromData"), col.numContrib()); + if (col.selection_bit(kIsTriggerTVX)) + histos.fill(HIST("MC_not_TF_ROF_borders/hNcontribAccFromData"), col.numContrib()); + + if (foundBcDiff == 0) { + histos.fill(HIST("MC_not_TF_ROF_borders/hNcontribColFromData_foundBcDiff0"), col.numContrib()); + if (col.selection_bit(kIsTriggerTVX)) + histos.fill(HIST("MC_not_TF_ROF_borders/hNcontribAccFromData_foundBcDiff0"), col.numContrib()); + } + } } } }