Skip to content

Commit a42be47

Browse files
[PWGLF] adding only primary in mc for pT dependent cut in DCA (#17142)
Co-authored-by: ALICE Action Bot <alibuild@cern.ch>
1 parent c36ea5e commit a42be47

1 file changed

Lines changed: 19 additions & 17 deletions

File tree

PWGLF/Tasks/Nuspex/nucleitpcpbpb.cxx

Lines changed: 19 additions & 17 deletions
Original file line numberDiff line numberDiff line change
@@ -918,9 +918,24 @@ struct NucleitpcPbPb {
918918
histomc.fill(HIST("hSpectramc"), particleAnti, collision.centFT0C(), ptReco, ptTOF);
919919

920920
// Basic track histograms
921-
histos.fill(HIST("dcaXY"), ptReco, track.dcaXY());
922-
histos.fill(HIST("dcaZ"), ptReco, track.dcaZ());
923-
histos.fill(HIST("Tpcsignal"), getRigidity(track) * track.sign(), track.tpcSignal());
921+
if (decayType == 0) {
922+
histos.fill(HIST("dcaXY"), ptReco, track.dcaXY());
923+
histos.fill(HIST("dcaZ"), ptReco, track.dcaZ());
924+
histos.fill(HIST("Tpcsignal"), getRigidity(track) * track.sign(), track.tpcSignal());
925+
}
926+
// Delta Pt histograms
927+
float ptGen = matchedMCParticle.pt();
928+
float deltaPt = ptReco - ptGen;
929+
930+
if (pdg == -particlePdgCodes.at(i) && decayType == 0) { // Anti-particle
931+
histomc.fill(HIST("histDeltaPtVsPtGenanti"), ptReco, deltaPt);
932+
histomc.fill(HIST("histPIDtrackanti"), ptReco, track.pidForTracking());
933+
}
934+
935+
if (pdg == particlePdgCodes.at(i) && decayType == 0) { // Particle
936+
histomc.fill(HIST("histDeltaPtVsPtGen"), ptReco, deltaPt);
937+
histomc.fill(HIST("histPIDtrack"), ptReco, track.pidForTracking());
938+
}
924939

925940
// Fill mass²/z² for MC - separate for particles and anti-particles
926941
if (track.hasTOF()) {
@@ -940,7 +955,7 @@ struct NucleitpcPbPb {
940955
}
941956
}
942957

943-
if (!skipHe4) {
958+
if (!skipHe4 && decayType == 0) {
944959
if (pdg > 0) {
945960
histomc.fill(HIST("hMassVsPtMC"), ptReco, massSquareOverChargeSquare, collision.centFT0C());
946961
} else {
@@ -949,19 +964,6 @@ struct NucleitpcPbPb {
949964
}
950965
}
951966
}
952-
953-
// Delta Pt histograms
954-
float ptGen = matchedMCParticle.pt();
955-
float deltaPt = ptReco - ptGen;
956-
957-
if (pdg == -particlePdgCodes.at(i)) { // Anti-particle
958-
histomc.fill(HIST("histDeltaPtVsPtGenanti"), ptReco, deltaPt);
959-
histomc.fill(HIST("histPIDtrackanti"), ptReco, track.pidForTracking());
960-
}
961-
if (pdg == particlePdgCodes.at(i)) { // Particle
962-
histomc.fill(HIST("histDeltaPtVsPtGen"), ptReco, deltaPt);
963-
histomc.fill(HIST("histPIDtrack"), ptReco, track.pidForTracking());
964-
}
965967
}
966968
}
967969
break; // Found the matching collision, break out of collision loop

0 commit comments

Comments
 (0)