Skip to content

Commit f85296f

Browse files
committed
feat: Add TPC PID cluster information for He3 and Pi
This commit adds the number of TPC PID clusters for the He3 and Pi tracks to the hypCand struct. This information is used to better characterize the particle identification of the tracks and improve the analysis of hypertriton candidates.
1 parent 2888347 commit f85296f

File tree

1 file changed

+8
-2
lines changed

1 file changed

+8
-2
lines changed

PWGLF/TableProducer/Nuspex/hyperRecoTask.cxx

Lines changed: 8 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -113,6 +113,8 @@ struct hyperCandidate {
113113
float massTOFHe3 = 0.f;
114114
uint8_t nTPCClustersHe3 = 0u;
115115
uint8_t nTPCClustersPi = 0u;
116+
uint8_t nTPCpidClusHe3 = 0u;
117+
uint8_t nTPCpidClusPi = 0u;
116118
uint32_t clusterSizeITSHe3 = 0u;
117119
uint32_t clusterSizeITSPi = 0u;
118120

@@ -405,8 +407,10 @@ struct hyperRecoTask {
405407
hypCand.nSigmaHe3 = computeNSigmaHe3(heTrack);
406408
hypCand.nTPCClustersHe3 = heTrack.tpcNClsFound();
407409
hypCand.tpcSignalHe3 = heTrack.tpcSignal();
410+
hypCand.nTPCpidClusHe3 = (int16_t)heTrack.tpcNClsFindable() - heTrack.tpcNClsFindableMinusPID();
408411
hypCand.clusterSizeITSHe3 = heTrack.itsClusterSizes();
409412
hypCand.nTPCClustersPi = piTrack.tpcNClsFound();
413+
hypCand.nTPCpidClusPi = (int16_t)piTrack.tpcNClsFindable() - piTrack.tpcNClsFindableMinusPID();
410414
hypCand.tpcSignalPi = piTrack.tpcSignal();
411415
hypCand.tpcChi2He3 = heTrack.tpcChi2NCl();
412416
hypCand.clusterSizeITSPi = piTrack.itsClusterSizes();
@@ -682,6 +686,7 @@ struct hyperRecoTask {
682686
hypCand.decVtx[0], hypCand.decVtx[1], hypCand.decVtx[2],
683687
hypCand.dcaV0dau, hypCand.he3DCAXY, hypCand.piDCAXY,
684688
hypCand.nSigmaHe3, hypCand.nTPCClustersHe3, hypCand.nTPCClustersPi,
689+
hypCand.nTPCpidClusHe3, hypCand.nTPCpidClusPi,
685690
hypCand.momHe3TPC, hypCand.momPiTPC, hypCand.tpcSignalHe3, hypCand.tpcSignalPi, hypCand.tpcChi2He3,
686691
hypCand.massTOFHe3,
687692
hypCand.clusterSizeITSHe3, hypCand.clusterSizeITSPi, hypCand.flags, trackedHypClSize);
@@ -716,6 +721,7 @@ struct hyperRecoTask {
716721
hypCand.decVtx[0], hypCand.decVtx[1], hypCand.decVtx[2],
717722
hypCand.dcaV0dau, hypCand.he3DCAXY, hypCand.piDCAXY,
718723
hypCand.nSigmaHe3, hypCand.nTPCClustersHe3, hypCand.nTPCClustersPi,
724+
hypCand.nTPCpidClusHe3, hypCand.nTPCpidClusPi,
719725
hypCand.momHe3TPC, hypCand.momPiTPC, hypCand.tpcSignalHe3, hypCand.tpcSignalPi, hypCand.tpcChi2He3,
720726
hypCand.massTOFHe3,
721727
hypCand.clusterSizeITSHe3, hypCand.clusterSizeITSPi, hypCand.flags, trackedHypClSize);
@@ -750,7 +756,7 @@ struct hyperRecoTask {
750756
hypCand.recoPtPi(), hypCand.recoPhiPi(), hypCand.recoEtaPi(),
751757
hypCand.decVtx[0], hypCand.decVtx[1], hypCand.decVtx[2],
752758
hypCand.dcaV0dau, hypCand.he3DCAXY, hypCand.piDCAXY,
753-
hypCand.nSigmaHe3, hypCand.nTPCClustersHe3, hypCand.nTPCClustersPi,
759+
hypCand.nSigmaHe3, hypCand.nTPCClustersHe3, hypCand.nTPCClustersPi, hypCand.nTPCpidClusHe3, hypCand.nTPCpidClusPi,
754760
hypCand.momHe3TPC, hypCand.momPiTPC, hypCand.tpcSignalHe3, hypCand.tpcSignalPi, hypCand.tpcChi2He3,
755761
hypCand.massTOFHe3,
756762
hypCand.clusterSizeITSHe3, hypCand.clusterSizeITSPi, hypCand.flags, trackedHypClSize,
@@ -825,7 +831,7 @@ struct hyperRecoTask {
825831
-1, -1, -1,
826832
-1, -1, -1,
827833
-1, -1, -1,
828-
-1, -1, -1, -1, -1, -1,
834+
-1, -1, -1, -1, -1, -1, 0, 0,
829835
-1, -1, -1, false,
830836
chargeFactor * hypCand.genPt(), hypCand.genPhi(), hypCand.genEta(), hypCand.genPtHe3(),
831837
hypCand.gDecVtx[0], hypCand.gDecVtx[1], hypCand.gDecVtx[2],

0 commit comments

Comments
 (0)