From 94c119412d68ab1776e0c2e1ac6c2ffe08e5d2e9 Mon Sep 17 00:00:00 2001 From: Francesco Mazzaschi Date: Fri, 5 Sep 2025 10:21:37 +0200 Subject: [PATCH] Add He3 fake cluster information --- PWGLF/DataModel/LFHypernucleiTables.h | 4 +++- PWGLF/TableProducer/Nuspex/hyperRecoTask.cxx | 11 +++++++---- 2 files changed, 10 insertions(+), 5 deletions(-) diff --git a/PWGLF/DataModel/LFHypernucleiTables.h b/PWGLF/DataModel/LFHypernucleiTables.h index f0f37bf4d86..58adbe10ed6 100644 --- a/PWGLF/DataModel/LFHypernucleiTables.h +++ b/PWGLF/DataModel/LFHypernucleiTables.h @@ -14,8 +14,8 @@ /// \brief Slim hypernuclei tables /// -#include "Framework/AnalysisDataModel.h" #include "Framework/ASoAHelpers.h" +#include "Framework/AnalysisDataModel.h" #ifndef PWGLF_DATAMODEL_LFHYPERNUCLEITABLES_H_ #define PWGLF_DATAMODEL_LFHYPERNUCLEITABLES_H_ @@ -80,6 +80,7 @@ DECLARE_SOA_COLUMN(GenXDecVtx, genXDecVtx, float); // Decay v DECLARE_SOA_COLUMN(GenYDecVtx, genYDecVtx, float); // Decay vertex of the candidate (y direction) DECLARE_SOA_COLUMN(GenZDecVtx, genZDecVtx, float); // Decay vertex of the candidate (z direction) DECLARE_SOA_COLUMN(IsReco, isReco, bool); // bool: true for reco +DECLARE_SOA_COLUMN(IsFakeHeOnITSLayer, isFakeHeOnITSLayer, uint8_t); // uint8_t: bit map for fake He on ITS layers DECLARE_SOA_COLUMN(IsSignal, isSignal, bool); // bool: true for signal DECLARE_SOA_COLUMN(IsRecoMCCollision, isRecoMCCollision, bool); // bool: true for reco MC collision DECLARE_SOA_COLUMN(IsSurvEvSel, isSurvEvSel, bool); // bool: true for survived event selection @@ -143,6 +144,7 @@ DECLARE_SOA_TABLE(MCHypCands, "AOD", "MCHYPCANDS", hyperrec::GenYDecVtx, hyperrec::GenZDecVtx, hyperrec::IsReco, + hyperrec::IsFakeHeOnITSLayer, hyperrec::IsSignal, hyperrec::IsRecoMCCollision, hyperrec::IsSurvEvSel); diff --git a/PWGLF/TableProducer/Nuspex/hyperRecoTask.cxx b/PWGLF/TableProducer/Nuspex/hyperRecoTask.cxx index e93828980b1..06f5611872a 100644 --- a/PWGLF/TableProducer/Nuspex/hyperRecoTask.cxx +++ b/PWGLF/TableProducer/Nuspex/hyperRecoTask.cxx @@ -124,8 +124,10 @@ struct hyperCandidate { unsigned int collisionID = 0; bool isMatter = false; - bool isSignal = false; // true MC signal - bool isReco = false; // true if the candidate is actually reconstructed + bool isSignal = false; // true MC signal + bool isReco = false; // true if the candidate is actually reconstructed + uint8_t isFakeHeOnITSLayer = 0u; // bit map for fake He on ITS layers + bool isRecoMCCollision = false; // true if the corresponding MC collision has been reconstructed bool isSurvEvSelection = false; // true if the corresponding event passed the event selection int pdgCode = 0; // PDG code of the hypernucleus @@ -648,6 +650,7 @@ struct hyperRecoTask { hypCand.gDecVtx[i] = secVtx[i] - primVtx[i]; } hypCand.isSignal = true; + hypCand.isFakeHeOnITSLayer = mcLabHe.mcMask() & 0x7F; // check if any of the first 7 bits is set hypCand.pdgCode = heMother.pdgCode(); hypCand.isRecoMCCollision = recoCollisionIds[heMother.mcCollisionId()] > 0; hypCand.isSurvEvSelection = isSurvEvSelCollision[heMother.mcCollisionId()]; @@ -766,7 +769,7 @@ struct hyperRecoTask { hypCand.clusterSizeITSHe3, hypCand.clusterSizeITSPi, hypCand.flags, trackedHypClSize, chargeFactor * hypCand.genPt(), hypCand.genPhi(), hypCand.genEta(), hypCand.genPtHe3(), hypCand.gDecVtx[0], hypCand.gDecVtx[1], hypCand.gDecVtx[2], - hypCand.isReco, hypCand.isSignal, hypCand.isRecoMCCollision, hypCand.isSurvEvSelection); + hypCand.isReco, hypCand.isFakeHeOnITSLayer, hypCand.isSignal, hypCand.isRecoMCCollision, hypCand.isSurvEvSelection); } // now we fill only the signal candidates that were not reconstructed @@ -839,7 +842,7 @@ struct hyperRecoTask { -1, -1, -1, false, chargeFactor * hypCand.genPt(), hypCand.genPhi(), hypCand.genEta(), hypCand.genPtHe3(), hypCand.gDecVtx[0], hypCand.gDecVtx[1], hypCand.gDecVtx[2], - hypCand.isReco, hypCand.isSignal, hypCand.isRecoMCCollision, hypCand.isSurvEvSelection); + hypCand.isReco, -1, hypCand.isSignal, hypCand.isRecoMCCollision, hypCand.isSurvEvSelection); } } PROCESS_SWITCH(hyperRecoTask, processMC, "MC analysis", false);