Skip to content

Commit 9013641

Browse files
committed
add momentum cuts to charged tracks
1 parent d59ce93 commit 9013641

File tree

1 file changed

+4
-4
lines changed

1 file changed

+4
-4
lines changed

PWGLF/TableProducer/Nuspex/ebyeMaker.cxx

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -571,8 +571,8 @@ struct EbyeMaker {
571571
continue;
572572
}
573573
histos.fill(HIST("QA/tpcSignal"), track.tpcInnerParam(), track.tpcSignal());
574-
575-
nTracksColl++;
574+
if (trackPt > ptMin[0] && trackPt < ptMax[0])
575+
nTracksColl++;
576576

577577
for (int iP{0}; iP < kNpart; ++iP) {
578578
if (trackPt < ptMin[iP] || trackPt > ptMax[iP]) {
@@ -853,7 +853,8 @@ struct EbyeMaker {
853853
if ((((mcPart.flags() & 0x8) || (mcPart.flags() & 0x2)) && (doprocessMcRun2 || doprocessMiniMcRun2)) || ((mcPart.flags() & 0x1) && !doprocessMiniMcRun2))
854854
continue;
855855
auto pdgCode = mcPart.pdgCode();
856-
if ((std::abs(pdgCode) == PDG_t::kPiPlus || std::abs(pdgCode) == PDG_t::kElectron || std::abs(pdgCode) == PDG_t::kMuonMinus || std::abs(pdgCode) == PDG_t::kKPlus || std::abs(pdgCode) == PDG_t::kProton) && mcPart.isPhysicalPrimary())
856+
auto genPt = std::hypot(mcPart.px(), mcPart.py());
857+
if ((std::abs(pdgCode) == PDG_t::kPiPlus || std::abs(pdgCode) == PDG_t::kElectron || std::abs(pdgCode) == PDG_t::kMuonMinus || std::abs(pdgCode) == PDG_t::kKPlus || std::abs(pdgCode) == PDG_t::kProton) && mcPart.isPhysicalPrimary() && genPt > ptMin[0] && genPt < ptMax[0])
857858
nChPartGen++;
858859
if (std::abs(pdgCode) == PDG_t::kLambda0) {
859860
if (!mcPart.isPhysicalPrimary() && !mcPart.has_mothers())
@@ -868,7 +869,6 @@ struct EbyeMaker {
868869
if (!foundPr) {
869870
continue;
870871
}
871-
auto genPt = std::hypot(mcPart.px(), mcPart.py());
872872
CandidateV0 candV0;
873873
candV0.genpt = genPt;
874874
candV0.geneta = mcPart.eta();

0 commit comments

Comments
 (0)