@@ -332,6 +332,8 @@ struct NucleitpcPbPb {
332332 // NEW: Add DCAxy vs pT histograms for MC
333333 histomc.add (" DCAxy_vs_pT_transport" , " DCA_{xy} vs p_{T} (Transport);p_{T} (GeV/c);DCA_{xy} (cm)" ,
334334 {HistType::kTH3F , {ptAxis, axisDCA, axisCent}});
335+ histomc.add (" DCAxy_vs_pT_primary" , " DCA_{xy} vs p_{T} (primary);p_{T} (GeV/c);DCA_{xy} (cm)" ,
336+ {HistType::kTH3F , {ptAxis, axisDCA, axisCent}});
335337 histomc.add (" DCAxy_vs_pT_weakdecay" , " DCA_{xy} vs p_{T} (Weak Decay);p_{T} (GeV/c);DCA_{xy} (cm)" ,
336338 {HistType::kTH3F , {ptAxis, axisDCA, axisCent}});
337339 histomc.add (" DCAxy_vs_pT_total" , " DCA_{xy} vs p_{T} (Total);p_{T} (GeV/c);DCA_{xy} (cm)" ,
@@ -1117,6 +1119,9 @@ struct NucleitpcPbPb {
11171119
11181120 if (pdg == particlePdgCodes.at (i)) { // He3
11191121 histomc.fill (HIST (" DCAxy_vs_pT_total" ), ptDCA, track.dcaXY (), collision.centFT0C ());
1122+ if (matchedMCParticle.isPhysicalPrimary ()) {
1123+ histomc.fill (HIST (" DCAxy_vs_pT_primary" ), ptDCA, track.dcaXY (), collision.centFT0C ());
1124+ }
11201125 if (decayType == 2 ) { // Transport/Material
11211126 histomc.fill (HIST (" DCAxy_vs_pT_transport" ), ptDCA, track.dcaXY (), collision.centFT0C ());
11221127 } else if (decayType == 1 ) { // Weak decay (including HF)
0 commit comments