Skip to content

Commit e06b310

Browse files
ArvindiitiArvindiiti
committed
[PWGJE] Fix event selection bits and added histograms for MC particles
1 parent cde8c99 commit e06b310

File tree

1 file changed

+83
-50
lines changed

1 file changed

+83
-50
lines changed

‎PWGJE/Tasks/nucleiInJets.cxx‎

Lines changed: 83 additions & 50 deletions
Original file line numberDiff line numberDiff line change
@@ -117,6 +117,7 @@ struct nucleiInJets {
117117
Configurable<int> centralityType{"centralityType", 0, "0: FT0M, 1: FT0C, 2: FV0A"};
118118
Configurable<std::vector<int>> cfgOccupancyRange{"cfgOccupancyRange", {0, 1000}, "Occupancy selection"};
119119
Configurable<bool> useOccupancy{"useOccupancy", true, "Events with custom occupancy selection"};
120+
Configurable<bool> useEtaSelForEffDen{"useEtaSelForEffDen", false, "eta selection for gen particles"};
120121

121122
Configurable<double> cfgtrkMinPt{"cfgtrkMinPt", 0.15, "set track min pT"};
122123
Configurable<double> cfgtrkMaxEta{"cfgtrkMaxEta", 0.8, "set track max Eta"};
@@ -158,6 +159,7 @@ struct nucleiInJets {
158159
Configurable<bool> useTOFVeto{"useTOFVeto", false, "true: use TOF veto, false: no TOF veto"};
159160
Configurable<bool> isRequireHitsInITSLayers{"isRequireHitsInITSLayers", true, "true: at least one hit in the its inner layes"};
160161
Configurable<bool> useMcC{"useMcC", true, "use mcC"};
162+
Configurable<bool> useRapidityCutForPID{"useRapidityCutForPID", false, "true: use rapidity cut for PID, false: no rapidity cut for PID"};
161163

162164
Configurable<bool> addpik{"addpik", true, "add pion and kaon hist"};
163165
ConfigurableAxis binsDCA{"binsDCA", {400, -1.f, 1.f}, ""};
@@ -1380,8 +1382,8 @@ struct nucleiInJets {
13801382
jetHist.fill(HIST("hNEventsInc"), 0.5);
13811383

13821384
bool isSel8 = jetderiveddatautilities::selectCollision(coll, jetderiveddatautilities::initialiseEventSelectionBits("sel8"));
1383-
bool isSelNoSameBunchPileup = jetderiveddatautilities::selectCollision(coll, jetderiveddatautilities::initialiseEventSelectionBits("selNoSameBunchPileup"));
1384-
bool isSelIsGoodZvtxFT0vsPV = jetderiveddatautilities::selectCollision(coll, jetderiveddatautilities::initialiseEventSelectionBits("selIsGoodZvtxFT0vsPV"));
1385+
bool isSelNoSameBunchPileup = jetderiveddatautilities::selectCollision(coll, jetderiveddatautilities::initialiseEventSelectionBits("NoSameBunchPileup"));
1386+
bool isSelIsGoodZvtxFT0vsPV = jetderiveddatautilities::selectCollision(coll, jetderiveddatautilities::initialiseEventSelectionBits("IsGoodZvtxFT0vsPV"));
13851387

13861388
if (sel8Coll && !isSel8)
13871389
return;
@@ -1421,113 +1423,137 @@ struct nucleiInJets {
14211423
if (!isTrackSelected(trk)) {
14221424
continue;
14231425
}
1426+
1427+
auto rapidityData = [&](float m2z) {
1428+
const float rap = trk.rapidity(m2z);
1429+
return rap > std::abs(cfgtrkMaxRap);
1430+
};
1431+
1432+
auto prRapidityWithinRange = rapidityData(o2::constants::physics::MassProton);
1433+
auto deRapidityWithinRange = rapidityData(o2::constants::physics::MassDeuteron);
1434+
14241435
if (std::fabs(trk.eta()) > cfgtrkMaxEta)
14251436
continue;
14261437
if (trk.sign() > 0) { // particle info
14271438
if (useTOFNsigmaPreSel && trk.hasTOF()) {
1428-
if (std::abs(trk.tofNSigmaPr()) < cfgnTPCPIDPrTOF)
1439+
if (std::abs(trk.tofNSigmaPr()) < cfgnTPCPIDPrTOF && (!useRapidityCutForPID || prRapidityWithinRange)) {
14291440
jetHist.fill(HIST("tracksInc/proton/h3PtVsProtonNSigmaTPCVsPt"), trk.pt(), trk.tpcNSigmaPr(), centrality);
1430-
if (std::abs(trk.tofNSigmaDe()) < cfgnTPCPIDDeTOF)
1441+
}
1442+
if (std::abs(trk.tofNSigmaDe()) < cfgnTPCPIDDeTOF && (!useRapidityCutForPID || deRapidityWithinRange)) {
14311443
jetHist.fill(HIST("tracksInc/deuteron/h3PtVsDeuteronNSigmaTPCVsPt"), trk.pt(), trk.tpcNSigmaDe(), centrality);
1444+
}
14321445
} else if (!useTOFNsigmaPreSel && !useTOFVeto) {
1433-
jetHist.fill(HIST("tracksInc/proton/h3PtVsProtonNSigmaTPCVsPt"), trk.pt(), trk.tpcNSigmaPr(), centrality);
1434-
jetHist.fill(HIST("tracksInc/deuteron/h3PtVsDeuteronNSigmaTPCVsPt"), trk.pt(), trk.tpcNSigmaDe(), centrality);
1435-
1446+
if (!useRapidityCutForPID || prRapidityWithinRange) {
1447+
jetHist.fill(HIST("tracksInc/proton/h3PtVsProtonNSigmaTPCVsPt"), trk.pt(), trk.tpcNSigmaPr(), centrality);
1448+
}
1449+
if (!useRapidityCutForPID || deRapidityWithinRange) {
1450+
jetHist.fill(HIST("tracksInc/deuteron/h3PtVsDeuteronNSigmaTPCVsPt"), trk.pt(), trk.tpcNSigmaDe(), centrality);
1451+
}
14361452
} else if (!useTOFNsigmaPreSel && useTOFVeto) {
14371453
if (trk.hasTOF()) {
14381454
if (std::abs(trk.tofNSigmaPr()) < cfgnTPCPIDPrTOF) {
1439-
jetHist.fill(HIST("tracksInc/proton/h3PtVsProtonNSigmaTPCVsPt"), trk.pt(), trk.tpcNSigmaPr(), centrality);
1455+
if (!useRapidityCutForPID || prRapidityWithinRange)
1456+
jetHist.fill(HIST("tracksInc/proton/h3PtVsProtonNSigmaTPCVsPt"), trk.pt(), trk.tpcNSigmaPr(), centrality);
14401457
}
14411458
} else {
1442-
jetHist.fill(HIST("tracksInc/proton/h3PtVsProtonNSigmaTPCVsPt"), trk.pt(), trk.tpcNSigmaPr(), centrality);
1459+
if (!useRapidityCutForPID || prRapidityWithinRange)
1460+
jetHist.fill(HIST("tracksInc/proton/h3PtVsProtonNSigmaTPCVsPt"), trk.pt(), trk.tpcNSigmaPr(), centrality);
14431461
}
14441462
if (trk.hasTOF()) {
14451463
if (std::abs(trk.tofNSigmaDe()) < cfgnTPCPIDDeTOF) {
1446-
jetHist.fill(HIST("tracksInc/deuteron/h3PtVsDeuteronNSigmaTPCVsPt"), trk.pt(), trk.tpcNSigmaDe(), centrality);
1464+
if (!useRapidityCutForPID || deRapidityWithinRange)
1465+
jetHist.fill(HIST("tracksInc/deuteron/h3PtVsDeuteronNSigmaTPCVsPt"), trk.pt(), trk.tpcNSigmaDe(), centrality);
14471466
}
14481467
} else {
1449-
jetHist.fill(HIST("tracksInc/deuteron/h3PtVsDeuteronNSigmaTPCVsPt"), trk.pt(), trk.tpcNSigmaDe(), centrality);
1468+
if (!useRapidityCutForPID || deRapidityWithinRange)
1469+
jetHist.fill(HIST("tracksInc/deuteron/h3PtVsDeuteronNSigmaTPCVsPt"), trk.pt(), trk.tpcNSigmaDe(), centrality);
14501470
}
14511471
}
14521472
if (addTOFplots && trk.hasTOF()) {
14531473
float massTOF = trk.p() * std::sqrt(1.f / (trk.beta() * trk.beta()) - 1.f);
14541474
if (!useTPCpreSel) {
1455-
jetHist.fill(HIST("tracksInc/proton/h2TOFmassProtonVsPt"), massTOF, trk.pt(), centrality);
1456-
jetHist.fill(HIST("tracksInc/proton/h2TOFmass2ProtonVsPt"), massTOF * massTOF - MassProton * MassProton, trk.pt(), centrality);
1457-
jetHist.fill(HIST("tracksInc/proton/h2TofNsigmaProtonVsPt"), trk.tofNSigmaPr(), trk.pt(), centrality);
1458-
1459-
jetHist.fill(HIST("tracksInc/deuteron/h2TOFmassDeuteronVsPt"), massTOF, trk.pt(), centrality);
1460-
jetHist.fill(HIST("tracksInc/deuteron/h2TOFmass2DeuteronVsPt"), massTOF * massTOF - MassDeuteron * MassDeuteron, trk.pt(), centrality);
1461-
jetHist.fill(HIST("tracksInc/deuteron/h2TofNsigmaDeuteronVsPt"), trk.tofNSigmaDe(), trk.pt(), centrality);
1475+
if (!useRapidityCutForPID || prRapidityWithinRange) {
1476+
jetHist.fill(HIST("tracksInc/proton/h2TOFmassProtonVsPt"), massTOF, trk.pt(), centrality);
1477+
jetHist.fill(HIST("tracksInc/proton/h2TOFmass2ProtonVsPt"), massTOF * massTOF - MassProton * MassProton, trk.pt(), centrality);
1478+
jetHist.fill(HIST("tracksInc/proton/h2TofNsigmaProtonVsPt"), trk.tofNSigmaPr(), trk.pt(), centrality);
1479+
}
1480+
if (!useRapidityCutForPID || deRapidityWithinRange) {
1481+
jetHist.fill(HIST("tracksInc/deuteron/h2TOFmassDeuteronVsPt"), massTOF, trk.pt(), centrality);
1482+
jetHist.fill(HIST("tracksInc/deuteron/h2TOFmass2DeuteronVsPt"), massTOF * massTOF - MassDeuteron * MassDeuteron, trk.pt(), centrality);
1483+
jetHist.fill(HIST("tracksInc/deuteron/h2TofNsigmaDeuteronVsPt"), trk.tofNSigmaDe(), trk.pt(), centrality);
1484+
}
14621485
} else {
1463-
if (std::abs(trk.tpcNSigmaPr()) < cfgnTPCPIDPr) {
1486+
if (std::abs(trk.tpcNSigmaPr()) < cfgnTPCPIDPr && (!useRapidityCutForPID || prRapidityWithinRange)) {
14641487
jetHist.fill(HIST("tracksInc/proton/h2TOFmassProtonVsPt"), massTOF, trk.pt(), centrality);
14651488
jetHist.fill(HIST("tracksInc/proton/h2TOFmass2ProtonVsPt"), massTOF * massTOF - MassProton * MassProton, trk.pt(), centrality);
14661489
jetHist.fill(HIST("tracksInc/proton/h2TofNsigmaProtonVsPt"), trk.tofNSigmaPr(), trk.pt(), centrality);
14671490
}
1468-
if (std::abs(trk.tpcNSigmaDe()) < cfgnTPCPIDDe) {
1491+
if (std::abs(trk.tpcNSigmaDe()) < cfgnTPCPIDDe && (!useRapidityCutForPID || deRapidityWithinRange)) {
14691492
jetHist.fill(HIST("tracksInc/deuteron/h2TOFmassDeuteronVsPt"), massTOF, trk.pt(), centrality);
14701493
jetHist.fill(HIST("tracksInc/deuteron/h2TOFmass2DeuteronVsPt"), massTOF * massTOF - MassDeuteron * MassDeuteron, trk.pt(), centrality);
14711494
jetHist.fill(HIST("tracksInc/deuteron/h2TofNsigmaDeuteronVsPt"), trk.tofNSigmaDe(), trk.pt(), centrality);
14721495
}
14731496
}
14741497
}
14751498

1476-
if (cEnableProtonQA && std::abs(trk.tpcNSigmaPr()) < cfgnTPCPIDPr) {
1499+
if (cEnableProtonQA && std::abs(trk.tpcNSigmaPr()) < cfgnTPCPIDPr && (!useRapidityCutForPID || prRapidityWithinRange)) {
14771500
jetHist.fill(HIST("tracksInc/proton/dca/after/hDCAxyVsPtProton"), trk.dcaXY(), trk.pt(), centrality);
14781501
jetHist.fill(HIST("tracksInc/proton/dca/after/hDCAzVsPtProton"), trk.dcaZ(), trk.pt(), centrality);
14791502
}
1480-
if (cEnableDeuteronQA && std::abs(trk.tpcNSigmaDe()) < cfgnTPCPIDDe) {
1503+
if (cEnableDeuteronQA && std::abs(trk.tpcNSigmaDe()) < cfgnTPCPIDDe && (!useRapidityCutForPID || deRapidityWithinRange)) {
14811504
jetHist.fill(HIST("tracksInc/deuteron/dca/after/hDCAxyVsPtDeuteron"), trk.dcaXY(), trk.pt(), centrality);
14821505
jetHist.fill(HIST("tracksInc/deuteron/dca/after/hDCAzVsPtDeuteron"), trk.dcaZ(), trk.pt(), centrality);
14831506
}
14841507

14851508
} else { // anti-particle info
14861509
if (useTOFNsigmaPreSel && trk.hasTOF()) {
1487-
if (std::abs(trk.tofNSigmaPr()) < cfgnTPCPIDPrTOF)
1510+
if (std::abs(trk.tofNSigmaPr()) < cfgnTPCPIDPrTOF && (!useRapidityCutForPID || prRapidityWithinRange))
14881511
jetHist.fill(HIST("tracksInc/antiProton/h3PtVsantiProtonNSigmaTPCVsPt"), trk.pt(), trk.tpcNSigmaPr(), centrality);
1489-
if (std::abs(trk.tofNSigmaDe()) < cfgnTPCPIDDeTOF)
1512+
if (std::abs(trk.tofNSigmaDe()) < cfgnTPCPIDDeTOF && (!useRapidityCutForPID || deRapidityWithinRange))
14901513
jetHist.fill(HIST("tracksInc/antiDeuteron/h3PtVsantiDeuteronNSigmaTPCVsPt"), trk.pt(), trk.tpcNSigmaDe(), centrality);
14911514
} else if (!useTOFNsigmaPreSel && !useTOFVeto) {
1492-
jetHist.fill(HIST("tracksInc/antiProton/h3PtVsantiProtonNSigmaTPCVsPt"), trk.pt(), trk.tpcNSigmaPr(), centrality);
1493-
jetHist.fill(HIST("tracksInc/antiDeuteron/h3PtVsantiDeuteronNSigmaTPCVsPt"), trk.pt(), trk.tpcNSigmaDe(), centrality);
1515+
if (!useRapidityCutForPID || prRapidityWithinRange)
1516+
jetHist.fill(HIST("tracksInc/antiProton/h3PtVsantiProtonNSigmaTPCVsPt"), trk.pt(), trk.tpcNSigmaPr(), centrality);
1517+
if (!useRapidityCutForPID || deRapidityWithinRange)
1518+
jetHist.fill(HIST("tracksInc/antiDeuteron/h3PtVsantiDeuteronNSigmaTPCVsPt"), trk.pt(), trk.tpcNSigmaDe(), centrality);
14941519
} else if (!useTOFNsigmaPreSel && useTOFVeto) {
14951520
if (trk.hasTOF()) {
1496-
if (std::abs(trk.tofNSigmaPr()) < cfgnTPCPIDPrTOF) {
1521+
if (std::abs(trk.tofNSigmaPr()) < cfgnTPCPIDPrTOF && (!useRapidityCutForPID || prRapidityWithinRange)) {
14971522
jetHist.fill(HIST("tracksInc/antiProton/h3PtVsantiProtonNSigmaTPCVsPt"), trk.pt(), trk.tpcNSigmaPr(), centrality);
14981523
}
14991524
} else {
1500-
jetHist.fill(HIST("tracksInc/antiProton/h3PtVsantiProtonNSigmaTPCVsPt"), trk.pt(), trk.tpcNSigmaPr(), centrality);
1525+
if (!useRapidityCutForPID || prRapidityWithinRange)
1526+
jetHist.fill(HIST("tracksInc/antiProton/h3PtVsantiProtonNSigmaTPCVsPt"), trk.pt(), trk.tpcNSigmaPr(), centrality);
15011527
}
15021528
if (trk.hasTOF()) {
1503-
if (std::abs(trk.tofNSigmaDe()) < cfgnTPCPIDDeTOF) {
1529+
if (std::abs(trk.tofNSigmaDe()) < cfgnTPCPIDDeTOF && (!useRapidityCutForPID || deRapidityWithinRange)) {
15041530
jetHist.fill(HIST("tracksInc/antiDeuteron/h3PtVsantiDeuteronNSigmaTPCVsPt"), trk.pt(), trk.tpcNSigmaDe(), centrality);
15051531
}
15061532
} else {
1507-
jetHist.fill(HIST("tracksInc/antiDeuteron/h3PtVsantiDeuteronNSigmaTPCVsPt"), trk.pt(), trk.tpcNSigmaDe(), centrality);
1533+
if (!useRapidityCutForPID || deRapidityWithinRange)
1534+
jetHist.fill(HIST("tracksInc/antiDeuteron/h3PtVsantiDeuteronNSigmaTPCVsPt"), trk.pt(), trk.tpcNSigmaDe(), centrality);
15081535
}
15091536
}
15101537
if (addTOFplots && trk.hasTOF()) {
15111538
float massTOF = trk.p() * std::sqrt(1.f / (trk.beta() * trk.beta()) - 1.f);
15121539
if (!useTPCpreSel) {
1513-
jetHist.fill(HIST("tracksInc/antiProton/h2TOFmassantiProtonVsPt"), massTOF, trk.pt(), centrality);
1514-
jetHist.fill(HIST("tracksInc/antiProton/h2TOFmass2antiProtonVsPt"), massTOF * massTOF - MassProton * MassProton, trk.pt(), centrality);
1515-
jetHist.fill(HIST("tracksInc/antiProton/h2TofNsigmaantiProtonVsPt"), trk.tofNSigmaPr(), trk.pt(), centrality);
1516-
1517-
jetHist.fill(HIST("tracksInc/antiDeuteron/h2TOFmassantiDeuteronVsPt"), massTOF, trk.pt(), centrality);
1518-
jetHist.fill(HIST("tracksInc/antiDeuteron/h2TOFmass2antiDeuteronVsPt"), massTOF * massTOF - MassDeuteron * MassDeuteron, trk.pt(), centrality);
1519-
jetHist.fill(HIST("tracksInc/antiDeuteron/h2TofNsigmaantiDeuteronVsPt"), trk.tofNSigmaDe(), trk.pt(), centrality);
1520-
} else {
1521-
if (std::abs(trk.tpcNSigmaPr()) < cfgnTPCPIDPr) {
1540+
if (!useRapidityCutForPID || prRapidityWithinRange) {
15221541
jetHist.fill(HIST("tracksInc/antiProton/h2TOFmassantiProtonVsPt"), massTOF, trk.pt(), centrality);
15231542
jetHist.fill(HIST("tracksInc/antiProton/h2TOFmass2antiProtonVsPt"), massTOF * massTOF - MassProton * MassProton, trk.pt(), centrality);
15241543
jetHist.fill(HIST("tracksInc/antiProton/h2TofNsigmaantiProtonVsPt"), trk.tofNSigmaPr(), trk.pt(), centrality);
15251544
}
1526-
if (std::abs(trk.tpcNSigmaDe()) < cfgnTPCPIDDe) {
1545+
if (!useRapidityCutForPID || deRapidityWithinRange) {
15271546
jetHist.fill(HIST("tracksInc/antiDeuteron/h2TOFmassantiDeuteronVsPt"), massTOF, trk.pt(), centrality);
15281547
jetHist.fill(HIST("tracksInc/antiDeuteron/h2TOFmass2antiDeuteronVsPt"), massTOF * massTOF - MassDeuteron * MassDeuteron, trk.pt(), centrality);
15291548
jetHist.fill(HIST("tracksInc/antiDeuteron/h2TofNsigmaantiDeuteronVsPt"), trk.tofNSigmaDe(), trk.pt(), centrality);
1530-
} else {
1549+
}
1550+
} else {
1551+
if (std::abs(trk.tpcNSigmaPr()) < cfgnTPCPIDPr && (!useRapidityCutForPID || prRapidityWithinRange)) {
1552+
jetHist.fill(HIST("tracksInc/antiProton/h2TOFmassantiProtonVsPt"), massTOF, trk.pt(), centrality);
1553+
jetHist.fill(HIST("tracksInc/antiProton/h2TOFmass2antiProtonVsPt"), massTOF * massTOF - MassProton * MassProton, trk.pt(), centrality);
1554+
jetHist.fill(HIST("tracksInc/antiProton/h2TofNsigmaantiProtonVsPt"), trk.tofNSigmaPr(), trk.pt(), centrality);
1555+
}
1556+
if (std::abs(trk.tpcNSigmaDe()) < cfgnTPCPIDDe && (!useRapidityCutForPID || deRapidityWithinRange)) {
15311557
jetHist.fill(HIST("tracksInc/antiDeuteron/h2TOFmassantiDeuteronVsPt"), massTOF, trk.pt(), centrality);
15321558
jetHist.fill(HIST("tracksInc/antiDeuteron/h2TOFmass2antiDeuteronVsPt"), massTOF * massTOF - MassDeuteron * MassDeuteron, trk.pt(), centrality);
15331559
jetHist.fill(HIST("tracksInc/antiDeuteron/h2TofNsigmaantiDeuteronVsPt"), trk.tofNSigmaDe(), trk.pt(), centrality);
@@ -2011,8 +2037,8 @@ struct nucleiInJets {
20112037
jetHist.fill(HIST("recInc/eventStat"), 0.5);
20122038

20132039
bool isSel8 = jetderiveddatautilities::selectCollision(coll, jetderiveddatautilities::initialiseEventSelectionBits("sel8"));
2014-
bool isSelNoSameBunchPileup = jetderiveddatautilities::selectCollision(coll, jetderiveddatautilities::initialiseEventSelectionBits("selNoSameBunchPileup"));
2015-
bool isSelIsGoodZvtxFT0vsPV = jetderiveddatautilities::selectCollision(coll, jetderiveddatautilities::initialiseEventSelectionBits("selIsGoodZvtxFT0vsPV"));
2040+
bool isSelNoSameBunchPileup = jetderiveddatautilities::selectCollision(coll, jetderiveddatautilities::initialiseEventSelectionBits("NoSameBunchPileup"));
2041+
bool isSelIsGoodZvtxFT0vsPV = jetderiveddatautilities::selectCollision(coll, jetderiveddatautilities::initialiseEventSelectionBits("IsGoodZvtxFT0vsPV"));
20162042

20172043
if (sel8Coll && !isSel8)
20182044
return;
@@ -2065,13 +2091,20 @@ struct nucleiInJets {
20652091
// require mc getProcess to get Decay and Material secondaries
20662092
int particleOriginType = 0;
20672093
auto isMcPrimary = false;
2068-
// auto isMcTransport = false; // auto isMcSecondaryFromMaterial = false; auto isMcSecondaryFromWeakDecay = false;
2069-
if (mcTrack.isPhysicalPrimary()) {
2070-
isMcPrimary = true;
2094+
auto isProdByGen = false;
2095+
auto isFromWeakDecay = false;
2096+
2097+
isMcPrimary = mcTrack.isPhysicalPrimary();
2098+
isProdByGen = mcTrack.producedByGenerator();
2099+
isFromWeakDecay = mcTrack.getProcess() == TMCProcess::kPDecay;
2100+
2101+
if (isMcPrimary) {
20712102
particleOriginType = 1;
2072-
} else if (mcTrack.getGenStatusCode() == -1) {
2073-
// isMcTransport = true;
2074-
particleOriginType = 2;
2103+
} else if (!isProdByGen) {
2104+
particleOriginType = 2; // from transport
2105+
if (isFromWeakDecay) {
2106+
particleOriginType = 3; // from weak decay
2107+
}
20752108
}
20762109

20772110
// Fill DCAxy histograms
@@ -2171,7 +2204,7 @@ struct nucleiInJets {
21712204
for (const auto& mcParticle : mcParticles_per_coll) {
21722205
if (!mcParticle.isPhysicalPrimary())
21732206
continue;
2174-
if (std::fabs(mcParticle.eta()) > cfgtrkMaxEta)
2207+
if (std::fabs(mcParticle.eta()) > cfgtrkMaxEta && useEtaSelForEffDen)
21752208
continue;
21762209
if (std::fabs(mcParticle.y()) > cfgtrkMaxRap)
21772210
continue;

0 commit comments

Comments
 (0)