Skip to content

[PWGLF/NuSpEx] Saved tpcNClsCrossedRows inside the LFLnnTable + minor changes #15507

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
maciacco merged 1 commit into from
Mar 24, 2026
+26 −15
Merged

[PWGLF/NuSpEx] Saved tpcNClsCrossedRows inside the LFLnnTable + minor changes #15507

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
6 changes: 4 additions & 2 deletions PWGLF/DataModel/LFLnnTables.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -53,6 +53,8 @@ DECLARE_SOA_COLUMN(CosPA, cosPA, double); // Cos
DECLARE_SOA_COLUMN(NSigma3H, nSigma3H, float); // Number of sigmas of the 3H daughter
DECLARE_SOA_COLUMN(NTPCclus3H, nTPCclus3H, uint8_t); // Number of TPC clusters of the 3H daughter
DECLARE_SOA_COLUMN(NTPCclusPi, nTPCclusPi, uint8_t); // Number of TPC clusters of the Pi daughter
DECLARE_SOA_COLUMN(NTPCNClsCrossedRows3H, nTPCclsCrossedRows3H, uint8_t); // Number of TPC clusters crossed rows of the 3H daughter
DECLARE_SOA_COLUMN(NTPCNClsCrossedRowsPi, nTPCclsCrossedRowsPi, uint8_t); // Number of TPC clusters crossed rows of the Pi daughter
DECLARE_SOA_COLUMN(TPCsignal3H, tpcSignal3H, uint16_t); // TPC signal of the 3H daughter
DECLARE_SOA_COLUMN(TPCsignalPi, tpcSignalPi, uint16_t); // TPC signal of the Pi daughter
DECLARE_SOA_COLUMN(Flags, flags, uint8_t); // Flags for PID in tracking (bits [0, 3] for negative daughter, [4,7] for positive daughter)
Expand Down Expand Up @@ -88,7 +90,7 @@ DECLARE_SOA_TABLE(DataLnnCands, "AOD", "LNNCANDS",
lnnrec::PtPi, lnnrec::PhiPi, lnnrec::EtaPi,
lnnrec::XDecVtx, lnnrec::YDecVtx, lnnrec::ZDecVtx,
lnnrec::DcaV0Daug, lnnrec::Dca3H, lnnrec::DcaPi,
lnnrec::NSigma3H, lnnrec::NTPCclus3H, lnnrec::NTPCclusPi,
lnnrec::NSigma3H, lnnrec::NTPCclus3H, lnnrec::NTPCclusPi, lnnrec::NTPCNClsCrossedRows3H, lnnrec::NTPCNClsCrossedRowsPi,
lnnrec::TPCmom3H, lnnrec::TPCmomPi, lnnrec::TPCsignal3H, lnnrec::TPCsignalPi,
lnnrec::MassTrTOF, lnnrec::TPCchi3H,
lnnrec::ITSclusterSizes3H, lnnrec::ITSclusterSizesPi,
Expand All @@ -105,7 +107,7 @@ DECLARE_SOA_TABLE(MCLnnCands, "AOD", "MCLNNCANDS",
lnnrec::XDecVtx, lnnrec::YDecVtx, lnnrec::ZDecVtx,
lnnrec::DcaV0Daug, lnnrec::Dca3H, lnnrec::DcaPi,
lnnrec::NSigma3H, lnnrec::NTPCclus3H, lnnrec::NTPCclusPi,
lnnrec::TPCmom3H, lnnrec::TPCmomPi, lnnrec::TPCsignal3H, lnnrec::TPCsignalPi,
lnnrec::TPCmom3H, lnnrec::TPCmomPi, lnnrec::TPCsignal3H, lnnrec::TPCsignalPi, lnnrec::NTPCNClsCrossedRows3H, lnnrec::NTPCNClsCrossedRowsPi,
lnnrec::MassTrTOF, lnnrec::TPCchi3H,
lnnrec::ITSclusterSizes3H, lnnrec::ITSclusterSizesPi,
lnnrec::Flags,
Expand Down
35 changes: 22 additions & 13 deletions PWGLF/TableProducer/Nuspex/lnnRecoTask.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -168,6 +168,8 @@ struct LnnCandidate {
uint16_t tpcSignalPi = 0u;
uint8_t nTPCClusters3H = 0u;
uint8_t nTPCClustersPi = 0u;
uint8_t nTPCClustersCrossedRows3H = 0u;
uint8_t nTPCClustersCrossedRowsPi = 0u;
uint32_t clusterSizeITS3H = 0u;
uint32_t clusterSizeITSPi = 0u;
bool isMatter = false;
Expand Down Expand Up @@ -481,6 +483,7 @@ struct lnnRecoTask {
auto& h3track = lnnCand.isMatter ? posTrack : negTrack;
auto& pitrack = lnnCand.isMatter ? negTrack : posTrack;
auto& h3Rigidity = lnnCand.isMatter ? posRigidity : negRigidity;
auto& piRigidity = lnnCand.isMatter ? negRigidity : posRigidity;

// fill QA track histogram studies to check track signal before selections
h2FT0CnClusTPCtoTrBfSel->Fill(collision.centFT0C(), h3track.tpcNClsFound());
Expand Down Expand Up @@ -514,16 +517,20 @@ struct lnnRecoTask {
continue;
}

lnnCand.tpcChi3H = lnnCand.isMatter ? h3track.tpcChi2NCl() : negTrack.tpcChi2NCl();
lnnCand.tpcChi3H = h3track.tpcChi2NCl();
lnnCand.nSigma3H = lnnCand.isMatter ? nSigmaTPCpos : nSigmaTPCneg;
lnnCand.nTPCClusters3H = lnnCand.isMatter ? h3track.tpcNClsFound() : negTrack.tpcNClsFound();
lnnCand.tpcSignal3H = lnnCand.isMatter ? h3track.tpcSignal() : negTrack.tpcSignal();
lnnCand.clusterSizeITS3H = lnnCand.isMatter ? h3track.itsClusterSizes() : negTrack.itsClusterSizes();
lnnCand.nTPCClustersPi = !lnnCand.isMatter ? h3track.tpcNClsFound() : negTrack.tpcNClsFound();
lnnCand.tpcSignalPi = !lnnCand.isMatter ? h3track.tpcSignal() : negTrack.tpcSignal();
lnnCand.clusterSizeITSPi = !lnnCand.isMatter ? h3track.itsClusterSizes() : negTrack.itsClusterSizes();
lnnCand.mom3HTPC = lnnCand.isMatter ? posRigidity : negRigidity;
lnnCand.momPiTPC = !lnnCand.isMatter ? posRigidity : negRigidity;
lnnCand.nTPCClusters3H = h3track.tpcNClsFound();
lnnCand.nTPCClustersCrossedRows3H = h3track.tpcNClsCrossedRows();
lnnCand.tpcSignal3H = h3track.tpcSignal();
lnnCand.clusterSizeITS3H = h3track.itsClusterSizes();

lnnCand.nTPCClustersPi = pitrack.tpcNClsFound();
lnnCand.nTPCClustersCrossedRowsPi = pitrack.tpcNClsCrossedRows();
lnnCand.tpcSignalPi = pitrack.tpcSignal();
lnnCand.clusterSizeITSPi = pitrack.itsClusterSizes();

lnnCand.mom3HTPC = h3Rigidity;
lnnCand.momPiTPC = piRigidity;

lnnCand.flags |= lnnCand.isMatter ? static_cast<uint8_t>((posTrack.pidForTracking() & 0xF) << 4) : static_cast<uint8_t>((negTrack.pidForTracking() & 0xF) << 4);
lnnCand.flags |= lnnCand.isMatter ? static_cast<uint8_t>(negTrack.pidForTracking() & 0xF) : static_cast<uint8_t>(posTrack.pidForTracking() & 0xF);
Expand Down Expand Up @@ -633,7 +640,9 @@ struct lnnRecoTask {
// Fill 2D map after all selections
hdEdx3HSel->Fill(chargeFactor * lnnCand.mom3HTPC, h3track.tpcSignal());
hdEdx3HPosTrack->Fill(lnnCand.mom3HTPC, h3track.tpcSignal());
hdEdx3HNegTrack->Fill(-lnnCand.mom3HTPC, h3track.tpcSignal());
if (!lnnCand.isMatter) {
hdEdx3HNegTrack->Fill(-lnnCand.mom3HTPC, h3track.tpcSignal());
}
hNsigma3HSel->Fill(chargeFactor * lnnCand.mom3HTPC, lnnCand.nSigma3H);
if (h3track.hasTOF()) {
h3HSignalPtTOF->Fill(chargeFactor * h3track.pt(), beta);
Expand Down Expand Up @@ -788,7 +797,7 @@ struct lnnRecoTask {
lnnCand.recoPtPi(), lnnCand.recoPhiPi(), lnnCand.recoEtaPi(),
lnnCand.decVtx[0], lnnCand.decVtx[1], lnnCand.decVtx[2],
lnnCand.dcaV0dau, lnnCand.h3DCAXY, lnnCand.piDCAXY,
lnnCand.nSigma3H, lnnCand.nTPCClusters3H, lnnCand.nTPCClustersPi,
lnnCand.nSigma3H, lnnCand.nTPCClusters3H, lnnCand.nTPCClustersPi, lnnCand.nTPCClustersCrossedRows3H, lnnCand.nTPCClustersCrossedRowsPi,
lnnCand.mom3HTPC, lnnCand.momPiTPC, lnnCand.tpcSignal3H, lnnCand.tpcSignalPi,
lnnCand.mass2TrTOF, lnnCand.tpcChi3H,
lnnCand.clusterSizeITS3H, lnnCand.clusterSizeITSPi, lnnCand.flags);
Expand Down Expand Up @@ -879,7 +888,7 @@ struct lnnRecoTask {
lnnCand.recoPtPi(), lnnCand.recoPhiPi(), lnnCand.recoEtaPi(),
lnnCand.decVtx[0], lnnCand.decVtx[1], lnnCand.decVtx[2],
lnnCand.dcaV0dau, lnnCand.h3DCAXY, lnnCand.piDCAXY,
lnnCand.nSigma3H, lnnCand.nTPCClusters3H, lnnCand.nTPCClustersPi,
lnnCand.nSigma3H, lnnCand.nTPCClusters3H, lnnCand.nTPCClustersPi, lnnCand.nTPCClustersCrossedRows3H, lnnCand.nTPCClustersCrossedRowsPi,
lnnCand.mom3HTPC, lnnCand.momPiTPC, lnnCand.tpcSignal3H, lnnCand.tpcSignalPi,
lnnCand.mass2TrTOF, lnnCand.tpcChi3H,
lnnCand.clusterSizeITS3H, lnnCand.clusterSizeITSPi, lnnCand.flags,
Expand Down Expand Up @@ -966,7 +975,7 @@ struct lnnRecoTask {
-1, -1, -1,
-1, -1, -1,
-1, -1, -1,
-1, -1, -1,
-1, -1, -1, -1, -1,
-1, -1, -1, -1,
-1, -1,
-1, -1, -1,
Expand Down
Loading