Skip to content

Commit 98163be

Browse files
fcatalan92gconesab
authored andcommitted
Small fix to heavy-quark rapidity check
1 parent 8130153 commit 98163be

1 file changed

Lines changed: 10 additions & 4 deletions

File tree

MC/config/PWGHF/external/generator/generator_pythia8_gaptriggered_hf.C

Lines changed: 10 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -88,12 +88,18 @@ protected:
8888
auto daughterList = event[iPart].daughterList();
8989
bool hasQ = false, hasQbar = false, atSelectedY = false;
9090
for (auto iDau : daughterList) {
91-
if (event[iDau].id() == mQuarkPdg)
91+
if (event[iDau].id() == mQuarkPdg) {
9292
hasQ = true;
93-
if (event[iDau].id() == -mQuarkPdg)
93+
if ((event[iDau].y() > mQuarkRapidityMin) && (event[iDau].y() < mQuarkRapidityMax)) {
94+
atSelectedY = true;
95+
}
96+
}
97+
if (event[iDau].id() == -mQuarkPdg) {
9498
hasQbar = true;
95-
if ((event[iDau].y() > mQuarkRapidityMin) && (event[iDau].y() < mQuarkRapidityMax))
96-
atSelectedY = true;
99+
if ((event[iDau].y() > mQuarkRapidityMin) && (event[iDau].y() < mQuarkRapidityMax)) {
100+
atSelectedY = true;
101+
}
102+
}
97103
}
98104
if (hasQ && hasQbar && atSelectedY) {
99105
isGoodAtPartonLevel = true;

0 commit comments

Comments
 (0)